Articles | Volume 15, issue 12
Research article
23 Jun 2022
Research article |  | 23 Jun 2022

An optimal estimation-based retrieval of upper atmospheric oxygen airglow and temperature from SCIAMACHY limb observations

Kang Sun, Mahdi Yousefi, Christopher Chan Miller, Kelly Chance, Gonzalo González Abad, Iouli E. Gordon, Xiong Liu, Ewan O'Sullivan, Christopher E. Sioris, and Steven C. Wofsy

An optimal estimation-based algorithm is developed to retrieve the number density of excited oxygen (O2) molecules that generate airglow emissions near 0.76 µm (b1Σg+ or A band) and 1.27 µm (a1Δg or 1Δ band) in the upper atmosphere. Both oxygen bands are important for the remote sensing of greenhouse gases. The algorithm is applied to the limb spectra observed by the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) instrument in both the nominal (tangent heights below  90 km) and mesosphere–lower thermosphere (MLT) modes (tangent heights spanning 50–150 km). The number densities of emitting O2 in the a1Δg band are retrieved in an altitude range of 25–100 km near-daily in 2010, providing a climatology of O2 a1Δg-band airglow emission. This climatology will help disentangle the airglow from backscattered light in nadir remote sensing of the a1Δg band. The global monthly distributions of the vertical column density of emitting O2 in a1Δg state show mainly latitudinal dependence without other discernible geographical patterns. Temperature profiles are retrieved simultaneously from the spectral shapes of the a1Δg-band airglow emission in the nominal limb mode (valid altitude range of 40–100 km) and from both a1Δg- and b1Σg+-band airglow emissions in the MLT mode (valid range of 60–105 km). The temperature retrievals from both airglow bands are consistent internally and in agreement with independent observations from the Atmospheric Chemistry Experiment Fourier transform spectrometer (ACE-FTS) and the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), with the absolute mean bias near or below 5 K and root mean squared error (RMSE) near or below 10 K. The retrieved emitting O2 number density and temperature provide a unique dataset for the remote sensing of greenhouse gases and constraining the chemical and physical processes in the upper atmosphere.

1 Introduction

This study of upper atmospheric airglow from oxygen (O2) is driven by the need to measure O2 simultaneously with methane (CH4) and carbon dioxide (CO2) in satellite remote sensing. Methane and CO2 are two of the most important anthropogenic greenhouse gases, but the spatiotemporal variations in their sources and sinks are poorly understood, leading to significant uncertainties in projections of future climate trends (Miller et al.2007; Turner et al.2019; Friedlingstein et al.2020; Saunois et al.2020). Methane, in particular, has a global warming potential of about 56–105 times higher than that of CO2 for a 20-year time period (Howarth2014). Reducing methane emissions is among the most impactful actions that can be taken to reduce the global warming, which requires the identification of emission sources with a high degree of accuracy. Spaceborne observations offer a powerful tool in quantifying the spatiotemporal distributions of methane and CO2 and inferring their sources and sinks due to the extensive spatial and temporal coverage of satellites (Eldering et al.2017; Turner et al.2018; Lorente et al.2021). To separate the sources and sinks of methane and CO2 from variations in surface pressure and specific humidity, the abundances of CO2 and methane observed by satellites are usually represented by the column-averaged dry mole fractions (XCO2 and XCH4). Simultaneous observation of O2 is often necessary for deriving XCO2 and XCH4 and accounting for contamination of aerosol and clouds (Butz et al.2011).

The two widely used O2 absorption bands in greenhouse gas remote sensing are the b1Σg+X3Σg- band near 0.76 µm (O2 A band hereafter) and the a1ΔgX3Σg- band near 1.27 µm (O2 1Δ band hereafter). The O2 A band is commonly used in existing and planned spaceborne CO2 and methane sensors (Bovensmann et al.1999; Crisp et al.2017; Veefkind et al.2012; Moore et al.2018; Buchwitz et al.2013). The O2 absorption feature in the A band is much stronger than the shortwave infrared CO2 and methane bands and separated by significant spectral gaps, which challenges the spectral fitting and radiative transfer modeling. The O2 A band also overlaps with strong terrestrial solar-induced fluorescence that provides valuable information on plant photosynthesis but perturbs the O2 absorption features and, if not properly accounted for, leads to systematically biased greenhouse gas retrievals (Frankenberg et al.2012).

The O2 1Δ band plays an instrumental role in ground-based greenhouse gas remote sensing (Wunch et al.2011; Fu et al.2014; Frey et al.2019) and is adopted by the MicroCarb (Bertaux et al.2020) and MethaneSAT (Staebell et al.2021) satellite missions. Both O2 A and 1Δ bands feature upper-atmospheric airglow emissions (Zarboo et al.2018) that overlap with the backscattered solar light containing O2 absorption signals. The A-band airglow, which peaks in the mesopause region, has previously been neglected in nadir remote sensing due to its low intensity (Sioris2003), but there is a lack of quantitative assessment of how the A-band airglow impacts greenhouse gas retrieval. The 1Δ-band airglow emitted from the upper stratosphere and mesosphere is much stronger due to photochemically generated O2 in a1Δg state from ozone photolysis and was the main reason for choosing the A band over the 1Δ band by spaceborne sensors (Kuang et al.2002). However, recent advances in modeling the O2 1Δ band airglow spectra enable the disentangling of the airglow from backscattered light (Sun et al.2018a; Bertaux et al.2020) and open up the opportunity for spaceborne nadir remote sensing of the O2 1Δ band.

In order to account for the airglow contributions in the nadir-observed O2 spectra, it is crucial to have an accurate understanding of the spatial, temporal, and spectral distribution of airglow emissions. Nevertheless, observation-based studies of O2 A- and 1Δ-band airglow are rare and generally lack the information needed for nadir spaceborne remote sensing of greenhouse gases. Zarboo et al. (2018) retrieved volume emission rates (VERs) for the O2 A- and 1Δ-band airglow over the vertical range of 50–150 km using the mesosphere and lower thermosphere (MLT) limb observation mode of the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) instrument. These MLT observations could not capture a significant portion of the O2 1Δ-band airglow, which peaks near the stratopause, and the linear inversion applied by Zarboo et al. (2018) is subject to systematic biases in the retrieved airglow VERs below 90 km for the O2 A band and below 60 km for the O2 1Δ band due to O2 self-absorption (Sun et al.2018a). Li et al. (2020) conducted high vertical and along-track resolution retrieval of the O2 1Δ-band airglow for the purpose of retrieving mesospheric ozone using the Optical Spectrograph and InfraRed Imaging System (OSIRIS) limb radiance, which could not reveal the necessary horizontal distributions and spectral shape of the airglow. Sun et al. (2018a) and Bertaux et al. (2020) fit SCIAMACHY limb spectra in the O2 1Δ band with onion-peeling algorithms and retrieved airglow from selected orbits, although the spectral variation in airglow emission lines due to temperature was not fully incorporated.

The distribution of temperature in the atmospheric regions where airglow is emitted is important for reconstructing the emission spectra. The spectral band shape has been used to retrieve upper atmospheric temperature using CO2 (Boone et al.2005; Marshall et al.2011; García-Comas et al.2014), O2 A-band absorption (Nowlan et al.2007), and the O2 A-band airglow (Sheese et al.2010; Yang et al.2021). Using the O2 1Δ-band airglow spectra, Sun et al. (2018a) was able to retrieve mesospheric temperature above 60 km, limited by excessive uncertainties below it. The mass spectrometer incoherent scatter (MSIS) model provides temperature estimates throughout the atmosphere and is often used when observations are unavailable (Picone et al.2002).

To this end, we develop an optimal estimation-based algorithm to retrieve the airglow emissions and upper atmospheric temperature from limb-viewed radiance spectra. Compared with the onion-peeling algorithm in Sun et al. (2018a) that fits limb spectra from high to low tangent heights progressively, this algorithm combines spectra over a range of tangent heights and simultaneously retrieves the vertical profiles of local emissions and temperature in a consistent manner. The use of Bayesian inversion enables the incorporation of a priori knowledge, balancing of measurement error and prior error, and detailed posterior error analysis, including the averaging kernel matrix and degrees of freedom for signal (DOFS; Rodgers2000; Brasseur and Jacob2017). The airglow emission spectra are simulated based on the spectral model for the O2 1Δ band proposed by Sun et al. (2018a), which we demonstrate can be extended to the O2 A band with a simple generalization. We compare the retrieved temperatures in the upper stratosphere, mesosphere, and lower thermosphere with independent measurements. The algorithm is applied to 1 year of SCIAMACHY limb observations, including the MLT mode, to construct a climatology of O2 airglow and upper atmospheric temperature at 10:00 local solar time (LST).

2 Data

The SCIAMACHY limb spectra and MSIS model outputs are necessary for the retrieval algorithm. The retrieved temperature profiles from SCIAMACHY are compared with observations from Atmospheric Chemistry Experiment Fourier transform spectroscopy (ACE-FTS) and Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) instruments.


The SCIAMACHY instrument is an eight-channel grating spectrometer that measures radiation that is backscattered, reflected, transmitted, or emitted by the Earth's atmosphere and surface in limb and nadir geometries from 240 to 2380 nm. The instrument was launched on board the Envisat satellite, which was operational on a sun-synchronous orbit with an Equator-crossing time in the descending node of 10:00 LST from March 2002 until April 2012. In this study, only limb scattering measurements in spectral channels 4 (597–789 nm) and 6 (990–1750 nm) were used. In the nominal limb mode, SCIAMACHY observed the atmosphere from the surface up to 93 km in 2010. It also observed the MLT region, covering 50–150 km altitudes in the limb view over 2 d in each month, from July 2008 until April 2012. We conduct spectral and radiometric calibrations of SCIAMACHY spectra using version 3.2.6 of the SciaL1C command line tool software package. We then use a customized algorithm to convert the Level 1 data to the NetCDF4 file format. In contrast to previous studies (Bender et al.2017; Sun et al.2018a; Zarboo et al.2018; Bertaux et al.2020), where SCIAMACHY spectra were averaged across the track, we keep the eight across-track positions separate to enhance the horizontal space resolution. Limb observations over different altitudes at the same across-track position (usually there are 15 observing altitudes per across-track position for both nominal and MLT modes) are grouped as a single vertical sounding in order to retrieve atmospheric profiles. The latitude and longitude of the tangent point near 50 km altitude is used to represent the location of the entire vertical sounding, which gives a location ambiguity of less than  20 km. The average vertical distance between adjacent tangent heights within the same across-track position is about 6.6 km for both nominal and MLT modes. We have noticed that the eight vertical soundings across the track can always be grouped into four pairs. The soundings within each pair are close in latitude and longitude but zigzag in tangent heights. Therefore, the eight vertical soundings can be combined into four vertical soundings with a doubled vertical sampling of 3.3 km. We retrieve profiles at eight native across-track positions and only pair the retrieved profiles into four when calculating the column number density of emitting O2.

To cover the A band, radiance spectra in the 759 to 772 nm range are obtained from channel 4. In addition to airglow emissions, the observed radiance spectra contain photons from Rayleigh scattering and multiple scattering by the atmosphere and the surface. The scattering signal can be approximated by limb views at maximum tangent heights and contains the O2 absorption feature for the A band. No O2 absorption is observed for the 1Δ band due to much weaker atmospheric scattering and lower O2 absorption. To account for the scattered light in the A band, a background signal consisting of averaged A-band spectra at high tangent heights (130–150 km) from the same sounding is subtracted from each limb spectrum (Zarboo et al.2018). Before each subtraction, the background signal is scaled to match its out-of-band radiance with the out-of-band radiance of the limb spectrum to be corrected for. This step assures the out-of-band radiance centers at zero after correction. This correction assumes the spectral shape of scattered light is the same for all tangent heights and may leads to systematic errors at low tangent heights from the A band, where airglow emission is low, and the scattering path may differ significantly from the thermosphere. Following Zarboo et al. (2018), we consider 750–759 and 767–780 nm to be out of band. The O2 1Δ-band spectra are taken from 1240 to 1300 nm in channel 6 with two bad pixels being removed. The multiple scattering is considered negligible in the O2 1Δ band, and only a linear background fitted using the out-of-band radiance at 1210–1240 and 1300–1340 nm is subtracted from each limb spectrum.

