- Articles & preprints
- Submission
- Policies
- Peer review
- Editorial board
- About
- EGU publications
- Manuscript tracking

Journal cover
Journal topic
**Atmospheric Measurement Techniques**
An interactive open-access journal of the European Geosciences Union

Journal topic

- Articles & preprints
- Submission
- Policies
- Peer review
- Editorial board
- About
- EGU publications
- Manuscript tracking

- Articles & preprints
- Submission
- Policies
- Peer review
- Editorial board
- About
- EGU publications
- Manuscript tracking

- Abstract
- Introduction
- Measurement principle
- Zeeman effect modeling
- Measurement and retrieval setting
- Retrieval errors
- Conclusions
- Appendix A: Spectroscopic parameters
- Appendix B: Matrix exponential
- Code availability
- Author contributions
- Competing interests
- Acknowledgements
- Financial support
- Review statement
- References

**Research article**
20 Jan 2020

**Research article** | 20 Jan 2020

Potential for the measurement of mesosphere and lower thermosphere (MLT) wind, temperature, density and geomagnetic field with Superconducting Submillimeter-Wave Limb-Emission Sounder 2 (SMILES-2)

Potential for the measurement of mesosphere and lower thermosphere (MLT) wind, temperature, density and geomagnetic field with Superconducting Submillimeter-Wave Limb-Emission Sounder 2 (SMILES-2)
Potential for the measurement of mesosphere and lower thermosphere (MLT) wind, temperature,...
Philippe Baron et al.

^{1}National Institute of Information and Communications Technology, Koganei, Tokyo 184-8795, Japan^{2}Center for Global Environmental Research, National Institute for Environmental Studies, Tsukuba, Ibaraki 305-8506, Japan^{3}Planets, Max Planck Institute for Solar System Research, Göttingen, Lower Saxony, Germany^{4}Department of Earth and Planetary Science, Kyushu University, Fukuoka 812-8581, Japan^{5}Center for Environmental Remote Sensing, Chiba University, Chiba-shi, Japan^{6}Department of Space, Earth and Environment, Chalmers University of Technology, 41296 Gothenburg, Sweden^{7}Institute for Space-Earth Environmental Research, Nagoya University, Nagoya Aichi 464-8601, Japan^{8}Ionosphere Research Unit, University of Oulu, Oulu, Finland^{9}National Institute of Polar Research, Tachikawa-shi, Tokyo 190-8518, Japan^{10}Division of Science, Kyoto Sangyo University, Kyoto 603-8555, Japan^{11}Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan^{12}Research Institute for Sustainable Humanosphere, Kyoto University, Kyoto 611-0011, Japan^{13}Japan Aerospace Exploration Agency, Sagamihara, Kanagawa 252-5210, Japan

^{1}National Institute of Information and Communications Technology, Koganei, Tokyo 184-8795, Japan^{2}Center for Global Environmental Research, National Institute for Environmental Studies, Tsukuba, Ibaraki 305-8506, Japan^{3}Planets, Max Planck Institute for Solar System Research, Göttingen, Lower Saxony, Germany^{4}Department of Earth and Planetary Science, Kyushu University, Fukuoka 812-8581, Japan^{5}Center for Environmental Remote Sensing, Chiba University, Chiba-shi, Japan^{6}Department of Space, Earth and Environment, Chalmers University of Technology, 41296 Gothenburg, Sweden^{7}Institute for Space-Earth Environmental Research, Nagoya University, Nagoya Aichi 464-8601, Japan^{8}Ionosphere Research Unit, University of Oulu, Oulu, Finland^{9}National Institute of Polar Research, Tachikawa-shi, Tokyo 190-8518, Japan^{10}Division of Science, Kyoto Sangyo University, Kyoto 603-8555, Japan^{11}Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan^{12}Research Institute for Sustainable Humanosphere, Kyoto University, Kyoto 611-0011, Japan^{13}Japan Aerospace Exploration Agency, Sagamihara, Kanagawa 252-5210, Japan

**Correspondence**: Philippe Baron (baron@nict.go.jp)

**Correspondence**: Philippe Baron (baron@nict.go.jp)

Abstract

Back to toptop
Submillimeter-Wave Limb-Emission Sounder 2 (SMILES-2) is a satellite mission proposed in Japan
to probe the middle and upper atmosphere (20–160 km). The main instrument is composed of 4 K cooled radiometers operating near 0.7 and 2 THz.
It could measure the diurnal changes of the horizontal wind above 30 km, temperature
above 20 km, ground-state atomic oxygen above 90 km and atmospheric density near the mesopause, as well as abundance of about 15 chemical species.
In this study we have conducted simulations to assess the wind, temperature and density retrieval performance in the mesosphere and lower thermosphere (60–110 km) using the radiometer at
760 GHz. It contains lines of
water vapor (H_{2}O), molecular oxygen (O_{2}) and nitric oxide (NO) that are the strongest signals measured with SMILES-2 at these altitudes. The Zeeman effect on the O_{2} line due to the geomagnetic field (** B**) is considered; otherwise, the retrieval errors would be underestimated by a factor of 2 above 90 km. The optimal configuration for the radiometer’s polarization is found to be vertical linear. Considering a retrieval vertical resolution of 2.5 km, the line-of-sight wind is retrieved with a precision of
2–5 m s

Download & links

How to cite

Back to top
top
How to cite.

Baron, P., Ochiai, S., Dupuy, E., Larsson, R., Liu, H., Manago, N., Murtagh, D., Oyama, S., Sagawa, H., Saito, A., Sakazaki, T., Shiotani, M., and Suzuki, M.: Potential for the measurement of mesosphere and lower thermosphere (MLT) wind, temperature, density and geomagnetic field with Superconducting Submillimeter-Wave Limb-Emission Sounder 2 (SMILES-2), Atmos. Meas. Tech., 13, 219–237, https://doi.org/10.5194/amt-13-219-2020, 2020.

1 Introduction

Back to toptop
The mesosphere and lower thermosphere (MLT) is a transitional region (60–110 km) between
atmospheric layers with very different characteristics, namely the
stratosphere (15–60 km) and the thermosphere (90–400 km) (Smith, 2012; Shiotani et al., 2019).
In the stratosphere, O_{3} controls the chemical and radiative processes; hence, it also regulates the temperature and the dynamics.
In the thermosphere, the chemistry and the radiative balance
are mainly controlled by the oxygen atoms. In this region, wind and temperature exhibit large
diurnal variations and are strongly influenced by tides generated in the lower
atmosphere. The thermosphere is also the region of interactions between the ionized (plasma) and neutral atmosphere.

The mean physical characteristics of the MLT (wind, temperature and density)
are primarily established by energy transferred from the troposphere via small-scale gravity waves (GWs) (Fritts and Alexander, 2003; Tsuda, 2014). Hence, the MLT state deviates significantly from the
radiative equilibrium, as illustrated by the occurrence of the coldest point of the Earth system (≈150 K) in the summer polar
mesopause.
Waves with planetary scales also contribute to the upper atmosphere climate (general circulation) through their momentum and energy
transport/deposition (Forbes et al., 2006; Pancheva and Mukhtarov, 2011). In particular, tides that are
mainly driven by diurnally varying diabatic heating in the troposphere and the
stratosphere propagate upward, with their amplitude reaching a maximum in the
MLT (Chapman and Lindzen, 1970; Sakazaki et al., 2015).
Hence, the MLT plays a key role in connecting the lower and upper atmosphere and also in
linking both hemispheres (Xu et al., 2009; Karlsson and Becker, 2016). Furthermore, the increase in anthropogenic CO_{2} is responsible for
a cooling of 1–3 K decade^{−1} in the MLT that has been measured since the early 1990s (Beig, 2011).

The processes behind these phenomena are still not well quantified. The difficulty arises from the nonlinear interactions between the GWs, tides, planetary waves, the background wind and the electromagnetic field (Sato et al., 2018; Immel et al., 2006). The system is further complicated by the interconnections between the dynamics and highly variable chemical species, as well as the very different temporal and spatial scales of these processes. Observations of the MLT, in particular of wind, temperature and density, are therefore essential to further our understanding of this region (Smith, 2012).

Continuous measurements of temperature and wind are performed from ground-based stations using lidars (Steinbrecht et al., 2009; Baumgarten, 2010), radars (Jacobi et al., 2015; Tsutsumi et al., 2017) and, up to 70 km, millimeter radiometers (Rüfenacht et al., 2014). Density was recently monitored using meteor radars (Yi et al., 2019), but measurements remain scarce. Satellite observations of the MLT have also been performed for several decades. The missions currently in operation and capable of measuring at these altitudes are listed in Table 1. Temperature is measured with various techniques and spectral domains (Schwartz et al., 2006; Sica et al., 2008; Sheese et al., 2010; Christensen et al., 2015; Eastes et al., 2017; Englert et al., 2017), but discrepancies larger than 10 K can be found between these measurements above 80 km (García-Comas et al., 2014). Baron et al. (2013) and Shepherd (2015) described the past and current wind measurements from space. Currently only TIDI and MLS (and soon MIGHTI) are capable of measuring MLT winds but with a poor sensitivity below 80 km (Niciejewski et al., 2006; Wu et al., 2008; Englert et al., 2017), and MLS, which is equipped with a single antenna, can only measure one component of the wind vector (it was not designed for wind measurement).

In the future, we clearly risk a lack of satellite observations since all the current
missions (except ICON) have already exceeded their theoretical lifetime. Sweden is preparing two Innosat-based missions that are of interest for the study of the MLT (Table 1).
The MATS mission aims at characterizing
the 3-D structure of the GWs near 90–100 km using the oxygen A-band
emission and the ultraviolet light scattered by noctilucent clouds (Gumbel et al., 2018).
Information on temperature will also be retrieved. The other
mission is SIW, a sub-millimeter limb
sounder that will measure horizontal wind, temperature and trace gases up to about
80 km (Baron et al., 2018). The MATS and SIW missions will be operational for 2 years between 2019 and 2021 and between 2023 and 2025, respectively. Other projects have not been selected yet and remain uncertain.
For example, Wu et al. (2016) proposed a THz limb sounder (TLS) to measure the atomic oxygen
line at 2 THz. Such an instrument could fly together with a new version of
SABER (Mlynczak and Yee, 2017). The European Space Agency (ESA) is studying a limb sounder
operating between 0.8 and 4 THz for the retrieval of the abundance of chemical
species such as atomic oxygen (O) or the hydroxyl radical (OH) (Gerber et al., 2013). TALIS, a limb sounder using similar spectral bands to Aura MLS, is being
studied in China (Wang et al., 2019). Kaufmann et al. (2018) described a concept for a limb
sounder on board a CubeSat to measure temperature with high horizontal resolution using the
molecular oxygen (O_{2}) A-band infra-red emission.

