Introduction

Ultraviolet imaging observation is a useful tool for diagnosing the state of the ionosphere. In particular, measurements of emissions associated with \({\mathrm O}^+\) ions, which are the dominant species in the lower ionosphere, provide valuable information on spatio-temporal variations of the ionosphere. The most prominent emission from \({\mathrm O}^+\) is the 83.4 extreme ultraviolet (EUV) emission (e.g., Meier 1991) and several observations of this emission have been conducted (e.g., McCoy et al. 1985; Stephan et al. 2019). The observations of airglow at 91.1 nm and 135.6 nm are used for investigating \({\mathrm O}^+\) distributions, especially on the nightside (e.g., Dymond et al. 2004; Sagawa et al. 2005; Immel et al. 2006; Dymond et al. 2017c; Eastes et al. 2019).

The International Space Station–ionosphere-mesosphere-atmosphere Plasmasphere cameras (ISS-IMAP) mission operated two kinds of imagers on board the International Space Station (ISS). One was the Visible and near-Infrared Spectral Imager (VISI) (Sakanoi et al. 2011) and the other was the Extreme Ultraviolet Imagers (EUVI) (Yoshikawa et al. 2011). EUVI consisted of two telescopes, EUVI-A and EUVI-B. EUVI-A was designed to observe the 30.4 nm radiation from \({\mathrm He}^+\) and EUVI-B was intended to observe the 83.4 nm radiation from \({\mathrm O}^+\). As illustrated in Uji et al. (2014) and Hozumi et al. (2017), these two telescopes are directed toward the aft direction of the ISS orbit to conduct limb observations of the Earth’s ionosphere. The EUVI telescopes thus acquired information on vertical and horizontal structures of the ionosphere and plasmasphere.

The signals taken from EUVI-A were likely to be 30.4 nm EUV due to resonant scattering by \({\mathrm He}^+\) ions (Hozumi et al. 2017). However, the source of the signals from EUVI-B has not been verified. Hozumi (2013) reported that EUVI-B sometimes detected evident EUV signals in the umbra of the Earth, though these nighttime signals were not satisfactorily explained. Indeed, modeling of the EUVI-B measurements is not easy. Since 83.4 nm EUV undergoes multiple scattering by \({\mathrm O}^+\) ions in the ionosphere (e.g., Stephan , 2016; Geddes et al. , 2016), the relation with the \({\mathrm O}^+\) density profile along the line of sight would be complicated. In addition, it should be noted that the wavelength range covered by EUVI-B was not narrow enough. Since EUVI-B covered the wavelength range from about 70 nm to 110 nm (Yoshikawa et al. 2011; Uji et al. 2014), it did not exclude the contribution of other nightglow emissions, such as the 91.1 nm, 98.9 nm, and 102.6 nm emissions. In particular, since the 91.1 nm emission due to recombination between \({\mathrm O}^{+}\) ions and electrons is one of the prominent emissions on the nightside (Chakrabarti et al. 1984; Meier 1991; Feldman et al. , 1992), it is necessary to take into account its contribution.

The purpose of this paper is to clarify the source of the signals observed by EUVI-B to facilitate the use of the EUVI-B data. Since the EUVI-A and EUVI-B telescopes were operated when the ISS orbited in the umbra of the Earth (Uji et al. 2014), the EUV flux observed on the nightside is discussed. We conduct a Monte Carlo simulation to examine the effect of multiple scattering of the 83.4 nm EUV. We also examine the effect of the 91.1 nm EUV. The contribution of each effect is then quantitatively evaluated.

Observation

The operation of ISS-IMAP started in October 2012. The altitude of the ISS is about 400 km and the orbital inclination is \(51^{\circ }\). The ISS-IMAP telescopes thus conducted imaging observations of the ionosphere at mid and low latitudes. The field of view of the imager EUVI-B was \(13.2^{\circ }\times 13.2^{\circ }\). Figure 1 shows the EUV images taken from EUVI-B during 21:08:26–21:09:25 UT and 21:18:38–21:19:37 UT on December 26, 2012. In each panel, the orange contours indicate the altitude at the tangential point of the line of sight for each pixel. The altitudes from about 100–400 km were covered by the field of view of this imager. Each pixel shows the count of photons. The expected count C for \(4\pi I\) Rayleighs of EUV flux was

