Deriving column-integrated thermospheric temperature with the N2 Lyman–Birge–Hopfield (2,0) band
- Ann and H. J. Smead Department of Aerospace Engineering Sciences, University of Colorado Boulder, Boulder, CO, USA
Correspondence: Clayton Cantrall (email@example.com)
This paper presents a new technique to derive thermospheric temperature from space-based disk observations of far ultraviolet airglow. The technique, guided by findings from principal component analysis of synthetic daytime Lyman–Birge–Hopfield (LBH) disk emissions, uses a ratio of the emissions in two spectral channels that together span the LBH (2,0) band to determine the change in band shape with respect to a change in the rotational temperature of N2. The two-channel-ratio approach limits representativeness and measurement error by only requiring measurement of the relative magnitudes between two spectral channels and not radiometrically calibrated intensities, simplifying the forward model from a full radiative transfer model to a vibrational–rotational band model. It is shown that the derived temperature should be interpreted as a column-integrated property as opposed to a temperature at a specified altitude without utilization of a priori information of the thermospheric temperature profile. The two-channel-ratio approach is demonstrated using NASA GOLD Level 1C disk emission data for the period of 2–8 November 2018 during which a moderate geomagnetic storm has occurred. Due to the lack of independent thermospheric temperature observations, the efficacy of the approach is validated through comparisons of the column-integrated temperature derived from GOLD Level 1C data with the GOLD Level 2 temperature product as well as temperatures from first principle and empirical models. The storm-time thermospheric response manifested in the column-integrated temperature is also shown to corroborate well with hemispherically integrated Joule heating rates, ESA SWARM mass density at 460 km, and GOLD Level 2 column ratio.
Remote sensing of Earth's far ultraviolet (FUV) airglow from space provides important insights into the energetics, dynamics, and composition of the upper atmosphere (Meier, 1991; Paxton et al., 2017). The N2 Lyman–Birge–Hopfield (LBH) bands (∼127–280 nm) are prominent daytime FUV airglow features that emanate from the lower to middle thermosphere (∼120–200 km). Currently operating instruments measuring the LBH bands include the Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics (TIMED) satellite's Global Ultraviolet Imager (GUVI) launched in 2001 (Christensen et al., 2003), the Defense Meteorological Satellite Program's (DMSP) Special Sensor Ultraviolet Spectrographic Imager (SSUSI) launched in 2003 (Paxton et al., 2002), the Global-scale Observations of the Limb and Disk (GOLD) launched in 2018 (McClintock et al., 2020a), and the Ionospheric Connection Explorer's Far UltraViolet imaging spectrograph (FUV) launched in 2019 (Mende et al., 2017).
The utility of the LBH bands for probing thermospheric temperature was demonstrated by Aksnes et al. (2006) with limb observations by the Advanced Research and Global Observation Satellite's (ARGOS) High-resolution Ionospheric and Thermospheric Spectrograph (HITS) instrument. Eastes et al. (2008) subsequently showed that disk observations of LBH bands could be used for global monitoring of thermospheric temperature. These authors fit LBH laboratory spectra to observed emissions using an optimal estimation routine with varying parameters such as the N2 rotational temperature, population rates of each vibrational band, the NI 149.3 nm line emission intensity, O2 photoabsorption, and background emission rates. GOLD is the first mission to provide a Level 2 data product of thermospheric temperature (TDISK) using LBH disk emissions between ∼132–162 nm with a similar retrieval implementation (Eastes et al., 2017). Thermospheric temperatures have also been derived from TIMED GUVI observations (Zhang et al., 2019) using an intensity ratio between the (0,0) band and (1,0) band that the authors found to be quasi-linearly dependent on the N2 rotational temperature. The authors attributed the temperatures to the altitude at the peak of the LBH contribution function (∼155 km) based on radiative transfer calculations.
This paper presents a new technique to derive thermospheric temperature from spectrographic measurements of FUV airglow. The technique, unlike in past work, uses the ratio of two spectral channels that span a single LBH band to determine the change in band shape with respect to a change in the rotational temperature of N2. Section 2 provides background and exploration of the LBH temperature signal with principal component analysis (PCA) to motivate the new technique. Section 3 details the implementation and provides a discussion on the error sources and a rationale behind our interpretation of the derived temperature as a column-integrated property that we refer to as column-integrated thermospheric temperature, Tci. Section 4 presents the demonstrative results of applying the technique to GOLD Level 1C radiance data for the period of 2–8 November 2018, during which a moderate geomagnetic storm event occurred. The derived thermospheric temperatures are compared to the GOLD Level 2 version 3 TDISK data product over the same period. Due to the lack of independent remotely sensed or in situ temperature measurements in the lower to middle thermosphere, the derived column-integrated temperatures are also compared to (1) synthetically generated column-integrated temperatures from model simulations by NOAA's Whole Atmosphere Model (WAM) (Akmaev, 2011) and the Naval Research Laboratory Mass Spectrometer and Incoherent Scatter Radar Extended (NRLMSISE-00) (Picone et al., 2002) and (2) observations of other thermospheric states, including the GOLD Level 2 data product (Correira et al., 2018) and mass density by ESA's SWARM constellation (Astafyeva et al., 2017), as well as hemispherically integrated Joule heating rates estimated from SuperDARN and ground-based magnetometer data by using the Assimilative Mapping of Geospace Observations (AMGeO) (AMGeO Collaboration, 2019; Matsuo, 2020).
The thermospheric temperature signal exists in the rotational structure of the N2 LBH bands. In the case of thermospheric N2, the rotational temperature is equivalent to the ambient neutral temperature (Aksnes et al., 2006). To motivate the new approach for extracting this signal from the LBH (2,0) band, this section presents results from PCA performed on simulated LBH emissions. Synthetic LBH emission data are generated by forward modeling WAM simulation results for the period of 2–8 November 2018. WAM simulation experiments are executed with solar and geomagnetic forcing conditions specified according to the actual values of the F10.7, Kp, hemispheric power indices, solar wind velocities and densities, and interplanetary magnetic fields. Section 2.1 discusses forward modeling of LBH emissions, and Sect. 2.2 presents the PCA results.
2.1 Forward modeling of LBH emissions
The forward model used to produce synthetic LBH emissions is built with the Global Airglow Model (GLOW) and a radiative transfer model (Solomon, 2017). GLOW computes LBH volume emission rates as a function of altitude, which are input into the radiative transfer model to produce line-of-sight emissions of the LBH band system. The most important component of the forward model for the purposes of deriving thermospheric temperatures is the LBH vibrational–rotational band model (Budzien et al., 2001). The band model is a lookup table of laboratory spectra that specifies, for a given temperature, a unique spectrum for the upper vibrational states v–9 of N2. In the current implementation of the forward model, the v–9 vibrational population rates are those provided in Ajello et al. (2020), which are based on GOLD observations and are held constant. The population rate distribution can vary with the energy distribution of the electron flux in addition to variation in excitation sources other than direct excitation such as radiative cascade and collision-induced electronic transition (Ajello et al., 2020; Eastes, 2000a, b; Ajello and Shemansky, 1985). Ajello and Shemansky (1985) state that excitation thresholding should be included in airglow models to accurately reproduce LBH band intensity. However, as discussed in the following section, absolute band intensity is not needed to extract the N2 rotational temperature.
2.2 PCA of simulated LBH emissions
PCA is a data reduction technique that is useful for identifying the dominant orthogonal modes of variability from data. PCA is applied here using eigenvalue decomposition of a sample covariance matrix, Sλλ, of LBH emissions, , at wavelengths, λ, computed from aggregated data sets of simulated emissions of the LBH band system during 2–8 November 2018 for a total of samples.
is the mean LBH spectrum of the N samples. The useful results of PCA for this investigation are a set of eigenvectors (principal components), v, that describe the mode of variability in the LBH band system, with associated eigenvalues, σ. Suppose that v is an orthonormal set of spatiotemporally invariant basis and spatiotemporally dependent coefficients, c, represent the amplitude of the mode for each disk emission sample at a given time, ti, and location, ri, then can be expressed:
where is the residual after subtracting the mean and the sum of n weighted modes from . The total variance of c matches σ2 for that mode.
Figure 1 shows the mean of simulated LBH radiance, , between 138–162 nm generated with a spectral pixel size of 0.04 nm and a spectral resolution of 0.19 nm full width at half maximum (FWHM). The first two leading modes of variability in the spectrum, vB and vT, scaled by their eigenvalues or total standard deviations, σB and σT, are also shown. The leading mode vB is identified as the overall scaling of the LBH intensity. The value of suggests that this mode accounts for 98.3 % of the total variability in the simulated LBH spectra. The second leading mode vT is identified as the temperature signal. According to the value of this secondary mode accounts for 1.6 % of the total variability in the simulated LBH spectra. The correlation coefficient, R, between time-dependent coefficients for this temperature mode cT and the simulated WAM temperatures at 155 km altitude over the course of 2–8 November 2018 is 0.71. Together these two principal components account for 99.9 % of the variability in the simulated LBH spectra, suggesting the LBH system is highly compressible.
Figure 2 focuses on the LBH (2,0) band identified in Fig. 1. Compared to the other LBH bands the (2,0) band is relatively bright and is isolated from the even brighter atomic oxygen emissions at 130.4 and 135.6 nm and the atomic nitrogen emission at 149.3 nm (not shown in Fig. 1). The temperature signal in LBH emissions is apparent in the morphological shape of vT displayed in Fig. 2. As the rotational temperature, Tr, of N2 increases there is an effective skewing of the LBH (2,0) band to longer wavelengths. The close inspection of vT indicates that LBH (2,0) band emissions at wavelengths above 138.56 nm are positively correlated with temperature while those below are negatively correlated. Emissions at 138.56 nm are not affected by temperature variability and thus have zero amplitude in vT. This observation substantiates an approach of binning the LBH (2,0) band into two channels using 138.56 nm as a boundary to preserve how the temperature signal manifests in the LBH emission's morphological shape. Channel A is defined as the sum of all wavelengths negatively correlated with temperature (), and channel B contains all wavelengths positively correlated with temperature (). The two-channel ratio, , is a linear function of temperature. A similar two-channel-ratio approach was adopted in Cantrall et al. (2019) for testing the feasibility of assimilating GOLD Level 1C data into the WAM, but a justification of such an approach was not provided.
This section details the derivation of column-integrated thermospheric temperature, Tci, from the N2 rotational structure observed in top-of-atmosphere LBH emissions using the ratio of two channels that together span the LBH (2,0) band as motivated in Sect. 2. Section 3.1 explains the step-by-step procedure, followed by a discussion on potential error sources of Tci in Sect. 3.2 and analysis in Sect. 3.3 that support the interpretation of Tci as a column-integrated temperature rather than a temperature attributed to a specific altitude.
The procedure to determine Tci using the two-channel ratio consists of four steps as follows:
generate a set of synthetic LBH (2,0) bands at the instrument's spectral pixel size for a range of temperature using the vibrational–rotational band model (Budzien et al., 2001);
apply an instrument model on each synthetic band to account for the instrument's spectral resolution and spectral registration;
bin each band into channels A and B and fit the ratio, , to temperature by least squares;
compute the ratio, , from the observed LBH (2,0) band and determine Tci by regressing the observed ratio on the predetermined relationship between the ratio and temperature.
The two-channel ratio has a number of benefits; most importantly, it can limit the impact of the following uncertainties: (1) uncertainty associated with LBH excitation and extinction processes that affect the absolute intensity of each band and (2) uncertainty associated with instrument performance variations across the LBH band system. This technique to derive Tci only requires measurement of the relative magnitudes between two spectral channels (two spectral channels of size 0.5 nm) and a vibrational–rotational band model to map temperature to measurements. Measurement of a fully resolved, radiometrically calibrated LBH band system is not required and neither is a forward model to produce absolute LBH intensity.
3.2 Sources of error
There are two categories of error associated with determining physical parameters from observations: measurement error and representativeness error. The measurement error is the error associated with the measuring device, while the representativeness error is the difference between the observation and the physical model's representation of the observation (Rodgers, 2000). There are two dominant sources of systematic measurement error in Tci stemming from variations in the instrument's spectral registration and resolution. Figure 3 shows the error in Tci as a function of the error in the modeled spectral registration and the error in the modeled spectral resolution. It is apparent in Fig. 3 that a significant temperature error of about 50 K (5 %–10 %) can occur if the errors exceed a hundredth of a nanometer level for the spectral registration and a tenth of a nanometer level for spectral resolution. A discussion on mitigating these two sources of systematic measurement error when deriving Tci from GOLD data is provided in Sect. 4.1.
The predominant source of random measurement error that determines the precision in Tci is shot noise. The Tci random measurement error due to shot noise is quantified using Monte Carlo samples of simulated Tci derivations considering the instrument performance (McClintock et al., 2020a, b). Particle background counts are at times an additional random noise source. For the case study with GOLD data, the particle backgrounds were low as indicated by the “High_Background” flag in the Level 1C data, and therefore this error source is not considered. The statistics of background counts and the associated temperature errors should be quantified for the general application of this technique to any period.
Sources of representativeness error in deriving Tci are those that cause relative differences in the channel intensity other than temperature that are not captured in the vibrational–rotational band model. Photoabsorption by O2 is one source to consider. There is a 1.5 % difference in the mean O2 absorption cross section between the two channels that corresponds to a negligible difference in transmittance along the line of sight considering the O2 absorption cross-section variation with temperature. Another source of representativeness error associated with the (2,0) band is due to the overlap of the bright (2,0) transition and the weak (5,2) transition. Inaccurate specification of the v and v vibrational population rates cause a slight change in shape of the band with respect to the observations that could be interpreted as a change in the rotational temperature. Figure 8 in Ajello et al. (2020) provides the v–6 population rates and their uncertainties. These uncertainties are used to determine the associated error in the derived temperatures using the (2,0) band due to inaccurate specification of the v and v population rates. It is important to note that this representativeness error does not exist if the (1,1) or (2,3) bands are used in the derivation instead of the (2,0) band because the (1,1) and (2,3) bands are isolated from other LBH bands. However, these bands are also much weaker and suffer from significantly larger random error due to shot noise. Figure 4 shows the total random measurement error and representativeness error in Tci using the (2,0) band. The representativeness error is a function of temperature and can range from 15 K at Tci=400 to 48 K at Tci=1200 K. Random measurement error from shot noise is a function of the (2,0) band intensity with values of 30 and 70 K for photon counts of 1500 and 250, respectively.
3.3 Interpretation of column-integrated temperature
Interpretation of column-integrated temperature, Tci, is addressed using synthetic LBH disk emission observations generated by forward modeling WAM simulation results. The column-integrated temperature computed from synthetic observations is hereafter denoted as to contrast to computed from GOLD LBH disk emission data that is introduced later. To examine if can be attributed to a certain pressure we compare the WAM pressure level with the temperature that most closely matches , denoted as , to the pressure level at the peak of the LBH contribution function, pτ=1, where the LBH optical depth, τ, is unity. and pτ=1 are computed over the entire simulation period of 2–8 November 2018.
The LBH contribution function peak, pτ=1, changes with solar zenith angle (SZA) and observing zenith angle (OZA) as shown in Fig. 5. pτ=1 decreases in pressure (increases in altitude) for increases in SZA and OZA with a stronger dependence on SZA. Removing the OZA dependence, Fig. 6 shows there is a clear difference in pτ=1 and in their respective dependences on SZA ( ranges – hPa and pτ=1 ranges – when SZA ranges 5–70∘). The weaker SZA dependence of can be explained by the FWHM of the contribution function that spans ∼60 km at low SZA and ∼90 km for high SZA (Laskar et al., 2021). The contribution function acts as an averaging kernel for temperature over these large vertical widths that tends to reduce the SZA effect. The net result is derived temperatures that are generally hotter than temperatures at pτ=1 () for low SZA and temperatures that are generally cooler than temperatures at pτ=1 () for high SZA. Figure 6 also shows variability in (up to hPa or ∼10 km for the simulation conditions) at a given SZA that reflects considerable variability in the vertical temperature structure within the width of the contribution function given varying forcing conditions.
Figure 6 reinforces that Tci derived from the procedural steps specified in Sect. 3.1 is a column-integrated quantity, containing information from a larger altitude range of the lower-middle thermosphere than just at pτ=1. Perhaps, Tci can be justified to be attributed to pτ=1 when measurement and representativeness errors exceed the gap between Tci and the temperatures at pτ=1 at a given SZA and OZA. In general, specific pressure or altitude attribution of Tci requires additional a priori knowledge of the thermospheric temperature profile.
The two-channel-ratio approach to derive the column-integrated temperature is demonstrated using NASA GOLD Level 1C disk FUV emission data, , for the period of 2–8 November 2018 during which a moderate geomagnetic storm (Kp=5.8, nT) has occurred. Due to the lack of independent thermospheric temperature observations, the efficacy of this approach is validated through comparisons with GOLD Level 2 version 3 temperature product (TDISK). and TDISK are equivalent variables only differing in their approach. is also compared to two-channel column-integrated temperatures computed from the synthetic observations by forward modeling the WAM simulations () as described in Sect. 2, along with two-channel column-integrated temperatures computed from synthetic observations by forward modeling NRLMSISE-00 (TMSIS). The approach is further corroborated through comparisons of the storm-time changes of to hemispherically integrated Joule heating rates (QJH) estimated from SuperDARN and ground-based magnetometer data using AMGeO, ESA SWARM mass density measurements at 460 km () based on calculation from precise orbit determinations using the Global Positioning System receivers on the spacecraft, and GOLD Level 2 version 3 product (). Section 4.1 provides a description of the GOLD LBH Level 1C disk emission data used in the derivation; Sect. 4.2 presents results comparing with TDISK, TMSIS, and ; and Sect. 4.3 presents results comparing the storm-time response of with QJH, , and . Table 1 defines each of the variable symbols introduced above.
4.1 GOLD LBH disk emission data
GOLD observes the daytime FUV airglow from ∼134–162 nm on Earth's disk between 06:00 and 23:00 universal time (UT) from geostationary orbit at 47.5∘ W longitude. A full disk image is produced every ∼30 min at a spatial resolution of 125×125 km by alternating between scans of the Northern Hemisphere and Southern Hemisphere. The GOLD Level 1C radiance data with a spectral pixel size of 0.04 nm are used to derive in this study. The GOLD Level 1C data are spatially binned 2×2 (250×250 km spatial resolution) to improve the SNR by a factor of 2. Prior to deriving , efforts were made to reduce the impact of systematic biases that are present in version 3 of the GOLD Level 1C data product. Variations in spectral resolution along the GOLD detector are identified with the FWHM of the OI 135.6 doublet through fitting a 2-Gaussian distribution. Variations in the spectral registration are identified by differencing the modeled peak wavelength given the fitted OI 135.6 doublet FWHM by the peak wavelength determined by fitting a log-normal distribution to the (2,0) band. Note that the degradation of the detector due to the strength of the OI 135.6 doublet can cause errors in the spectral resolution estimate, but significant degradation had not occurred by 2–8 November 2018. Corrections for spectral registration and resolution are incorporated into Step 2 of the Tci algorithm (see Sect. 3).
4.2 Comparing to TDISK, TMSIS, and
Figure 7 displays along with TDISK, TMSIS, and over Earth's disk viewed by GOLD from 3–7 November 2018 at 15:00 UT, noon local time (LT) at the center of the disk (0∘ N, 47.5∘ W). A moderate geomagnetic storm commenced the evening of 4 November and lasted through 5 November (Gan et al., 2020). Figure 8 shows the mean bias difference (MBD) of from TDISK, TMSIS, and as a function of longitude (considering latitudes between ±10∘) and latitude (considering all longitudes viewed by GOLD) for 2–8 November 2018 at 15:00 UT. During this period the temperatures derived from observations (i.e., and TDISK) exhibit globally similar temperature amplitudes and display a similar morphological temperature response to geomagnetic activity over the disk. Note that there is slight banding near the Equator in and TDISK where the Southern Hemisphere and Northern Hemisphere scans meet that is likely due to systematic errors at the top and bottom edge of the detector that were not completely corrected. and TDISK show strong agreement near the center of the disk with an MBD less than 15 K (1 %–3 %) increasing to a maximum of ∼40 K (4 %–8 %) near the disk edge. The slope of –TDISK with respect to latitude and longitude indicates has a stronger south–north and west–east temperature gradient than TDISK. There is also agreement in the temperature morphology over the disk between and prior to the storm, but the storm-time response simulated by WAM, as manifest in , shows considerably higher temperatures in the mid- and high latitudes and a longer post-storm recovery time in comparison to and TDISK. – displays a similar west–east slope to –TDISK except for the region just west of the sub-solar point (−80 to −50∘ longitude) where is ∼25 K cooler than . TMSIS and show agreement in the temperature morphology over the disk, but TMSIS displays cooler temperatures, particularly just west of the subsolar point at low and mid-latitudes (up to 60 K), and a stronger west–east temperature gradient.
The and TDISK comparison is expanded in Fig. 9 to include all times in the range 07:00–22:00 UT for the period of 2–8 November 2018. Figure 9 shows that and TDISK have different behavior with the viewing conditions determined by SZA and OZA. increases with both SZA and OZA with a stronger trend for SZA. TDISK increases with OZA but remains relatively uniform with SZA, even decreasing slightly for SZA > 25∘. There are two likely explanations for the dependence of the derived temperatures on viewing conditions: (1) the derived temperatures reflect real temperature changes with viewing conditions because the contribution function is peaking at different pressures (Fig. 5); (2) the derived temperatures reflect temperature biases with viewing conditions because the LBH emission intensity is changing. Intensity decreases with increasing SZA due to reduced LBH excitation but increases with increasing OZA due to a larger air mass along the line of sight. To test which explanation best describes the dependence of and TDISK on viewing conditions, Fig. 9 is correlated to the pressure at the peak of the LBH contribution function, pτ=1, (Fig. 5) and to the mean LBH intensity measured by GOLD over the same period as a function of SZA and OZA. TDISK is weakly correlated () with pτ=1 and strongly correlated (R=0.72) with LBH intensity. In contrast, is strongly correlated () with pτ=1 and weak-moderately correlated () with LBH intensity. The stronger correlation between and pτ=1 compared to TDISK and pτ=1 and weaker correlation between and LBH intensity compared to TDISK and LBH intensity over this analysis period is suggestive that is more sensitive to real temperature changes as the probed pressures change with viewing conditions and less susceptible to biases due to a change in LBH intensity with viewing conditions. This is attributed to the fact that derivation does not require measurement of a fully resolved, radiometrically calibrated LBH band system nor a forward model to produce absolute LBH intensity. There are likely still biases in with LBH intensity as indicated by the weak-moderate correlation (), particularly at low intensities (high SZA) where shot noise can lead to positive biases up to 15 K in the two-channel-ratio approach.
4.3 Storm-time response
Figure 10 displays the response to the geomagnetic storm in , QJH, , and . , , and are shown as percent differences from the quiet-time conditions on 2 November 2018. The global temporal evolution of these variables is in good agreement with each other and consistent with known storm-time responses of thermospheric variables (e.g., Fuller-Rowell et al., 1994). A rise of magnetospheric energy influx as suggested by QJH leads to increased temperatures and upwelling of heavy molecular-rich air in the high and mid-latitudes as indicated by depletions of (−40 % near 50∘ latitude and −20 % near −50∘ latitude) and enhancements of (∼20 % near ±50∘ latitude) and (∼250 % near ±50∘ latitude). Enhancements of (20 %–30 % near 30∘ latitude) in the low latitudes suggest a subsequent development of downwelling following the pole-to-Equator global circulation in response to the storm-time Joule heating rise. Global thermospheric expansion is also apparent on 5 November as suggested by an increase in and over all latitudes. Note that the first detection of the temperature change was on the evening of 4 November when Joule heating rates have started to increase but are still relatively low (<50 GW). The post-storm recovery times are also in good agreement and appear to be on the order of 2–3 d.
A new technique to derive thermospheric temperature from space-based disk observations of FUV airglow is presented. The technique uses a ratio of the emissions in two spectral channels that together span the Lyman–Birge–Hopfield (LBH) (2,0) band to determine the change in band shape with respect to a change in the rotational temperature of N2. While this study focused on the LBH (2,0) band to derive thermospheric temperature, the described technique can be applied to any LBH band or combination of bands. The derived temperature from this technique is shown to be a column-integrated property referred to as column-integrated thermospheric temperature, Tci. Tci should not be attributed to the peak of the LBH contribution function without consideration of the viewing conditions and Tci derivation uncertainty. The definition of column-integrated thermospheric temperatures and other parameters used for comparison in the paper is given in Table 1. Specific findings of this work are as follows.
The LBH spectrum quantified with PCA of synthetic daytime LBH disk emission data is found to be highly compressible (two principal components explain 99.9 % of the variability). Analysis of the secondary principal component mode, which characterizes how the LBH temperature signal manifests as the change in band shape, substantiates the approach to bin an LBH spectral band into two channels such that the temperature-induced band shape change is best preserved. The study has shown that thermospheric temperatures can be derived from an observed two-channel ratio by using a precomputed relationship of the ratio to temperature from an LBH vibrational–rotational band model. In this two-channel-ratio approach, representativeness errors originating from forward modeling are reduced because radiometrically calibrated LBH band intensities are not required in the derivation procedure, and negative impacts of systematic measurement errors, stemming from variations across the band system in the instrument's spectral registration and resolution, are reduced because a fully resolved LBH band system is not required.
The derived temperature from the two-channel approach can have significant systematic biases of about 50 K (5 %–10 %) if the spectral registration and resolution are not known to the hundredth of a nanometer level and tenth of a nanometer level, respectively, as shown in Fig. 3. In addition to these known sources of systematic biases, there is intrinsic random error in Tci due primarily to shot noise and representativeness error due to misspecification of the and population rates in the vibrational–rotational band model. The random measurement error is estimated to be 20–60 K (3 %–9 %), and the representativeness error is estimated to be 15–30 K (2 %–5 %) for the case study with GOLD L1C data.
For the period of 2–8 November 2018 during which a moderate geomagnetic storm has occurred, the temperatures derived from observations (i.e., and TDISK) exhibit globally similar temperatures. is in good agreement with and TDISK at low latitudes but exhibits considerably higher temperatures at mid- and high latitudes during the storm response. TMSIS exhibits globally cooler temperatures to the observations. However, there are clear differences between and TDISK with respect to viewing conditions. There is stronger correlation between and pτ=1 () compared to TDISK and pτ=1 () and weaker correlation between and LBH intensity () compared to TDISK and LBH intensity (R=0.72) over the analysis period. These differences highlight a potential benefit of the two-channel-ratio approach to reduce the representativeness error by measurement of the relative intensities between two channels that only requires a vibrational–rotational band model for the forward model instead of a full radiative transfer model. The temporal evolution of global Tci corroborates well with temporal changes of hemispherically integrated Joule heating rates QJH, SWARM mass density at 460 km , and GOLD , which is consistent with known storm-time responses of thermospheric variables.
The lookup table for the two-channel ratio versus N2 rotational temperature considering the GOLD spectral registration and resolution variation along the detector used to derive column-integrated temperatures and the resulting column-integrated temperatures for the period of 2–8 November 2018 presented in this paper are available at https://doi.org/10.17605/OSF.IO/KHNQ7 (Cantrall and Matsuo, 2021). GOLD L1C and L2 data can be accessed at the GOLD Science Data Center (http://gold.cs.ucf.edu/search/; NASA, 2021a) and at NASA's Space Physics Data Facility (https://spdf.gsfc.nasa.gov; NASA, 2021b). The code for NOAA's WAM model is available at https://github.com/NOAA-SWPC/WAM (NOAA-SWPC, 2021). The NRLMSISE-00 neutral atmosphere model is available from the NASA CCMC at https://ccmc.gsfc.nasa.gov/modelweb/models/nrlmsise00.php (CCMC, 2021). The Python interface for the NRLMSISE-00 neutral atmosphere model is available at https://github.com/st-bender/pynrlmsise00 (Bender, 2021). Near-Earth solar wind data are provided by the Goddard Space Flight Center Space Physics Data Facility and are available at https://omniweb.gsfc.nasa.gov/ (NASA, 2021c). The density measurements (L2 DNSxPOD data product) from Swarm can be obtained at https://earth.esa.int/web/guest/swarm/data-access (ESA, 2021) upon registration. AMGeO is an open-source software available from https://amgeo.colorado.edu (AMGeO, 2021) upon registration. SuperMAG ground magnetometer data are available at https://supermag.jhuapl.edu/ (SuperMAG, 2021). SuperDARN radar data are available at http://vt.superdarn.org (VT, 2021).
CC developed the presented technique and performed the analyses. TM contributed AMGeO, determined the validation approach, and provided interpretation of the analyses. CC and TM prepared the manuscript.
The contact author has declared that neither they nor their co-author has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors acknowledge Wlliam E. McClintock for his assistance with the GOLD data, Stanley Solomon for his assistance with the use of GLOW model, Adam Kubaryk for his assistance with the use of WAM, and Liam Kilcommons for his assistance with the use of AMGeO.
AMGeO is supported by the NSF EarthCube awards ICER 1928403, ICER 1928327, and ICER 1928358. The authors acknowledge the use of SuperDARN data. SuperDARN is a collection of radars funded by national scientific funding agencies of Australia, Canada, China, France, Italy, Japan, Norway, South Africa, the United Kingdom, and the United States of America. For SuperMAG data we are grateful for INTERMAGNET, Alan Thomson; CARISMA, Ian Mann; CANMOS, Geomagnetism Unit of the Geological Survey of Canada; the S-RAMP Database, Kiyohumi Yumoto and Kazuo Shiokawa; the SPIDR database; AARI, Oleg Troshichev; the MACCS program, Mark Engebretson; GIMA; MEASURE, UCLA IGPP and Florida Institute of Technology; SAMBA, Eftyhia Zesta; 210 Chain, Kiyohumi Yumoto; SAMNET, Farideh Honary; IMAGE, Liisa Juusola; Finnish Meteorological Institute, Liisa Juusola; Sodankylä Geophysical Observatory, Tero Raita; UiT the Arctic University of Norway, Tromsø Geophysical Observatory, Magnar G. Johnsen; GFZ German Research Centre For Geosciences, Jürgen Matzka; Institute of Geophysics, Polish Academy of Sciences, Anne Neska and Jan Reda; Polar Geophysical Institute, Alexander Yahnin and Yarolav Sakharov; Geological Survey of Sweden, Gerhard Schwarz; Swedish Institute of Space Physics, Masatoshi Yamauchi; AUTUMN, Martin Connors; DTU Space, Thom Edwards and Anna Willer; South Pole and McMurdo Magnetometer, Louis J. Lanzarotti and Alan T. Weatherwax; ICESTAR; RAPIDMAG; British Artarctic Survey; McMac, Peter Chi; BGS, Susan Macmillan; Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation (IZMIRAN); MFGI, Balazs Heilig; Institute of Geophysics, Polish Academy of Sciences, Anne Neska and Jan Reda; University of L'Aquila, Massimo Vellante; BCMT, Vincent Lesur and Aude Chambodut; data obtained in cooperation with Geoscience Australia, Andrew Lewis; AALPIP, co-principal investigators Bob Clauer and Michael Hartinger; SuperMAG, Jesper W. Gjerloev; and data obtained in cooperation with the Australian Bureau of Meteorology, Richard Marshall.
This research has been supported by the National Aeronautics and Space Administration (grant no. 80NSSC19K1432) and the National Science Foundation (grant no. AGS-1848544).
This paper was edited by Jorge Luis Chau and reviewed by two anonymous referees.
Ajello, J. and Shemansky, D.: A reexamination of important N2 cross sections by electron impact with application to the dayglow: The Lyman–Birge–Hopfield Band System and NI (119.99 nm), J. Geophys. Res.-Atmos., 90, 9845–9861, https://doi.org/10.1029/JA090iA10p09845, 1985.
Ajello, J. M., Evans, J. S., Veibell, V., Malone, C. P., Holsclaw, G. M., Hoskins, A. C., Lee, R. A., McClintock, W. E., Aryal, A., Eastes, R. W., and Schneider, N.: The UV spectrum of the Lyman–Birge–Hopfield band system of N2 induced by cascading from electron impact, J. Geophys. Res.-Space, 125, e2019JA027546, https://doi.org/10.1029/2019JA027546, 2020.
Akmaev, R. A.: Whole atmosphere modeling: Connecting terrestrial and space weather, Rev. Geophys., 49, RG4004, https://doi.org/10.1029/2011RG000364, 2011.
Aksnes, A., Eastes, R., Budzien, S., and Dymond, K.: Neutral temperatures in the lower thermosphere from N2 Lyman–Birge–Hopfield (LBH) band profiles, Geophys. Res. Lett., 33, L15103, https://doi.org/10.1029/2006GL026255, 2006.
AMGeO: AMGeO Maps, available at: https://amgeo.colorado.edu, last access: 25 October 2021.
AMGeO Collaboration: A Collaborative Data Science Platform for the Geospace Community: Assimilative Mapping of Geospace Observations (AMGeO) v1.0.0, Zenodo, https://doi.org/10.5281/zenodo.3564914, 2019.
Astafyeva, E., Zakharenkova, I., Huba, J. D., Doornbos, E., and van den Ijssel, J.: Global Ionospheric and thermospheric effects of the June 2015 geomagnetic disturbances: Multi-instrumental observations and modeling, J. Geophys. Res.-Space, 122, 11716–11742, https://doi.org/10.1002/2017JA024174, 2017.
Bender, S.: pynrlmsise00, GitHub [code], https://github.com/st-bender/pynrlmsise00, last access: 25 October 2021.
Budzien, S. A., Fortna, C. B., Dymond, K. F., Thonnard, S. E., Nicholas, A. C., McCoy, R. P., and Thomas, R. J.: Thermospheric temperature derived from ARGOS observations of N2 Lyman–Birge–Hopfield emission, in: Eos Trans. AGU, 82, Spring Meet. Suppl., Abstract SA62A-03, 2001.
Cantrall, C. and Matsuo, T.: Deriving column-integrated thermospheric temperature with the N2 Lyman-Birge-Hopfield (2,0) band, OSF [data set], https://doi.org/10.17605/OSF.IO/KHNQ7, 2021.
Cantrall, C. E., Matsuo, T., and Solomon, S. C.: Upper atmosphere radiance data assimilation: A feasibility study for GOLD far ultraviolet observations, J. Geophys. Res.-Space, 124, 8154–8164, https://doi.org/10.1029/2019JA026910, 2019.
CCMC: NRLMSISE-00 Atmosphere Model, CCMC [code], https://ccmc.gsfc.nasa.gov/modelweb/models/nrlmsise00.php, last access: 27 October 2021.
Christensen, A. B., Paxton, L.J., Avery, S., Craven, J., Crowley, G., Humm, D. C., Kil, H., Meier, R. R., Meng, C.-I., Morrison, D., Orgorzalek, B. B., Straus, P., Strickland, D. J., Swenson, R. M., Walterscheid, R. L., Wolven, B., and Zhang, Y.: Initial observations with the Global Ultraviolet Imager (GUVI) in the NASA TIMED satellite mission, J. Geophys. Res.-Space, 108, 1451, https://doi.org/10.1029/2003JA009918, 2003.
Correira, J, Evans, J. S., Krywonos, A., Lumpe, J. D., Codrescu, M., Veibell, V., McClintock, B. E., and Eastes, R. W.: Global-scale Observations of Limb and Disk (GOLD): Overview of and QEV science data products, in: 2018 AGU Fall Meeting, Washington, DC, SA21A-3169, 2018.
Eastes, R. W.: Emissions from the N2 Lyman–Birge–Hopfield bands in the Earth's atmosphere, Phys. Chem. Earth Pt. C, 25, 523– 527, https://doi.org/10.1016/S1464-1917(00)00069-6, 2000a.
Eastes, R. W.: Modeling the N2 Lyman-Birge-Hopfield bands in the dayglow: Including radiative and collisional cascading between the singlet states, J. Geophys. Res.-Space, 105, 18557–18573, https://doi.org/10.1029/1999JA000378, 2000b.
Eastes, R. W., McClintock, W. E., Codrescu, M. V., Aksnes, A., Anderson, D. N., Andersson, L., Baker, D., Burns, A., Budzien, S. Daniell, R., Dymond, K., Eparvier, F., Harvey, J., Immel, T., Krywonos, A., Lankton, M., Lump, J., Prolss, G., Richmond, A., and Woods, T.: Global-Scale Observations of the Limb and Disk (GOLD): New Observing Capabilities for the Ionosphere-Thermosphere, in: American Geophysical Union Geophysical Monograph Series, American Geophysical Union, Washington, DC, 319–326, https://doi.org/10.1029/181GM29, 2008.
Eastes, R. W., McClintock, W. E., Burns, A. G., Anderson, D. N., Andersson, L., Codrescu, M., Correira, J. T., Daniell, R. E., England, S. L., Evans, J. S., Harvey, J., Krywonos, A., Lumpe, J. D., Richmond, A. D., Rusch, D. W., Siegmund, O., Solomon, S. C., Strickland, D. J., Woods, T. N., Aksnes, A., Budzien, S. A., Dymond, K. F., Eparvier, F. G., Martinis, C. R., and Oberheide, J.: The Global-Scale Observations of the Limb and Disk (GOLD) Mission, Space Sci. Rev., 212, 383–408, https://doi.org/10.1007/s11214-017-0392-2, 2017.
ESA: Swarm Data, available at: https://earth.esa.int/web/guest/swarm/data-access, last access: 25 October 2021.
Fuller-Rowell, T. J., Codrescu, M. V., Moffett, R. J., and Quegan, S.: Response of the thermosphere and ionosphere to geomagnetic storms, J. Geophys. Res., 99, 3893–3914, https://doi.org/10.1029/93JA02015, 1994.
Gan, Q., Eastes, R. W., Burns, A. G., Wang, W., Qian, L., Solomon, S. C., Codrescu, M. V., McInerney, J., and McClintock, W. E.: First synoptic observations of geomagnetic storm effects on the global-scale OI 135.6-nm dayglow in the thermosphere by the GOLD mission, Geophys. Res. Lett., 47, e2019GL085400, https://doi.org/10.1029/2019GL085400, 2020.
Laskar, F. I., Pedatella, N. M., Codrescu, M. V., Eastes, R. W., Evans, J. S., Burns, A. G., and McClintock, W. E.: Impact of GOLD retrieved thermospheric temperatures on a whole atmosphere data assimilation model. J. Geophys. Res.-Space, 126, e2020JA028646, https://doi.org/10.1029/2020JA028646, 2021.
Matsuo T.: Recent Progress on Inverse and Data Assimilation Procedure for High-Latitude Ionospheric Electrodynamics, in: Ionospheric Multi-Spacecraft Analysis Tools, ISSI Scientific Report Series, vol. 17, edited by: Dunlop, M. and Lühr, H., Springer, Cham, https://doi.org/10.1007/978-3-030-26732-2_10, 2020.
McClintock, W. E., Eastes, R. W., Hoskins, A. C., Siegmund, O. H. W., McPhate, J. B., Krywonos, A., Solomon, S. C., and Burns, A. G.: Global-scale Observations of the Limb and Disk mission implementation: 1. Instrument design and early flight performance, J. Geophys. Res.-Space, 125, e2020JA027797, https://doi.org/10.1029/2020JA027797, 2020a.
McClintock, W. E., Eastes, R. W., Beland, S., Bryant, K. B., Burns, A. G., Correira, J., Daniell, R. E., Evans, J. S., Harper, C. S., Karan, D. K., Krywonos, A., Lumpe, J. D., Plummer, T. M., Solomon, S. C., Vanier, B. A., and Veibel, V.: Global-scale Observations of the Limb and Disk mission implementation: 2. Observations, data pipeline, and Level 1 data products, J. Geophys. Res.-Space, 125, e2020JA027809, https://doi.org/10.1029/2020JA027809, 2020b.
Meier, R. R.: Ultraviolet spectroscopy and remote sensing of the upper atmosphere, Space Sci. Rev., 58, 1–185, https://doi.org/10.1007/BF01206000, 1991.
Mende, S. B., Frey, H. U., Rider, K., Chou, C., Harris, S., Siegmund, O., England, S., Wilkins, C., Craig, W., Immel, T., Turin, P., Darling, N., Loicq, J., Blain, P., Syrstad, E., Thompson, B., Burt, R., Champagne, J., Sevilla, P., and Ellis, S.: The Far Ultra-Violet Imager on the Icon Mission, Space Sci. Rev., 212, 655–696, https://doi.org/10.1007/s11214-017-0386-0, 2017.
NASA: GOLD, available at: http://gold.cs.ucf.edu/search/ (last access: 25 October 2021), 2021a.
NASA: NASA's Space Physics Data Facility (SPDF), available at: https://spdf.gsfc.nasa.gov (last access: 25 October 2021), 2021b.
NASA: OMNIWeb Plus, available at: https://omniweb.gsfc.nasa.gov/ (last access: 25 October 2021), 2021c.
NOAA-SWPC: WAM, GitHub [code], https://github.com/NOAA-SWPC/WAM, last access: 25 October 2021
Paxton, L. J., Morrison, D., Kil, H., Zhang, Y., Ogorzalek, B. S., and Meng, C.: Validation of remote sensing products produced by the Special Sensor Ultraviolet Scanning Imager (SSUSI): A far UV imaging spectrograph on DMSP F-16, in: SPIE Optical Spectroscopic Techniques and Instrumentation for Atmospheric and Space Research IV, vol. 4485, edited by: Larar, A. M. and Mlynczak, M. G., 338–348, 2002.
Paxton, L. J., Schaefer, R. K., Zhang, Y., and Kil, H.: Far ultraviolet instrument technology, J. Geophys. Res.-Space, 122, 2706–2733, https://doi.org/10.1002/2016JA023578, 2017.
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, 1468, https://doi.org/10.1029/2002JA009430, 2002.
Rodgers, C. D.: Inverse methods for atmospheric sounding: Theory and practice, World Scientific, Singapore, 2000.
Solomon, S. C.: Global modeling of thermospheric airglow in the far ultraviolet, J. Geophys. Res.-Space, 122, 7834–7848, https://doi.org/10.1002/2017JA024314, 2017.
SuperMAG: https://supermag.jhuapl.edu/, last acces: 25 October 2021.
VT: http://vt.superdarn.org, last access: 25 October 2021.
Zhang, Y., Paxton, L. J., and Schaefer, R. K.: Deriving thermospheric temperature from observations by the global ultraviolet imager on the thermosphere ionosphere mesosphere energetics and dynamics satellite, J. Geophys. Res.-Space, 124, 5848–5856, https://doi.org/10.1029/2018JA026379, 2019.