Superconducting Submillimeter-Wave Limb-Emission Sounder 2 (SMILES-2) is a middle and upper atmospheric satellite mission proposed to the Japan Aerospace Exploration Agency (JAXA) (Ochiai et al., 2017, 2019; Shiotani et al., 2019). If selected, it will be launched around 2026 on a JAXA M-class satellite. The objectives are to provide geophysical information with unprecedented precision and altitude coverage such as the temperature between 15 and 160 km, horizontal wind between 30 and 160 km, atmospheric density up to 110 km, ground state of atomic oxygen between 90 and 160 km, and more than 15 trace gases' abundance (Baron et al., 2019a, b). The proposed satellite will be equipped with two antennas for the limb measurement of horizontal winds, and three radiometers near 0.7 and 2 THz cooled at 4 K, a technology successfully tested with JEM/SMILES (Kikuchi et al., 2010). With a precessing orbit and the high receiver precision, it will be possible to retrieve diurnal variations of very weak signals, as demonstrated with JEM/SMILES (Sakazaki et al., 2013; Khosravi et al., 2013).

In this study we discuss the potential for SMILES-2 to measure the main characteristics of the neutral MLT, namely wind, temperature and atmospheric density.
An essential source of information is the O_{2} transition at 773.8 GHz.
As a magnetic dipole, O_{2} is subject to the Zeeman effect induced by the Earth's magnetic field (** B**). Special care is taken to properly include this effect in the simulations in order to correctly assess the measurement performance. Retrieval errors induced by uncertainties on

2 Measurement principle

Back to toptop
The observation characteristics are summarized in Table 2. The atmospheric limb is scanned from about 20 to 180 km. Scans are performed alternatively with two antennas looking at perpendicular directions to each other.
Both antennas can probe the same atmospheric column with a 7 min delay (Fig. 1), allowing us to derive the 2-D
horizontal winds. The same method will be used for SIW and more information is given in Baron et al. (2018). The limb geometry provides a high vertical resolution of 2–3 km, and the zonal and meridional samplings at the Equator are about 20^{∘} (2200 km) and 6^{∘} (650 km), respectively.

The orbit precesses with a period of about 3 months.
The satellite orientation is reversed after every half precession cycle in order to keep the solar panels properly illuminated and the radiative-cooling panels in the shadow side. The latitude coverage is between 50^{∘} S and 80^{∘} N or 80^{∘} S and 50^{∘} N depending on the satellite orientation.
At low and middle latitudes, the same latitude is observed twice per orbit, with LT differences close to 12 h. Hence, gathering the observations between each maneuver allows us to piece together the complete diurnal cycle of the retrieved parameters.

Three spectral bands near 638, 763 GHz and 2 THz are measured simultaneously (Ochiai et al., 2018). The band at 638 GHz contains a strong stratospheric and lower mesospheric
signal from ozone (O_{3}). This band is the same as that selected for
SIW and its main characteristics are described in Baron et al. (2018). Two THz bands are measured alternatively: one contains OH lines and the
second one an O line (Ochiai et al., 2017; Baron et al., 2019a). The O line is used to retrieve between 90 and 160 km, the abundance of O in its ground state, wind and temperature (Baron et al., 2015, 2019b).

^{∗} Estimated for a tangent height of 80 km including the antenna field of view and the scan velocity.

The 763 GHz band (Table 3) is the band considered in this study. It contains lines of water vapor
(H_{2}O) at 752.03 GHz and O_{2} at 773.84 GHz (Fig. 2) that provide a strong signal in the MLT.
It also contains other molecular lines, weaker but still suitable for our study: nitric oxide (NO, 751.67–752.00 and 773.02–773.05 GHz), O_{3} (754.46 and 776.66 GHz) and carbon monoxide isotopologue (^{13}CO) at 771.183 GHz. The bands have changed compared to those originally described by Ochiai et al. (2017), a change motivated to reduce the power consumption. In the new setting, the CO line is about 50 times weaker than that previously selected.

Most of the lines in the spectral bands are emitted by chemical species
in their ground state under local thermodynamic equilibrium.
The molecular abundance and the temperature are retrieved from the amplitude of the
lines. Their Doppler shift (2.5 kHz for 1 m s^{−1}) is used to retrieve the line-of-sight (LOS) wind.
The atmospheric
density is derived from the O_{2} abundance considering that the volume mixing
ratio of O_{2} is well known below 110 km (Schwartz et al., 2006).

Above about 70 km, the lines are broadened by the random molecular motions, i.e., Doppler broadening, and they do not carry direct information on the pressure (Appendix A). Consequently, the density of the molecule can be retrieved and not the volume-mixing ratio (VMR) as in the lower altitudes.

Molecular oxygen is a magnetic dipole that interacts with ** B**. It is subject to the Zeeman effect (Lenoir, 1968) and the selected spectroscopic transition is split into

In this study, we consider LOS tangent heights between 60 and 110 km. They are provided as input for the inversion algorithm; therefore they must be known before inverting the spectra. Height registration for a complete scan is calculated differently in the lower part of the scan and in the range of interest (between ∼20 and 60 and between 60 and 110 km, respectively).

Between 20 and 60 km, an approach similar to that used for Aura MLS (Schwartz et al., 2006) can be used. The
LOS tangent pressure and atmospheric temperature would be retrieved simultaneously from the O_{2} line near 763 GHz and from O_{3} lines in the 638 GHz band.
The height of the pressure levels would then be derived from the hydrostatic equilibrium equation.
The resulting precisions are estimated to be better than 1 % and 75 m for the LOS tangent pressure and height, respectively (Baron et al., 2019b).

In the altitude range of interest (>60 km), the LOS tangent heights are inferred from the extrapolation of those calculated previously for the lower altitudes and attitude data from the star-trackers and GPS on board the satellite. Based on JEM/SMILES results, the expected precision on the retrieved LOS tangent heights will be 100 m or better (Ochiai et al., 2013).

3 Zeeman effect modeling

Back to toptop
The Zeeman effect on atmospheric molecular-oxygen lines has been extensively studied (Lenoir, 1968; Pardo et al., 1995; Schwartz et al., 2006; Larsson et al., 2014; Navas-Guzmán et al., 2015). In this study, we describe the polarized radiance with Stokes vectors as in Landi Degl'Innocenti and Landolfi (2004) (e.g., Eq. 1.32), Larsson et al. (2014) and Steiner et al. (2016). The magnetic field characteristics (amplitude and orientation angles with respect to the LOS) are defined at the LOS tangent height (Fig. 4) and are assumed to be constant over the LOS. This approximation is the same as that used by Yee et al. (2017), and it is justified since most of the retrieved information comes from a thin altitude range around the tangent point.

The interaction between the radiation and the atmosphere are described by the 4×4 absorption matrix **K**:

$$\begin{array}{}\text{(1)}& \mathbf{K}={k}_{\mathrm{a}}\mathbf{I}+{\mathbf{K}}_{\mathrm{o}},\end{array}$$

where **I** is the identity matrix, *k*_{a} is the scalar absorption coefficient and **K**_{o} is a matrix with off-diagonal components:

$$\begin{array}{}\text{(2)}& {\mathbf{K}}_{\mathrm{o}}=\left[\begin{array}{cccc}\mathrm{0}& q& u& v\\ q& \mathrm{0}& {v}^{\prime}& -{u}^{\prime}\\ u& -{v}^{\prime}& \mathrm{0}& {q}^{\prime}\\ v& {u}^{\prime}& -{q}^{\prime}& \mathrm{0}\end{array}\right].\end{array}$$

The scalar absorption coefficient is computed using a line-by-line model and the Zeeman effect is only applied on the O_{2} transition:

$$\begin{array}{}\text{(3)}& \begin{array}{rl}& {k}_{\mathrm{a}}(\mathit{\nu},z)=\\ & \sum _{M,t}\phantom{\rule{0.125em}{0ex}}{n}_{M}\left(z\right)\phantom{\rule{0.125em}{0ex}}{S}_{t}\left(z\right)\phantom{\rule{0.125em}{0ex}}F(\mathit{\nu},{\mathit{\nu}}_{t},{\mathrm{\Gamma}}_{M,z})\\ & +{\displaystyle \frac{{n}_{{\mathrm{O}}_{\mathrm{2}}}\left(z\right)\phantom{\rule{0.125em}{0ex}}{S}_{x}\left(z\right)}{\mathrm{2}}}\phantom{\rule{0.125em}{0ex}}\left({\mathrm{sin}}^{\mathrm{2}}\right(\mathit{\theta})\sum _{\mathit{\pi}}\left[{s}_{\mathit{\pi}}\phantom{\rule{0.125em}{0ex}}F(\mathit{\nu},{\mathit{\nu}}_{\mathit{\pi}},{\mathrm{\Gamma}}_{{\mathrm{O}}_{\mathrm{2}},z})\phantom{\rule{0.125em}{0ex}}\right]\\ & +\left(\mathrm{1}+{\mathrm{cos}}^{\mathrm{2}}\left(\mathit{\theta}\right)\right)\sum _{{\mathit{\sigma}}^{+},{\mathit{\sigma}}^{-}}\left[{\displaystyle \frac{{s}_{\mathit{\sigma}}}{\mathrm{2}}}\phantom{\rule{0.125em}{0ex}}F(\mathit{\nu},{\mathit{\nu}}_{\mathit{\sigma}},{\mathrm{\Gamma}}_{{\mathrm{O}}_{\mathrm{2}},z})\phantom{\rule{0.125em}{0ex}}\right]),\end{array}\end{array}$$

where *ν* is the frequency, *z* is the altitude, *t* denotes a spectroscopic transition of the species *M* that is not affected by the geomagnetic field, *n*_{M} (${n}_{{\mathrm{O}}_{\mathrm{2}}}$) is the number density of *M* (O_{2}), *S*_{t} is the line strength, *F* is the Voigt function (Schreier et al., 2014; Larsson et al., 2014) and Γ_{t,z} represents the parameters related to the linewidth (Appendix A). The angle *θ* is the inclination angle of the magnetic field with respect to the LOS (Fig. 4a).

The frequencies *ν*_{σ,π} (Hz) are those of the Zeeman components (Fig. 3). They are dependent on the magnetic field (Larsson et al., 2014):

$$\begin{array}{}\text{(4)}& {\mathit{\nu}}_{\mathit{\sigma},\mathit{\pi}}=-{\displaystyle \frac{{\mathit{\mu}}_{\mathrm{b}}}{{h}_{\mathrm{p}}}}\phantom{\rule{0.125em}{0ex}}\left|\mathit{B}\right|{g}_{\mathrm{s}}{\mathit{\beta}}_{\mathrm{m}}=\mathrm{2.80209}\times {\mathrm{10}}^{\mathrm{10}}\phantom{\rule{0.125em}{0ex}}\left|\mathit{B}\right|{\mathit{\beta}}_{\mathrm{m}},\end{array}$$

where *g*_{s}=2.002064, *μ*_{b} is the Bohr magneton ($\mathrm{9.27401}\times {\mathrm{10}}^{-\mathrm{24}}$ J T^{−1}), *h*_{p} is the Planck constant ($\mathrm{6.62618}\times {\mathrm{10}}^{-\mathrm{34}}$ m^{2} kg s^{−1}) and