$$\begin{aligned} C= \int 4\pi I(t) S\,dt\simeq 4\pi IST, \end{aligned}$$
(1)

where S is the total sensitivity and T is the exposure time. S was equal to 0.0019 [counts/s/Rayleigh/bin] for EUVI-B at 83.4 nm (Uji et al. 2014). The exposure time T was 60 s here. The expected count for 100 Rayleighs is thus about 11.4 for the 83.4 nm EUV. Figure 2 shows the position of the ISS from 21:02 UT to 21:25 UT on December 26, 2012. Since each telescope was directed toward the backward direction of the ISS orbit as described above, Fig. 1 captured the northward (equatorward) view from the ISS location. The ISS was located near the equator around 21:08 UT, and it was located at mid latitude in the southern hemisphere around 21:18 UT.

In Fig. 3, the blue line indicates the time profile of the EUV intensity from 21:02 UT to 21:25 UT on December 26, 2012, which corresponds to the period indicated in Fig. 2. The EUV intensity is shown by the averaged count over the \(20\times 20\) pixels surrounded by the red square in each panel of Fig. 1. The error bars in this figure shows the \(\pm 2\sigma \) intervals where \(\sigma \) is the standard errors of the averaged count. Note that the observation was not necessarily taken every 60 s while the exposure time was 60 s. For reference, we plot the total electron content (TEC) value from IONEX (IOnosphere Map Exchange) at the ISS location with the orange line, and that at the tangential point of the line of sight with the green line. The averaged EUV count showed a bimodal profile, likely corresponding to the electron density structure indicated by the TEC observation. Figure 2 shows the position of the ISS for the corresponding period. The trough in the bimodal structure of the TEC was observed when the ISS crossed the dip equator, and thus it was regarded as the equatorial anomaly (e.g., Appleton 1946). The EUV observation thus corresponds well to the electron density structure. In the following, we work to reproduce this observational result.

Fig. 1
figure 1

EUV image taken from EUVI-B during 21:08:26–21:09:25 UT (a) and 21:18:38–21:19:37 UT (b) on December 26, 2012. The red square in each panel indicates the region where the EUV count is averaged to obtain Fig. 3. The orange contours in each panel indicate the altitude (in \({\mathrm km}\)) at the tangential point along the line of sight

Fig. 2
figure 2

Position of the ISS for the period from 21:02 UT to 21:25 UT on December 26, 2012. The tip of an arrow indicates the tangential point at each time. The blue solid line on the map indicates the dip equator. The TEC value is shown with the color contour

Fig. 3
figure 3

Time profile of the averaged EUV count over the pixels surrounded by the red square in Fig. 1 (blue) during the period from 21:02 UT to 21:25 UT on December 26, 2012. The total electron content (TEC) value at the ISS location (orange) and that at the tangential point along the line of sight (green) are also plotted for reference. The horizontal axis indicates the universal time (UT) and the local time (LT) at the position of the ISS

Modeling of scattering

The EUV at 83.4 nm is attributed to resonant scattering by \({\mathrm O}^{+}\) ions. The EUV flux due to the scattering, \(F_s {\mathrm [photons/(m^{2}\cdot str\cdot s)]}\), is obtained by the following line-of-sight integral:

$$\begin{aligned} F_s =\int p(\nu ')e^{-\tau (\nu ', \ell )}\left[ \frac{1}{4\pi }\int {\bar{\sigma }}_s(\nu ) N({\varvec{r}}(\ell ))\Phi (\nu , {\varvec{r}}(\ell ), \omega )\,{\mathrm d}\omega \,{\mathrm d}\nu \right] \, {\mathrm d}\ell \,{\mathrm d}\nu ', \end{aligned}$$
(2)