2.2 MSIS model

The temperature, pressure, and the ground state O2 number densities are sampled from the NRLMSISE-00 model, a MSIS model developed at the U.S. Naval Research Laboratory (Picone et al.2002) using a Python package by Hirsch and Kastinen (2021). The MSIS temperature profile is used as a priori values in the optimal estimation. The pressure profile is fixed in the algorithm, whereas the relative changes in O2 number densities from MSIS values are retrieved, which affects the O2 self-absorption.


The ACE-FTS instrument was launched aboard the Canadian SCISAT-1 satellite in August 2003 and is still active. It operates from 750 to 4400 cm−1 (2.2 to 13.3 µm) at a high spectral resolution (0.02 cm−1). A sequence of atmospheric transmission spectra in the limb geometry is recorded using the Sun as an infrared source during sunrise and sunset (i.e., solar occultation). The temperature and pressure retrievals from the ACE-FTS were performed at 15–125 km by fitting a set of spectral windows containing CO2 lines with the inverse of temperature as a free parameter (Boone et al.2020). The temperature retrieval process was done by dividing the atmosphere into two altitude regimes with a crossover at about 50 km. The temperature retrieval was performed assuming a fixed CO2 mixing ratio up to about 50 km, while above it an empirical function was used to describe CO2 mixing ratio in order to force the results to exhibit smooth behavior. We use the temperature profiles from version 4.1 of ACE-FTS Level 2 data for the comparison of temperature retrieved from SCIAMACHY airglow. The ACE-FTS temperature profiles have been extensively used in the validation of profiles from other instruments such as SOFIE (Solar Occultation for Ice Experiment; Marshall et al.2011), MIPAS (García-Comas et al.2014), and OSIRIS (Sheese et al.2012).


The MIPAS instrument was a Fourier transform spectrometer for the detection of limb emission spectra in the middle and upper atmosphere on board the Envisat satellite since March 2002. It observed a wide spectral interval throughout the mid-infrared with a high spectral resolution of 0.025 cm−1. Operating in a wavelength range from 4.15 to 14.6 µm, MIPAS detected and spectrally resolved a large number of emission features of atmospheric minor constituents playing a major role in atmospheric chemistry. The retrieval of temperature is done from measurements of the CO2 atmospheric radiance at 15 µm for each MIPAS single limb scan. Because MIPAS was on the same platform as SCIAMACHY, there were abundant colocated observations, enabling monthly and zonally resolved temperature intercomparison. The measurement modes of MIPAS used in this study include the nominal measurement mode with an altitude coverage of roughly 6–70 km, the middle atmosphere (MA) mode covering 18–102 km, the upper atmospheric (UA) mode covering 42–172 km, and the noctilucent cloud (NLC) mode covering 39–102 km. The nominal measurement mode makes up the bulk of MIPAS measurements, whereas the MA and UA modes were available at least every 10 d, and the NLC mode only happened on a few days in 2010. We use the nominal temperature profiles from version 8 of MIPAS Level 2 data retrieved by the European Space Agency (ESA; Dinelli et al.2021). Version 8 data from the other modes are obtained through the Institute of Meteorology and Climate Research in cooperation with the Instituto de Astrofísica de Andalucía (IMK/IAA) retrieval algorithm (García-Comas et al.2012; Kiefer et al.2021). The typical total errors are 0.5–2 K below 70 km and 2–7 K above (for MA, UA, and NLC modes). The typical vertical resolutions in the comparison range of this study are 3–7 km.

3 Methodology

We assume the atmosphere as homogeneous layers and calculate airglow emission and O2 self-absorption using HITRAN (high-resolution transmission molecular absorption database; Rothman2021). The local emission and absorption from these layers are calculated and integrated along the line-of-sight to simulate the limb-viewed radiance. A Bayesian inversion is applied to retrieve airglow emission and temperature profiles by minimizing the difference between simulated and observed limb radiance spectra with a priori regularization.

3.1 Airglow emission from a single layer

An atmospheric layer bounded by two tangent heights of SCIAMACHY limb observations is the basic spatial resolving unit of this study. We create an additional layer above the outermost tangent height by assuming a layer thickness equal to the average difference between adjacent tangent heights. The atmospheric properties are assumed to be uniform within a layer, and the layer height is represented by its middle altitude when sampling meteorological parameters from the MSIS model and comparing with other observations. The excited O2 molecules that give rise to A and 1Δ band airglow are denoted as O2(b1Σg+) and O2(a1Δg), respectively, and generalized as O2. Their number densities ([O2]) within the layer are proportional to the VERs via the band-integrated Einstein A coefficients a as follows:

(1) VER = O 2 a .

Here VER is measured by photons per cubic centimeter per second (hereafter photons cm−3 s−1), [O2] is measured by molecules per cubic centimeter (hereafter molec. cm−3), and Einstein A coefficient is 2.27×10-4 s−1 for the 1Δ band (Sun et al.2018a) and 0.08693 s−1 for the A band (Long et al.2010). It is worth pointing out that [O2] is the number density that directly contributes to the airglow at the wavelengths of the O2 A band or 1Δ band and does not necessarily include all excited O2 in the b1Σg+ or a1Δg state. For example, the emission channel b1Σg+a1Δg, which corresponds to the so-called Noxon band at 1.91 µm (Noxon1961), is ignored. Nevertheless, this does not change the derived parameters.

The absorption and emission spectra in a layer can be resolved by either wavelength (λ) or wavenumber (ν=1/λ). We will use the subscript λ to denote a parameter resolved by wavelength and the subscript ν to denote a parameter resolved by wavenumber. We calculate the monochromatic absorption cross section of O2 using the Python library HITRAN Application Programming Interface (HAPI; Kochanov et al.2016) at a grid space of 0.0002 nm for the O2 A band and 0.001 nm for the 1Δ band, assuming the Voigt line shape profile. Note the monochromatic absorption cross section in centimeters squared per molecule (hereafter cm2 molec.−1) is the same at wavelength and wavenumber coordinates, i.e., σλ=σν, so we may simplify it as σ. We have tested both HITRAN2016 (Gordon et al.2017) and HITRAN2020 (Gordon et al.2022) spectroscopic databases and provide results using HITRAN2020 unless otherwise noted. The Jacobian of σ with respect to temperature T, σ/T, is calculated by finite difference of 0.01 K using HAPI.

The monochromatic VER (i.e., emissivity) for airglow in wavenumber space is denoted as εν with a unit of photons per cubic centimeter per second per centimeter inverse (hereafter photons cm−3 s−1 (cm−1)−1). Integrating εν across the band should recover the VER, as follows:

(2) VER = ν ε ν d ν .

According to the spectroscopic airglow model described in Sun et al. (2018a) and Bertaux et al. (2020), the O2 airglow emission spectrum εν is related to its absorption spectrum σ by the following:

(3) ε ν = c 0 σ ν 2 exp ( c 2 ν / T ) - 1 ,

where c0 is a scaling constant, and c2 is the second radiation constant in Planck's law with a value of 1.4387769 cm K. Combining Eqs. (1), (2), and (3), we can solve for c0 as follows:

(4) c 0 = [ O 2 ] a ν σ ν 2 exp ( c 2 ν / T ) - 1 d ν .

Plugging in Eq. (4) into Eq. (3), we obtain the following:

(5) ε ν = [ O 2 ] a t ν ν t ν d ν ,


(6) t ν = σ ν 2 exp ( c 2 ν / T ) - 1

is an intermediate function to facilitate Jacobian derivation in the following step.

The SCIAMACHY retrieval works in wavelength space, so it is more convenient to define airglow emissivity as a function of wavelength (i.e., using ελ instead of εν). The unit of ελ is photons per cubic centimeter per second per nanometer (hereafter photons cm−3 s−1 nm−1). Integrating ελ across the band in wavelength space should equally give the VER, as follows:

(7) VER = λ ε λ d λ .

ελ and εν are, hence, related via the following:

(8) ε λ = ε ν ν λ ,

and similar relationships apply for the Jacobians of ελ and εν. The Jacobian of airglow emissivity with respect to [O2] is as follows:

(9) ε λ [ O 2 ] = ε ν [ O 2 ] ν λ = a ν t ν λ ν t ν d ν .

The Jacobian of ελ with respect to temperature is more complicated and can be derived by differentiating Eq. (5) as follows:

(10) ε λ T = ε ν T ν λ = [ O 2 ] a t ν T ν t ν d ν - t ν ν t ν T d ν ν t ν d ν 2 ν λ .

Here, we leverage the following relationship:

(11) ν t ν d ν T = ν t ν T d ν .

Equation (10) is completed by plugging in the following:

(12) t ν T = σ T ν 2 exp ( c 2 ν / T ) - 1 + σ c 2 ν 3 exp ( c 2 ν / T ) T exp ( c 2 ν / T ) - 1 2 ,

which is the derivative of Eq. (6) with respect to T.

3.2 Forward model for limb-viewed airglow spectra

Assuming that we include N limb observations in a vertical sounding, limited by an altitude range, the number of atmospheric layers will also be N. We use index i to denote each limb view in a collection of limb observations, and hence, i ranges from 1 to N. Figure 1 illustrates the line of sight at tangent layer i=2 when there are N=7 layers (and equivalently, tangent heights) under consideration. We use index j to denote layers penetrated by the line of sight of the limb view i, and j ranges from i to N. The line of sight slices through those layers twice under the assumption of homogeneous layers. To uniquely identify each segment of the line of sight, another index l is introduced to denote those segments, from near to far, as shown in Fig. 1. Namely, a segment is a slant path of the line of sight through a layer. The segment index l ranges from 1 to 2Ni. A layer j can be uniquely identified for any segment l, as follows:

(13) j = l - N + i - 1 2 + i + 1 2 .

However, a given layer j will correspond to two segments (i.e., each layer is penetrated twice), as follows:

(14) l { N + j - i , N - j + 1 } .

Figure 1Illustration of the airglow limb-viewing line of sight. The observer is located at the left end of the horizontal line. The layer thicknesses are greatly exaggerated relative to the Earth's radius for visualization. The dashed vertical line indicates the tangent point and the thickness of each layer.


We denote the native resolution radiance spectrum at limb view i as rλ,i, where the subscript λ implies it is resolved in wavelength space. Its unit is photons per centimeter squared per second per nanometer per steradian (hereafter photons cm−2 s−1 nm−1 sr−1). To model this radiance, we need to sum emissions of all segments along the line of sight while considering absorption by downstream ground state O2 molecules, as follows:

(15) r λ , i = l = 1 2 N - i L i j ε λ , j 4 π exp - τ ̃ λ , i j - l = 1 l - 1 τ λ , i j .

Here the outermost summation loops from the nearest (l=1) to the farthest (l=2N-i) segment, and each segment l maps uniquely to a layer j through Eq. (13). Lij is the path length of the tangent view i through layer j in units of centimeters. Lijελ,j/(4π) is the unattenuated radiance emitted from such a segment. If there were no absorption, then the exponential term in Eq. (15) would be absent, and the observed limb radiance is simply the sum of radiances of all segments.

The exponential term in Eq. (15) represents the O2 self-absorption between the emitting segment l and the observer. The summation within the exponential term gives the total absorption optical depth of the downstream segments ranging from l=l-1, the segment immediately downstream, to l=1, the segment closest to the observer. Each downstream segment l again maps uniquely to a layer j, whose optical depth is given by the following:

(16) τ λ , i j = [ O 2 ] j σ λ , j L i j .

Here [O2]j is the O2 number density (in ground state, not to be confused with the emitting O2 or O2) in layer j, and σλ,j is the monochromatic absorption cross section of O2 at layer j (the subscript λ is added back to emphasize that σ is resolved in wavelength space).

The τ̃λ,ij term in Eq. (15) represents the self-absorption happening within the emitting segment. It has been neglected in the onion-peeling algorithm in Sun et al. (2018a), leading to large errors in spectral fitting at lower tangent heights where self-absorption is significant. We treat it here as an effective optical depth and as if all emitters are concentrated at one end of the segment. The following can be proven:

(17) τ ̃ λ , i j = - ln 1 - exp - [ O 2 ] j σ λ , j L i j [ O 2 ] j σ λ , j L i j .

Moreover, such an effective optical depth of the emitting segment approaches half of its optical depth as a transmitting-only segment (i.e., ([O2]jσλ,jLij)/2) when it is optically thin. The proof is given in Appendix A.

Equation (15) serves as the forward model for limb-viewed radiance, and the Jacobians, with respect to parameters to be retrieved (i.e., the state vector elements), can be obtained by finite difference. To reduce the computing cost, we derive the analytical Jacobians when possible. The derivations of analytical Jacobians with respect to temperature, emitting O2 number density, and the relative change in the ground state O2 number density are provided in Appendix B. The analytical Jacobians are also validated with finite difference. In addition, we include a squeeze factor of the SCIAMACHY instrument line shape (ILS) and a wavelength shift in the state vector. Before each inversion iteration, the radiance and Jacobians at all tangent heights are concatenated, convolved with the ILS, and sampled at the shifted SCIAMACHY wavelength grid. The Jacobians with respect to all state vector elements are then assembled into the full Jacobian matrix to be used in the inversion.

3.3 Optimal estimation-based inversion

After background subtraction, the limb-viewed radiance from tangent heights within the retrieval range is concatenated as the measurement vector r. The forward model described in the previous section maps the state vector x into the measurement space and is denoted in the following as f(x):

(18) r = f ( x ) + e ,

where e is the measurement error representing SCIAMACHY detector noise and any other inadequacy of the retrieval system in reproducing the limb-observed radiance. We denote the error covariance matrix of e as Sr and assume no correlation between individual measurements. The diagonal elements of Sr are calculated as the sum of a scaled radiance and the square of the readout noise. The scaled radiance term approximates the detector shot noise and model–data discrepancy, whose variances are assumed to scale linearly with radiance. The scaling factors applied to the radiance are tuned to balance the model–data discrepancy and retrieval–prior discrepancy. Their values are chosen to be 5×108 for the O2 1Δ band and 1×107 for the O2 A band. The readout noise is approximated by the standard deviation of the out-of-band radiance for each limb-viewed spectrum.

The retrieval algorithm minimizes the cost function as follows:

(19) χ 2 = r - f ( x ) T S r - 1 r - f ( x ) + x - x a T S a - 1 x - x a ,

where xa is the prior vector, and Sa is the prior error covariance matrix. A smaller prior error leads to heavier influence by the prior values and vice versa. For each retrieval, we first conduct a linear inversion of the spectra and use the vertical mean value of the inverted [O2] profile as the prior values for the [O2] profile. Note that the linear inversion ignores self-absorption and leads to inaccurate profile shape (Sun et al.2018a), but the purpose here is only to obtain a very rough estimate of [O2]. The corresponding prior error is set to be 100 times the prior, which is a constant value for all altitudes. This effectively gives no prior constraint to the [O2] profile and assures its information all comes from observations through near-unity DOFS of retrieved [O2]. The a priori temperature profile is sampled from the MSIS model. Retrieval of temperature becomes difficult for the O2 1Δ band in the stratosphere because of strong self-absorption and interference from upper layers. To stabilize the temperature retrieval, we apply a tighter prior error of 10 K below 50 km and relax to 30 K above 50 km. A logistic function with a steepness scale of 2.5 km is applied to smooth the transition at 50 km. Above 90 km, which is only observed in the MLT mode, the temperature prior error is further relaxed to 60 K to account for the potentially large observation–prior difference in the thermosphere. In the state vector, we also include changes in the O2 profile relative to the a priori from MSIS, which is equivalent to retrieving the natural log of [O2]. The prior values are set to be zero, i.e., no change to [O2] from MSIS, and the prior errors are set at 50 %. The O2 number density in the fitting is not constrained with temperature and pressure through the ideal gas law or hydrostatic balance. This lends some freedom to the O2 number density in deviating from the MSIS model and slightly enhances the goodness of fit in layers with large self-absorption. The DOFS of the retrieved ln ([O2]) also serves as a good qualitative indicator for the emergence of O2 self-absorption when going from high to low altitudes. The prior errors of each individual profile ([O2], temperature, and ln ([O2])) are correlated internally with a length scale of 1 scale height.

The forward model is nonlinear, so the cost function has to be minimized iteratively. We adopt the Levenberg–Marquardt modification of the Gauss–Newton method. On each iteration i, we solve for the state vector update dxi+1, as follows:

(20) d x i + 1 = ( 1 + γ ) S a - 1 + K i T S r - 1 K i - 1 K i T S r - 1 ( r - f ( x i ) ) - S a - 1 ( x i - x a ) ,

where Ki is the full Jacobian matrix (i.e., r/x) at iteration i, and the state vector is initialized by the a priori, i.e., x0=xa. γ is the Levenberg–Marquardt parameter that helps stabilize the iteration compared with the standard Gauss–Newton method. We initialize and update the value of γ and determine convergence following the OCO-2 retrieval algorithm (Crisp et al.2021).

After convergence is achieved, the posterior error covariance matrix can be calculated as follows:

(21) S ^ = K T S r - 1 K + S a - 1 - 1 ,

where K is the Jacobian matrix in the final iteration, and the averaging kernel matrix is as follows:

(22) A = S ^ K T S r - 1 K .

In this study, the DOFS of each state vector element refers to the corresponding diagonal element in the averaging kernel (Liu et al.2010).

4 Results

The algorithm described in the previous section is applied to retrieve O2 1Δ-band airglow from SCIAMACHY nominal limb observations and both 1Δ- and A-band airglow from SCIAMACHY MLT mode limb observations. For nominal limb observations, the measurement data are limited to tangent heights above 25 km and below 100 km (usually the uppermost tangent height for nominal limb scans is at  90 km). When the solar zenith angle (SZA) is above 70, the lower height bound is gradually increased to 40 km as the airglow vanishes at lower altitude. The A-band airglow retrieval is not attempted using the nominal limb observations due to limited coverage.

4.1 Retrieval demonstration using individual vertical soundings

Figure 2 demonstrates the spectral fitting of a vertical sounding of SCIAMACHY O2 1Δ-band nominal limb spectra on 3 January 2010 at 28.0 N, 99.5 E. In total, 10 limb observations (i.e., N=10) with tangent heights ranging from 28.4 to 87.4 km are included. Figure 2a shows the concatenated radiance spectra, and Fig. 2b shows the fitting residuals. The predicted 95 % confidence intervals, approximated by twice the measurement error, are also overlaid with the residuals. The measurement errors at the 1Δ band are dominated by the portion that scales with radiance, whereas readout noise that is constant across the band is relatively insignificant, except at the highest two tangent heights. This highlights the importance of appropriately considering the portion of measurement errors that changes with radiance. The measurement error would be grossly underestimated if only the out-of-band variance were considered. The goodness of fit is indicated by the χ2 value of 1.74 at convergence. Retrieval using the HITRAN2016 line list gives a slightly higher χ2 of 1.79 (result not shown), indicating a more accurate O2 1Δ-band spectroscopy in HITRAN2020. Indeed, the latest edition of the database makes use of very accurate experimental work that was recently carried out in Grenoble (Konefał et al.2020; Tran et al.2020) and at National Institute of Standards and Technology (NIST; Fleurbaey et al.2021). These measurements allowed improved spectroscopic models that efficiently decouple contributions of the magnetic dipole and much weaker but not negligible electric quadrupole transitions. The theoretical background is provided in Gordon et al. (2010) and Mishra et al. (2011).

Figure 2(a) SCIAMACHY limb-viewed spectra concatenated over 10 tangent heights within one vertical sounding (black), simulated concatenated spectra using the prior values of the state vector (blue), and simulated concatenated spectra using the posterior values of the state vector (red). Each limb-viewed 1Δ-band spectrum contains 77 data points, so the measurement vector length is 770. Altitudes of each tangent height are labeled above the corresponding spectra. The radiance unit is photons cm−2 s−1 nm−1 sr−1. (b) Residuals from the fit in panel (a) and the 95 % confidence intervals of residuals are shown.


The retrieved atmospheric profiles from the vertical sounding in Fig. 2 are shown in Fig. 3. The retrieved O2(a1Δg) number density is displayed in Fig. 3a, with error bars indicating twice the posterior error or approximately the 95 % confidence interval. As demonstrated in previous studies, O2 1Δ-band airglow concentrates in the lower mesosphere and upper stratosphere and peaks above the stratopause. The retrieved [O2(a1Δg)] becomes increasingly uncertain down to the stratosphere because (1) the lower the layer, the smaller the number of tangent views that can detect them; (2) those limb views include the contributions by emissions from all layers above, and (3) O2 self-absorption becomes significant. The nominal limb observations that stop at  90 km capture most of the upper tail of O2 1Δ-band airglow, as supported by comparison to the results of the MLT mode retrieval shown in Figs. 4 and 5 that extend beyond 120 km. The DOFS for the retrieved [O2(a1Δg)] is effectively unity at all retrieval altitudes, as shown in Fig. 3c. The retrieved temperature profile with posterior uncertainty is shown in Fig. 3b together with MSIS temperature that is used for prior values (green) and a colocated MIPAS temperature profile (black). The MIPAS observation is located at 28.7 N, 99.1 E with a 91 km spatial separation and a 15 min temporal separation from SCIAMACHY. In the mesosphere, the retrieved SCIAMACHY temperature deviates from prior values and matches closely with independent MIPAS observations. This lends confidence to the spectral modeling of the O2 1Δ-band airglow and the forward model for limb-viewed spectra. The DOFS for the retrieved temperature is close to unity above 50 km and quickly drops below 50 km (Fig. 3c), indicating heavier prior influence for the stratospheric temperature. This is by design as it is challenging to extract observational information for the stratosphere due to O2 self-absorption, whereas the prior stratospheric temperature is more trustworthy. Also shown in Fig. 3c is the DOFS for the relative changes in the O2 number density or, equivalently, ln ([O2]). The influence of O2 self-absorption is felt by the retrieval below 80 km and becomes significant below 60 km. Below 40 km, the DOFS drops due to the general loss of observational constraint. The DOFS for [O2(a1Δg)] does not appear to drop because its prior constraint is much lower.

Figure 3(a) Retrieved [O2(a1Δg)] profile from the spectral fitting in Fig. 2. (b) Retrieved temperature profile simultaneously with the [O2(a1Δg)] (pink), colocated MSIS model temperature (green), and colocated MIPAS observation (black). Error bars indicate twice the posterior error in panels (a) and (b). (c) DOFS for the retrieved [O2(a1Δg)], temperature, and ln ([O2]) profiles. Note that the profiles are plotted at the layer center altitude that is  3.3 km above each tangent height.


Figure 4Similar to Fig. 2 but using 10 limb views in an MLT vertical sounding.


Figure 5Similar to Fig. 3 but using an MLT vertical sounding.