$$\begin{array}{}\text{(5)}& {\displaystyle}\begin{array}{rl}& {\mathit{\beta}}_{\mathrm{m}}=\left({\displaystyle \frac{{J}_{u}({J}_{u}+\mathrm{1})+S(S+\mathrm{1})-{N}_{u}({N}_{u}+\mathrm{1})}{\mathrm{2}({J}_{u}+\mathrm{1}){J}_{u}}}\right.\\ & \left.{m}_{u}-{\displaystyle \frac{{J}_{l}({J}_{l}+\mathrm{1})+S(S+\mathrm{1})-{N}_{l}({N}_{l}+\mathrm{1})}{\mathrm{2}({J}_{l}+\mathrm{1}){J}_{l}}}\phantom{\rule{0.125em}{0ex}}{m}_{l}\right).\end{array}\end{array}$$

Here, the lower scripts *u* and *l* denote the upper and lower levels of the transition, respectively, and *N*, *J*, *S* and *m* are quantum numbers associated with the angular momentum, the spin, the total momentum *N*+*S* and the projection of *J* on the *B* axis.

The coefficients of **K**_{o} are derived from Landi Degl'Innocenti and Landolfi (2004) (Eq. 5.36):

$$\begin{array}{}\text{(6)}& \begin{array}{rl}& q={\displaystyle \frac{{\mathrm{sin}}^{\mathrm{2}}\left(\mathit{\theta}\right)\mathrm{cos}\left(\mathrm{2}\mathit{\varphi}\right)}{\mathrm{2}}}\\ & (\sum _{\mathit{\pi}}\left[{s}_{\mathit{\pi}}\phantom{\rule{0.125em}{0ex}}F(\mathit{\nu},{\mathit{\nu}}_{\mathit{\pi}},{\mathrm{\Gamma}}_{x,z})\phantom{\rule{0.125em}{0ex}}\right]-\sum _{{\mathit{\sigma}}^{\pm}}\left[{\displaystyle \frac{{s}_{\mathit{\sigma}}}{\mathrm{2}}}\phantom{\rule{0.125em}{0ex}}F(\mathit{\nu},{\mathit{\nu}}_{\mathit{\sigma}},{\mathrm{\Gamma}}_{{\mathrm{O}}_{\mathrm{2}},z})\phantom{\rule{0.125em}{0ex}}\right]),\\ & u={\displaystyle \frac{{\mathrm{sin}}^{\mathrm{2}}\left(\mathit{\theta}\right)\mathrm{sin}\left(\mathrm{2}\mathit{\varphi}\right)}{\mathrm{2}}}\\ & (\sum _{\mathit{\pi}}\left[{s}_{\mathit{\pi}}\phantom{\rule{0.125em}{0ex}}F(\mathit{\nu},{\mathit{\nu}}_{\mathit{\pi}},{\mathrm{\Gamma}}_{x,z})\phantom{\rule{0.125em}{0ex}}\right]-\sum _{{\mathit{\sigma}}^{\pm}}\left[{\displaystyle \frac{{s}_{\mathit{\sigma}}}{\mathrm{2}}}\phantom{\rule{0.125em}{0ex}}F(\mathit{\nu},{\mathit{\nu}}_{\mathit{\sigma}},{\mathrm{\Gamma}}_{{\mathrm{O}}_{\mathrm{2}},z})\phantom{\rule{0.125em}{0ex}}\right]),\\ & v=\mathrm{cos}\left(\mathit{\theta}\right)\\ & (\sum _{{\mathit{\sigma}}^{\pm}}\pm \phantom{\rule{0.125em}{0ex}}{\displaystyle \frac{{s}_{{\mathit{\sigma}}^{\pm}}}{\mathrm{2}}}\phantom{\rule{0.125em}{0ex}}F(\mathit{\nu},{\mathit{\nu}}_{{\mathit{\sigma}}^{\pm}},{\mathrm{\Gamma}}_{{\mathrm{O}}_{\mathrm{2}},z}\left)\right).\end{array}\end{array}$$

The parameters *u*^{′}, *v*^{′} and *q*^{′} are computed by replacing the term *F* with *F*^{′},
the dispersive part of the complex Voigt function (see Appendix A and Schreier et al., 2014).

The LOS is divided in narrow ranges of size d*s* (typically 5 km long) in which the atmospheric parameters are
considered constant. The change of the polarized radiance passing through an homogeneous range is derived from a matrix equation which is similar
to the scalar radiative transfer one used for a nonpolarized radiation (Semel and López, 1999):

$$\begin{array}{}\text{(7)}& {\mathit{b}}_{\mathrm{a}}(s+\mathrm{d}s)=(\mathbf{I}-\mathbf{\Lambda}(s,s+\mathrm{d}s\left)\right)\cdot {\mathit{b}}_{\mathrm{p}}\left(s\right)+\mathbf{\Lambda}(s,s+\mathrm{d}s)\cdot {\mathit{b}}_{\mathrm{a}}\left(s\right),\end{array}$$

where *b*_{a}(*s*) is the Stokes vector at the position *s* on the LOS (the frequency dependence is omitted), “⋅” is the matrix multiplication operator, ${\mathit{b}}_{\mathrm{p}}\left(s\right)=\left[P\right(s),\mathrm{0},\mathrm{0},\mathrm{0}{]}^{T}$ describes the nonpolarized source function between *s* and *s*+d*s*, *P*(*s*) is the Planck function, and $\mathbf{\Lambda}(s,s+\mathrm{d}s)$ is 4×4 evolution operator matrix defined as follows:

$$\begin{array}{}\text{(8)}& \mathbf{\Lambda}(s,s+\mathrm{d}s)=\mathrm{exp}(-\mathbf{K}(s\left)\phantom{\rule{0.125em}{0ex}}\mathrm{d}s\right).\end{array}$$

The integration over the LOS is performed by applying the scalar equation given by Urban et al. (2004) to Stokes parameters:

$$\begin{array}{}\text{(9)}& \begin{array}{rl}& {\mathit{b}}_{\mathrm{a}}\left(\mathrm{sat}\right)=\sum _{i=\mathrm{0}}^{N-\mathrm{1}}\mathbf{\Lambda}(i+\mathrm{1},\mathrm{sat})\cdot (\mathbf{I}-\mathbf{\Lambda}(\mathrm{0},i+\mathrm{1})\cdot \mathbf{\Lambda}(\mathrm{0},i+\mathrm{1}\left)\right)\\ & \cdot \left({\mathit{b}}_{\mathrm{p}}\right(i)-{\mathit{b}}_{\mathrm{p}}(i+\mathrm{1}\left)\right)+(\mathbf{I}-\mathbf{\Lambda}(\mathrm{0},\mathrm{sat})\cdot \mathbf{\Lambda}(\mathrm{0},\mathrm{sat}\left)\right)\\ & \cdot {\mathit{b}}_{\mathrm{p}}\left(N\right),\end{array}\end{array}$$

where *b*_{a}(sat) is the Stokes vector representing the radiation state at the antenna position,
*i* is the index of the level at *s*_{i} (*i*=0 for the tangent point) and *N* is the number of levels above the tangent point. The cosmic background radiation is neglected.
We use the relationship $\mathbf{\Lambda}(i,j)=\mathbf{\Lambda}(k,j)\cdot \mathbf{\Lambda}(i,k)$
with $i<k<j$ (the two matrices on the right-hand side of the equality do not commute).

4 Measurement and retrieval setting

Back to toptop
The measured radiance for antenna *a* (*a*=1 or 2) at the elevation angle *θ* and the IF *ν* is as follows:

$$\begin{array}{}\text{(10)}& {y}_{\mathit{\theta},\mathit{\nu}}^{a}={\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}\left({\mathit{R}}_{\mathit{\theta},\mathit{\nu}}^{\text{u}}\phantom{\rule{0.25em}{0ex}}\textcolor{red}{}\phantom{\rule{0.25em}{0ex}}{\mathit{y}}^{\text{a,u}}({\mathit{\nu}}_{\text{LO}}+\mathit{\nu})+{\mathit{R}}_{\mathit{\theta},\mathit{\nu}}^{\text{l}}\phantom{\rule{0.25em}{0ex}}\textcolor{red}{}\phantom{\rule{0.25em}{0ex}}{\mathit{y}}^{\text{a,l}}({\mathit{\nu}}_{\text{LO}}-\mathit{\nu})\right),\end{array}$$

where *y*^{a,u} and *y*^{a,l} are the atmospheric specific intensities in the upper and lower sidebands around the local oscillator frequency *ν*_{LO}, *R*_{θ,ν} represents the antenna and spectrometer functions, and is the convolution operator (Baron et al., 2018). A simple case with a constant upper and lower sideband ratio is considered.
The Zeeman model is only used within a
bandwidth of 200 MHz encompassing the O_{2} line (upper sideband). Outside this range, the nonpolarized radiative transfer model described in Baron et al. (2018) is used.
In order to transform the Stokes vector (Eq. 9) to the specific intensity associated with the radiometer's polarization, we first rotate the vector from the atmospheric frame to the detector frame as follows:

$$\begin{array}{}\text{(11)}& {\mathit{b}}_{\mathrm{d}}={\mathrm{M}}_{\mathrm{r}}\left({\mathit{\alpha}}_{\mathrm{d}}\right)\cdot {\mathit{b}}_{\mathrm{a}},\end{array}$$

where *b*_{d} is the Stokes vector in the instrument frame and M_{r}(*α*_{d}) is the Mueller matrix for a rotation *α*_{d}:

$$\begin{array}{}\text{(12)}& {\mathrm{M}}_{\mathrm{r}}\left({\mathit{\alpha}}_{\mathrm{d}}\right)=\left[\begin{array}{cccc}\mathrm{1}& \mathrm{0}& \mathrm{0}& \mathrm{0}\\ \mathrm{0}& \mathrm{cos}\left(\mathrm{2}\phantom{\rule{0.125em}{0ex}}{\mathit{\alpha}}_{\mathrm{d}}\right)& \mathrm{sin}\left(\mathrm{2}\phantom{\rule{0.125em}{0ex}}{\mathit{\alpha}}_{\mathrm{d}}\right)& \mathrm{0}\\ \mathrm{0}& -\mathrm{sin}\left(\mathrm{2}\phantom{\rule{0.125em}{0ex}}{\mathit{\alpha}}_{\mathrm{d}}\right)& \mathrm{cos}\left(\mathrm{2}\phantom{\rule{0.125em}{0ex}}{\mathit{\alpha}}_{\mathrm{d}}\right)& \mathrm{0}\\ \mathrm{0}& \mathrm{0}& \mathrm{0}& \mathrm{1}\end{array}\right].\end{array}$$

The specific intensity *y* corresponding to the detector polarization is

