the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Reconstruction of precipitating electrons and threedimensional structure of a pulsating auroral patch from monochromatic auroral images obtained from multiple observation points
Mizuki Fukizawa
Takeshi Sakanoi
Yoshimasa Tanaka
Yasunobu Ogawa
Keisuke Hosokawa
Björn Gustavsson
Kirsti Kauristie
Alexander Kozlovsky
Tero Raita
Urban Brändström
Tima Sergienko
In recent years, aurora observation networks using highsensitivity cameras have been developed in the polar regions. These networks allow dimmer auroras, such as pulsating auroras (PsAs), to be observed with a high signaltonoise ratio. We reconstructed the horizontal distribution of precipitating electrons using computed tomography with monochromatic PsA images obtained from three observation points. The threedimensional distribution of the volume emission rate (VER) of the PsA was also reconstructed. The characteristic energy of the reconstructed precipitating electron flux ranged from 6 to 23 keV, and the peak altitude of the reconstructed VER ranged from 90 to 104 km. We evaluated the results using a model aurora and compared the model's electron density with the observed one. The electron density was reconstructed correctly to some extent, even after a decrease in PsA intensity. These results suggest that the horizontal distribution of precipitating electrons associated with PsAs can be effectively reconstructed from groundbased optical observations.
Please read the corrigendum first before continuing.

Notice on corrigendum
The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.

Article
(4876 KB)

The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.
 Article
(4876 KB)  Corrigendum
 BibTeX
 EndNote