The same O2 1Δ-band retrieval is applied to an MLT mode vertical sounding on 4 January 2010 at 55.8 N, 92.0 E. The spectral fitting is shown in Fig. 4. The lowest tangent height of this particular MLT sounding is at 57 km, and the upper limit is set to below 120 km, above which the airglow is negligible. The retrieval χ2 of 0.70 (0.74 if using HITRAN2016) is significantly lower than the nominal vertical sounding, indicating a better fitting quality. This is consistent with the fact that most challenges of the 1Δ-band retrieval are in the stratosphere, which is not observed in the MLT mode.

Figure 5 shows profiles retrieved from the spectral fit in Fig. 4. The portions below 90 km are qualitatively similar to the nominal mode retrieval in Fig. 3, although the MLT mode reveals further information at the mesopause region. There is a minor peak of [O2(a1Δg)] at near 90 km, whose maximum is  20 times lower than the major peak near the stratopause (comparing Figs. 5a and 3a). The nominal vertical soundings miss the upper part of the minor peak. Figure 5b compares retrieved MLT temperature with MSIS temperature (green) and colocated ACE-FTS temperature (black). The ACE-FTS observation is located at 55.5 N, 92.7 E, with a spatial separation of 56 km and a temporal separation of 1:57 h from SCIAMACHY. The retrieved temperature from SCIAMACHY 1Δ-band airglow agrees well with ACE-FTS in the upper mesosphere but becomes more uncertain above 95 km, where the 1Δ-band airglow becomes too dim to adequately reflect temperature from its spectral shape. The behaviors of DOFS in Fig. 5c are consistent with the nominal mode retrieval where they overlap and show the loss of observational constraint on temperature above 95 km due to fading airglow.

Figure 6Similar to Fig. 4 but showing the O2 A-band fitting using 12 limb views from the same MLT vertical sounding. Each limb-viewed A-band spectrum contains 62 data points, so the entire measurement vector has 744 elements.


The O2 A-band retrieval is applied to the same MLT mode vertical sounding and shown in Fig. 6. The upper limit of the tangent height is increased to 130 km as the A-band airglow extends further into the thermosphere. The airglow radiance observed at the A band is 2 orders of magnitude lower than the 1Δ band and closer to the readout noise level, which appears to dominate the total measurement uncertainty. The residual does not show significantly larger variance at high radiance, in contrast to the 1Δ band. Overall, the A-band airglow spectra shapes are well captured. The χ2 value for this fit is 0.77, with little difference between HITRAN2016 and HITRAN2020.

Figure 7 shows retrieved profiles from the A-band spectral fitting in Fig. 6. The emitting O2 number density ([O2(b1Σg+)]) can be retrieved with high precision above  90 km. The retrieval error quickly grows below the mesopause region (below 80 km) due to O2 self-absorption and loss of observational constraint. The temperature profiles from both A-band and 1Δ-band airglow are shown together with MSIS and ACE-FTS temperature profiles in Fig. 7b. Because the A-band airglow extends further into the thermosphere, the A-band-derived temperature shows less uncertainty than the 1Δ-band-derived temperature above the mesopause. Below the mesopause, the 1Δ-band-derived temperature shows better precision. Both SCIAMACHY-based temperature profiles capture the double mesopause feature observed by ACE-FTS and agree well between each other at 80–95 km, where both airglow are strong and less interfered. These indicate that the airglow spectral model can be reliably applied to both O2 A and 1Δ bands.

Figure 7Similar to Fig. 5 but in the A band. The A-band-derived emission of the O2 number density and temperature are shown in blue. In panel (b), the 1Δ-band-derived temperature, the same one shown in Fig. 5b, is overlaid in pink.


4.2 Climatology of O2 1Δ-band airglow and temperature retrieved from nominal limb mode observations

The [O2(a1Δg)] profiles are retrieved at eight vertical soundings across track. A numerical artifact appears when integrating the [O2(a1Δg)] profiles vertically, as four vertical soundings are defined at staggered altitudes compared to the other four. Therefore, we interweave the [O2(a1Δg)] profiles retrieved from the adjacent two across-track positions into one single profile and integrate the interwoven profile to obtain the column number density of O2(a1Δg). As a result, the O2(a1Δg) column number density is defined at four across-track positions instead of eight. The latitude and longitude of the paired across-track positions are averaged from the ones of the original vertical soundings. Figure 8 exemplifies the column number density of O2(a1Δg) calculated from 1 d of retrieval on 3 January 2010 using the SCIAMACHY nominal limb observations. Out of the 14 orbits in Fig. 8, one (orbit no. 41015) is highlighted to show the spatial extent of a single orbit. Some missing sounding locations are identifiable and due to failure to converge. These missing points are generally located at high latitudes and high solar zenith angles. In these transition regions between daytime and nighttime, the horizontal variation in airglow intensity is significant, which violates the homogeneous layer assumption for the retrieval algorithm. Retrieved data are often available for part of the ascending phase of the orbit at the summer hemisphere (most valid data are in the descending phase), leading to some repeated observations at the same latitude, although at different SZAs and potentially in the nighttime. To eliminate such a latitudinal ambiguity and nighttime data, we remove the ascending portion when averaging over multiple days and limit the SZA to within 100. The temperature sounding locations of MIPAS on the same day are shown in Fig. 8 as small gray and black symbols. The MIPAS temperature will be used for comparison in Sect. 5.

Figure 8Large symbols are colored by the O2(a1Δg) column number density retrieved from SCIAMACHY nominal limb orbits (nos. 41010–41023) on 3 January 2010. A single SCIAMACHY orbit (no. 41015) is highlighted in opaque colors with a black edge. Small gray symbols mark the temperature sounding locations from MIPAS that share orbits with SCIAMACHY. A single MIPAS orbit (no. 41023) is highlighted in black.

Figure 9Monthly average of O2(a1Δg) column number density retrieved from SCIAMACHY nominal limb orbits.

Figure 10Vertical–latitudinal distribution of [O2(a1Δg)]. Each data point in the retrieved profiles is binned to 3D grid boxes with sizes of 3×3×3.2 km and then averaged zonally.


We then apply the retrieval algorithm to the O2 1Δ-band spectra from all SCIAMACHY nominal limb observations over the year 2010. The use of HITRAN2020 improves the goodness of fit in mid- to high latitude but not in the tropical region, especially over land. This is not fully understood, as HITRAN2020 clearly reduces some structural residuals through more accurate line intensities. Figure 9 illustrates the global maps of monthly averaged O2(a1Δg) column number density in 2010. The sounding locations are binned to a 3 × 3 grid using the drop-in-the-box method (Sun et al.2018b). Although the single-day coverage is sparse (Fig. 8), continuous coverage is achievable at a monthly scale. In general, the value in each grid box is averaged from 5–10 sounding profiles. The maximum abundance of emitting O2 follows the amount of solar radiation closely and generally colocates with the subsolar latitude. The zonal variation is insignificant, and no other spatial pattern of O2 1Δ-band airglow can be convincingly identified. Figure 10 displays the vertical and latitudinal distribution of O2 1Δ-band airglow, represented by the zonal mean of retrieved [O2(a1Δg)] profiles for each month in 2010. The seasonal shift of airglow following the subsolar latitude is consistent with the spatial distribution of O2(a1Δg) column number density in Fig. 9. The vertical peak locations of [O2(a1Δg)] at 45–50 km does not appear to change significantly over the low and midlatitudes or at different times of the year. The spatial and temporal distributions of O2 1Δ-band airglow are in line with previous reports using SCIAMACHY and OSIRIS limb observations (Wiensz2005; Bertaux et al.2020; Li et al.2020), while providing a more complete characterization.

The O2 1Δ-band airglow is generated by several mechanisms (Zarboo et al.2018; He et al.2019; Bertaux et al.2020), and the most import mechanism of O2(a1Δg) production is the solar UV photolysis of ozone. The primary ozone layer exists in the stratosphere, resulting from the photolysis via the Herzberg continuum. A secondary ozone layer, generated by photolysis via the Schumann–Runge continuum, has been observed in the mesopause region (Smith and Marsh2005; Li et al.2020). The O2 1Δ-band airglow also features two peaks. The declining ozone concentration and increasing UV with altitude give rise to the primary peak of [O2(a1Δg)] at 45–50 km at low and midlatitudes. Figure 11 gives a closer look at the dependency of [O2(a1Δg)] profiles on the SZA. The retrieved [O2(a1Δg)] profiles north of the subsolar latitude over January 2010 are averaged over 2 SZA bins. The [O2(a1Δg)] profile shape shows little variation for SZA < 50, whereas the primary airglow peak shifts to higher altitudes as the SZA further increases. At the maximum observed SZA of 96 (left-most line in Fig. 11), the primary peak occurs at around 60 km. This is because, at large SZAs, the solar UV radiation no longer penetrates into the stratosphere and is mostly absorbed by ozone. The secondary [O2(a1Δg)] peak corresponding to the mesopause ozone layer is vaguely identifiable at 80–100 km and does not vary much over the SZA range shown in Fig. 11.

Figure 11Binned profiles of [O2(a1Δg)] for SZAs between 32–96 at 2 intervals during January 2010. Only data north of the subsolar latitude are used.


Figure 12Retrieved temperature using the 1Δ-band airglow observed by nominal limb orbits. Only temperature points with DOFS greater than 0.5 are included in the averaging, and only grid boxes with more than five temperature points are shown.


Vertical and latitudinal distributions of temperature retrieved from O2 1Δ airglow in the nominal limb mode is shown in Fig. 12 for each month in 2010. Only temperature retrievals with DOFS greater than 0.5 are included in the averaging, and only grid cells with more than five averaged measurements are shown, which limits the temperature profiles to above  40 km. The retrieved temperature profiles peak at the stratopause region (40–60 km) due to the ozone absorption of the solar UV radiation. The temperature minima in the mesopause region (80–100 km) are also largely captured. Temperature in the mesopause region is determined by a combination of radiative effects, including the absorptive heating of solar UV radiation by O2 and ozone (Mlynczak and Solomon1993) and radiative cooling of CO2 infrared emissions (Rodgers et al.1992), chemical effects of exothermic reactions of odd hydrogen and odd oxygen (Mlynczak and Solomon1993), and the dynamics of stratosphere–mesosphere, including turbulence heating (Fritts and Vanzandt1993; Lübken et al.1993) and the vertical transfer of heat by the upwind and gravity waves over the summer mesosphere. As shown in Fig. 12, the mesopause temperature is cold in polar summer and relatively warm in polar winter, driven by dynamical effects. In the polar summer, stronger radiative heating from ozone and O2 leads to upward air motion in the mesosphere, which in turn causes an adiabatic expansion of rising air and cools the mesopause. The upward motion that drives the cooling process is amplified by breaking upward gravity waves in the polar summer caused by the stratospheric easterly winds. An opposite process in the winter pole causes an adiabatic compression that results in a warm mesosphere (Björn1984; Smith2004). Another interesting feature in Fig. 12 is the implausibly high temperature in the northern winter at a high latitude and large SZA (generally larger than 90). A closer examination of the individual fits at those sounding locations reveals large residuals at lower tangent heights than 60 km, abnormally high χ2, and slow (if ever) convergence. The same region also has a large warm bias relative to ACE-FTS and MIPAS (to be shown in Sect. 5). We believe those are retrieval artifacts potentially due to the horizontal gradient of airglow intensity at sunrise or sunset conditions.

5 Comparison of temperature profiles

Here we focus on comparing the retrieved temperature profiles with other instruments, as temperature is the most challenging part of the state vector to be retrieved. Reliable temperature retrievals will support the validity of the algorithm and product. An accurate upper atmospheric temperature climatology will also help better simulate the airglow spectral shapes.

5.1 Internal temperature comparison using MLT mode orbits