$$\begin{array}{}\text{(13)}& y={\mathit{b}}_{\mathrm{d}}\left[\mathrm{1}\right]+m\phantom{\rule{0.25em}{0ex}}{\mathit{b}}_{\mathrm{d}}\left[n\right],\end{array}$$

where *b*_{d}[*n*] is the *n*th component of the Stokes vector and (*m*,*n*) is
$(-\mathrm{1},\mathrm{1})$, (1,1), $(-\mathrm{1},\mathrm{2})$, (1,2), $(-\mathrm{1},\mathrm{3})$ and (1,3) for horizontal,
vertical, +45, $-\mathrm{45}{}^{\circ}$, right and left circular
polarizations, respectively.

Figure 5 shows simulated spectra of the O_{2} line over the
Equator and at 80^{∘} N when the satellite is moving toward north
(ascending orbit branch). The tangent height is 100 km and the atmospheric
conditions are representative of the Northern Hemisphere in
wintertime (Baron et al., 2018). The magnetic field characteristics are zonal means inferred from a quiet solar day (Fig. 6).
Spectra are shown for different radiometer's polarizations.
Over the Equator, ** B** is along
the meridional direction and clear differences are seen between the
radiances measured with both antennas, except if the detector has a
vertical polarization. In that case, the radiometer detects only the

Over the polar region, the spectra measured by both antennas
are very similar since the vector ** B** is almost vertical and
perpendicular to both LOSs (Fig. 6). Only the Zeeman components

The geomagnetic field may exhibit rapid temporal and spatial variations that can be as large as hundreds of nanoteslas (Doumbia et al., 2007; Yee et al., 2017). Such variations will be difficult to take into account when processing the data and may lead to retrieval errors with the same magnitude as those induced by the measurement noise.

Such errors are mitigated by retrieving the three components of ** B** simultaneously with other atmospheric parameters.
It is done by using the scans of the same atmospheric column measured with the two antennas (Fig. 4).
The measurement vector

$$\begin{array}{}\text{(14)}& {\mathit{y}}^{T}=\left[{\mathit{y}}^{{a}_{\mathrm{1}}},\phantom{\rule{0.25em}{0ex}}{\mathit{y}}^{{a}_{\mathrm{2}}}\right],\end{array}$$

where the superscripts *a*_{1} and *a*_{2} denote that the parameters are associated with the antennas 1 and 2, respectively.
The vector ** x** describing the retrieved parameters contains the profiles of the chemical species having the most significant features in the MLT spectra, namely O

$$\begin{array}{}\text{(15)}& \begin{array}{rl}{\mathit{x}}^{T}=& \left[{\mathit{x}}_{{\mathrm{O}}_{\mathrm{2}}}^{{a}_{\mathrm{1}}},\mathrm{\dots},{\mathit{x}}_{\mathrm{T}}^{{a}_{\mathrm{1}}},{\mathit{x}}_{\mathrm{LW}}^{{a}_{\mathrm{1}}},\right.\\ & {\mathit{x}}_{{\mathrm{O}}_{\mathrm{2}}}^{{a}_{\mathrm{2}}},\mathrm{\dots},{\mathit{x}}_{\mathrm{T}}^{{a}_{\mathrm{2}}},{\mathit{x}}_{\mathrm{LW}}^{{a}_{\mathrm{2}}},\\ & \left.{\mathit{x}}_{\mathrm{Bw}},{\mathit{x}}_{\mathrm{Bu}},{\mathit{x}}_{\mathrm{Bv}}\right],\end{array}\end{array}$$

where *x*_{Bw},
*x*_{Bu} and *x*_{Bv} are the profiles of the vertical, zonal and meridional components of ** B**. The abundance and temperature profiles are retrieved for each antenna in order to account for differences between both scan locations. This is a similar approach to that used by Hagen et al. (2018) for the measurement of winds with the ground-based radiometer WIRA.

The retrieval error induced by the measurement noise is (Rodgers, 2000)

$$\begin{array}{}\text{(16)}& {\mathit{\u03f5}}_{n}^{\mathrm{2}}=\mathrm{diag}\left\{{\left({\mathbf{K}}^{T}{\mathbf{S}}_{y}^{-\mathrm{1}}\mathbf{K}+{\mathbf{U}}^{-\mathrm{1}}\right)}^{-\mathrm{1}}\right\},\end{array}$$

where $\mathbf{K}=\frac{\mathrm{d}\mathit{y}}{\mathrm{d}\mathit{x}}$ is the Jacobian matrix of the retrieved parameters ** x** and

$$\begin{array}{}\text{(17)}& {\mathbf{S}}_{y,i,i}={\displaystyle \frac{{\left({T}_{\text{sys}}+{y}_{i}\right)}^{\mathrm{2}}}{\mathit{\delta}\mathit{\nu}\phantom{\rule{0.125em}{0ex}}\mathit{\delta}t}},\end{array}$$

and ${\mathbf{S}}_{y,i,i}$ is the noise induced variance on the *i*th component of the measurement vector ** y**,

The radiative transfer model computes the Jacobian ${\mathbf{K}}_{B}=\partial {\mathit{y}}^{{a}_{i}}/\partial {\mathit{x}}_{B}$ with respect to antenna-*i* frame
($\mathit{\{}{\mathit{x}}_{\mathit{i}},{\mathit{y}}_{\mathit{i}},{\mathit{z}}_{\mathit{i}}\mathit{\}}$ in Fig. 4a). The matrix **K**_{B} is then computed in the atmospheric frame (Fig. 4):

$$\begin{array}{}\text{(18)}& {\displaystyle \frac{\partial {\mathit{y}}^{{a}_{i}}}{\partial {\mathit{B}}_{q}}}=\sum _{k=\mathit{\{}{x}_{i},{y}_{i},{z}_{i}\mathit{\}}}{\displaystyle \frac{\partial {\mathit{y}}^{{a}_{i}}}{\partial {\mathit{B}}_{k}}}\phantom{\rule{0.25em}{0ex}}{\displaystyle \frac{\partial {\mathit{B}}_{k}}{\partial {\mathit{B}}_{q}}},\end{array}$$

where $q=\mathit{\{}\mathit{u},\mathit{v},\mathit{w}\mathit{\}}$ denote the atmospheric frame axes, and

$$\begin{array}{}\text{(19)}& \begin{array}{rl}& {\mathit{B}}_{{x}_{i}}={\mathit{B}}_{w},\\ & {\mathit{B}}_{{y}_{i}}=\mathrm{cos}\left({\mathrm{\Phi}}_{i}\right)\phantom{\rule{0.125em}{0ex}}{\mathit{B}}_{u}+\mathrm{sin}\left({\mathrm{\Phi}}_{i}\right)\phantom{\rule{0.125em}{0ex}}{\mathit{B}}_{v},\\ & {\mathit{B}}_{{z}_{i}}=-\mathrm{sin}\left({\mathrm{\Phi}}_{i}\right)\phantom{\rule{0.125em}{0ex}}{\mathit{B}}_{u}+\mathrm{cos}\left({\mathrm{\Phi}}_{i}\right)\phantom{\rule{0.125em}{0ex}}{\mathit{B}}_{v},\end{array}\end{array}$$

where Φ_{i} is the angle between the antenna-*i* LOS and the meridional direction (Fig. 4a).

5 Retrieval errors