where \({\varvec{r}}\) denotes the location, \(\ell \) is the distance of \({\varvec{r}}\) from the observation site, \(p(\nu ')\) is the probability density function for the frequency of the scattered EUV \(\nu '\), \(\tau \) is the optical depth, \({\bar{\sigma }}_s\) is the scattering cross section, N is the \({\mathrm O}^{+}\) number density, and \(\Phi \) is the photon flux to be scattered. The photon flux is written as a function of the frequency \(\nu \), the location \({\varvec{r}}\), and the direction \(\omega \). We hereinafter assume that the scattering is isotropic. We take the \({\mathrm O}^{+}\) number density N from the International Reference Ionosphere (IRI) model (Bilitza et al. 2014). The optical depth \(\tau \) is associated with absorption due to neutrals and it is given by combining the line-of-sight integral for each species:

$$\begin{aligned} \tau (\nu , \ell ) =\int _0^{\ell }{\bar{\sigma }}_s(\nu )N({\varvec{r}}(s))\,{\mathrm d}s+\sum _i\int _0^{\ell }\alpha _{i}n_i({\varvec{r}}(s))\,{\mathrm d}s , \end{aligned}$$
(3)

where i denotes absorbing species and \(\alpha _i\) and \(n_i\) are the absorbing cross section and number density for each of the species, respectively. In this study, the absorption due to \({\mathrm N}_2\), \({\mathrm O}\), and \({\mathrm O}_2\) is taken into account. The number densities for these absorbing species are given by the Naval Research Laboratory Mass Spectrometer and Incoherent Scatter Radar Exosphere (NRLMSISE-00) model (Picone et al. 2002), and the absorbing cross sections are taken to be \(10.1\times 10^{-22}\,{\mathrm m}^2\), \(3.9\times 10^{-22}\,{\mathrm m}^2\), and \(10.4\times 10^{-22}\,{\mathrm m}^2\) for \({\mathrm N}_2\), \({\mathrm O}\), and \({\mathrm O}_2\), respectively, according to Link et al. (1994). The probability density function of \(\nu '\) is given as a Gaussian distribution:

$$\begin{aligned} p(\nu ')=\frac{1}{\sqrt{2\pi \Delta _\nu ^2}} \exp \left( -\frac{(\nu ' - \nu _0)^2}{2\Delta _\nu ^2}\right) , \end{aligned}$$
(4)

according to Eq. (14) in Appendix 1.

The scattering cross section \({\bar{\sigma }}_s\) should be dependent on the EUV frequency \(\nu \) because of the thermal motion of the scattering \({\mathrm O}^{+}\). We obtain the scattering cross section \({\bar{\sigma }}_s\) by the following equation

$$\begin{aligned} {\bar{\sigma }}_s(\nu )=\frac{e^2}{4\varepsilon _0 mc\sqrt{2\pi \Delta _\nu ^2}}f_{12} \exp \left( -\frac{(\nu -\nu _0 )^2}{2\Delta _\nu ^2}\right) , \end{aligned}$$
(5)

as shown in Appendix1, where e is the elementary charge, \(\varepsilon _0\) is the vacuum permittivity, m is the mass of a scattering ion (\({\mathrm O}^+\)), c is the light speed, \(f_{12}\) is the oscillator strength, \(\nu _0\) is the wave frequency at the center of the \(83.4\,{\mathrm nm}\) resonance line, and \(\Delta _\nu \) is the scale of Doppler broadening given by Eq. (15). The oscillator strength is taken at \(f_{12}=0.122\) according to Meier (1991). In obtaining \(\Delta _\nu \), the ion temperature is taken from the IRI model. The photon flux \(\Phi \) in Eq. (2) depends on the direction \(\omega \). It is demanding to obtain \(\Phi \) as a function of \(\omega \) over the line of sight for computing Eq. (2). However, since we assume isotropic scattering, it is not necessary to resolve the dependence of the photon flux on the direction \(\omega \), and instead the number density of photons at \({\varvec{r}}\) is sufficient to obtain \(F_s\). Denoting the number density of photons (per volume and frequency) as \(\rho \), Eq. (2) can be rewritten as