For the MLT mode, it is possible to intercompare temperature retrievals at the same observing locations using both the O2 1Δ- and A-band airglow. Figure 13 compares temperature profiles retrieved by both bands in an MLT orbit (orbit no. 41252 on 19 January 2010) and colored by the altitude of retrieved layer. Figure 13a shows all data with DOFS larger than 0.1 over the full vertical range of profiles (roughly 50–120 km). The reddish points above 300 K are generally above 100 km and in the thermosphere. The other points below 300 K can be further classified into the mesopause region (80–100 km) and the mesosphere (below 80 km). Figure 13b magnifies the mesopause region with a much higher DOFS threshold of 0.8 and shows a tight correlation (Pearson correlation coefficient r=0.91) between 1Δ- and A-band temperatures. As illustrated in the retrieved profiles in Figs. 5 and 7, 80–100 km is the sweet spot where both 1Δ- and A-band airglow are strong and generally free from self-absorption. Although the 1Δ- and A-band retrievals use the same MSIS temperature as the prior, the information in the mesopause temperature is mostly from observations with high DOFS values. Out of the 504 temperature pairs with DOFS larger than 0.1 for both bands, 449 pairs have DOFS larger than 0.8. Below the mesopause region, the A-band-derived temperature becomes unreliable due to interference from self-absorption and loss of observational constraint. Over the full vertical range, the differences between those two temperatures are within the 95 % confidence intervals labeled by the horizontal and vertical error bars in Fig. 13, which indicates that the posterior error generated by the optimal-estimation-based algorithm adequately characterizes the uncertainty.

Figure 13(a) Temperatures retrieved from both 1Δ and A bands in an MLT mode orbit on 19 January 2010. Gray error bars indicate the 95 % confidence intervals approximated by twice the posterior error standard deviation. Data are filtered by DOFS greater than 0.1 for both bands and colored by altitude. The red dashed line is 1 : 1. Panel (b) is similar to panel (a) but only shows data between 80 and 100 km. Data with DOFS greater than 0.1 but less than 0.8 are in gray (N=504). Data with DOFS greater than 0.8 are colored by altitude (N=449). The Pearson correlation coefficient (r) and root mean squared error (RMSE) are also labeled.


Figure 14Comparisons between temperatures retrieved from 1Δ- and A-band airglow over the mesopause region. Each panel includes all MLT orbits in each month of 2010.


Figure 14 extends the temperature comparison between the O2 1Δ and A bands over the mesopause region to all MLT orbits in each month of 2010. The MLT mode observations are conducted on 2 d in each month, so there are  30 MLT orbits per month and 352 MLT orbits in 2010. The total number of data pairs, Pearson correlation coefficient, RMSE, and mean bias are labeled in the corresponding panel of each month. The mesopause temperature spans broader ranges in both polar winter and summer (January and December, as well as June and July) than in other months due to the dynamics-driven cold summer polar mesopause. Overall, the correlation coefficients of mesopause temperatures from the two bands ranges from 0.75 to 0.90, with the A-band temperature colder than the 1Δ-band temperature by 5–8 K. The low bias in the A-band temperature is likely caused by errors propagated from the lower altitude, where retrieving temperature from A-band airglow becomes challenging due to strong self-absorption. The RMSE between the two bands is 10–12 K.

5.2 MLT and nominal limb temperature comparison with ACE-FTS

The ACE-FTS temperature retrieval covers the vertical ranges of airglow-derived temperature in both MLT and nominal modes of SCIAMACHY and can be compared with all temperature retrievals in this study. With coincidence criteria of 500 km and 2 h, we identify  170 SCIAMACHY MLT vertical soundings and  1800 nominal vertical soundings that colocate with ACE-FTS soundings. The colocated soundings are sparse and skew toward mid- to high latitudes (70–35 S and 42–70 N). Figure 15a and b display the mean bias and RMSE within 3.3 km vertical bins between the temperature profiles retrieved from 1Δ (pink) and A (blue) bands relative to the colocated temperature profiles from ACE-FTS. The statistics for the prior temperature from MSIS are shown in green. The retrieved temperature data are limited to DOFS greater than 0.5. Compared with the prior MSIS temperature, both 1Δ and A-band-derived posterior temperatures demonstrate smaller RMSE relative to ACE-FTS, except at the high end (above 100 km for the A band and above 95 km for the 1Δ band), where the airglow weakens, and at the low end (below 72 km) for the A band, where the O2 self-absorption strongly reduces observational constraints (see Fig. 7). The mean bias for the 1Δ-band temperature stays generally less than ±5 K below 90 km. In the mesopause region (80–100 km), the A-band temperature outperforms the 1Δ-band temperature with lower absolute mean bias and RMSE relative to ACE-FTS. However, the quality of A-band temperature quickly degrades below 80 km. The relative mean biases are consistent with the internal airglow temperature comparisons in Sect. 5.1. To merge these two temperature retrievals from the MLT observations, one may use the 1Δ-band-derived temperature below 80 km, use the A-band-derived temperature above 90 km, and average the two with linearly varying weights between 80 and 90 km.

Figure 15The mean bias (a) and RMSE (b) of 1Δ band (pink) and A band (blue) temperature retrieved from the MLT orbits relative to the colocated ACE-FTS profiles are shown. Panels (c) and (d) show similar statistics between the 1Δ-band temperature retrieved from nominal limb orbits and the colocated ACE-FTS profiles. Retrievals using HITRAN2016 spectroscopic data are shown as dotted lines, and retrievals using HITRAN2020 are shown as solid lines with the same color.


A similar comparison is made between the 1Δ-band-derived temperature in the nominal limb orbits and ACE-FTS in Fig. 15c–d. The overall shapes of the mean biases of the 1Δ-band temperature are consistent between Fig. 15a and c, with the absolute mean bias in the nominal limb orbits slightly lower. The 1Δ-band-derived temperature agrees well with the ACE-FTS temperature (mean bias <±5 K and significant reduction of RMSE relative to the prior) from 55 to 90 km. The substantial warm bias relative to ACE-FTS seen below 55 km originates from the retrieval artifact in northern high latitude winter and is overly amplified because the SCIAMACHY vs. ACE-FTS colocations coincidentally cluster in northern high latitude winter. A more clear picture will be depicted in the latitudinally resolved comparison with MIPAS in the following subsection. Results using HITRAN2016 spectroscopic data are shown as dotted lines in Fig. 15, while results using HITRAN2020 are in solid lines. Overall, the 1Δ-band retrievals become warmer and show more pronounced artifacts near the northern high latitude stratopause when using HITRAN2020.

5.3 Nominal limb temperature comparison with MIPAS

Both MIPAS and SCIAMACHY were aboard the Envisat platform. Consequently, the MIPAS temperature soundings are spatiotemporally close to the SCIAMACHY soundings for all available orbits (see an example of 1 d SCIAMACHY and MIPAS orbits in Fig. 8). Here we focus on the MIPAS comparison to 1Δ-band-derived temperature for the nominal limb orbits. Figure 16 shows the latitudinally and vertically resolved mean bias between O2 1Δ-band-derived temperature and MIPAS nominal mode temperature. The mean bias is calculated by subtracting the interpolated, same-orbit MIPAS temperature from each SCIAMACHY temperature profile and then regridding the resultant bias profile in the same way that is shown in Fig. 12. To avoid interpolation artifacts, we use the extended profiles in MIPAS Level 2 data, which fill the space above the highest retrieval level with a seasonally and diurnally varying climatology. As a result, the comparison should be limited to below the top MIPAS nominal tangent height at 70 km. The mean bias is in general within ± 5 K in the mesosphere and upper stratosphere throughout the year, and the airglow-derived temperature tends to be warmer than MIPAS. García-Comas et al. (2014) compared MIPAS with a range of satellite and ground-based temperature observations and found that MIPAS temperature differs from others by 2 K at 50–80 km in spring, autumn, and winter at all latitudes and in summer at low- to midlatitudes. Differences between MIPAS and the other instruments in the summer high latitudes are typically smaller than 2 K at 50–65 km and 5 K at 65–80 km. MIPAS in general shows colder temperatures in the mid-mesosphere.

Figure 16Monthly mean biases between temperature retrieved from O2 1Δ airglow in nominal limb orbits and the MIPAS nominal mode temperature. MIPAS temperature in the same orbit is first interpolated to the latitude and altitude of SCIAMACHY layers. The longitudinal difference is neglected. The comparison should be limited to below 70 km, which is the maximum height for MIPAS nominal mode.

Figure 17Monthly RMSE between temperature retrieved from O2 1Δ airglow in nominal limb orbits and the MIPAS nominal mode temperature product.

Figure 18Similar to Fig. 16 but using MIPAS MA, UA, and NLC observation modes instead of the nominal mode.

A consistent cold bias up to 10 K exists in the Southern Hemisphere stratopause region. As the DOFS for 1Δ-band-derived temperature drops quickly below 50 km (Fig. 3), this bias stems mainly from a cold bias of the MSIS prior relative to MIPAS. The implausibly warm winter polar stratopause seen in Fig. 12 is more prominent in Fig. 16, with warm biases of over 20 K compared with MIPAS. Although it is clearly a retrieval artifact, we choose to keep those soundings because the temperature retrieval in the upper portion is still valid, and the retrieved emitting O2 number densities are valuable for information about the airglow distributions at large SZA. Figure 17 shows the RMSE between the 1Δ-band temperature and MIPAS temperature. The RMSE values are largely below 10 K throughout the upper stratosphere and mesosphere with significant portions below 5 K. The exceptions are the aforementioned southern stratopause and northern polar winter stratopause biases. The fact that the RMSE values are similar to the absolute mean bias values indicates low variances in the retrieved 1Δ-band temperature.

In addition, we compare the O2 1Δ-band-derived temperature with the MIPAS IMK/IAA temperature product for the MA, UA, and NLC modes, as shown in Fig. 18. The numbers of coincidence with SCIAMACHY for these modes are  20 % of the MIPAS nominal mode, but they provide coverage above 70 km through the top of SCIAMACHY nominal mode retrieval. The SCIAMACHY–MIPAS temperature difference is consistent with the MIPAS nominal mode, as in Fig. 16. The absolute temperature difference in the 70–100 km vertical range is generally within ±5 K, except the summer mesopause at northern high latitudes.

6 Conclusions and discussion

We develop an optimal-estimation-based algorithm to retrieve O2 1Δ- and A-band airglow, as well as temperature ranging from the upper stratosphere to lower thermosphere, using limb-viewed spectra by the SCIAMACHY instrument. The algorithm is applied to nominal limb orbits that are available near-daily to retrieve the 1Δ-band airglow and MLT limb orbits that appear on 2 d per month to retrieve both 1Δ- and A-band airglow in the year 2010. We demonstrate the monthly climatology for O2 1Δ-band airglow and temperature retrieved from nominal limb observations, which will provide crucial information for the consideration of airglow for future remote sensing of the O2 1Δ band. The global monthly distributions of the vertical column density of emitting O2 in the a1Δg state show mainly latitudinal dependence without other discernible geographical patterns. The O2 1Δ-band-derived temperature agrees well with MIPAS Level 2 temperature with mean bias generally near or within ±5 K and RMSE below 10 K. Notable discrepancies include the stratospheric cold bias attributable to a priori influence, strong warm bias at winter polar stratopause region likely due to the horizontal gradient of the airglow intensity that violates the homogeneous layer assumption for the retrieval algorithm, and disagreements at the summer mesopause in northern high latitudes. The reliable retrieval of temperature from airglow indicates that we can confidently reproduce the spectral shape of airglow emissions.