Back to toptop
Figure 8 shows the retrieval errors on the
atmospheric density, temperature, LOS wind and the main chemical species at three latitudes (50^{∘} S,
Equator, 80^{∘} N).
For the instrumental setting, we considered a radiometer with a linear vertical polarization and the forward-looking antenna (antenna-1). The vertical resolution
of the retrieved profiles is 2.5 km for the main parameters (temperature, LOS wind, H_{2}O and O_{3}), 5 km for NO and 20 km for the components of ** B**.
Errors are computed for the same winter (DJF) climatology described in the previous section. The corresponding atmospheric state includes a stable polar vortex and does not show any NO enhancement due to energetic particle
precipitation. The results at the Equator and the Southern Hemisphere (SH)
mid-latitudes (50

The results for the full band are compared with those computed for the inversion of a 200 MHz band containing only the O_{2} line.
The purpose is to isolate and characterize the contribution of the H_{2}O, O_{3} and NO spectral lines to the retrieval of MLT parameters, in terms of altitude range and impact on the retrieval errors.
Latitudinal differences are induced by the mean meridional circulation (from the summer pole to the winter pole). In the winter
hemisphere, it is
responsible for an increase in NO and a decrease in H_{2}O,
especially over the polar region.
The largest sensitivity to NO is found in the
upper part of the MLT. The precision is better than 10 % above 95 km at 50^{∘} S and above 78 km over the winter polar region (NH in this study).
A precision of 10 % or better is achieved above 95 km at 50^{∘} S and 78 km in the winter polar region.
The sensitivity to H_{2}O decreases with increasing altitude, more sharply
above 90 km.
The precision is better than 1 % up to 75 km in the SH and 65 km in the NH polar region.

The relative error in O_{3} retrieval is ∼1 % around 60 km and strongly
increases with increasing altitude and outside of the polar night, because of the
daytime photo-dissociation of O_{3}.

The achieved precision of the atmospheric density (or O_{2}) profile is better than 5 % up to about 95 km at all latitudes.
Above 90 km, the signal intensity drops significantly and errors quickly increase, up to 20 % at 110 km.
Outside of the 70–90 km range, there are significant differences between the
error profiles calculated for the full- and narrow-band inversions. This shows that
spectral lines from other molecular species also have an impact on the O_{2}
retrievals.
This impact probably occurs through the temperature retrieval. For instance, over the winter polar region, the strong NO signal significantly improves the temperature retrievals and thus indirectly improves the O_{2} abundance retrieval.
Similarly, including H_{2}O and O_{3} lines leads to an improvement of the O_{2} retrieval quality below 70 km.

For all latitudes, the temperature retrieval error is better than 5 K below 90 km and 30 K at 110 km.
The O_{2} line is the main source of information on the temperature near 90 km.

The LOS wind, a key product for SMILES-2, is retrieved with a precision of 2–4 m s^{−1} up to 90 km. Above this altitude, the retrieval errors strongly increase, up to 20 m s^{−1} or more at 110 km.
The O_{2} line is the main source of information on the LOS wind above 70 km. Over the polar region and above 100 km, spectral lines of NO contribute
significantly to the LOS wind retrievals.

Figures 9 and 10 show the achieved retrieval precisions for temperature and LOS wind, at altitudes between 80 and 110 km and for different polarization settings.
Results are shown within the latitude range 50^{∘} S–80^{∘} N, for both antennas and for both the ascending and descending orbit branches. The results obtained with ** B**=0 are also presented.

For atmospheric temperature and below 90 km, the Zeeman effect has a negligible impact on the retrieval errors. Differences can be seen only at high latitudes, where the decrease in the H_{2}O abundance explains the larger impact of the O_{2} line on the retrieval. In terms of LOS wind retrieval, the Zeeman effect is negligible below 80 km. Above 90 km, the approximation ** B**=0 leads to a significant underestimation of the retrieval errors, with differences of up to a factor of 2. This clearly shows that the retrieval errors depend on the radiometer polarization, the LOS orientation and on the characteristics of the magnetic field.

The best overall precision is found for a radiometer with a linear vertical polarization. For instance, at NH high latitudes, the LOS wind retrieval error at 99 km is 6 m s^{−1} using a linear vertical polarization, but degrades to about 10 m s^{−1} for other polarization settings. Furthermore, using the linear vertical polarization yields homogeneous results for different observation geometries: we could not find significant differences between ascending and descending orbits or between the two antennas.

Figure 11 shows the retrieval errors on the three components of ** B** at 85 km and 105 km (vertical resolution of 20 km).
The results strongly depend on the radiometer's polarization. Best performance is achieved with a $\pm \mathrm{45}{}^{\circ}$ linear polarization. Errors are clearly smaller when the retrieved component is aligned with the background magnetic field: the error on

Contrary to the results shown in Sect. 5.1 where it was the optimal configuration, the linear vertical polarization yields a worse retrieval performance for ** B**. In this case, the retrieval errors on

Our results show that the sensitivity of the SMILES-2 instrument is high enough to potentially measure the electrojet-induced variations of ** B** at high latitudes even under quiet sun conditions, provided that the data are properly averaged.
Yee et al. (2017) used the Zeeman effect on the AURA/MLS O

Perturbations of the geomagnetic field near the Equator (30 and 80 nT for the surface vertical and horizontal components of ** B**) are much smaller than the retrieval precision (Doumbia et al., 2007). Therefore, extracting interesting information on the equatorial jet will be more challenging, and a receiver with a slant polarization could be necessary.

6 Conclusions

Back to toptop
This analysis demonstrates the potential of SMILES-2 for the measurement of the temperature, atmospheric density and LOS wind in the MLT (60–110 km).
The retrieval precision was assessed, focusing on the SMILES-2 band at 760 GHz, the most suitable for such measurements.
Special care was taken to properly include the Zeeman effect on the O_{2} line.
Our results showed that neglecting it could lead to underestimating the retrieval errors by a factor of up to 2 above
90 km.
Because the O_{2} line is polarized, the radiometer's polarization configuration had to be investigated. We found that the optimal configuration was vertically linear. The LOS wind is retrieved with a precision of
2–5 m s^{−1} up to 90 km (30 m s^{−1} at 110 km) and
a vertical resolution of 2.5 km. Temperature and atmospheric density are
retrieved with a precision better than 5 K
(30 K) and 7 % (20 %) up to 90 km (110 km), respectively.
The achieved precision of the wind measurements, a key product for SMILES-2, is
comparable to the requirements for the new ICON mission (Englert et al., 2017).
However, unlike optical sensors, SMILES-2 can acquire high-precision measurements during day and night, and at all latitudes, even during auroral events. The low noise level achieved by the 4 K super-cooled radiometers is essential to achieve good performance above 90 km, where sensitivity becomes critical due to significantly weaker signals.

The retrieval of the geomagnetic field using the O_{2} line was also discussed.
We showed that valuable information on the horizontal and vertical components of ** B** could be determined directly near the E-region auroral electrojets.
Yee et al. (2017) highlighted the need for such observations since, currently, only measurements from the ground or from low-orbit satellites near 400 km are available.
Yee et al. (2017) proposed a CubeSat constellation, with the purpose of measuring the O

The final instrumental setup is still under discussion. In terms of possible instrumental developments, the spectral bandwidth of the 763 GHz band might be reduced in the definitive configuration of SMILES-2. Narrowing the bandwidth by a factor of 2 (while ensuring a correct adjustment of the LO frequency) would cause minimal degradation of the measurement performance, limited to altitudes below about 40 km.

Future work to improve MLT retrievals will include the two other SMILES-2 bands.
Indeed, the atomic oxygen line at 2 THz contains temperature and wind information above 100 km. This line can help us to improve the wind retrieval precision to 10 m s^{−1} at 110 km (Baron et al., 2019b). In the 638 GHz band, a strong signal from O_{3} will be measured
below about 70 km in daytime and 90 km in nighttime.
Furthermore, new parameters for the Zeeman model became recently available (Larsson et al., 2019). Applying the updated parameters should induce a change of the O_{2} and O line intensities, of up to a few percent.
The Zeeman effect on other spectral lines, OH, NO and ClO, should also be studied.

Appendix A: Spectroscopic parameters

Back to toptop
The spectroscopic parameters are taken from the HITRAN database (Rothman et al., 2009). The line strength at the temperature *T* is as follows:

$$\begin{array}{}\text{(A1)}& \begin{array}{rl}& S\left(T\right)={\displaystyle \frac{{C}_{\mathrm{H}}}{{r}_{\mathrm{iso}}}}{S}_{\mathrm{H}}\left({T}_{\mathrm{0}}\right){\displaystyle \frac{{e}^{-{C}_{\mathrm{E}}\phantom{\rule{0.125em}{0ex}}{E}_{\mathrm{L}}/{k}_{\mathrm{b}}T}}{{e}^{-{C}_{\mathrm{E}}\phantom{\rule{0.125em}{0ex}}{E}_{\mathrm{L}}/{k}_{\mathrm{b}}{T}_{\mathrm{0}}}}}\\ & \left({\displaystyle \frac{\mathrm{1}-{e}^{-{C}_{\mathrm{E}}\phantom{\rule{0.125em}{0ex}}{\stackrel{\mathrm{\u203e}}{\mathit{\nu}}}_{\mathrm{0}}/{k}_{\mathrm{b}}T}}{\mathrm{1}-{e}^{-{C}_{\mathrm{E}}\phantom{\rule{0.125em}{0ex}}{\stackrel{\mathrm{\u203e}}{\mathit{\nu}}}_{\mathrm{0}}/{k}_{\mathrm{b}}{T}_{\mathrm{0}}}}}\right){\displaystyle \frac{Q\left({T}_{\mathrm{0}}\right)}{Q\left(T\right)}}\left(\mathrm{Hz}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{\mathrm{molecule}}^{-\mathrm{1}}\right),\end{array}\end{array}$$

where ${k}_{\mathrm{b}}=\mathrm{1.380662}\times {\mathrm{10}}^{-\mathrm{23}}$ J K^{−1} is the Boltzmann constant, ${\stackrel{\mathrm{\u203e}}{\mathit{\nu}}}_{\mathrm{0}}$ (cm^{−1}) is the transition wavenumber, *S*_{H}(*T*_{0}) is the HITRAN line strength (cm^{−1} cm^{2} molecule^{−1}), *T*_{0}=296 K and *E*_{L} (cm^{−1}) is the lowest energy of the transition. The partition function *Q* is calculated from tabulated values between 120 and 500 K, a range that encompasses the temperatures found between 50 and 130 km (*Q*(296)=215.77). The constants *C*_{E}=10^{2} *h*_{p} *c* and ${C}_{\mathrm{H}}={\mathrm{10}}^{-\mathrm{2}}\phantom{\rule{0.125em}{0ex}}c$ allow the conversion of the HITRAN units to the International System (SI) ones. The isotopic ratio *r*_{iso} is taken away from *S*_{H} and added to the density profile instead. The Table A1 shows parameters of the main O_{2} millimeter lines.

Above the altitude of about 70 km, the real part of the Voigt function *F* (Eq. 3) is close to the Gauss function that describes lines broadened by random molecular velocities (Doppler broadening):

$$\begin{array}{}\text{(A2)}& F\left(\mathit{\nu}\right)={\displaystyle \frac{\mathrm{1}}{\mathrm{\Delta}{\mathit{\nu}}_{\mathrm{d}}}}{\left({\displaystyle \frac{\mathrm{ln}\mathrm{2}}{\mathit{\pi}}}\right)}^{{\scriptscriptstyle \frac{\mathrm{1}}{\mathrm{2}}}}{e}^{-\mathrm{ln}\mathrm{2}{\left({\scriptscriptstyle \frac{\mathit{\nu}-{\mathit{\nu}}_{\mathrm{0}}}{\mathrm{\Delta}{\mathit{\nu}}_{\mathrm{d}}}}\right)}^{\mathrm{2}}}\left({\text{Hz}}^{-\mathrm{1}}\right)\end{array}$$

with

$$\begin{array}{}\text{(A3)}& \mathrm{\Delta}{\mathit{\nu}}_{\mathrm{d}}={\displaystyle \frac{{\mathit{\nu}}_{\mathrm{0}}}{c}}{\left({\displaystyle \frac{\mathrm{2}\mathrm{ln}\mathrm{2}\phantom{\rule{0.125em}{0ex}}R\phantom{\rule{0.125em}{0ex}}T}{M}}\right)}^{{\scriptscriptstyle \frac{\mathrm{1}}{\mathrm{2}}}}\left(\text{Hz}\right),\end{array}$$

and Δ*ν*_{d} is the Doppler broadening parameter (i.e., the half-width at half-maximum, HWHM) of *F*, *ν*_{0} is the frequency of the transition, $c=\mathrm{2.997924}\times {\mathrm{10}}^{\mathrm{8}}$ m s^{−1} is the speed of light in vacuum, *R*=8.31446 J K^{−1} mol^{−1} is
the gas constant and *M* is the molar mass (0.031980 kg mol^{−1} for O_{2}). At 80 km, Δ*ν*_{d} is about 0.6–0.7 MHz for the O_{2} line at 773 GHz, while the pressure broadening HWHM is only 0.01–0.02 MHz.

The dispersion profile used for the calculation of the coefficient *q*^{′}, *u*^{′} and *v*^{′} (Eq. 2) is given by the following:

$$\begin{array}{}\text{(A4)}& {F}^{\prime}\left(\mathit{\nu}\right)=\sqrt{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}F\left(\mathit{\nu}\right)\phantom{\rule{0.125em}{0ex}}\mathrm{erfi}\left(v\right)\left({\mathrm{Hz}}^{-\mathrm{1}}\right),\end{array}$$

with $v=\mathrm{ln}\mathrm{2}\left(\frac{\mathit{\nu}-{\mathit{\nu}}_{\mathrm{0}}}{\mathrm{\Delta}{\mathit{\nu}}_{\mathrm{d}}}\right)$ and $\mathrm{erfi}\left(v\right)=\sqrt{\mathrm{2}/\mathit{\pi}}{\int}_{\mathrm{0}}^{v}\phantom{\rule{0.125em}{0ex}}\mathrm{exp}\left({t}^{\mathrm{2}}\right)\phantom{\rule{0.125em}{0ex}}d$t the imaginary error function (Eq. 5.54 in Landi Degl'Innocenti and Landolfi, 2004).

Appendix B: Matrix exponential

Back to toptop
The computation of the matrix exponential in Eq. (1) is the performance bottleneck in our implementation of the radiative
transfer solver if we use a general algorithm.
A significantly faster algorithm has been implemented using the symmetry in
**K**_{o} (Eq. 2). The evolution operator **Λ** (Eq. 8) is written as
$\mathrm{exp}(-{k}_{\mathrm{a}}\phantom{\rule{0.125em}{0ex}}\mathrm{d}s)\mathrm{exp}\left({\stackrel{\mathrm{\u0303}}{\mathbf{K}}}_{\mathrm{o}}\right),$
with *k*_{a} the scalar absorption coefficient (Eq. 3) and ${\stackrel{\mathrm{\u0303}}{\mathbf{K}}}_{\mathrm{o}}=-{\mathbf{K}}_{\mathrm{o}}\phantom{\rule{0.125em}{0ex}}\mathrm{d}s$. The Cayley–Hamilton theorem is used to compute $\mathrm{exp}\left({\stackrel{\mathrm{\u0303}}{\mathbf{K}}}_{\mathrm{o}}\right)$:

$$\begin{array}{}\text{(B1)}& \mathrm{exp}\left({\stackrel{\mathrm{\u0303}}{\mathbf{K}}}_{\mathrm{o}}\right)=\sum _{k=\mathrm{0}}^{\mathrm{3}}{\mathit{\kappa}}_{k}{\stackrel{\mathrm{\u0303}}{\mathbf{K}}}_{\mathrm{o}}^{k},\end{array}$$

where ${\stackrel{\mathrm{\u0303}}{\mathbf{K}}}_{\mathrm{o}}^{\mathrm{0}}$ is the identity matrix.
The coefficient *κ*_{k} are derived using the four eigenvalues of ${\stackrel{\mathrm{\u0303}}{\mathbf{K}}}_{\mathrm{o}}$:

$$\left\{\begin{array}{l}{\mathrm{exp}}^{{\mathit{\lambda}}_{\mathrm{1}}}={\mathit{\kappa}}_{\mathrm{0}}\phantom{\rule{0.125em}{0ex}}+{\mathit{\kappa}}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}+{\mathit{\kappa}}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}+{\mathit{\kappa}}_{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{3}},\\ {\mathrm{exp}}^{-{\mathit{\lambda}}_{\mathrm{1}}}={\mathit{\kappa}}_{\mathrm{0}}\phantom{\rule{0.125em}{0ex}}-{\mathit{\kappa}}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}+{\mathit{\kappa}}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}-{\mathit{\kappa}}_{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{3}},\\ {\mathrm{exp}}^{j\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{2}}}={\mathit{\kappa}}_{\mathrm{0}}\phantom{\rule{0.125em}{0ex}}+j\phantom{\rule{0.125em}{0ex}}{\mathit{\kappa}}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}-{\mathit{\kappa}}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}-j\phantom{\rule{0.125em}{0ex}}{\mathit{\kappa}}_{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{3}},\\ {\mathrm{exp}}^{-j\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{2}}}={\mathit{\kappa}}_{\mathrm{0}}\phantom{\rule{0.125em}{0ex}}-j\phantom{\rule{0.125em}{0ex}}{\mathit{\kappa}}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}-{\mathit{\kappa}}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}+j\phantom{\rule{0.125em}{0ex}}{\mathit{\kappa}}_{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{3}},\end{array}\right.$$