$$\begin{aligned} \begin{aligned} F_s&=\int p(\nu ')e^{-\tau (\nu ', \ell )}\left[ \frac{1}{4\pi }\int {\bar{\sigma }}_s(\nu ) N({\varvec{r}}(\ell ))c\rho (\nu , {\varvec{r}}(\ell ))\,{\mathrm d}\nu \right] \,{\mathrm d}\ell \,{\mathrm d}\nu ' \\&=\frac{c}{4\pi }\int p(\nu ')e^{-\tau (\nu ', \ell )}N({\varvec{r}}(\ell ))\left[ \int {\bar{\sigma }}_s(\nu )\rho (\nu , {\varvec{r}}(\ell ))\,{\mathrm d}\nu \right] \,{\mathrm d}\ell \,{\mathrm d}\nu ' \\&=\frac{c}{4\pi }\int p(\nu ')e^{-\tau (\nu ', \ell )}N({\varvec{r}}(\ell ))A\,{\mathrm d}\ell \,{\mathrm d}\nu ', \end{aligned} \end{aligned}$$
(6)

where we denote

$$\begin{aligned} A=\int {\bar{\sigma }}_s(\nu )\rho (\nu , {\varvec{r}}(\ell ))\,{\mathrm d}\nu . \end{aligned}$$
(7)

While the integral with \({\mathrm d}\ell \) and \({\mathrm d}\nu '\) was computed using the trapezoidal rule, we estimated the integral A by the Monte Carlo simulation in which the effect of multiple-scattering was considered. In this Monte Carlo simulation, we computed the propagation of a large number of photons from the EUV seeds on the sunlit ionosphere, and estimated the spatial distribution of the photon density resulting from the scattering by \({\mathrm O}^+\) ions. The details of the Monte Carlo simulation are described in Appendix 2. As explained in the appendix, our Monte Carlo model would tend to overestimate the 83.4 nm flux for the target event on 26 December 2012 because of the choice of the parameter of the incident solar flux. The result of the Monte Carlo simulation would thus provide an upper-limit estimate of the 83.4 nm EUV flux.

Modeling of \({\mathrm O}^{+}\) recombination

The EUV at 91.1 nm is emitted as a result of recombination of \({\mathrm O}^{+}\) ions and electrons. The EUV flux due to the \({\mathrm O}^{+}\) recombination is obtained by the following equation:

$$\begin{aligned} F_{\mathrm r}=\frac{1}{4\pi }\int e^{-\tau (\ell )} \lambda n_{\mathrm e}({\varvec{r}}(\ell ))\,n_{{\mathrm O}^+}({\varvec{r}}(\ell ))\,{\mathrm d}\ell \end{aligned}$$
(8)

where \(n_{\mathrm e}\) is the electron density and \(n_{{\mathrm O}^+}\) is the \({\mathrm O}^{+}\) ion density. The electron density \(n_{\mathrm e}\) and the \({\mathrm O}^{+}\) density \(n_{{\mathrm O}^+}\) are taken from the IRI model; that is, we assume the same \({\mathrm O}^{+}\) density distribution as for the 83.4 nm scattering model. The factor \(\lambda \) is the radiative recombination coefficient and \(\lambda =3.5\times 10^{23}\times (1160/T)\,{\mathrm m^3/s}\) (Meléndez-Alvira et al. 1999) where T is the electron temperature in Kelvin. The optical depth \(\tau \) can be obtained according to Eq. (3). The absorbing species taken into account are \({\mathrm N}_2\), \({\mathrm O}\), and \({\mathrm O}_2\). The absorption coefficients for 91.1 nm EUV are set to \(14.5\times 10^{-22}\,{\mathrm m}^2\), \(3.93\times 10^{-22}\,{\mathrm m}^2\), and \(15.34\times 10^{-22}\,{\mathrm m}^2\) for \({\mathrm N}_2\), \({\mathrm O}\), and \({\mathrm O}_2\), respectively, according to Dymond et al. (2017c), which uses the values provided by Kirby et al. (1979) and Conway (1988). The number densities for the absorbing species are given according to the NRLMSISE-00 model.

Result

Figure 4 shows the model EUV time profile to be observed by EUVI-B for the 83.4 nm EUV as a result of resonant scattering (red line) and for the 91.1 nm EUV as a result of recombination (green line). In this figure, the expected count was obtained by the line-of-sight integral taken in the direction of the center of the square region over which the averaged EUV count was calculated in Fig. 3 (i.e., the region surrounded by the red square in each panel of Fig. 1). According to Fig. 3, we simulated the EUV measurement during the period 21:00–21:30 UT on December 26, 2012 with the expected count value for one minute which can be obtained with Eq. (1). The total sensitivity S was 0.0019 at 83.4 nm and 0.0018 at 91.1 nm. The simulation result indicated that the contribution of the 91.1 nm EUV dominated that of the 83.4 nm EUV during the period of interest when the ISS orbited on the nightside. In this figure, the mean count observed by EUVI-B, which was shown in Fig. 3, is overplotted (blue line). The expected count value for the simulated 91.1 nm EUV was quantitatively comparable to the observed averaged count although the simulated EUV intensity was a little weaker than the observed intensity for the second peak around 21:18 UT. On the other hand, the EUV intensity at 83.4 nm calculated with the Monte Carlo model was much less than the observed value. As explained in Appendix 2, our Monte Carlo simulation would gives an upper-limit estimate of the 83.4 nm EUV for this event. This result therefore suggests that the 83.4 nm EUV could not explain the observed EUV count in the umbra of the Earth. The simulation result for the 91.1 nm EUV also roughly reproduced the bimodal structure seen in the averaged count in Fig. 3. In contrast, Fig. 4 shows that the 83.4 nm EUV should only be discernable when the ISS orbited near the dusk. When the ISS crossed the equator, the 83.4 nm EUV was not detectable in this event because the ISS was located in the deep umbra. We therefore conclude that the observed signal from EUVI-B in Fig. 3 can be attributed to the 91.1 nm EUV. On the righthand side of Fig. 4, we show the scale in Rayleighs which is obtained by assuming that the observed signal is entirely due to the 91.1 nm EUV. The estimated photon flux values might contain systematic errors due to the uncertainty of the sensor sensitivity. However, since the sensor uncertainty cannot be quantified, the error bars in Fig. 4 just represents the \(\pm 2\sigma \) intervals for the standard errors of the estimated flux.

Figure 5 shows the simulated images of the 91.1 nm EUV for EUVI-B at 21:08:56 UT and 21:19:08 UT on December 26, 2012. The color scale for each panel was adjusted in accord with Fig. 1. There are some discrepancies between Figs. 1 and 5. For example, while Fig. 1a suggests the source of the EUV emission shifted toward higher altitudes around 21:09 UT, the simulated emissions from the lower ionosphere were not as weak as the observed emissions. Figure 6 display the differences of the observed images in Fig. 1 from the simulated images in Fig. 5. The blue colors indicate that the observed EUV intensity exceeded the simulated intensity. Figure 6a suggests that the simulation underestimated the emissions from the lower altitudes. However, the EUV flux from the higher altitudes were well reproduced. Figure 6b suggests that the model underestimated the observed EUV signals around 21:19 UT as also suggested in Fig. 4. These discrepancies could be caused by variation of local conditions compared to the climatology represented in the models our simulation relied on. We determined the electron and \({\mathrm O}^+\) densities by the IRI model. The densities of the absorbing neutrals were given by the NRLMSISE-00 model. Although these empirical models provide the average state for a given condition, they are not necessarily highly accurate for each particular event. In general, spatio-temporal variations tend to be smoothed in empirical models. That would be a reason why the dip of the EUV flux observed around 21:14 UT was less clear in the simulation result. Figure 7 compares the IONEX TEC values shown in Fig. 3 with the reconstructed TEC values derived from the IRI model. The IONEX TEC values at the ISS location and the tangential point are plotted with the thick orange and green lines, respectively. The IRI TEC values at the ISS location and the tangential point are plotted with the pale orange and green lines, respectively. There are visible differences between the IONEX TEC and the IRI TEC, which would introduce errors in our model. The trough in the bimodal structure is less clear in the IRI TEC. It should also be noted that the discrepancies might partly be caused by other nightglow emissions. Considering the wavelength range covered by EUVI-B (70–110 nm), EUVI-B can also detect 98.9 nm and 102.6 nm emissions. The 98.9 nm and 102.6 nm emissions are typically weaker than the 83.4 nm emission by one order (Abreu et al. 1984; Chakrabarti et al. 1984; Meier 1991; Dymond et al. 2017b). Moreover, the sensitivity of EUVI-B was 0.0011 at 98.9 nm and 0.00047 at 102.6 nm, while it was 0.0018 at 91.1 nm. The contributions of the 98.9 nm and 102.6 nm EUV would therefore be too small to explain the entire difference. However, those emissions might, in part, contribute to the difference of the observation from the model prediction which assumes that the EUVI-B signals are totally attributed to the 91.1 nm EUV. At any rate, rough patterns in the observations were successfully reproduced by the simulation for the 91.1 nm EUV based on Eq. (8).

Fig. 4
figure 4

Simulated EUV intensity for EUVI-B for 21:00–21:30 UT on December 26, 2012. The red line indicates the time profile for 83.4 nm EUV and the green line indicates the time profile for 91.1 nm EUV. The blue line shows the mean count observed by EUVI-B which was shown in Fig. 3 for comparison. The scale for converting the expected count into the photon flux in Rayleighs is also shown on the righthand side of this figure

Fig. 5
figure 5

Simulated images of 91.1 nm EUV for EUVI-B at 21:08:56 UT (a) and 21:19:08 UT (b) on December 26, 2012

Fig. 6
figure 6

Difference of the EUV images shown in Fig. 1 from the simulated images in Fig. 5 for the periods 21:08:26–21:09:25 UT (a) and 21:18:38–21:19:37 UT (b) on December 26, 2012

Fig. 7
figure 7

Comparison of the IONEX TEC values shown in Fig. 3 with the TEC values reconstructed from the IRI model. The IONEX TEC values at the ISS location and the tangential point along the line of sight are plotted with the thick orange and green lines, respectively. The IRI TEC values at the ISS location and the tangential point are plotted with the pale orange and green lines, respectively

Discussion

As shown in the previous section, the EUV signals observed by EUVI-B in the umbra of the Earth can most likely be attributed to the 91.1 nm EUV. Although some discrepancies were seen in detailed structures between the simulated signals at 91.1 nm and the observed signals from EUVI-B, the observations were mostly reproduced by the model for the 91.1 nm EUV described in “Modeling of O+ recombination” section. Figure 8 is a histogram of the observed photon flux to show how often EUV signals were observed by EUVI-B for each local time interval in nighttime. This histogram displays the frequency distribution of the photon flux, which was estimated by assuming that EUVI-B observed the 91.1 nm EUV, for the period from 21 December 2012 to 3 March 2013. In order to focus on low-latitude observations, the data taken when the ISS latitude was higher than \(20^{\circ }\) (\(20^{\circ }\)N and \(20^{\circ }\)S) were excluded from this histogram. This figure shows that EUVI-B frequently observed EUV signals especially from 20:00 to 0:00 local time. This result is encouraging for the use of IMAP/EUVI-B data for studies of electron and \({\mathrm O}^+\) density distributions in the ionosphere. Various efforts have been made to obtain the \({\mathrm O}^+\) density distribution from remote optical measurements (e.g., Dymond et al. 1997; Qin et al. 2015). However, it is complicated to make an estimation from the 83.4 nm EUV data because of multiple scattering (e.g., Geddes et al. 2016). In contrast, estimating the \({\mathrm O}^+\) density from the 91.1 nm EUV would be much easier because the relationship between the observed signals and the \({\mathrm O}^+\) density is relatively simple. In the F-region, since the electron density is approximately equal to the \({\mathrm O}^+\) density, Eq. (8) can be approximated as

$$\begin{aligned} F_{\mathrm r}\approx \frac{1}{4\pi }\int e^{-\tau (\ell )} \lambda n_{{\mathrm O}^+}^2\,{\mathrm d}\ell . \end{aligned}$$
(9)

If the densities of the absorbing neutrals are known and thus the optical depth \(\tau \) is known, only the \({\mathrm O}^+\) density \(n_{{\mathrm O}^+}\) is unknown. The problem for obtaining \(n_{{\mathrm O}^+}\) thus becomes simple. Indeed, Dymond et al. (2017a, 2017c) demonstrated the reconstruction of the \({\mathrm O}^+\) density distribution based on EUV data at 91.1 nm. Multiple images from different viewpoints according to the orbital motion can be used to estimate the three-dimensional ion distribution (e.g., Nakano et al. 2014). The EUVI-B data could thus be used to reconstruct the three-dimensional \({\mathrm O}^+\) distribution in the vicinity of the ISS orbit.

Fig. 8
figure 8

Two-dimensional histogram of the estimated EUV flux when the ISS latitude was less than \(20^{\circ }\) for each local time interval in nighttime

When the EUVI-B signals are treated as the 91.1 nm EUV, the contribution of the 83.4 nm EUV should be regarded as noise. Although the simulation result in Fig. 4 suggested that the contribution of the 83.4 nm EUV was negligible for most of the period of interest, it would be meaningful to evaluate how deeply the 83.4 nm EUV penetrates the umbra of the Earth. We computed the omnidirectional flux of the 83.4 nm EUV at an altitude of 400 km using the Monte Carlo model described in “Modeling of scattering” section, where the omnidirectional flux is the flux averaged over \(4\pi \) steradian. Figure 9 shows the estimated omnidirectional flux in Rayleighs for the nightside. The sunlit region at 400 km altitude was masked in this figure to highlight the penetration into the umbra. Since the EUV flux on the nightside would not be isotropic, the omnidirectional flux would not match the EUV flux observed by an imager. Figure 9, however, would provide a rough estimate of the contribution of the 83.4 nm EUV. The estimated omnidirectional flux could be discerned only near the terminator, and even there it was less than 1 Rayleigh. Therefore, it would be reasonable to assume the contribution of the 83.4 nm EUV to the EUVI-B measurements is negligible.

When using the EUVI-B data for estimating the \({\mathrm O}^+\) density, we should take care of other factors which can affect the estimation. One issue is the uncertainty of the sensitivity value S. Although the uncertainty of S is unknown and it has not been quantified, it might introduce some bias to the estimate. We should also take into consideration the possible contributions of the 98.9 nm and 102.6 nm emissions. Although those emissions are inferred to be small as described above, their effects might cause some errors when estimating the \({\mathrm O}^+\) density. However, these problems can be calibrated by comparing other independent ionospheric data such as the TEC data.

Fig. 9
figure 9

Estimated omnidirectional flux of the 83.4 nm EUV at an altitude of 400 km on the nightside (in Rayleighs). The sunlit region at 400 km altitude is masked and only the flux in the umbra of the Earth is displayed

Summary

In order to explain the EUV signals observed from IMAP/EUVI-B on board the ISS on the nightside, we conducted two numerical modelings. One is a Monte Carlo simulation to examine the effect of the 83.4 nm emission due to the resonant scattering by \({\mathrm O}^+\) ions. The other is a forward model to examine the effect of the 91.1 nm emission due to the recombination between electrons and \({\mathrm O}^+\) ions. The simulation result suggests that the EUV signals observed from EUVI-B during nighttime were mainly the 91.1 nm EUV. The contribution of the 83.4 nm emission was typically negligible during nighttime. Since reconstructing the \({\mathrm O}^+\) density distribution from the 91.1 nm EUV would be relatively easy, it is expected that the EUV signals observed from IMAP/EUVI-B will be useful for studies of electron and \({\mathrm O}^+\) density distributions in the ionosphere.