The successful retrieval of the O2 A-band airglow and the associated temperature demonstrates the generalizability of the airglow spectral model and the algorithm. The A-band- and 1Δ-band-derived temperature profiles are in good agreement in the mesopause region, with the A-band-derived temperature lower by 5–8 K. Intercomparisons with ACE-FTS suggest that the A-band-derived temperature is of higher quality above  90 km due to the sharper decline of 1Δ-band airglow with altitude. At lower altitudes, the 1Δ-band may provide more reliable temperature retrieval. Driven by the need for nadir observations of the 1Δ band by MethaneSAT, this study focuses on the 1Δ-band retrieval, whereas the synergy between the 1Δ- and A-band airglow may be further explored by future studies to improve understanding of the chemistry, radiation, and dynamics of the MLT region. It would also be interesting to include the Noxon band (Noxon1961) at 1.91 µm in future analyses. This band corresponds to the emission from the b1Σg+ to the a1Δg excited electronic states and is another avenue for depopulation of the photochemically produced oxygen molecules in the b1Σg+ state. It is allowed through the electric quadrupole mechanism (Gordon et al.2010).

The SCIAMACHY instrument made 10 years of near-daily nominal limb observations and 4 years of MLT-mode limb observations. It is a natural future step to process the entire SCIAMACHY lifetime to generate long-term records of O2 airglow emissions and upper atmospheric temperature. Comparison with MIPAS indicates that the current a priori temperature from the MSIS model may be inadequate, so it may be desirable to use more accurate reanalyzed temperature profiles, where available, and stitch the MSIS temperature above. A recent development improved the MSIS model especially at the MLT regions (Emmert et al.2021), and future adoption of the new MSIS profiles will improve the a priori estimates. In contrast to previous studies using the SCIAMACHY limb observations, we do not aggregate across the track, which leads to the options of keeping eight soundings across the track or interweaving each pair of adjacent soundings into four soundings across the track. The retrievals in this study are conducted at the native across-track resolution and, hence, a vertical resolution of 6.6 km. Future efforts will test the second option, which enhances the vertical resolution to 3.3 km, in line with previous SCIAMACHY limb studies (Bender et al.2017; Zarboo et al.2018).

Appendix A: Derivation of emitting layer optical depth

Radiance at one end of the emitting segment, with light path length L and uniform absorber/emitter concentrations, is as follows:

(A1) r λ = ε λ 4 π 0 L exp - [ O 2 ] σ λ ( L - L ) d L = ε λ L 4 π 1 - exp ( - [ O 2 ] σ λ L ) [ O 2 ] σ λ L .

The effective optical depth of the emitting segment (τ̃λ), assuming all emitters concentrated at the far end, should satisfy the following:

(A2) r λ = ε λ L 4 π exp - τ ̃ λ .

The optical depth of this segment as a transmitting layer is simply as follows:

(A3) τ λ = [ O 2 ] σ λ L .

Combining Eqs. (A1), (A2), and (A3) results in the following:

(A4) τ ̃ λ = - ln 1 - exp ( - τ λ ) τ λ .

This is equivalent to Eq. (17). At the optical thin limit, applying L'Hôpital's rule, leads to the following:

(A5) lim τ λ 0 τ ̃ λ τ λ = lim τ λ 0 - ln 1 - exp ( - τ λ ) τ λ τ λ = lim τ λ 0 τ λ exp ( - τ λ ) + exp ( - τ λ ) - 1 τ λ ( exp ( - τ λ ) - 1 ) = lim τ λ 0 - exp ( - τ λ ) + τ λ exp ( - τ λ ) - 2 exp ( - τ λ ) + τ λ exp ( - τ λ ) = 1 2 .

Namely, the effective optical depth of the emitting segment due to the self-absorption is half of its optical depth as a transmitting segment at optical thin limit, which makes intuitive sense. At the optical thick limit, the effective optical depth of emitting segment due to self-absorption approaches the natural logarithm of its optical depth as a transmitting segment, as follows:

(A6) lim τ λ τ ̃ λ = ln τ λ .
Appendix B: Jacobians of limb-viewed radiance with respect to temperature, [O2], and ln ([O2])

The Jacobian of radiance observed at limb view i with respect to the temperature of layer j (ji, otherwise the Jacobian values are all zero) is as follows:

(B1) r λ , i T j = l { N + j - i , N - j + 1 } L i j 4 π ε λ , j T j × exp - τ ̃ λ , i j - l = 1 l - 1 τ λ , i j + l { N + j - i , N - j + 1 } L i j ε λ , j 4 π - τ ̃ λ , i j T j × exp - τ ̃ λ , i j - l = 1 l - 1 τ λ , i j + l { N + j - i , N - j + 1 } l = l + 1 2 N - i L i j ε λ , j 4 π × - τ λ , i j T j exp - τ ̃ λ , i j - l ′′ = 1 l - 1 τ λ , i j ′′ .

Here the first term on the right-hand side reflects the sensitivity of airglow emission to the temperature of layer j (with two corresponding segments as in Eq. 14), the second term reflects the temperature sensitivity of self-absorption of layer j as an emitting layer, and the third term reflects the temperature sensitivity of layer j as a transmitting layer, where temperature alters its absorption from upstream emissions.

The temperature sensitivities of the transmitting optical depth (τλ/T) and its effective optical depth as an emitting layer (τ̃λ/T) for both layers can be related to the temperature sensitivity of the absorption cross section (σλ/T, obtained by finite difference using HAPI) by applying the chain rule to Eq. 16 and 17. We assume that the ground state O2 number density [O2] does not depend on the temperature and retrieval [O2] profile independently to account for errors in the a priori O2 profile. The temperature sensitivity of emissivity (ελ/T) is given by Eq. 10. The derived temperature Jacobian in Eq. B1 is confirmed with numerical finite difference derivative of Eq. 15 at the instrument resolution and sampling grid. The comparison is shown in Fig. B1. The analytical Jacobians are consistent with finite difference within 10−6 for all tangent heights.

The Jacobians of emitting O2 number density are derived similarly, but only involving the emitting segment, as follows:

(B2) r λ , i [ O 2 ] j = l { N + j - i , N - j + 1 } L i j 4 π ε λ , j [ O 2 ] j × exp - τ ̃ λ , i j - l = 1 l - 1 τ λ , i j ,

where ελ,j/[O2]j is the Jacobian of airglow emissivity with respect to [O2] for layer j and given by Eq. (9).

Besides, we also include relative changes to the ground state O2 number density in the state vector, which is equivalent to retrieving the natural logarithm of ground state O2 number density, ln ([O2]). The corresponding Jacobians are given by the following:

(B3) r λ , i ln ( [ O 2 ] j ) = [ O 2 ] j r λ , i [ O 2 ] j = [ O 2 ] j l { N + j - i , N - j + 1 } L i j ε λ , j 4 π × - τ ̃ λ , i j [ O 2 ] j exp - τ ̃ λ , i j - l = 1 l - 1 τ λ , i j + [ O 2 ] j l { N + j - i , N - j + 1 } l = l + 1 2 N - i L i j ε λ , j 4 π × - τ λ , i j [ O 2 ] j exp - τ ̃ λ , i j - l ′′ = 1 l - 1 τ λ , i j ′′ .

Here the first term on the right-hand side of Eq. (B3) reflects the sensitivity of the effective optical depth of layer j to its own O2 number density, and the second term reflects the sensitivity of its optical depth to its O2 number density as a transmitting layer. Both τλ/[O2]) and τ̃λ/[O2] can be readily calculated by differentiating Eqs. (16) and (17).

Figure B1Validation of analytical Jacobians of radiance with respect to temperature, based on Eq. (B1) (black lines), with finite difference (white dots). Each vertical column represents one limb view (fixed i) over all layers (different j). In total, 10 tangent views (i.e., N=10) are included using a single vertical sounding at the O2 1Δ band from SCIAMACHY orbit no. 41011 on 3 January 2010. Only Jacobians of layers above the tangent point are non-zero.


Code availability

The underlying software code can be accessed at (last access: 17 June 2022;, Sun2022a).

Data availability

The data for this paper can be accessed at (Sun2022b).

Author contributions

GGA and EO'S generated the SCIAMACHY Level 1 files. KS developed and implemented the forward model and retrieval algorithm with inputs from XL and CCM. IEG provided expertise in spectroscopy. CES provided expertise in upper atmospheric chemistry and observations and helped with scientific discussion and interpretation. XL, SCW, KC, and KS managed the project. KS and MY performed the temperature intercomparison and wrote the draft. All co-authors contributed to the editing of the paper.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We acknowledge MethaneSAT, LLC for supporting science and algorithms at the University at Buffalo and the Smithsonian Astrophysical Observatory. We thank Roman V. Kochanov, for the help with HAPI. The SCIAMACHY Level 1 data processing is supported by NASA Making Earth System Data Records for Use in Research Environments program. Some of the computations in this paper were conducted at the Center for Computational Research (UB CCR2022) at the University at Buffalo and the Smithsonian High Performance Cluster, Smithsonian Institution (SI HPC2022).

Financial support

This research has been supported by MethaneSAT, LLC, and NASA Earth Sciences Division (grant no. 80NSSC18M0091).

Review statement

This paper was edited by Christian von Savigny and reviewed by Miriam Sinnhuber and one anonymous referee.


Bender, S., Sinnhuber, M., Langowski, M., and Burrows, J. P.: Retrieval of nitric oxide in the mesosphere from SCIAMACHY nominal limb spectra, Atmos. Meas. Tech., 10, 209–220,, 2017. a, b

Bertaux, J.-L., Hauchecorne, A., Lefèvre, F., Bréon, F.-M., Blanot, L., Jouglet, D., Lafrique, P., and Akaev, P.: The use of the 1.27 µm O2 absorption band for greenhouse gas monitoring from space and application to MicroCarb, Atmos. Meas. Tech., 13, 3329–3374,, 2020. a, b, c, d, e, f, g

Björn, L.: The cold summer mesopause, Adv. Space Res., 4, 145–151,, 1984. a

Boone, C., Bernath, P., Cok, D., Jones, S., and Steffen, J.: Version 4 retrievals for the atmospheric chemistry experiment Fourier transform spectrometer (ACE-FTS) and imagers, J. Quant. Spectrosc. Ra., 247, 106939,, 2020. a

Boone, C. D., Nassar, R., Walker, K. A., Rochon, Y., McLeod, S. D., Rinsland, C. P., and Bernath, P. F.: Retrievals for the atmospheric chemistry experiment Fourier-transform spectrometer, Appl. Opt., 44, 7218–7231,, 2005. a

Bovensmann, H., Burrows, J. P., Buchwitz, M., Frerick, J., Noël, S., Rozanov, V. V., Chance, K. V., and Goede, A. P. H.: SCIAMACHY: Mission Objectives and Measurement Modes, J. Atmos. Sci., 56, 127–150,<0127:SMOAMM>2.0.CO;2, 1999. a

Brasseur, G. P. and Jacob, D. J.: Modeling of atmospheric chemistry, Cambridge University Press,, 2017. a

Buchwitz, M., Reuter, M., Bovensmann, H., Pillai, D., Heymann, J., Schneising, O., Rozanov, V., Krings, T., Burrows, J. P., Boesch, H., Gerbig, C., Meijer, Y., and Löscher, A.: Carbon Monitoring Satellite (CarbonSat): assessment of atmospheric CO2 and CH4 retrieval errors by error parameterization, Atmos. Meas. Tech., 6, 3477–3500,, 2013. a

Butz, A., Guerlet, S., Hasekamp, O., Schepers, D., Galli, A., Aben, I., Frankenberg, C., Hartmann, J.-M., Tran, H., Kuze, A., Keppel-Aleks, G., Toon, G., Wunch, D., Wennberg, P., Deutscher, N., Griffith, D., Macatangay, R., Messerschmidt, J., Notholt, J., and Warneke, T.: Toward accurate CO2 and CH4 observations from GOSAT, Geophys. Res. Lett., 38, L14812,, 2011. a