where *λ*_{1,2} is positive real-valued numbers that determine the four eigenvalues ±*λ*_{1} and ±*j**λ*_{2} of ${\stackrel{\mathrm{\u0303}}{\mathbf{K}}}_{\mathrm{o}}$.
This gives the following:

$$\begin{array}{}\text{(B2)}& \begin{array}{rl}& {\mathit{\kappa}}_{\mathrm{0}}={\displaystyle \frac{{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}\mathrm{cosh}\left({\mathit{\lambda}}_{\mathrm{1}}\right)\phantom{\rule{0.125em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}\mathrm{cos}\left({\mathit{\lambda}}_{\mathrm{2}}\right)}{{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{2}}+{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{2}}}},\\ & {\mathit{\kappa}}_{\mathrm{1}}={\displaystyle \frac{{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}\mathrm{sinh}\left({\mathit{\lambda}}_{\mathrm{1}}\right)/{\mathit{\lambda}}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}\mathrm{sin}\left({\mathit{\lambda}}_{\mathrm{2}}\right)/{\mathit{\lambda}}_{\mathrm{2}}}{{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{2}}+{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{2}}}},\\ & {\mathit{\kappa}}_{\mathrm{2}}={\displaystyle \frac{\mathrm{cosh}\left({\mathit{\lambda}}_{\mathrm{1}}\right)-\mathrm{cos}\left({\mathit{\lambda}}_{\mathrm{2}}\right)}{{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{2}}+{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{2}}}},\\ & {\mathit{\kappa}}_{\mathrm{3}}={\displaystyle \frac{\mathrm{sinh}\left({\mathit{\lambda}}_{\mathrm{1}}\right)/{\mathit{\lambda}}_{\mathrm{1}}-\mathrm{sin}\left({\mathit{\lambda}}_{\mathrm{2}}\right)/{\mathit{\lambda}}_{\mathrm{2}}}{{\mathit{\lambda}}_{\mathrm{1}}^{\mathrm{2}}+{\mathit{\lambda}}_{\mathrm{2}}^{\mathrm{2}}}}.\end{array}\end{array}$$

The eigenvalue parameters are ${\mathit{\lambda}}_{\mathrm{1}}=-\mathrm{d}s\phantom{\rule{0.125em}{0ex}}\sqrt{(A+B)/\mathrm{2}}$ and ${\mathit{\lambda}}_{\mathrm{2}}=-\mathrm{d}s\sqrt{(A-B)/\mathrm{2}}$, where