Aurora computed tomography (ACT) is a method for reconstructing the threedimensional (3D) volume emission rate (VER) of auroral emission based on monochromatic auroral images obtained from multiple observation points (e.g., Aso et al., 1990). The horizontal distribution of precipitating electron flux can be simultaneously obtained using ACT without rocket or satellite observations (Tanaka et al., 2011). Previous studies have applied ACT to bright and wellshaped discrete auroras, such as the quiet arc during the substorm growth phase and multiple auroral arcs (Aso et al., 1990, 1993, 1998; Frey et al., 1996; Nygrén et al., 1997; Tanaka et al., 2011). However, ACT has not been applied to pulsating auroras (PsAs).
A PsA is a type of diffuse aurora that appears as irregular patches showing quasiperiodic on–off switching of its intensity with a periodicity of ∼ 2–20 s (Yamamoto, 1988). The intensity is somewhat dimmer than that of a typical discrete aurora (some hundreds of rayleighs (R) up to tens of kilorayleighs at 557.7 nm; a few hundred rayleighs to ∼ 10 kR at 427.8 nm) (McEwen et al., 1981). It has been difficult to apply ACT to PsAs because the signaltonoise ratio (SNR) of PsA images is lower than those of discrete aurora images. However, remote operation of many highsensitivity cameras via the internet and an archive system capable of storing a massive amount of aurora data make it possible to observe PsAs with a high SNR.
The Magnetometers Ionospheric Radars Allsky Cameras Large Experiment (MIRACLE) network consists of nine allsky cameras (ASCs) located in the Fennoscandian region. Two of the ASCs with intensified chargecoupled devices (ICCDs) were replaced with cameras possessing newer electronmultiplying CCD (EMCCD) technology in 2007 (Sangalli et al., 2011). Ogawa et al. (2020a) developed a lowcost multiwavelength imaging system for aurora and airglow studies and installed Watec monochromatic imagers (WMIs) at several locations in the north and south polar regions. A WMI consists of a highly sensitive CCD camera made by Watec Co., Ltd (Japan). These cameras are suitable for studying very faint auroral structures such as PsAs. In this study, we attempted to use these highsensitivity cameras and ACT methods to reconstruct the 3D VER of a PsA patch and the horizontal distribution of precipitating electrons for the first time.
MIRACLE ASCs observed PsA patches from Kilpisjärvi (KIL; 69.05^{∘} N, 20.36^{∘} E), Abisko (ABK; 68.36^{∘} N, 18.82^{∘} E), and WMI ASCs at Skibotn (SKB; 69.35^{∘} N, 18.82^{∘} E) during the substorm recovery phase from 00:00 to 02:00 UT on 18 February 2018. These ASCs have a common field of view, as shown in Fig. 1b. The position of Tromsø (TRO; 69.58^{∘} N, 19.23^{∘} E), where the European Incoherent Scatter Scientific Association (EISCAT) radar operates, is also shown. We selected 427.8 nm auroral images in which PsA patches were detected at the EISCAT radar observation point. The reconstructed results were compared with the electron density observed using the EISCAT radar in Sect. 3.4. Figure 1a shows 427.8 nm auroral images obtained by the three ASCs from 00:53:30 to 00:53:42 UT. The temporal resolution of each ASCs was 2 s. A median filter of 3 × 3 pixels was applied to auroral images to improve the SNR. We also composited auroral images obtained from four WMI CCD cameras of the same type at SKB. The auroral images at SKB were calibrated to derive the absolute emission intensity using the optical facilities at the National Institute of Polar Research, Japan (Ogawa et al., 2020b), and have a time ambiguity of ∼ 1–2 s. We determined the time when the auroral image was obtained by aligning the temporal changes in the PsA patch, as shown in auroral images from ABK, KIL, and SKB. The ACT method used in this study is based on the method proposed by Tanaka et al. (2011). We adopted an oblique coordinate system with the origin (O) at coordinates of 69.4^{∘} N, 19.2^{∘} E. The x axis was antiparallel to the horizontal component of the geomagnetic field, the y axis was eastward, and the z axis was antiparallel to the geomagnetic field and perpendicular to the y axis (see Fig. 2 in Tanaka et al., 2011). The simulation region ranged from −75 to 75 km, from −100 to 100 km, and from 80 to 180 km for the x, y, and z axes, respectively. We set the energy (E) range to extend from 300 eV to 100 keV. The energy axis contains the information on the auroral emission altitude. This is because the higher the electron energy, the lower the stopping height of the precipitating electrons. This reconstruction region was divided linearly into n_{x} × n_{y} × n_{z} voxels along the x, y, and z axes and logarithmically into n_{E} bins in the E direction. We set the parameters n_{x}, n_{y}, n_{z}, and n_{E} to 75, 100, 50, and 50, respectively, corresponding to a spatial mesh size of 2 × 2 × 2 km. These parameters were selected so that each voxel has at least one lineofsight crossing from the pixels in the auroral images.
The differential flux of precipitating electrons was reconstructed by maximizing the posterior probability $P\left(\mathit{f}\stackrel{\mathrm{\u0303}}{\mathit{g}}\right)$, where f is a vector of the differential flux of precipitating electrons, and $\stackrel{\mathrm{\u0303}}{\mathit{g}}$ is a vector of gray levels at pixels in the auroral images obtained with ASCs. According to Bayes' theorem, the posterior probability $P\left(\mathit{f}\stackrel{\mathrm{\u0303}}{\mathit{g}}\right)$ is given by (Tanaka et al., 2011)
where g(f) is a vector of gray levels obtained by lineintegrating the VER in the lineofsight direction from each pixel (Eq. 8 in Tanaka et al., 2011). The VER was derived from model f using the aurora emission model (Eq. 3 in Tanaka et al., 2011). Σ^{−1} is the inverse covariance matrix, σ is the variance of ∇^{2}f, and the secondorder derivative of f is taken with respect to x, y, and E. The second term in Eq. (1) represents the smoothness of f in space and energy directions. We determined Σ^{−1} as the standard deviation calculated from each auroral image. To determine the standard deviation at each pixel in an auroral image, we calculated the mean value and standard deviation in the central 5 × 5 pixel regions for 120 images. Next, we derived a regression line between the mean value and the standard deviation. Finally, we converted the gray level at each pixel to the standard deviation using the equation of the derived regression line. To maximize the posterior probability, it is necessary to minimize the function
where λ, λ_{E}, and c_{j} are the socalled hyperparameters, which are constants corresponding to the weighting factors for the spatial (λ) and energy (λ_{E}) derivative terms and the correction factors for the relative sensitivity between cameras (c_{j}), respectively. The subscript j signifies the three observation points (ABK, KIL, and SKB). The parameter c_{SKB} was fixed at 1. The summation was conducted for the first term in Eq. (2) because c_{j} and ${\mathbf{\Sigma}}_{j}^{\mathrm{1}}$ were different for the three ASCs.
We carried out the change of variables f=exp (Σ) to take advantage of the nonnegative constraint on the differential flux f (i.e., f≥0). We then minimized the function $\mathit{\phi}\left(\mathit{x};\phantom{\rule{0.125em}{0ex}}\mathit{\lambda},\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{E}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{ABK}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{KIL}}\right)$ by implementing the Gauss–Newton algorithm with the initial value Σ^{(0)}=log (f^{(0)}), where f^{(0)}=10^{7} [m^{−2} s^{−1} eV^{−1}].
The hyperparameters were determined using the fivefold crossvalidation method (Stone, 1974). First, elements of the vector $\stackrel{\mathrm{\u0303}}{\mathit{g}}$ were divided into five subsets. Then, one was selected as a test set (${\stackrel{\mathrm{\u0303}}{\mathit{g}}}_{j}^{\mathrm{tes}})$ and the others as a training set (${\stackrel{\mathrm{\u0303}}{\mathit{g}}}_{j}^{\mathrm{tra}})$. We found the solution $\widehat{x}$ to minimize $\mathit{\phi}\left(\mathit{x};\phantom{\rule{0.125em}{0ex}}\mathit{\lambda},\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{E}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{ABK}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{KIL}}\right)$ using only the training set ${\stackrel{\mathrm{\u0303}}{\mathit{g}}}_{j}^{\mathrm{tra}}$ and then predicted the test set ${\mathit{g}}_{j}^{\mathrm{tes}}\left(\widehat{x}\right)$. We then calculated the sum of squares of the residuals between the test data and predicted
The crossvalidation score $\stackrel{\mathrm{\u203e}}{\mathit{\delta}}\left(\mathit{\lambda},\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{E}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{ABK}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{KIL}}\right)$ was calculated by averaging over five values of $\mathit{\delta}\left(\mathit{\lambda},\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{E}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{ABK}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{KIL}}\right)$, which were obtained by replacing the test set with one of the training sets in turn. The hyperparameters λ, λ_{E}, c_{ABK}, and c_{KIL} were determined by minimizing $\stackrel{\mathrm{\u203e}}{\mathit{\delta}}\left(\mathit{\lambda},\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{\mathrm{E}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{ABK}},\phantom{\rule{0.125em}{0ex}}{c}_{\mathrm{KIL}}\right)$ with a trialanderror method. In addition, the number of iterations for the Gauss–Newton algorithm was also simultaneously determined to be 200 to minimize $\stackrel{\mathrm{\u203e}}{\mathit{\delta}}$.
The PsA patches shown in Fig. 1a are embedded in the background diffuse auroral emission. We found that a horizontally uniform diffuse aurora causes ambiguity in the reconstruction result because the altitude of the uniform auroral structure cannot be determined from the singlewavelength images. Thus, we subtracted the background emission from the images prior to ACT reconstruction. We created the background emission image by assuming the same value for all voxels. The background VER was taken to be 75 cm^{−3} s^{−1}, corresponding to the spatially averaged observed background emission intensity. In this analysis, we assumed that the diffuse aurora showed a uniform VER in all voxels; however, the VER of the diffuse aurora depends generally on the altitude. We note that the altitude dependence of the VER did not affect the analysis result. This is because the horizontal distribution of the background emission intensity (i.e., diffuse aurora) in the auroral image was dependent only on the zenith angle θ (∝cos θ). This distribution is not dependent on the altitude distribution of the VER if the VER is horizontally uniform.
3.1 Reconstruction of a pulsating aurora patch model
We reconstructed a PsA patch model from pseudoauroral images to evaluate the analytical error of ACT before reconstructing the PsA patch from the observed auroral images. To create the pseudoauroral images, we prepared the horizontal distributions of the total energy (Q_{0}) and the characteristic energy (E_{c}). We then derived the 3D VER (L), as shown in Fig. 2a. The total energy was assumed to have a Gaussian shape in the horizontal directions with a maximum value of 6 mW m^{−2}. The energy distribution was considered to be a Maxwellian distribution with a uniform characteristic energy of 15 keV. Pseudoauroral images were obtained from L by solving the forward problem (Fig. 2b). We added random noise from a normal distribution with a mean value of 0 and the standard deviation determined from observed auroral images.
Figure 2c shows Q_{0}, E_{c}, and L reconstructed from the pseudoauroral images. The values of Q_{0} were calculated as ${Q}_{\mathrm{0}}={\sum}_{i}{E}_{i}f\left({E}_{i}\right)\left({E}_{i+\mathrm{1}}{E}_{i}\right)$. When we assume the energy distribution to be a Maxwellian distribution, the characteristic energy can be written as ${E}_{\mathrm{c}}=\frac{\mathrm{1}}{\mathrm{2}}\langle E\rangle =\frac{\mathrm{1}}{\mathrm{2}}\frac{{Q}_{\mathrm{0}}}{{\sum}_{i}f\left({E}_{i}\right)\left({E}_{i+\mathrm{1}}{E}_{i}\right)}$. We calculated the errors between the model and the result for Q_{0}, E_{c}, and L (Fig. 2d). The median values of the errors were −5 % for Q_{0}, −21 % for E_{c}, and −11 % for L. The northwestern part of Q_{0} was overestimated by at most 23 %, the edge part was underestimated by at most 29 %, and the central part was underestimated by ∼ 8 %. The central part of E_{c} was reconstructed correctly. In comparison, the edge part (especially the northwestern part) was underestimated by at most 56 %. The underestimation of E_{c} was caused by the overestimation of the emission altitude (Fig. 2d). The information regarding the PsA emission altitude is easily lost in obtaining the auroral image, as the structure of the PsA patch is vertically thin and horizontally wide. In addition, the SNR in the edge part is lower than in the central part because we assumed a Gaussian shape for the horizontal distribution of Q_{0}. These factors would tend to reduce the accuracy in the edge part. In this particular event, the optical observation points of ABK, SKB, and KIL were located in the southeastern direction from the PsA patch at the EISCAT radar observation point (TRO) (see Fig. 1). Therefore, an accurate reconstruction of the VER peak altitude was more difficult in the northwestern part than in other parts. This was because of the insufficient information on the northwestern part of the PsA patch owing to the bias in the optical observation point distribution.
Notably, the reconstructed results using the hyperparameters determined by the crossvalidation method revealed unexpected fine structures. To avoid this phenomenon, we set the lower limit of λ using a different method, namely by minimizing the sum of squares of the residuals between the model and the reconstructed result of Q_{0} and E_{c}. The lower limit on λ makes it challenging to reconstruct actual finescale structures in the patches.
3.2 Precipitating electrons
Figure 3a and b show Q_{0} and E_{c} as reconstructed from the observed auroral images (Fig. 1a). The maximum value of Q_{0} was ∼ 6 mW m^{−2}. The reconstructed E_{c} ranged from 6 to 23 keV. These energies are consistent with observation results from sounding rockets and lowaltitude satellites (e.g., McEwen et al., 1981; Miyoshi et al., 2015). The horizontal distribution of E_{c} was neither uniform nor stable in the patch during the pulsation. Thus, the ACT method is helpful for investigating PsAassociated temporal variations in the horizontal distribution of precipitating electrons without rocket or satellite observations. In particular, the southwestern part of E_{c} was enhanced at 00:53:38 UT. It should be noted that the edge and in the northwestern parts of E_{c} are expected to be underestimated owing to analytical error, as shown in Fig. 2d. These temporal variations indicate changes in the cyclotron resonance energy of whistlermode chorus waves during the pulsation. The chorus waves scatter electrons into a loss cone near the magnetic equator. The cyclotron resonance energy of chorus waves depends on the background magnetic field, electron density, and wave frequency (e.g., Kennel and Petschek, 1966). Therefore, the observed temporal variations indicate changes in the magnetospheric source region's background magnetic or plasma environment during the pulsation.
The characteristic energy can also be enhanced by the fieldaligned current (FAC). The absence of the FAC in the PsA patch has been observed because the PsA patch has no shear motion and the energy spectra of precipitating electrons observed by rockets and satellites did not show a fieldaligned acceleration, which are typically observed in discrete aurora (Davis, 1978). On the other hand, several studies have reported the FAC associated with PsA patches (Fujii et al., 1985; Gillies et al., 2015; Hosokawa and Ogawa, 2010). If the upward and downward FACs flow at the edge of the PsA path, the potential drop due to the upward FAC can accelerate precipitating electrons and enhance their characteristic energy (Sato et al., 2004; Shepherd and Fälthammar, 1980). However, we cannot suggest the existence of the FAC at the PsA patch's edge for this event because the total energy was not enhanced at the edge (Fig. 3). If the FAC exists, both total and characteristic energy should be increased. Multievent analysis is needed to examine the FAC in PsA patches. The 3D current structure in the PsA patch is beyond the scope of this study. Its reconstruction using ACT and EISCAT_3D radar is planned in the future.
3.3 Volume emission rate
The 3D distributions of the VERs derived from the reconstructed electron flux were obtained by solving the forward problem described in Sect. 2 (see Fig. 4a). Cross sections in the horizontal plane at an altitude of 94 km are shown in Fig. 4b. The peak altitude ranges from 90 to 104 km (Fig. 4c). The full width at half maximum is almost uniform with a median value of ∼ 20 km (Fig. 4e). The errors of the peak altitude and the altitude width in Fig. 4d and f were derived using the model and reconstructed VER shown in Fig. 2. The high peak altitude in the northwestern part is expected to be overestimated by at most 8 % owing to analytical error, as shown in Fig. 2d. For the altitude width, the most part is expected to be overestimated by ∼ 2 % (Fig. 4f). The reconstructed peak altitude and width are consistent with those determined in studies using stereoscopic observations or an incoherent scatter radar (Brown et al., 1976; Jones et al., 2009; Kataoka et al., 2016). StenbaekNielsen and Hallinan (1979) reported the existence of thin (<1 km vertical extent) PsA patches based on stereoscopic observations, but our results do not support their results.
The peak altitude of the PsA patch was also estimated using a different method (Appendix A). We projected the observed auroral images at altitudes ranging from 80 to 120 km with an interval of 2 km (Video A1, https://doi.org/10.5446/57558). The emission altitude was determined to be the altitude at which the sum of squares of the residuals between the two projected images reached a minimum value (Fig. A1). The estimated peak altitude range was 92 to 106 km from 00:53:30 to 00:53:40 UT (Fig. A2). These altitudes closely match those determined by ACT.
3.4 Electron density
The altitude profiles of VER at the EISCAT radar observation point shown in Fig. 4g were converted to the ionospheric electron density and compared with the actual data observed by the EISCAT radar. The continuity equation for the electron density can be written as
where n_{e} [m^{−3}] is the electron density, L [m^{−3} s^{−1}] is the VER, k is a positive constant for converting VER to the ionization rate (see Appendix B), and α_{eff} [m^{3} s^{−1}] is the effective recombination rate. We derived the electron density from the VER by solving Eq. (4) with the Runge–Kutta method. The initial value was derived from Eq. (4) under steadystate conditions (i.e., $\partial {n}_{\mathrm{e}}/\partial t=\mathrm{0})$ using reconstructed L at 00:53:36 UT. The VERs were interpolated linearly in time to use the Runge–Kutta method. The altitude profile of α_{eff} has been investigated by several studies using rocket and groundbased measurements. Vickrey et al. (1982) summarized many of these results and proposed the following best fit parameterization:
where z [km] is the altitude. Semeter and Kamalabadi (2005) used the effective recombination coefficients ${\mathit{\alpha}}_{{\mathrm{NO}}^{+}}$ and ${\mathit{\alpha}}_{{\mathrm{O}}_{\mathrm{2}}^{+}}$ for NO^{+} and O${}_{\mathrm{2}}^{+}$, respectively (Walls and Dunn, 1974), as upper and lower bounds on α_{eff}:
Here, T_{n} [K] is the neutral temperature. The red lines in Fig. 5 show the derived electron densities using these three recombination coefficients. We note that these values are underestimated compared with the electron densities observed by the EISCAT radar (black lines in Fig. 5). This underestimation probably comes from the background emission subtraction from the auroral images prior to ACT and from ambiguity in the effective recombination coefficients. The electron densities reconstructed from auroral images without background emission subtraction are shown as blue lines for reference in Fig. 5. The reconstruction results from the images, including background emission, approached the electron density profile observed with the EISCAT radar. We note that the electron density was reconstructed correctly to some extent, even after the auroral emission intensity decreased at 00:53:40 UT. This correct reconstruction considered the time change in the continuity equation. The electron density would seem to have rapidly decreased after 00:53:40 UT if the time change was not considered. This result suggested that the time change should be considered ($\mathrm{d}{n}_{\mathrm{e}}/\mathrm{d}t\ne \mathrm{0}$ in Eq. 4) when using the continuity equation to derive the electron densities associated with PsAs.
It should be noted that the electron density is underestimated at higher altitudes (> ∼ 140 km), even if the background emission is included. This underestimation can be improved by reconstructing lowenergy electron flux from auroral images of various wavelengths (e.g., 844.6 nm).
We applied the ACT method to a PsA patch for the first time and reconstructed the horizontal distribution of precipitating electrons from 427.8 nm auroral images obtained from three observation points. We improved the previously proposed ACT method by adding the following processes: the subtraction of the background diffuse aurora from the auroral images prior to ACT, the estimation of the relative sensitivity between ASCs, and the determination of the hyperparameters of the regularization term. The characteristic energies of the reconstructed electron fluxes (6 to 23 keV) and the peak altitudes of the reconstructed VERs (90 to 104 km) were consistent with those found in previous studies. We determined that the horizontal distribution of the characteristic energy was neither uniform nor stable in the patch during the pulsation, further underlining the shortcomings of rocket and satellite observations with respect to investigating PsAs. These spatial and temporal variations imply changes in the electron density and magnetic field in the magnetosphere. ACT error was evaluated using an auroral patch model. The characteristic energy of electron flux was correctly reconstructed in the center part of the patch but underestimated at the patch edge by at most 56 %. The reconstructed electron flux will be improved in future work by incorporating auroral images of various wavelengths.
Although we reconstructed the differential flux of precipitating electrons from auroral images using ACT, Tanaka et al. (2011) extended ACT to a method called generalizedaurora computed tomography (GACT). GACT uses multiinstrument data, such as ionospheric electron density from incoherent scatter radar, cosmic noise absorption from imaging riometers, and the auroral images. Tanaka et al. (2011) demonstrated that the incorporation of the ionospheric electron density from the EISCAT radar improved the accuracy of the reconstructed electron flux. Furthermore, 3D ionospheric observation by EISCAT_3D (http://www.eiscat3d.se, last access: 7 July 2022) is scheduled to begin in 2023. In the future, we will improve the reconstructed electron flux by conducting GACT analysis using electron density data from the EISCAT or EISCAT_3D radar.
Here, we estimate the peak altitude of the PsA patch using a different method to validate the results from ACT in Sect. 3.3. We projected the observed auroral images at altitudes from 80 to 120 km with an interval of 2 km. As an example, projected images for 00:53:36 UTC on 18 February 2018 are shown in Video A1 (https://doi.org/10.5446/57558). The emission altitude was determined to be the altitude at which the sum of squares of the residuals of the auroral intensity between the two projected images reached a minimum value (Fig. A1). The estimated peak altitudes from 00:53:30 to 00:53:40 UTC are shown in Fig. A2. These altitudes agreed with the results from ACT in Sect. 3.3.
In this section, we describe how to obtain the positive constant k(z) from Sect. 3.4. The ${\mathrm{N}}_{\mathrm{2}}^{+}$ (427.8 nm) emission is owing to the transition from ${\mathrm{N}}_{\mathrm{2}}^{+}\phantom{\rule{0.125em}{0ex}}{\left({B}^{\mathrm{2}}{\mathbf{\Sigma}}_{u}^{+}\right)}_{v=\mathrm{0}}$ to ${\mathrm{N}}_{\mathrm{2}}^{+}\phantom{\rule{0.125em}{0ex}}{\left({X}^{\mathrm{2}}{\mathbf{\Sigma}}_{g}^{+}\right)}_{v=\mathrm{1}}$. According to Sergienko and Ivanov (1993), the VER L(z) [m^{−3 }s^{−1}] is approximated by
where A_{0−1} is the Einstein coefficient for the transition from ${\mathrm{N}}_{\mathrm{2}}^{+}\phantom{\rule{0.125em}{0ex}}{\left({B}^{\mathrm{2}}{\mathbf{\Sigma}}_{u}^{+}\right)}_{v=\mathrm{0}}$ to ${\mathrm{N}}_{\mathrm{2}}^{+}\phantom{\rule{0.125em}{0ex}}{\left({X}^{\mathrm{2}}{\mathbf{\Sigma}}_{g}^{+}\right)}_{v=\mathrm{1}}$, w(z) [m^{−3} s^{−1}] is the production rate of ${\mathrm{N}}_{\mathrm{2}}^{+}\phantom{\rule{0.125em}{0ex}}\left({B}^{\mathrm{2}}{\mathbf{\Sigma}}_{u}^{+}\right)$, q_{0−0} is the Franck–Condon factor for the electronic transition from ${\mathrm{N}}_{\mathrm{2}}^{+}\phantom{\rule{0.125em}{0ex}}{\left({X}^{\mathrm{1}}{\mathbf{\Sigma}}_{g}^{+}\right)}_{v=\mathrm{0}}$ to ${\mathrm{N}}_{\mathrm{2}}^{+}\phantom{\rule{0.125em}{0ex}}{\left({B}^{\mathrm{2}}{\mathbf{\Sigma}}_{u}^{+}\right)}_{v=\mathrm{0}}$, p(z) is the probability that ε(z) excites N_{2}, ε(z) [eV m^{−3} s^{−1}] is the energy deposition rate, and Δε [eV] is the excitation energy cost of ${\mathrm{N}}_{\mathrm{2}}^{+}\phantom{\rule{0.125em}{0ex}}\left({B}^{\mathrm{2}}{\mathbf{\Sigma}}_{u}^{+}\right)$. The ionization rate due to the precipitating electrons q_{ion}(z) [m^{−3} s^{−1}] is given by
where Δε_{ion} [eV] is the energy used to produce an ion–electron pair. Substituting Eq. (B1) into Eq. (B2) gives
Therefore, the positive constant k(z) for converting VER to the ionization rate is
The parameters used for the calculation are summarized in Table B1.
The MIRACLE EMCCD camera data from ABK and KIL are available from the University of Oulu and the Finnish Meteorological Institute (2022, https://doi.org/10.23729/e86df44b8dad44e589f34bea0d3d1236, Raita and Kauristie, 2022). The auroral images obtained by four WMI CCD cameras at SKB can be obtained at http://pc115.seg20.nipr.ac.jp/www/optical/watec/skb/awi/rawdata/ (Ogawa, 2022a). The EISCAT data are available at http://pc115.seg20.nipr.ac.jp/www/eiscatdata/ (Ogawa, 2022b).
Video A1 is available at https://doi.org/10.5446/57558 (Fukizawa, 2022).
YT developed the aurora computed tomography method and code. YO conducted the EISCAT radar observation and prepared the ionospheric electron density data. KK, AK, and TR maintained the MIRACLE camera network and prepared the auroral images. MF analyzed the data prepared by the coauthors and prepared the manuscript with contributions from all coauthors. TaS, KH, BG, UB, and TiS contributed to the discussion and interpretation of the analysis results.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The first author is a research fellow (DC) of the Japan Society for the Promotion of Science (JSPS). This work was performed using facilities of the National Institute of Polar Research (NIPR). The study has been supported by JSPS KAKENHI (grant nos. JP17K05672, JP20J11829, and JP21H01152). EISCAT is an international association supported by research organizations in China (CRIRP), Finland (SA), Japan (NIPR), Norway (NFR), Sweden (VR), and the United Kingdom (UKRI). We thank Kellinsalmi Mirjam and CarlFredrik Enell for maintaining the MIRACLE camera network and data flow. The database construction for the imager data at Skibotn and the EISCAT radar data has been supported by the IUGONET (Interuniversity Upper atmosphere Global Observation NETwork) project (http://www.iugonet.org/, last access: 7 July 2022).
This research has been supported by the Japan Society for the Promotion of Science (grant nos. JP17K05672, JP20J11829, and JP21H01152).
This paper was edited by Dalia Buresova and reviewed by two anonymous referees.
Aso, T., Hashimoto, T., Abe, M., Ono, T., and Ejiri, M.: On the analysis of aurora stereo observations, J. Geomagn. Geoelectr., 42, 579–595, https://doi.org/10.5636/jgg.42.579, 1990.
Aso, T., Ejiri, M., Miyaoka, H., Hashimoto, T., Abu, T. Y., and Abe, M.: Aurora stereo observation in Iceland, Proc. NIPR Symp. Up. Atmos. Phys., 6, 1–14, 1993.
Aso, T., Ejiri, M., Urashima, A., Miyaoka, H., Steen, Å., Brändstrom, U., and Gustavsson, B.: First results of auroral tomography from ALISJapan multistation observations in March, 1995, Earth Planet. Space, 50, 81–86, https://doi.org/10.1186/BF03352088, 1998.
Brown, N. B., Davis, T. N., Hallinan, T. J., and StenbaekNielsen, H. C.: Altitude of pulsating aurora determined by a new instrumental thechnique, Geophys. Res. Lett., 3, 403–404, https://doi.org/10.1029/GL003i007p00403, 1976.
Davis, T. N.: Observed characteristics of auroral forms, Space Sci. Rev., 22, 77–113, https://doi.org/10.1007/BF00215814, 1978.
Frey, S., Frey, H. U., Carr, D. J., Bauer, O. H., and Haerendel, G.: Auroral emission profiles extracted from threedimensionally reconstructed arcs, J. Geophys. Res.Space, 101, 21731–21741, https://doi.org/10.1029/96ja01899, 1996.
Fujii, R., Oguti, T., and Yamamoto, T.: Relationships between pulsating auroras and fieldaligned electric currents, Mem. Natl. Inst. Polar Res. Spec. Issue, 36, 95–1003, 1985.
Fukizawa, M.: Movie A1, Logo TIB AVPortal [video], https://doi.org/10.5446/57558, 2022.
Gillies, D. M., Knudsen, D., Spanswick, E., Donovan, E., Burchill, J., and Patrick, M.: Swarm observations of fieldaligned currents associated with pulsating auroral patches, J. Geophys. Res. A, 120, 9484–9499, https://doi.org/10.1002/2015JA021416, 2015.
Hosokawa, K. and Ogawa, Y.: Pedersen current carried by electrons in auroral Dregion, Geophys. Res. Lett., 37, 1–5, https://doi.org/10.1029/2010GL044746, 2010.
Jones, A. V.: Aurora, D. Reidel Publishing Company, Dordrecht, Springer, ISBN 9789401020992, 1974.
Jones, S. L., Lessard, M. R., Fernandes, P. A., Lummerzheim, D., Semeter, J. L., Heinselman, C. J., Lynch, K. A., Michell, R. G., Kintner, P. M., StenbaekNielsen, H. C., and Asamura, K.: PFISR and ROPA observations of pulsating aurora, J. Atmos. Sol.Terr. Phys., 71, 708–716, https://doi.org/10.1016/j.jastp.2008.10.004, 2009.
Kataoka, R., Fukuda, Y., Uchida, H. A., Yamada, H., Miyoshi, Y., Ebihara, Y., Dahlgren, H., and Hampton, D.: Highspeed stereoscopy of aurora, Ann. Geophys., 34, 41–44, https://doi.org/10.5194/angeo34412016, 2016.
Kennel, C. F. and Petschek, H. E.: Limit on stably trapped particle fluxes, J. Geophys. Res., 71, 1–28, https://doi.org/10.1029/JZ071i001p00001, 1966.
McEwen, D. J., Yee, E., Whalen, B. A., and Yau, A. W.: Electron energy measurements in pulsating auroras, Can. J. Phys., 59, 1106–1115, https://doi.org/10.1139/p81146, 1981.
Miyoshi, Y., Saito, S., Seki, K., Nishiyama, T., Kataoka, R., Asamura, K., Katoh, Y., Ebihara, Y., Sakanoi, T., Hirahara, M., Oyama, S., Kurita, S., and Santolik, O.: Relation between energy spectra of pulsating aurora electrons and frequency spectra of whistlermode chorus waves, J. Geophys. Res.Space, 120, 7728–7736, https://doi.org/10.1002/2015JA021562, 2015.
Nygrén, T., Markkanen, M., Lehtinen, M., and Kaila, K.: Application of stochastic inversion in auroral tomography, Ann. Geophys., 14, 1124–1133, https://doi.org/10.1007/s0058599611241, 1997.
Ogawa, Y.: Watec observation database in NIPR [data set], http://pc115.seg20.nipr.ac.jp/www/optical/watec/skb/awi/rawdata/, last access: 6 June 2022a.
Ogawa, Y.: EISCAT database in NIPR, EISCAT [data set], http://pc115.seg20.nipr.ac.jp/www/eiscatdata/, last access: 27 February, 2022b.
Ogawa, Y., Tanaka, Y., Kadokura, A., Hosokawa, K., Ebihara, Y., Motoba, T., Gustavsson, B., Brändström, U., Sato, Y., Oyama, S., Ozaki, M., Raita, T., Sigernes, F., Nozawa, S., Shiokawa, K., Kosch, M., Kauristie, K., Hall, C., Suzuki, S., Miyoshi, Y., Gerrard, A., Miyaoka, H., and Fujii, R.: Development of lowcost multiwavelength imager system for studies of aurora and airglow, Polar Sci., 23, 100501, https://doi.org/10.1016/j.polar.2019.100501, 2020a.
Ogawa, Y., Kadokura, A., and Ejiri, M. K.: Optical calibration system of NIPR for aurora and airglow observations, Polar Sci., 26, 100570, https://doi.org/10.1016/j.polar.2020.100570, 2020b.
Picone, J. M., Hedin, A. E., Drob, D. P., and Aikin, A. C.: NRLMSISE00 empirical model of the atmosphere: Statistical comparisons and scientific issues, J. Geophys. Res.Space, 107, 1468, https://doi.org/10.1029/2002JA009430, 2002.
Raita, T. and Kauristie, K.: Subset of the MIRACLE emCCD allsky camera data at 427.8 nm from ABK and KIL on 18 Feb 2018 00–02 UT, Etsin [data set], https://doi.org/10.23729/e86df44b8dad44e589f34bea0d3d1236, 2022.
Sangalli, L., Partamies, N., Syrjäsuo, M., Enell, C. F., Kauristie, K., and Mäkinen, S.: Performance study of the new EMCCDbased allsky cameras for auroral imaging, Int. J. Remote Sens., 32, 2987–3003, https://doi.org/10.1080/01431161.2010.541505, 2011.
Sato, N., Wright, D. M., Carlson, C. W., Ebihara, Y., Sato, M., Saemundsson, T., Milan, S. E., and Lester, M.: Generation region of pulsating aurora obtained simultaneously by the FAST satellite and a SyowaIceland conjugate pair of observatories, J. Geophys. Res.Space, 109, 1–15, https://doi.org/10.1029/2004JA010419, 2004.
Semeter, J. and Kamalabadi, F.: Determination of primary electron spectra from incoherent scatter radar measurements of the auroral e region, Radio Sci., 40, RS2006, https://doi.org/10.1029/2004RS003042, 2005.
Sergienko, T. and Ivanov, V.: A new approach to calculate the excitation of atmospheric gases by auroral electron impact, Ann. Geophys. Eur. Geophys. Soc., 11, 717–727, 1993.
Shepherd, G. G. and Fälthammar, C.G.: Implications of extreme thinness of pulsating auroral structures, J. Geophys. Res.Space, 85, 217–218, https://doi.org/10.1029/ja085ia01p00217, 1980.
StenbaekNielsen, H. C. and Hallinan, T. J.: Pulsating auroras: Evidence for noncollisional thermalization of precipitating electrons, J. Geophys. Res., 84, 3257–3271, https://doi.org/10.1029/JA084iA07p03257, 1979.
Stone, M.: Crossvalidatory choice and assessment of statistical predictions (with discussion), J. R. Stat. Soc. Ser. B, 38, 102–102, https://doi.org/10.1111/j.25176161.1976.tb01573.x, 1974.
Tanaka, Y. M., Aso, T., Gustavsson, B., Tanabe, K., Ogawa, Y., Kadokura, A., Miyaoka, H., Sergienko, T., Brändström, U., and Sandahl, I.: Feasibility study on GeneralizedAurora Computed Tomography, Ann. Geophys., 29, 551–562, https://doi.org/10.5194/angeo295512011, 2011.
Vickrey, J. F., Vondrak, R. R., and Matthews, S. J.: Energy deposition by precipitating particles and Joule dissipation in the auroral ionosphere, J. Geophys. Res.Space, 87, 5184–5196, https://doi.org/10.1029/ja087ia07p05184, 1982.
Walls, F. L. and Dunn, G. H.: Measurement of total cross sections for electron recombination with NO^{+} and O${}_{\mathrm{2}}^{+}$ using ion storage techniques, J. Geophys. Res., 79, 1911–1915, https://doi.org/10.1029/ja079i013p01911, 1974.
Yamamoto, T.: On the temporal fluctuations of pulsating auroral luminosity, J. Geophys. Res., 93, 897–911, https://doi.org/10.1029/JA093iA02p00897, 1988.
 Abstract
 Introduction
 Data and methods
 Results and discussion
 Conclusions
 Appendix A: Estimation of the pulsating auroral emission peak altitude
 Appendix B: Derivation of k
 Data availability
 Video supplement
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.
 Article
(4876 KB)  Fulltext XML
 Abstract
 Introduction
 Data and methods
 Results and discussion
 Conclusions
 Appendix A: Estimation of the pulsating auroral emission peak altitude
 Appendix B: Derivation of k
 Data availability
 Video supplement
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References