Crisp, D., Pollock, H. R., Rosenberg, R., Chapsky, L., Lee, R. A. M., Oyafuso, F. A., Frankenberg, C., O'Dell, C. W., Bruegge, C. J., Doran, G. B., Eldering, A., Fisher, B. M., Fu, D., Gunson, M. R., Mandrake, L., Osterman, G. B., Schwandner, F. M., Sun, K., Taylor, T. E., Wennberg, P. O., and Wunch, D.: The on-orbit performance of the Orbiting Carbon Observatory-2 (OCO-2) instrument and its radiometrically calibrated products, Atmos. Meas. Tech., 10, 59–81,, 2017. a

Crisp, D., O’Dell, C., Eldering, A., Fisher, B., Oyafuso, F., Payne, V., Drouin, B., Toon, G., Laughner, J., Somkuti, P., McGarragh, G., Merrelli, A., Nelson, R., Gunson, M., Frankenberg, C., Osterman, G., Boesch, H., Brown, L., Castano, R., Christi, M. Connor, B., McDuffie, J. Miller, C., Natraj, V., O’Brien, D., I., P., Smyth, M., Thompson, D., and Granat, R.: Level 2 full physics algorithm theoretical basis document, OCO-2 Level 2 Full Physics Retrieval ATBD, Documents, last access: 9 September 2021. a

Dinelli, B. M., Raspollini, P., Gai, M., Sgheri, L., Ridolfi, M., Ceccherini, S., Barbara, F., Zoppetti, N., Castelli, E., Papandrea, E., Pettinari, P., Dehn, A., Dudhia, A., Kiefer, M., Piro, A., Flaud, J.-M., López-Puertas, M., Moore, D., Remedios, J., and Bianchini, M.: The ESA MIPAS/Envisat level2-v8 dataset: 10 years of measurements retrieved with ORM v8.22, Atmos. Meas. Tech., 14, 7975–7998,, 2021. a

Eldering, A., Wennberg, P. O., Crisp, D., Schimel, D. S., Gunson, M. R., Chatterjee, A., Liu, J., Schwandner, F. M., Sun, Y., O'Dell, C. W., Frankenberg, C., Taylor, T., Fisher, B., Osterman, G. B., Wunch, D., Hakkarainen, J., Tamminen, J., and Weir, B.: The Orbiting Carbon Observatory-2 early science investigations of regional carbon dioxide fluxes, Science, 358, eaam5745,, 2017. a

Emmert, J. T., Drob, D. P., Picone, J. M., Siskind, D. E., Jones Jr., M., Mlynczak, M. G., Bernath, P. F., Chu, X., Doornbos, E., Funke, B., Goncharenko, L. P., Hervig, M. E., Schwartz, M. J., Sheese, P. E., Vargas, F., Williams, B. P., and Yuan, T.: NRLMSIS 2.0: A Whole-Atmosphere Empirical Model of Temperature and Neutral Species Densities, Earth Space Sci., 8, e2020EA001321,, 2021. a

Fleurbaey, H., Reed, Z. D., Adkins, E. M., Long, D. A., and Hodges, J. T.: High accuracy spectroscopic parameters of the 1.27 µm band of O2 measured with comb-referenced, cavity ring-down spectroscopy, J. Quant. Spectrosc. Ra., 270, 107684,, 2021. a

Frankenberg, C., O'Dell, C., Guanter, L., and McDuffie, J.: Remote sensing of near-infrared chlorophyll fluorescence from space in scattering atmospheres: implications for its retrieval and interferences with atmospheric CO2 retrievals, Atmos. Meas. Tech., 5, 2081–2094,, 2012. a

Frey, M., Sha, M. K., Hase, F., Kiel, M., Blumenstock, T., Harig, R., Surawicz, G., Deutscher, N. M., Shiomi, K., Franklin, J. E., Bösch, H., Chen, J., Grutter, M., Ohyama, H., Sun, Y., Butz, A., Mengistu Tsidu, G., Ene, D., Wunch, D., Cao, Z., Garcia, O., Ramonet, M., Vogel, F., and Orphal, J.: Building the COllaborative Carbon Column Observing Network (COCCON): long-term stability and ensemble performance of the EM27/SUN Fourier transform spectrometer, Atmos. Meas. Tech., 12, 1513–1530,, 2019. a

Friedlingstein, P., O'Sullivan, M., Jones, M. W., Andrew, R. M., Hauck, J., Olsen, A., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Canadell, J. G., Ciais, P., Jackson, R. B., Alin, S., Aragão, L. E. O. C., Arneth, A., Arora, V., Bates, N. R., Becker, M., Benoit-Cattin, A., Bittig, H. C., Bopp, L., Bultan, S., Chandra, N., Chevallier, F., Chini, L. P., Evans, W., Florentie, L., Forster, P. M., Gasser, T., Gehlen, M., Gilfillan, D., Gkritzalis, T., Gregor, L., Gruber, N., Harris, I., Hartung, K., Haverd, V., Houghton, R. A., Ilyina, T., Jain, A. K., Joetzjer, E., Kadono, K., Kato, E., Kitidis, V., Korsbakken, J. I., Landschützer, P., Lefèvre, N., Lenton, A., Lienert, S., Liu, Z., Lombardozzi, D., Marland, G., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Niwa, Y., O'Brien, K., Ono, T., Palmer, P. I., Pierrot, D., Poulter, B., Resplandy, L., Robertson, E., Rödenbeck, C., Schwinger, J., Séférian, R., Skjelvan, I., Smith, A. J. P., Sutton, A. J., Tanhua, T., Tans, P. P., Tian, H., Tilbrook, B., van der Werf, G., Vuichard, N., Walker, A. P., Wanninkhof, R., Watson, A. J., Willis, D., Wiltshire, A. J., Yuan, W., Yue, X., and Zaehle, S.: Global Carbon Budget 2020, Earth Syst. Sci. Data, 12, 3269–3340,, 2020. a

Fritts, D. C. and Vanzandt, T. E.: Spectral Estimates of Gravity Wave Energy and Momentum Fluxes. Part I: Energy Dissipation, Acceleration, and Constraints, J. Atmos. Sci., 50, 3685–3694,<3685:SEOGWE>2.0.CO;2, 1993. a

Fu, D., Pongetti, T. J., Blavier, J.-F. L., Crawford, T. J., Manatt, K. S., Toon, G. C., Wong, K. W., and Sander, S. P.: Near-infrared remote sensing of Los Angeles trace gas distributions from a mountaintop site, Atmos. Meas. Tech., 7, 713–729,, 2014. a

García-Comas, M., Funke, B., López-Puertas, M., Bermejo-Pantaleón, D., Glatthor, N., von Clarmann, T., Stiller, G., Grabowski, U., Boone, C. D., French, W. J. R., Leblanc, T., López-González, M. J., and Schwartz, M. J.: On the quality of MIPAS kinetic temperature in the middle atmosphere, Atmos. Chem. Phys., 12, 6009–6039,, 2012. a

García-Comas, M., Funke, B., Gardini, A., López-Puertas, M., Jurado-Navarro, A., von Clarmann, T., Stiller, G., Kiefer, M., Boone, C. D., Leblanc, T., Marshall, B. T., Schwartz, M. J., and Sheese, P. E.: MIPAS temperature from the stratosphere to the lower thermosphere: Comparison of vM21 with ACE-FTS, MLS, OSIRIS, SABER, SOFIE and lidar measurements, Atmos. Meas. Tech., 7, 3633–3651,, 2014. a, b, c

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

Gordon, I., Rothman, L., Hargreaves, R., Hashemi, R., Karlovets, E., Skinner, F., Conway, E., Hill, C., Kochanov, R., Tan, Y., Wcisło, P., Finenko, A., Nelson, K., Bernath, P., Birk, M., Boudon, V., Campargue, A., Chance, K., Coustenis, A., Drouin, B., Flaud, J., Gamache, R., Hodges, J., Jacquemart, D., Mlawer, E., Nikitin, A., Perevalov, V., Rotger, M., Tennyson, J., Toon, G., Tran, H., Tyuterev, V., Adkins, E., Baker, A., Barbe, A., Canè, E., Császár, A., Dudaryonok, A., Egorov, O., Fleisher, A., Fleurbaey, H., Foltynowicz, A., Furtenbacher, T., Harrison, J., Hartmann, J., Horneman, V., Huang, X., Karman, T., Karns, J., Kassi, S., Kleiner, I., Kofman, V., Kwabia-Tchana, F., Lavrentieva, N., Lee, T., Long, D., Lukashevskaya, A., Lyulin, O., Makhnev, V., Matt, W., Massie, S., Melosso, M., Mikhailenko, S., Mondelain, D., Müller, H., Naumenko, O., Perrin, A., Polyansky, O., Raddaoui, E., Raston, P., Reed, Z., Rey, M., Richard, C., Tóbiás, R., Sadiek, I., Schwenke, D., Starikova, E., Sung, K., Tamassia, F., Tashkun, S., Vander Auwera, J., Vasilenko, I., Vigasin, A., Villanueva, G., Vispoel, B., Wagner, G., Yachmenev, A., and Yurchenko, S.: The HITRAN2020 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 277, 107949,, 2022. a

Gordon, I. E., Kassi, S., Campargue, A., and Toon, G. C.: First identification of the a1ΔgX3Σg- electric quadrupole transitions of oxygen in solar and laboratory spectra, J. Quant. Spectrosc. Ra., 111, 1174–1183,, 2010. a, b

He, W., Wu, K., Feng, Y., Fu, D., Chen, Z., and Li, F.: The Radiative Transfer Characteristics of the O2 Infrared Atmospheric Band in Limb-Viewing Geometry, Remote Sens., 11, 2702,, 2019. a

Hirsch, M. and Kastinen, D.: space-physics/msise00: Matlab tests, enhance build, Zenodo [code],, 2021. a

Howarth, R. W.: A bridge to nowhere: methane emissions and the greenhouse gas footprint of natural gas, Energy Sci. Eng., 2, 47–60,, 2014. a

Kiefer, M., von Clarmann, T., Funke, B., García-Comas, M., Glatthor, N., Grabowski, U., Kellmann, S., Kleinert, A., Laeng, A., Linden, A., López-Puertas, M., Marsh, D. R., and Stiller, G. P.: IMK/IAA MIPAS temperature retrieval version 8: nominal measurements, Atmos. Meas. Tech., 14, 4111–4138,, 2021. a

Kochanov, R. V., Gordon, I. E., Rothman, L. S., Wcisło, P., Hill, C., and Wilzewski, J. S.: HITRAN Application Programming Interface (HAPI): A comprehensive approach to working with spectroscopic data, J. Quant. Spectrosc. Ra., 177, 15–30, 2016. a

Konefał, M., Kassi, S., Mondelain, D., and Campargue, A.: High sensitivity spectroscopy of the O2 band at 1.27 µm: (I) pure O2 line parameters above 7920 cm−1, J. Quant. Spectrosc. Ra., 241, 106653,, 2020. a

Kuang, Z., Margolis, J., Toon, G., Crisp, D., and Yung, Y.: Spaceborne measurements of atmospheric CO2 by high-resolution NIR spectrometry of reflected sunlight: An introductory study, Geophys. Res. Lett., 29, 11-1–11-4,, 2002. a

Li, A., Roth, C. Z., Pérot, K., Christensen, O. M., Bourassa, A., Degenstein, D. A., and Murtagh, D. P.: Retrieval of daytime mesospheric ozone using OSIRIS observations of O2 (a1Δg) emission, Atmos. Meas. Tech., 13, 6215–6236,, 2020. a, b, c

Liu, X., Bhartia, P. K., Chance, K., Spurr, R. J. D., and Kurosu, T. P.: Ozone profile retrievals from the Ozone Monitoring Instrument, Atmos. Chem. Phys., 10, 2521–2537,, 2010. a