$$\begin{array}{rl}& A=\left[\right.\mathrm{8}\phantom{\rule{0.25em}{0ex}}(q\phantom{\rule{0.125em}{0ex}}{q}^{\prime}\phantom{\rule{0.125em}{0ex}}v\phantom{\rule{0.125em}{0ex}}{v}^{\prime}\phantom{\rule{0.25em}{0ex}}+\phantom{\rule{0.25em}{0ex}}u\phantom{\rule{0.125em}{0ex}}{u}^{\prime}\phantom{\rule{0.125em}{0ex}}q\phantom{\rule{0.125em}{0ex}}{q}^{\prime}\phantom{\rule{0.25em}{0ex}}+\phantom{\rule{0.25em}{0ex}}u\phantom{\rule{0.125em}{0ex}}{u}^{\prime}\phantom{\rule{0.125em}{0ex}}v\phantom{\rule{0.125em}{0ex}}{v}^{\prime})+{q}^{\mathrm{4}}\phantom{\rule{0.125em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{u}^{\mathrm{4}}\phantom{\rule{0.125em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{v}^{\mathrm{4}}\\ & +\phantom{\rule{0.125em}{0ex}}{{q}^{\prime}}^{\mathrm{4}}\phantom{\rule{0.125em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{{u}^{\prime}}^{\mathrm{4}}\phantom{\rule{0.125em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{{v}^{\prime}}^{\mathrm{4}}+\mathrm{2}\phantom{\rule{0.125em}{0ex}}\left(\right.{q}^{\mathrm{2}}({u}^{\mathrm{2}}-{{u}^{\prime}}^{\mathrm{2}}+{v}^{\mathrm{2}}-{{v}^{\prime}}^{\mathrm{2}}+{{q}^{\prime}}^{\mathrm{2}})\\ & +{u}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}({{u}^{\prime}}^{\mathrm{2}}+{v}^{\mathrm{2}}-{{v}^{\prime}}^{\mathrm{2}}-{{q}^{\prime}}^{\mathrm{2}})\\ & +{v}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}(-{{u}^{\prime}}^{\mathrm{2}}+{{v}^{\prime}}^{\mathrm{2}}-{{q}^{\prime}}^{\mathrm{2}})\\ & +{{v}^{\prime}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}({{u}^{\prime}}^{\mathrm{2}}+{{q}^{\prime}}^{\mathrm{2}})\\ & +{{u}^{\prime}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{{q}^{\prime}}^{\mathrm{2}}){}^{{\scriptscriptstyle \frac{\mathrm{1}}{\mathrm{2}}}}\end{array}$$

and

$$B={q}^{\mathrm{2}}-{{q}^{\prime}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{u}^{\mathrm{2}}-{{u}^{\prime}}^{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}+\phantom{\rule{0.125em}{0ex}}{v}^{\mathrm{2}}-{{v}^{\prime}}^{\mathrm{2}}.$$

Author contributions

Back to toptop
Author contributions.

PB designed the study, performed the simulations and wrote the paper. All co-authors provided valuable information. MS is the PI of the mission proposal.

Competing interests

Back to toptop
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements

Back to toptop
Acknowledgements.

Philippe Baron would like to thank Franz Schreier (German Aerospace Center, DLR) for providing the python implementation of the complex Voigt function used in GARLIC (see reference given for Eq. 3). We would like to thank Hugh C. Pumphrey and the two anonymous referees for their valuable comments that helped us to improve the paper.

Financial support

Back to toptop
Financial support.

SMILES-2 studies are supported by the strategic development research fund from the Institute of Space and Astronautical Science (ISAS)/JAXA. Huixin Liu acknowledges support by JSPS KAKENHI (grants no. 18H01270, 18H04446 and 17KK0095).

Review statement

Back to toptop
Review statement.

This paper was edited by Christian von Savigny and reviewed by Hugh C. Pumphrey and two anonymous referees.

References

Back to toptop
Baron, P., Murtagh, D. P., Urban, J., Sagawa, H., Ochiai, S., Kasai, Y., Kikuchi, K., Khosrawi, F., Körnich, H., Mizobuchi, S., Sagi, K., and Yasui, M.: Observation of horizontal winds in the middle-atmosphere between 30^{∘} S and 55^{∘} N during the northern winter 2009–2010, Atmos. Chem. Phys., 13, 6049–6064, https://doi.org/10.5194/acp-13-6049-2013, 2013. a

Baron, P., Manago, N., Ozeki, H., Irimajiri, Y., Murtagh, D., Uzawa, Y., Ochiai, S., Shiotani, M., and Suzuki, M.: Measurement of stratospheric and mesospheric winds with a submillimeter wave limb sounder: Results from JEM/SMILES and simulation study for SMILES-2, in: Sensors, Systems, and Next-Generation Satellites XIX, edited by: Meynart, R., Neeck, S. P., and Shimoda, H., International Society for Optics and Photonics, SPIE, 9639, 140–159, https://doi.org/10.1117/12.2194741, 2015. a

Baron, P., Murtagh, D., Eriksson, P., Mendrok, J., Ochiai, S., Pérot, K., Sagawa, H., and Suzuki, M.: Simulation study for the Stratospheric Inferred Winds (SIW) sub-millimeter limb sounder, Atmos. Meas. Tech., 11, 4545–4566, https://doi.org/10.5194/amt-11-4545-2018, 2018. a, b, c, d, e, f, g

Baron, P., Manago, N., Ochiai, S., and Suzuki, M.: SMILES-2 band selection study for chemical species, IGARSS 2019–2019, Int. Geosci. Remote Se., 7698–7701, https://doi.org/10.1109/IGARSS.2019.8897703, 2019a. a, b

Baron, P., Ochiai, S., Murtagh, D., Sagawa, H., Saito, A., Shiotani, M., and Suzuki, M.: Performance assessment of Super conducting Submillimeter-Wave Limb-Emission Sounder-2 (SMILES-2), IGARSS 2019–2019, Int. Geosci. Remote Se., 7556–7559, https://doi.org/10.1109/IGARSS.2019.8898496, 2019b. a, b, c, d

Baumgarten, G.: Doppler Rayleigh/Mie/Raman lidar for wind and temperature measurements in the middle atmosphere up to 80 km, Atmos. Meas. Tech., 3, 1509–1518, https://doi.org/10.5194/amt-3-1509-2010, 2010. a

Beig, G.: Long-term trends in the temperature of the mesosphere/lower thermosphere region: 1. Anthropogenic influences, J. Geophys. Res.-Space, 116, A00H11, https://doi.org/10.1029/2011JA016646, 2011. a

Chapman, S. and Lindzen, R. S.: Atmospheric tides, D. Reidel Publishing Company and Dordrecht-Holland, Springer Netherlands, 200 pp., https://doi.org/10.1007/978-94-010-3399-2, 1970. a

Christensen, O. M., Eriksson, P., Urban, J., Murtagh, D., Hultgren, K., and Gumbel, J.: Tomographic retrieval of water vapour and temperature around polar mesospheric clouds using Odin-SMR, Atmos. Meas. Tech., 8, 1981–1999, https://doi.org/10.5194/amt-8-1981-2015, 2015. a

Doumbia, V., Maute, A., and Richmond, A. D.: Simulation of equatorial electrojet magnetic effects with the thermosphere-ionosphere-electrodynamics general circulation model, J. Geophys. Res.-Space, 112, A09309, https://doi.org/10.1029/2007JA012308, 2007. a, b

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. a

Englert, C. R., Brown, C. M., Bach, B., Bach, E., Bach, K., Harlander, J. M., Seely, J. F., Marr, K. D., and Miller, I.: High-efficiency echelle gratings for MIGHTI, the spatial heterodyne interferometers for the ICON mission, Appl. Opt., 56, 2090–2098, https://doi.org/10.1364/AO.56.002090, 2017. a, b, c

Forbes, J. M., Russell, J., Miyahara, S., Zhang, X., Palo, S., Mlynczak, M., Mertens, C. J., and Hagan, M. E.: Troposphere-thermosphere tidal coupling as measured by the SABER instrument on TIMED during July–September 2002, J. Geophys. Res.-Space, 111, A10S06, https://doi.org/10.1029/2005JA011492, 2006. a

Fritts, D. C. and Alexander, M. J.: Gravity wave dynamics and effects in the middle atmosphere, Rev. Geophys., 41, 1003, https://doi.org/10.1029/2001RG000106, 2003. 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, https://doi.org/10.5194/amt-7-3633-2014, 2014. a

Gerber, D., Swinyard, B. M., Ellison, B. N., Plane, J. M. C., Feng, W., Navarathinam, N., Eves, S. J., Bird, R., Linfield, E. H., Davies, A. G., and Parkes, S.: LOCUS: Low cost upper atmosphere sounder, Proc. SPIE, 8889, 88891I, https://doi.org/10.1117/12.2028675, 2013. a

Gumbel, J., Megner, L., Christensen, O. M., Chang, S., Dillner, J., Ekebrand, T., Giono, G., Hammar, A., Hedin, J., Ivchenko, N., Karlsson, B., Kruse, M., Li, A., McCallion, S., Murtagh, D. P., Olentšenko, G., Pak, S., Park, W., Rouse, J., Stegman, J., and Witt, G.: The MATS Satellite Mission – Gravity Waves Studies by Mesospheric Airglow/Aerosol Tomography and Spectroscopy, Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2018-1162, in review, 2018. a

Hagen, J., Murk, A., Rüfenacht, R., Khaykin, S., Hauchecorne, A., and Kämpfer, N.: WIRA-C: a compact 142-GHz-radiometer for continuous middle-atmospheric wind measurements, Atmos. Meas. Tech., 11, 5007–5024, https://doi.org/10.5194/amt-11-5007-2018, 2018. a

Immel, T. J., Sagawa, E., England, S. L., Henderson, S. B., Hagan, M. E., Mende, S. B., Frey, H. U., Swenson, C. M., and Paxton, L. J.: Control of equatorial ionospheric morphology by atmospheric tides, Geophys. Res. Lett., 33, L15108, https://doi.org/10.1029/2006GL026161, 2006. a

Jacobi, C., Lilienthal, F., Geissler, C., and Krug, A.: Long-term variability
of mid-latitude mesosphere-lower thermosphere winds over Collm (51^{∘} N, 13^{∘} E),
J. Atmos. Sol.-Terr. Phy., 136, 174–186,
https://doi.org/10.1016/j.jastp.2015.05.006, 2015. a

Karlsson, B. and Becker, E.: How Does Interhemispheric Coupling Contribute to Cool Down the Summer Polar Mesosphere?, J. Climate, 29, 8807–8821, https://doi.org/10.1175/JCLI-D-16-0231.1, 2016. a

Kaufmann, M., Olschewski, F., Mantel, K., Solheim, B., Shepherd, G., Deiml, M., Liu, J., Song, R., Chen, Q., Wroblowski, O., Wei, D., Zhu, Y., Wagner, F., Loosen, F., Froehlich, D., Neubert, T., Rongen, H., Knieling, P., Toumpas, P., Shan, J., Tang, G., Koppmann, R., and Riese, M.: A highly miniaturized satellite payload based on a spatial heterodyne spectrometer for atmospheric temperature measurements in the mesosphere and lower thermosphere, Atmos. Meas. Tech., 11, 3861–3870, https://doi.org/10.5194/amt-11-3861-2018, 2018. a

Khosravi, M., Baron, P., Urban, J., Froidevaux, L., Jonsson, A. I., Kasai, Y., Kuribayashi, K., Mitsuda, C., Murtagh, D. P., Sagawa, H., Santee, M. L., Sato, T. O., Shiotani, M., Suzuki, M., von Clarmann, T., Walker, K. A., and Wang, S.: Diurnal variation of stratospheric and lower mesospheric HOCl, ClO and HO_{2} at the equator: comparison of 1-D model calculations with measurements by satellite instruments, Atmos. Chem. Phys., 13, 7587–7606, https://doi.org/10.5194/acp-13-7587-2013, 2013. a

Kikuchi, K., Nishibori, T., Ochiai, S., Ozeki, H., Irimajiri, Y., Kasai, Y., Koike, M., Manabe, T., Mizukoshi, K., Murayama, Y., Nagahama, T., Sano, T., Sato, R., Seta, M., Takahashi, C., Takayanagi, M., Masuko, H., Inatani, J., Suzuki, M., and Shiotani, M.: Overview and early results of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES), J. Geophys. Res., 115, D23306, https://doi.org/10.1029/2010JD014379, 2010. a

Landi Degl'Innocenti, M. and Landolfi, M.: Polarization in Spectral Lines, Astrophysics and Space Science Library, Springer Netherlands, https://books.google.co.jp/books?id=8sl2CkmZNWIC, 2004. a, b, c, d, e

Larsson, R., Ramstad, R., Mendrok, J., Buehler, S. A., and Kasai, Y.: A method for remote sensing of weak planetary magnetic fields: Simulated application to Mars, Geophys. Res. Lett., 40, 5014–5018, https://doi.org/10.1002/grl.50964, 2013. a

Larsson, R., Buehler, S. A., Eriksson, P., and Mendrok, J.: A treatment of the Zeeman effect using Stokes formalism and its implementation in the Atmospheric Radiative Transfer Simulator (ARTS), J. Quant. Spectrosc. Ra., 133, 445–453, https://doi.org/10.1016/j.jqsrt.2013.09.006, 2014. a, b, c, d, e

Larsson, R., Lankhaar, B., and Eriksson, P.: Updated Zeeman effect splitting coefficients for molecular oxygen in planetary applications, J. Quant. Spectrosc. Ra., 224, 431–438, https://doi.org/10.1016/j.jqsrt.2018.12.004, 2019. a

Lenoir, W. B.: Microwave spectrum of molecular oxygen in the mesosphere, J. Geophys. Res., 73, 361–376, https://doi.org/10.1029/JA073i001p00361, 1968. a, b

Mlynczak, M. G. and Yee, J. H.: LATTICE: The Lower Atmosphere-Thermosphere-Ionosphere Coupling Experiment, AGU Fall Meeting Abstracts, 1, 2017. a

Navas-Guzmán, F., Kämpfer, N., Murk, A., Larsson, R., Buehler, S. A., and Eriksson, P.: Zeeman effect in atmospheric O_{2} measured by ground-based microwave radiometry, Atmos. Meas. Tech., 8, 1863–1874, https://doi.org/10.5194/amt-8-1863-2015, 2015. a

Niciejewski, R., Wu, Q., Skinner, W., Gell, D., Cooper, M., Marshall, A., Killeen, T., Solomon, S., and Ortland, D.: TIMED Doppler Interferometer on the Thermosphere Ionosphere Mesosphere Energetics and Dynamics Satellite: Data product overview, J. Geophys. Res.-Space, 111, A11S90, https://doi.org/10.1029/2005JA011513, 2006. a

Ochiai, S., Kikuchi, K., Nishibori, T., Manabe, T., Ozeki, H., Mizobuchi, S., and Irimajiri, Y.: Receiver Performance of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on the International Space Station, IEEE T. Geosci. Remote, 51, 3791–3802, https://doi.org/10.1109/TGRS.2012.2227758, 2013. a

Ochiai, S., Baron, P., Nishibori, T., Irimajiri, Y., Uzawa, Y., Manabe, T., Maezawa, H., Mizuno, A., Nagahama, T., Sagawa, H., Suzuki, M., and Shiotani, M.: SMILES-2 mission for temperature, wind, and composition in the whole atmosphere, SOLA, 13A, 13–18, https://doi.org/10.2151/sola.13A-003, 2017. a, b, c, d

Ochiai, S., Baron, P., Irimajiri, Y., Uzawa, Y., Nishibori, T., Akinori, S., Suzuki, M., and Shiotani, M.: Superconducting Submillimeter-Wave Limb-Emission Sounder, SMILES-2, for middle and upper atmospheric study, IGARSS 2018–2018, IEEE T. Geosci. Remote, 9153–9156, https://doi.org/10.1109/IGARSS.2018.8518300, 2018. a

Ochiai, S., Baron, P., Irimajiri, Y., Nishibori, T., Hasegawa, T., Uzawa, Y., Maezawa, H., Manabe, T., Mizuno, A., Nagahama, T., Kimura, K., Suzuki, M., Saito, A., and Shiotani, M.: Conceptual study of Superconducting Submillimeter-Wave Limb-Emission Sounder-2 (SMILES-2) receiver, IGARSS 2019–2019, IEEE T. Geosci. Remote, 8792–8795, https://doi.org/10.1109/IGARSS.2019.8898693, 2019. a

Pancheva, D. and Mukhtarov, P.: Atmospheric Tides and Planetary Waves: Recent Progress Based on SABER/TIMED Temperature Measurements (2002–2007), Springer Netherlands, Dordrecht, 19–56, https://doi.org/10.1007/978-94-007-0326-1_2, 2011. a

Pardo, J., Pagani, L., Gerin, M., and Prigent, C.: Evidence of the zeeman
splitting in the $\mathrm{21}->\mathrm{01}$ rotational transition of the atmospheric
^{16}O^{18}O molecule from ground-based measurements, J. Quant. Spectrosc. Ra., 54, 931–943,
https://doi.org/10.1016/0022-4073(95)00129-9, 1995. a

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding: Theory and Practise, 2 of Series on Atmospheric, Oceanic and Planetary Physics, World Scientific, 256 pp., https://doi.org/10.1142/3171, 2000. a

Rothman, L., Gordon, I., Barbe, A., Benner, D., Bernath, P., Birk, M., Boudon, V., Brown, L., Campargue, A., Champion, J.-P., Chance, K., Coudert, L., Dana, V., Devi, V., Fally, S., Flaud, J.-M., Gamache, R., Goldman, A., Jacquemart, D., Kleiner, I., Lacome, N., Lafferty, W., Mandin, J.-Y., Massie, S., Mikhailenko, S., Miller, C., Moazzen-Ahmadi, N., Naumenko, O., Nikitin, A., Orphal, J., Perevalov, V., Perrin, A., Predoi-Cross, A., Rinsland, C., Rotger, M., Smith, M., Sung, K., Tashkun, S., Tennyson, J., Toth, R., Vandaele, A., and Auwera, J. V.: The HITRAN 2008 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 110, 533–572, https://doi.org/10.1016/j.jqsrt.2009.02.013, 2009. a

Rüfenacht, R., Murk, A., Kämpfer, N., Eriksson, P., and Buehler, S. A.: Middle-atmospheric zonal and meridional wind profiles from polar, tropical and midlatitudes with the ground-based microwave Doppler wind radiometer WIRA, Atmos. Meas. Tech., 7, 4491–4505, https://doi.org/10.5194/amt-7-4491-2014, 2014. a

Sakazaki, T., Fujiwara, M., Mitsuda, C., Imai, K., Manago, N., Naito, Y., Nakamura, T., Akiyoshi, H., Kinnison, D., Sano, T., Suzuki, M., and Shiotani, M.: Diurnal ozone variations in the stratosphere revealed in observations from the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on board the International Space Station (ISS), J. Geophys. Res.-Atmos., 118, 2991–3006, https://doi.org/10.1002/jgrd.50220, 2013. a

Sakazaki, T., Sato, K., Kawatani, Y., and Watanabe, S.: Three-dimensional structures of tropical nonmigrating tides in a high-vertical-resolution general circulation model, J. Geophys. Res.-Atmos., 120, 1759–1775, https://doi.org/10.1002/2014JD022464, 2015. a

Sato, K., Yasui, R., and Miyoshi, Y.: The Momentum Budget in the Stratosphere, Mesosphere, and Lower Thermosphere, Part I: Contributions of Different Wave Types and In Situ Generation of Rossby Waves, J. Atmos. Sci., 75, 3613–3633, https://doi.org/10.1175/JAS-D-17-0336.1, 2018. a

Schreier, F., Garcia, S. G., Hedelt, P., Hess, M., Mendrok, J., Vasquez, M., and Xu, J.: GARLIC – A general purpose atmospheric radiative transfer line-by-line infrared-microwave code: Implementation and evaluation, J. Quant. Spectrosc. Ra., 137, 29–50, https://doi.org/10.1016/j.jqsrt.2013.11.018, 2014. a, b

Schwartz, M. J., Read, W. G., and Snyder, W. V.: EOS MLS forward model polarized radiative transfer for Zeeman-split oxygen lines, IEEE T. Geosci. Remote, 44, 1182–1191, https://doi.org/10.1109/TGRS.2005.862267, 2006. a, b, c, d, e

Semel, M. and López, A.: Integration of the radiative transfer equation for polarized light: The exponential solution, A&A, 342, 201–211, 1999. 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 O_{2} A-band emission
spectra, Can. J. Phys., 88, 919–925, https://doi.org/10.1139/p10-093, 2010. a

Shepherd, G. G.: Development of wind measurement systems for future space missions, Acta Astronaut., 115, 206–217, https://doi.org/10.1016/j.actaastro.2015.05.015, 2015. a

Shiotani, M., Saito, A., Sakazaki, T., Ochiai, S., Baron, P., Nishibori, T., Suzuki, M., Abe, T., Maezawa, H., and Oyama, S.: A proposal for satellite observation of the whole atmosphere – Superconducting Submillimeter-Wave Limb-Emission Sounder-2 (SMILES-2), IGARSS 2019–2019, IEEE T. Geosci. Remote, 8788–8791, https://doi.org/10.1109/IGARSS.2019.8898423, 2019. a, b

Sica, R. J., Izawa, M. R. M., Walker, K. A., Boone, C., Petelina, S. V., Argall, P. S., Bernath, P., Burns, G. B., Catoire, V., Collins, R. L., Daffer, W. H., De Clercq, C., Fan, Z. Y., Firanski, B. J., French, W. J. R., Gerard, P., Gerding, M., Granville, J., Innis, J. L., Keckhut, P., Kerzenmacher, T., Klekociuk, A. R., Kyrö, E., Lambert, J. C., Llewellyn, E. J., Manney, G. L., McDermid, I. S., Mizutani, K., Murayama, Y., Piccolo, C., Raspollini, P., Ridolfi, M., Robert, C., Steinbrecht, W., Strawbridge, K. B., Strong, K., Stübi, R., and Thurairajah, B.: Validation of the Atmospheric Chemistry Experiment (ACE) version 2.2 temperature using ground-based and space-borne measurements, Atmos. Chem. Phys., 8, 35–62, https://doi.org/10.5194/acp-8-35-2008, 2008. a

Smith, A. K.: Global Dynamics of the MLT, Surv. Geophys., 33, 1177–1230, https://doi.org/10.1007/s10712-012-9196-9, 2012. a, b

Steinbrecht, W., McGee, T. J., Twigg, L. W., Claude, H., Schönenborn, F., Sumnicht, G. K., and Silbert, D.: Intercomparison of stratospheric ozone and temperature profiles during the October 2005 Hohenpeißenberg Ozone Profiling Experiment (HOPE), Atmos. Meas. Tech., 2, 125–145, https://doi.org/10.5194/amt-2-125-2009, 2009. a

Steiner, O., Züger, F., and Belluzzi, L.: Polarized radiative transfer in discontinuous media, A42, 586, 14 pp., https://doi.org/10.1051/0004-6361/201527158, 2016. a

Tsuda, T.: Characteristics of atmospheric gravity waves observed using the MU (Middle and Upper atmosphere) radar and GPS (Global Positioning System) radio occultation, P. Jpn. Acad., 90, 12–27, 2014. a

Tsutsumi, M., Sato, K., Sato, T., Kohma, M., Nakamura, T., Nishimura, K., and Tomikawa, Y.: Characteristics of Mesosphere Echoes over Antarctica Obtained Using PANSY and MF Radars, SOLA, 13A, 19–23, https://doi.org/10.2151/sola.13A-004, 2017. a

Urban, J., Baron, P., Lautié, N., Schneider, N., Dassas, K., Ricaud, P., and De La Noë, J.: MOLIERE (v5): a versatile forward- and inversion model for the millimeter and sub-millimeter wavelength range, J. Quant. Spectrosc. Ra. 83, 529–554, https://doi.org/10.1016/S0022-4073(03)001040-3, 2004. a

Wang, W., Wang, Z., and Duan, Y.: Performance evaluation of THz Atmospheric Limb Sounder (TALIS) of China, Atmos. Meas. Tech., 13, 13–38, https://doi.org/10.5194/amt-13-13-2020, 2020. a

Wu, D. L., Schwartz, M. J., Waters, J. W., Limpasuvan, V., Wu, Q. A., and Killeen, T. L.: Mesospheric doppler wind measurements from Aura Microwave Limb Sounder (MLS), Adv. Space Res., 42, 1246–1252, 2008. a

Wu, D. L., Yee, J.-H., Schlecht, E., Mehdi, I., Siles, J., and Drouin, B. J.: THz limb sounder (TLS) for lower thermospheric wind, oxygen density, and temperature, J. Geophys. Res.-Space, 121, 7301–7315, https://doi.org/10.1002/2015JA022314, 2016. a

Xu, X., Manson, A. H., Meek, C. E., Chshyolkova, T., Drummond, J. R., Hall, C. M., Riggin, D. M., and Hibbins, R. E.: Vertical and interhemispheric links in the stratosphere-mesosphere as revealed by the day-to-day variability of Aura-MLS temperature data, Ann. Geophys., 27, 3387–3409, https://doi.org/10.5194/angeo-27-3387-2009, 2009. a

Yamazaki, Y. and Maute, A.: Sq and EEJ – A Review on the Daily Variation of the Geomagnetic Field Caused by Ionospheric Dynamo Currents, Space Sci. Rev., 206, 299–405, https://doi.org/10.1007/s11214-016-0282-z, 2017. a

Yee, J. H., Gjerloev, J., Wu, D., and Schwartz, M. J.: First Application of the Zeeman Technique to Remotely Measure Auroral Electrojet Intensity From Space, Geophys. Res. Lett., 44, 10134–10139, https://doi.org/10.1002/2017GL074909, 2017. a, b, c, d, e, f

Yi, W., Xue, X., Reid, I. M., Murphy, D. J., Hall, C. M., Tsutsumi, M., Ning, B., Li, G., Vincent, R. A., Chen, J., Wu, J., Chen, T., and Dou, X.: Climatology of the mesopause relative density using a global distribution of meteor radars, Atmos. Chem. Phys., 19, 7567–7581, https://doi.org/10.5194/acp-19-7567-2019, 2019. a

Short summary

Submillimeter-Wave Limb-Emission Sounder 2 (SMILES-2) is a satellite mission proposed in Japan to probe the middle and upper atmosphere (20–160 km). The key products are wind, temperature and density. If selected, this mission could provide new insights into vertical coupling in the atmosphere and could help improve weather and climate models. We conducted simulation studies to assess the measurement performances in the altitude range 60–110 km, with a special focus on the geomagnetic effects.

Submillimeter-Wave Limb-Emission Sounder 2 (SMILES-2) is a satellite mission proposed in Japan...

Atmospheric Measurement Techniques

An interactive open-access journal of the European Geosciences Union