Long, D. A., Havey, D. K., Okumura, M., Miller, C. E., and Hodges, J. T.: O2 A-band line parameters to support atmospheric remote sensing, J. Quant. Spectrosc. Ra., 111, 2021–2036,, 2010. a

Lorente, A., Borsdorff, T., Butz, A., Hasekamp, O., aan de Brugh, J., Schneider, A., Wu, L., Hase, F., Kivi, R., Wunch, D., Pollard, D. F., Shiomi, K., Deutscher, N. M., Velazco, V. A., Roehl, C. M., Wennberg, P. O., Warneke, T., and Landgraf, J.: Methane retrieved from TROPOMI: improvement of the data product and validation of the first 2 years of measurements, Atmos. Meas. Tech., 14, 665–684,, 2021. a

Lübken, F.-J., Hillert, W., Lehmacher, G., and von Zahn, U.: Experiments revealing small impact of turbulence on the energy budget of the mesosphere and lower thermosphere, J. Geophys. Res.-Atmos., 98, 20369–20384,, 1993. a

Marshall, B. T., Deaver, L. E., Thompson, R. E., Gordley, L. L., McHugh, M. J., Hervig, M. E., and Russell III, J. M.: Retrieval of temperature and pressure using broadband solar occultation: SOFIE approach and results, Atmos. Meas. Tech., 4, 893–907,, 2011. a, b

Miller, C. E., Crisp, D., DeCola, P. L., Olsen, S. C., Randerson, J. T., Michalak, A. M., Alkhaled, A., Rayner, P., Jacob, D. J., Suntharalingam, P., Jones, D. B. A., Denning, A. S., Nicholls, M. E., Doney, S. C., Pawson, S., Boesch, H., Connor, B. J., Fung, I. Y., O'Brien, D., Salawitch, R. J., Sander, S. P., Sen, B., Tans, P., Toon, G. C., Wennberg, P. O., Wofsy, S. C., Yung, Y. L., and Law, R. M.​​​​​​​: Precision requirements for space-based XCO2 data, J. Geophys. Res.-Atmos., 112, D10314,, 2007. a

Mishra, A., Balasubramanian, T., and Shetty, B.: Generalized electric quadrupole branch linestrengths for the infrared atmospheric oxygen bands, J. Quant. Spectrosc. Ra., 112, 2303–2309,, 2011. a

Mlynczak, M. G. and Solomon, S.: A detailed evaluation of the heating efficiency in the middle atmosphere, J. Geophys. Res.-Atmos., 98, 10517–10541,, 1993. a, b

Moore III, B., Crowell, S. M. R., Rayner, P. J., Kumer, J., O'Dell, C. W., O'Brien, D., Utembe, S., Polonsky, I., Schimel, D., and Lemen, J.: The Potential of the Geostationary Carbon Cycle Observatory (GeoCarb) to Provide Multi-scale Constraints on the Carbon Cycle in the Americas, Front. Environ. Sci., 6, 109,, 2018. a

Nowlan, C., McElroy, C., and Drummond, J.: Measurements of the O2 A- and B-bands for determining temperature and pressure profiles from ACE–MAESTRO: Forward model and retrieval algorithm, J. Quant. Spectrosc. Ra., 108, 371–388,, 2007. a

Noxon, J. F.: Observation of the b1Σg+​​​​​​​–a1Δg transition in O2, Can. J. Phys., 39, 1110–1119,, 1961. a, b

Picone, J. M., Hedin, A. E., Drob, D. P., and Aikin, A. C.: NRLMSISE-00 empirical model of the atmosphere: Statistical comparisons and scientific issues, J. Geophys. Res.-Space, 107, SIA 15-1–SIA 15-16,, 2002. a, b

Rodgers, C. D.: Inverse methods for atmospheric sounding: theory and practice, vol. 2, World scientific,, 2000. a

Rodgers, C. D., Taylor, F. W., Muggeridge, A. H., López-Puertas, M., and López-Valverde, M. A.: Local thermodynamic equilibrium of carbon dioxide in the upper atmosphere, Geophys. Res. Lett., 19, 589–592,, 1992. a

Rothman, L. S.: History of the HITRAN Database, Nature Reviews Physics, 3, 302–304,, 2021. a

Saunois, M., Stavert, A. R., Poulter, B., Bousquet, P., Canadell, J. G., Jackson, R. B., Raymond, P. A., Dlugokencky, E. J., Houweling, S., Patra, P. K., Ciais, P., Arora, V. K., Bastviken, D., Bergamaschi, P., Blake, D. R., Brailsford, G., Bruhwiler, L., Carlson, K. M., Carrol, M., Castaldi, S., Chandra, N., Crevoisier, C., Crill, P. M., Covey, K., Curry, C. L., Etiope, G., Frankenberg, C., Gedney, N., Hegglin, M. I., Höglund-Isaksson, L., Hugelius, G., Ishizawa, M., Ito, A., Janssens-Maenhout, G., Jensen, K. M., Joos, F., Kleinen, T., Krummel, P. B., Langenfelds, R. L., Laruelle, G. G., Liu, L., Machida, T., Maksyutov, S., McDonald, K. C., McNorton, J., Miller, P. A., Melton, J. R., Morino, I., Müller, J., Murguia-Flores, F., Naik, V., Niwa, Y., Noce, S., O'Doherty, S., Parker, R. J., Peng, C., Peng, S., Peters, G. P., Prigent, C., Prinn, R., Ramonet, M., Regnier, P., Riley, W. J., Rosentreter, J. A., Segers, A., Simpson, I. J., Shi, H., Smith, S. J., Steele, L. P., Thornton, B. F., Tian, H., Tohjima, Y., Tubiello, F. N., Tsuruta, A., Viovy, N., Voulgarakis, A., Weber, T. S., van Weele, M., van der Werf, G. R., Weiss, R. F., Worthy, D., Wunch, D., Yin, Y., Yoshida, Y., Zhang, W., Zhang, Z., Zhao, Y., Zheng, B., Zhu, Q., Zhu, Q., and Zhuang, Q.: The Global Methane Budget 2000–2017, Earth Syst. Sci. Data, 12, 1561–1623,, 2020. a

Sheese, P. E., Llewellyn, E. J., Gattinger, R. L., Bourassa, A. E., Degenstein, D. A., Lloyd, N. D., and McDade, I. C.: Temperatures in the upper mesosphere and lower thermosphere from OSIRIS observations of O2 A-band emission spectra, Can. J. Phys., 88, 919–925,, 2010. a

Sheese, P. E., Strong, K., Llewellyn, E. J., Gattinger, R. L., Russell III, J. M., Boone, C. D., Hervig, M. E., Sica, R. J., and Bandoro, J.: Assessment of the quality of OSIRIS mesospheric temperatures using satellite and ground-based measurements, Atmos. Meas. Tech., 5, 2993–3006,, 2012. a

SI HPC: High Performance Computing Cluster, Research Computing, Smithsonian Research Computing​​​​​​​,, 2022. a

Sioris, C. E.: Impact of the dayglow and the Ring effect on the retrieval of surface pressure from the A and B bands of O2: application to Orbiting Carbon Observatory, Tech. rep., Internal Report, Harvard-Smithsonian Center for Astrophysics, Cambridge, MA, USA, (last access: 17 June 2022), 2003. a

Smith, A. K.: Physics and chemistry of the mesopause region, J. Atmos. Sol.-Terr. Phy., 66, 839–857,, 2004. a

Smith, A. K. and Marsh, D. R.: Processes that account for the ozone maximum at the mesopause, J. Geophys. Res.-Atmos., 110, D23305,, 2005. a

Staebell, C., Sun, K., Samra, J., Franklin, J., Chan Miller, C., Liu, X., Conway, E., Chance, K., Milligan, S., and Wofsy, S.: Spectral calibration of the MethaneAIR instrument, Atmos. Meas. Tech., 14, 3737–3753,, 2021. a

Sun, K.: Kang-Sun-CfA/Methane: for the AMT paper (1.1), Zenodo [code],, 2022a. a

Sun, K.: Level 2 data for SCIAMACHY airglow retrieval in 2010, V1, Harvard Dataverse [data set],, 2022b. a

Sun, K., Gordon, I. E., Sioris, C. E., Liu, X., Chance, K., and Wofsy, S. C.: Reevaluating the Use of O2 a1Δg Band in Spaceborne Remote Sensing of Greenhouse Gases, Geophys. Res. Lett., 45, 5779–5787,, 2018a. a, b, c, d, e, f, g, h, i, j, k

Sun, K., Zhu, L., Cady-Pereira, K., Chan Miller, C., Chance, K., Clarisse, L., Coheur, P.-F., González Abad, G., Huang, G., Liu, X., Van Damme, M., Yang, K., and Zondlo, M.: A physics-based approach to oversample multi-satellite, multispecies observations to a common grid, Atmos. Meas. Tech., 11, 6679–6701,, 2018b. a

Tran, D. D., Tran, H., Vasilchenko, S., Kassi, S., Campargue, A., and Mondelain, D.: High sensitivity spectroscopy of the O2 band at 1.27 µm: (II) air-broadened line profile parameters, J. Quant. Spectrosc. Ra., 240, 106673,, 2020. a

Turner, A. J., Jacob, D. J., Benmergui, J., Brandman, J., White, L., and Randles, C. A.: Assessing the capability of different satellite observing configurations to resolve the distribution of methane emissions at kilometer scales, Atmos. Chem. Phys., 18, 8265–8278,, 2018. a

Turner, A. J., Frankenberg, C., and Kort, E. A.: Interpreting contemporary trends in atmospheric methane, P. Natl. Acad. Sci. USA, 116, 2805–2813,, 2019. a

UB CCR: CCR Facility Description, Center for Computational Research (CCR),, last access: 11 January 2022.  a

Veefkind, J., Aben, I., McMullan, K., Förster, H., de Vries, J., Otter, G., Claas, J., Eskes, H., de Haan, J., Kleipool, Q., van Weele, M., Hasekamp, O., Hoogeveen, R., Landgraf, J., Snel, R., Tol, P., Ingmann, P., Voors, R., Kruizinga, B., Vink, R., Visser, H., and Levelt, P.: TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83,, 2012. a

Wiensz, J. T.: Ozone retrievals from the oxygen infrared channels of the OSIRIS infrared imager, Master's thesis, Library and Archives Canada, (last access: 17 June 2022), 2005. a

Wunch, D., Toon, G. C., Blavier, J.-F. L., Washenfelder, R. A., Notholt, J., Connor, B. J., Griffith, D. W. T., Sherlock, V., and Wennberg, P. O.: The Total Carbon Column Observing Network, Philos. T. Roy. Soc. A, 369, 2087–2112,, 2011. a

Yang, X.-J., Wang, H.-M., Li, Y.-F., Wang, Y.-M., and Hu, X.-Q.: Temperature in the Near Space From the Emission Spectra of Oxygen A Band, Spectrosc. Spect. Anal., 41, 5, (last access: 4 May 2022)​​​​​​​, 2021. a

Zarboo, A., Bender, S., Burrows, J. P., Orphal, J., and Sinnhuber, M.: Retrieval of O2(1Σ) and O2(1Δ) volume emission rates in the mesosphere and lower thermosphere using SCIAMACHY MLT limb scans, Atmos. Meas. Tech., 11, 473–487,, 2018. a, b, c, d, e, f, g, h

Short summary
This study of upper atmospheric airglow from oxygen is motivated by the need to measure oxygen simultaneously with methane and CO2 in satellite remote sensing. We provide an accurate understanding of the spatial, temporal, and spectral distribution of airglow emissions, which will help in the satellite remote sensing of greenhouse gases and constraining the chemical and physical processes in the upper atmosphere.