![pageSearch](/themes/hestia/images/page-search.webp)
Paige DeSorbo’s Chic 30th Birthday Included a Crop Top, Maxi Skirt, and Pink Balloons
Apr 27, 2023Best Fanny Packs for Men and Women 2023
Jul 02, 2023Five Charles Oliveira Fights To Watch Before UFC 289
May 24, 2023Road to Double Gold: Pessanha Becomes First Woman With All Subs At Worlds
Jun 09, 2023Seat belt violation leads to four felony arrests for drugs
Jul 18, 2023Localized mesospheric ozone destruction corresponding to isolated proton aurora coming from Earth’s radiation belt
![May 04, 2023](/themes/hestia/images/news-details-icon1.webp)
Scientific Reports volume 12, Article number: 16300 (2022) Cite this article
3007 Accesses
1 Citations
178 Altmetric
Metrics details
Relativistic electron precipitation (REP) from the Earth's radiation belt plays an important role in mesospheric ozone loss as a connection between space weather and the climate system. However, the rapid (tens of minutes) destruction of mesospheric ozone directly caused by REP has remained poorly understood due to the difficulty of recognizing its location and duration. Here we show a compelling rapid correspondence between localized REP and ozone destruction during a specific auroral phenomenon, the called an isolated proton aurora (IPA). The IPA from the Earth's radiation belt becomes an important spatial and temporal proxy of REP, distinct from other auroral phenomena, and allowing visualizing micro-ozone holes. We found ozone destruction of as much as 10–60% within 1.5 h of the initiation of IPA. Electromagnetic ion cyclotron waves in the oxygen ion band observed as the driver of REP likely affect through resonance with mainly ultra-relativistic (> 2 mega-electron-volts) energy electrons. The rapid REP impact demonstrates its crucial role and direct effect on regulating the atmospheric chemical balance.
Atmospheric effects due to energetic particle precipitation (EPP) at hundreds of kilo-electron-volts (keV) to mega-electron-volts (MeV) have significant atmospheric chemical impacts in the mesosphere (50 to 80 km) and the upper stratosphere (~ 50 km) below the thermosphere (100 to 200 km)1,2,3,4. EPP is one of the major sources of the catalytic destruction of ozone in the polar region (magnetic latitudes > 55°) due to the production of EPP-driven odd nitrogen (NOx) and odd hydrogen (HOx)5,6,7,8. Plasma particle energy is important for determining the altitude of ionization in the atmosphere9. EPP-NOx interactions affect the vertical transport of EPP-driven NOx in the lower thermosphere (its generation altitude of ~ 80 km) down to stratospheric altitude as a so-called indirect effects10. As a result, the polar vortex effectively transports NOx in the polar region, then NOx transport plays a significant role in mesospheric ozone loss by 10 to 20% on a time scale of several months to decades11. In contrast, the local production of the EPP-driven NOx and HOx directly contributes to ozone destruction at their production altitudes9,12. EPP-HOx interactions, particularly, can rapidly occur with a shorter time period of hours, because of the lifetime of a few hours of the HOx family4. The ozone density in the mesosphere is much smaller than that in the stratosphere, but the mesospheric ozone and atmospheric ionizations may play an important role for the global climate system via the chemistry and transport processes13,14. Solar proton events, which bring a strong enhancement in energetic proton flux (> 10 MeV) caused by an active solar eruption, are a major source of EPP, and such solar proton precipitation plays an important role in ozone destruction5,6,7 and abnormal electron density enhancement in the whole polar cap region15. Given that effects from solar proton events are global, we are prompted to ask a new question: can EPP direct impacts on atmospheric chemistry be observed in a localized area with a short duration? In principle, EPP can show clear localization like aurora at specific latitudes, longitudes, and appearance periods. Thus, if EPP directly and rapidly acts as a major driver for the mesospheric ozone destruction as the same with solar proton events, localized ozone destruction associated with localized EPP should be observable. Simulation studies of the direct effect in for specific EPP events (pulsating aurora (> 200-keV electrons)16 and relativistic electron microbursts (> 1-MeV electrons)17) have predicted up to ~ 20% ozone destruction in the mesosphere, which can be equivalent to the effects created by solar proton events18. Such simulation studies predict ozone loss associated with spatially and temporally localized EPP events. However, to date, no observational evidence has verified such localized ozone loss, because identifying any localized, short-duration EPP events in observations is difficult.
In this study, we overcome this difficulty using a specific auroral phenomenon, the so-called isolated proton aurora (IPA) associated with Pc1 geomagnetic pulsations19. Most auroras seen at magnetic latitudes usually from 65° to 70° (i.e., in the "auroral oval") are caused by the particle precipitation of low-energy electrons (keV to tens of keV range), which cannot penetrate to the mesosphere, but the IPA at subauroral latitudes (magnetic latitudes from 55° to 60°), although mostly energized by protons, is accompanied by relativistic electron precipitation (REP) from the outer Earth's radiation belt20,21, which is a torus-shaped zone of relativistic energetic electrons around the Earth's magnetosphere22. The IPA is caused by cyclotron resonant interactions of ring current (tens to a few hundred keV) protons23,24 and radiation belt (1- to 10-MeV) electrons25,26 with electromagnetic ion cyclotron (EMIC) waves27,28,29, which are coherent electromagnetic wave emissions in the Earth's magnetosphere and are observed as Pc1 geomagnetic pulsations of frequency a few Hz on the ground. The precipitating protons cause the IPA to emit with characteristic hydrogen optical emissions at subauroral latitudes, making it a proxy for the spatial distribution of REP21. The clear temporal correlation between IPA and EMIC/Pc1 waves is further useful as a proxy for the temporal localization of REP30,31. Our results demonstrate that EMIC-driven IPA visualizes the localization of the mesospheric ozone loss, like a "needle hole" in the ozone layer, due to precipitation of the Earth's radiation belt electrons. Quantifying the direct effects of EMIC-driven REP on the atmospheric chemistry and confirming that it accounts for ozone accumulation in such a needle hole contributes to a better understanding of the Sun-magnetosphere-climate connection.
Figure 1 shows a schematic illustration of localized mesospheric ozone loss associated with an IPA. The IPA optical emission serves as an ionospheric indicator of the concentrated downpour of REP. The MeV electrons and lower energy ring current protons are scattered through EMIC wave-particle interactions at the magnetic equator32 so that some of them precipitate following their geomagnetic field lines; then the IPA from the protons and accompanying ozone loss due to the MeV electron precipitation (i.e., REP) can be seen at the same location at different altitudes. We utilize innovative satellite remote sensing for identifying IPA and ozone profiles (see "Methods"). The satellite remote sensing of IPA and ozone can allow large spatial coverage regardless of surface weather conditions. The IPA at subauroral latitudes is separated from the equatorward boundary of the auroral oval. In satellite-based limb observations of ozone, this separation is enough to distinguish IPA effects from those of other auroral phenomena near an ozone observation tangent point, which is a point where the line of sight of a sounding ray is closest to the Earth's surface. The related EMIC waves in the magnetosphere propagate to the ground along the magnetic field line. The EMIC/Pc1 waves are observed at magnetic conjugate ground stations in both the northern and southern hemispheres33. The temporal information of IPA is supplemented by a ground observation of related EMIC/Pc1 waves at a fixed ground station. Here we show two prime examples of localized ozone loss related to the IPA.
Effects of EMIC wave-particle interactions at different altitudes. Both EMIC-driven protons (about tens of kilo-electron-volts) and relativistic electrons precipitate into the upper atmosphere along the same geomagnetic field lines from the magnetosphere. The localized IPA optical emission due to the protons and the mesospheric ozone destruction by relativistic electron precipitation (REP) can be seen on the same field line at different altitudes. The Earth image is obtained from NASA, https://earthobservatory.nasa.gov/features/BlueMarble/BlueMarble_2002.php. This illustration created with Adobe Illustrator version 26.2.1 (https://www.adobe.com/).
Figure 2a and b show a coordinated observation of IPA [Lyman alpha (121.6 nm) hydrogen emissions] by the Special Sensor Ultraviolet Spectrographic Imager (SSUSI)34 on board the Defense Meteorological Satellite Program (DMSP) satellites and the mesospheric (altitude of 52 km) ozone profiles (stars) by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER)35 on board the Thermosphere Ionosphere Mesosphere Energetics Dynamics (TIMED) satellite from 21:30 to 25:20 UT (13 to 17 magnetic local time (MLT) in the dusk sector), June 22, 2015, during the early main phase of a large geomagnetic storm (minimum Dst index = − 121 nT). A localized REP is seen, as opposed to the effects of solar proton precipitation, which would be observed in the whole polar cap region15. The related REP (triangles) is measured by the vertical component of the Radiation Belt Monitor (RBM-Z) of the Monitor of All-sky X-ray Image (MAXI)36 on board the International Space Station (ISS). The global image is plotted in dipole coordinates using the coefficients from the 13th International Geomagnetic Reference Field model37. The coordinated IPA and the ozone observation shows a spatial correspondence between the localized IPA and the ozone needle hole. The IPA is seen at a magnetic latitude of ~ 57°, far from the auroral oval boundary at around 62° north latitude. The spatial scale of the IPA is 390 km in the latitudinal direction and 1340 km in the longitudinal direction. The transverse size across the geomagnetic field line of the EMIC wave-particle interaction region is 0.3 to 0.5 Earth radii (Re) in the radial direction, as determined by mapping to the magnetic equator using an empirical geomagnetic field model38. The IPA and REP distributions at the magnetic equator were projected from the L = 4.3 to 4.8, MLT = 16 to 17 region. RBM-Z also shows a perfect spatial agreement between the REP and IPA in the longitudinal direction, reflecting the increasing count rate coming from measuring loss cone electrons as REP. Figure 2c shows the altitude profile of ozone measured by the SABER along tangent points of orbit 73,365 of the TIMED satellite just crossing the IPA location. Figure 2d shows the ozone concentration difference between the ozone profiles crossing the IPA and the monthly mean value, with the standard deviation. The monthly mean near the IPA location is calculated from the SABER observations in magnetic latitudes from 55°N to 60°N and longitudes from 15°W to 60°W during 22 to 25 UT between May 23 and June 21, 2015. A significant ozone loss during orbit 73,365 crossing the IPA location during event numbers 34 and 35 is identified at the altitude from 47 to 66 km in the upper stratosphere. A positive bump of ozone around an altitude of 58 km during event number 35 (blue curve) would be caused by the effects of vertical wind shear due to the localized REP39. The average ozone loss at the altitude from 47 to 66 km during event numbers 34 and 35 shows a decrease by 11% from the monthly mean value and by 17% from the ozone measurements in the vicinity of the IPA, but outside the IPA in the TIMED orbits 73,364 and 73,366 (see Supplementary Fig. S1). The related EMIC/Pc1 waves in the oxygen ion (O +) band were observed on the ground station at Athabasca (ATH, L = 4.5)40, Canada, as shown in Fig. 2e. The EMIC/Pc1 waves were generated at 22:05 UT, and then TIMED/SABER observed the ozone needle hole below the IPA at 23:36 UT during TIMED orbit 73,365. Thus, the ozone variation below the IPA can decrease by up to 17% more in 1.5 h after the generation of EMIC/Pc1 waves with respect to the surrounding region of the IPA.
Localized ozone loss associated with the EMIC-driven IPA during a geomagnetic storm. (a) Global map of the mesospheric ozone (SABER, stars), REP (RMB-Z, triangles \(\ge 90\) counts), and aurora (superposition acquired from 19:44 to 25:46 UT by DMSP18 and 16/SSUSI). (b) Magnified map (aurora taken at 23:45 UT by DMSP16/SSUSI). Black arrows indicate the orbit direction of the ISS and TIMED satellite, with the orbit number and the time information in parentheses giving the measurement time for TIMED/SABER. The black dot indicates the location of Athabasca (ATH), Canada. Maps were generated by MATLAB R2020b, https://www.mathworks.com. (c) Tangent altitude profile of ozone during orbit 73,365 through the IPA. (d) Difference in the mesospheric ozone profiles. "Needle hole" region highlighted. (e) EMIC/Pc1 wave activity on June 22, 2015, at ATH. The black dotted line is the O + gyrofrequency at the magnetic equator estimated from an empirical geomagnetic field model38. Red lines indicate the time each instrument passed the magnetic latitude of the IPA (~ 57°). Numbers in parentheses denote the orbit number of each satellite.
The crucial role of EMIC-driven REP in decreasing mesospheric ozone is further shown in Fig. 3 for a different, non-storm event during 20:30 to 24:45 UT (14 to 18 MLT in the dusk sector), on August 11, 2014. The IPA arc was distributed with a narrow latitude range of 360 km and a wide longitudinal width of 70° (Fig. 3a and b). The equatorial transverse size of the EMIC wave-particle interaction region is 0.8 to 1.0 Re in the radial direction using an empirical geomagnetic field model41. The IPA and REP distributions at the magnetic equator were projected from the L = 4.7 to 5.7, MLT = 20 to 22 region. Figure 3 has the same format as Fig. 2, but the REP was measured with the E4 channel of the Medium Energy Proton/Electron Detector (MEPED)42 aboard the Polar-orbiting Operational Environmental Satellite (POES) constellation, and the related EMIC/Pc1 waves in the O + band (Fig. 3e) were observed on the ground station at Fort Churchill (FCHU, L = 7.2)43, Canada, which is a geomagnetic conjugate point for the IPA observed in Antarctica. After weak EMIC/Pc1 waves were generated at 21:15 UT, POES 15 observed the related REP with a 40-degree wide longitudinal range at 21:24 UT, POES 18 at 23:32 UT, and MetOP02 at 23:57 UT. TIMED/SABER observed the ozone loss during orbits 68,682 and 68,683, which crossed the IPA at different longitudes. A dense mesospheric ozone cloud is seen around altitudes from 62 to 75 km, but a clear destruction of the ozone cloud is seen at event numbers 74 and 75 during the TIMED/SABER orbit 68,683, when it crossed the IPA region (Fig. 3c). In Fig. 3d, the monthly mean of ozone near the IPA location is calculated from the SABER observations at magnetic latitudes from 65°S to 70°S and longitudes from 0°E to 45°E during 21:30 to 24:30 UT between July 11 and August 9, 2014. The ozone loss at altitudes from 64 to 72 km observed during event numbers 74, 75, and 78 of TIMED/SABER orbits 68,683 and 68,682 decreased by 51% from the monthly mean value. The ozone loss in event numbers 74 and 75 of TIMED/SABER orbit 68,683 was 61% from the ozone measurements by event numbers 73 and 76 of the same orbit 68,683, but they are observed outside the IPA (see Supplementary Fig. S2). The ozone variation below the IPA is observed to more rapidly decrease ~ 1 h after the generation of EMIC/Pc1 waves compared with that in the surrounding region of the IPA. The clear spatial agreement between the IPA and the mesospheric ozone destruction strongly supports the idea that the EMIC-induced REP directly and rapidly affects the ozone decreases via direct ionization in the mesosphere.
Localized ozone loss associated with EMIC-driven IPA in the absence of a storm. (a) Global map of the mesospheric ozone (SABER, stars), REP (MEPED, triangles \(\ge 4\) counts), and aurora (superposition acquired from 22:18 to 27:32 UT by DMSP19/SSUSI). The black dot indicates the magnetic conjugate point of Fort Churchill (FCHU), Canada. (b) Magnified map (aurora taken at 24:07 UT by DMSP19/SSUSI). Black arrows indicate the orbit direction of the TIMED and POES satellites and the time information in parentheses is the measurement time for TIMED/SABER. Maps were generated by MATLAB R2020b, https://www.mathworks.com. (c) Tangent altitude profile of ozone by orbit 68,683 through the IPA. (d) Difference in the mesospheric ozone profiles. e EMIC/Pc1 wave activity on August 11, 2014, at FCHU. The black dotted line is the O + gyrofrequency at the magnetic equator estimated from an empirical geomagnetic field model41. Red lines indicate the time each instrument passed the same magnetic latitude of the IPA (~ 63°). Numbers in parentheses denote the orbit number of each satellite.
The frequency bands of EMIC/Pc1 waves are characterized by related ion species (e.g., proton (H +), Helium (He +), and Oxygen (O +) etc.), and the frequency band is important for determining the resonance plasma energy25. EMIC/Pc1 waves along with the IPA are usually observed in the helium ion (He +) band20,21,30,31, but EMIC/Pc1 waves in both events studied were in the O + band. We evaluated the equatorial pitch angle diffusion rates25 (see "Methods") and the minimum resonant energy of electrons by the observed EMIC/Pc1 waves to confirm the validation of the EMIC-driven REP. The minimum resonant energy \(E_{{{\text{min}}}}\) of electrons interacting with EMIC waves is written as25
where \({\Omega }_{e}\) is the electron gyrofrequency; \(\omega\) is the angular frequency of EMIC waves; n is the refractive index of EMIC waves, which depends on the number density of each ion; and \(E_{0}\) is the rest mass energy of an electron (0.511 MeV). Figure 4 shows the equatorial pitch angle diffusion rate (panels a and b) and the minimum resonant energy (panels c and d) for an EMIC wave packet. We used the global core plasma model (GCPM)44 2.2 for the electron and ion densities, and empirical geomagnetic field models38,41 with the observed date and location. We assumed parallel plane wave propagation of EMIC waves at the equator to simplify the calculation. Oblique propagation of EMIC waves can reduce the efficiency of REP, but the dependence on the wave normal angles is minor in the case of EMIC waves in the O + band45, so our assumption using parallel propagation is adequate. The equatorial diffusion rates under the observed conditions are given in an extremely high energy range > 10 MeV, but the diffusion rates in the case of a typical low geomagnetic field strength (170 nT and 110 nT for the 2015 and 2014 events, respectively) and a high cold electron density (250 cm−3 and 100 cm−3) are given in ultra-relativistic energies > 2 MeV over a wide pitch angle range (see Fig. 4a, b). EMIC waves in the H + and He + bands cannot typically resonate with ultra-relativistic electrons in such a wide pitch angle range at the equator25. Such EMIC waves mainly affect the pitch angle scattering of a few MeV electrons25,46, but the O + band EMIC waves effectively enhance the loss of ultra-relativistic (> 2 MeV) electrons. The wave frequency46 and hot ion density47,48 are also quite important for the loss of radiation belt electrons. This essential pitch angle scattering of ultra-relativistic electrons by EMIC waves is similar to that reported in previous studies49,50. Furthermore, non-resonant electrons at lower energies of up to ~ 100 keV, which are below the resonance cutoff of MeV energies in Fig. 4c and d, can still precipitate due to non-resonant interactions by EMIC wave packets with narrow edges51. The effects of non-resonant electrons are not taken into account the calculated diffusion rates, and such ~ 100 keV electron precipitation can still produce important atmospheric impacts because of a much larger population of ~ 100 keV electrons than that for ultra-relativistic electrons49,50.
Equatorial pitch angle diffusion and the minimum resonant energy of electrons for observed conditions. (a, b) Equatorial pitch angle diffusion as a function of equatorial pitch angle. Each curve is plotted every 1 MeV. Solid curves are the calculation results using observed conditions and dotted curves are for a low geomagnetic field strength and a high electron density. (c, d) Minimum resonant energy as a function of the EMIC wave frequency normalized by the cyclotron frequency of O + ions \({\Omega }_{\mathrm{O}+}\). The yellow rectangle indicates the observed frequency range of the EMIC/Pc1 waves, where \({B}_{0}\) is the geomagnetic field intensity at the magnetic equator and \({N}_{e}\) is the electron density. Red bold curves are the calculation results using observed conditions. Black solid and dotted curves are for the typical values of \({B}_{0}\) and \({N}_{e}\).
Since the observed amplitudes of EMIC/Pc1 waves in both events are close to several % of the background geomagnetic field at the equator, the nonlinear pitch angle scattering26,52,53 can significantly contribute to the observed REP events. The typical values of the quasi-linear diffusion rates of the radiation belt electrons range from several hours to a day, but the nonlinear pitch angle scattering by large amplitude EMIC waves having rising-tone frequency structures can cause more rapid (< 1 min) precipitation26. In the nonlinear effects, the radiation belt electrons having high pitch angles are guided to lower pitch angles by nonlinear wave trapping from the large amplitude EMIC waves53, but it is not sufficient to push them into the loss cone. Then, strong radiation belt electron precipitation is caused by a nonlinear scattering process at low pitch angles by the EMIC waves without nonlinear wave trapping53. The combination between nonlinear wave trapping and the scattering at low pitch angle by nonlinear effects of large amplitude EMIC waves leads to effective and rapid precipitation of radiation belt electrons into the atmosphere, with a timescale on the order of seconds, like radiation belt electron microbursts53. Previous observational studies on the relativistic electron microbursts in a co-located IPA54 and on 1-Hz IPA luminous modulations associated with rapid temporal variations of cosmic noise absorption55 support the existence of fast REP on the order of seconds by nonlinear effects. Therefore, the nonlinear effects by large amplitude EMIC waves can strongly contribute to the rapid REP events.
In this study, a clear spatial correspondence between mesospheric ozone loss and localized REP from the Earth's radiation belt was identified using the IPA during both geomagnetic storm and non-storm conditions. The unprecedented clear correspondence between the EMIC-driven IPA location and the localized mesospheric ozone loss was locally shown on a narrow spatial scale of ~ 400 km in the latitudinal direction. The observed IPA shows a wide longitudinal range over 1000 km, with the effects of the mesospheric jet, up to 60 m s−1 in mainly an east to west direction56, too small to be important for transport of ozone in the studied events. The mesospheric ozone can quickly respond within 1 to 1.5 h after the EMIC-driven REP. Another candidate driver of REP is whistler mode chorus waves57,58, but chorus waves cannot scatter the energetic protons causing the IPA. The direct link between IPA and mesospheric ozone destruction is extremely useful for quantitatively understanding of the atmospheric impacts on the loss of the radiation belt electrons by EMIC waves without the effects of chorus waves. Precipitation of MeV electrons was observed by the MAXI and MEPED instruments during both events as a concentrated downpour of electrons, but the observed EMIC/Pc1 waves in the O + band can resonate with ultra-relativistic electrons beyond the observed range of a few mega-electron-volts. The observed REP events are expected to widely range from the energy covered by MAXI/RBM-Z and MEPED to ultra-relativistic energies. Ultra-relativistic electrons directly penetrate to the mesosphere and the upper stratosphere and create chemical changes in this altitude range9. On the other hand, the flux of ultra-relativistic electrons over 10 MeV would be extremely small compared to that of relativistic electrons49,50. We have two possibilities to interpret the observed mesospheric ozone destruction. One is the direct effect by the ultra-relativistic (> 2 MeV) electron precipitation, which can strongly interact with EMIC/Pc1 waves in the O + band as shown in Fig. 4. The stopping height of the ultra-relativistic electrons is in the stratosphere below the mesosphere9, but the ionization rate of the atmosphere up to the stopping altitude is almost the same as that for the relativistic electrons9. Thus, the observed mesospheric ozone destruction can be caused by the effects of the ionization by the ultra-relativistic electrons above the stopping height. The other scenario is the direct effect by precipitating non-resonant electrons in the lower relativistic and/or ~ 100 keV energies51. These lower energy electrons as compared to the smaller numbers of ultra-relativistic electrons have a lower fractional scattering efficiency with the EMIC/Pc1 waves, but a much larger population. The lower fractional scattering of electrons of lower energy acting on a much larger population can produce a significant atmospheric impact. If the later scenario is the major reason, the rapid ozone destruction events would be observed during EMIC waves in other ion bands, not only for the O + band. On the other hand, a simulation study showed not only an effective REP, but also a nonlinear blocking of REP by large amplitude EMIC waves59. The precipitation blocking was effective at limited low pitch angles less than 10 deg. The precipitating flux of REP can be determined by a balance between the precipitation effects guiding particles into the loss cone in a wide pitch angle range and the precipitation blocking at low pitch angles. Our study should motivate future studies using combined wave and ozone data with incorporation with these precipitation and precipitation blocking effects through wave-particle interaction (quasi-linear, nonlinear, and non-resonance etc.) processes. Identifying the most important EMIC-driven precipitating electron energy to have the major atmospheric impacts remains an open question, so multi-coupled magnetosphere-ionosphere-atmosphere simulations using the quantitative flux measurements of radiation belt electrons into the IPA are also suggested so that further modeling and simulations can investigate the conditions under which immediate ozone destruction processes are created by various wave-particle interaction processes and other atmospheric effects. A previous simulation study suggested a weak mesospheric ozone destruction up to ~ 10% caused by EMIC-driven REP60, while the present observations show a greater 10 to 60% ozone destruction, which could have an impact similar to that from other EPP phenomena (pulsating aurora16 and microbursts17). The accumulated impacts of needle holes in the ozone layer by the IPA cannot be ignored when considering overall ozone changes in the mesosphere.
The EMIC wave-particle interactions shed light on the direct chemical effects of relativistic electrons (REP) on the Earth's atmospheric processes as an additional EPP source, along with other similar phenomena (e.g., solar proton events, pulsating aurora, and microbursts). Further study is needed to estimate the global effects of the EMIC-driven IPA from individual localized events. The type of IPA events shown here can be excited in different local time sectors at subauroral latitudes61. Unfortunately, satellite observations are limited due to orbit and viewing method. It remains unclear how the ozone loss by the IPA occurs in different local time sectors. If the continuous global distribution of IPAs is captured by a ground-based network (such as PWING (study of dynamical variation of Particles and Waves in the INner magnetosphere using Ground-based network observations)40) and satellite constellation62 observations covering a whole longitudinal range at subauroral latitudes, then quantifying the localized processes of REP-driven chemical impacts in the mesosphere could contribute to an improvement of global circulation models for the middle atmosphere63,64. EMIC waves in Jupiter's magnetosphere similarly scatter energetic heavy ions (sulfur and oxygen) to its ionosphere in a similar manner to the case for Earth, causing similar IPA emissions in the Jovian polar regions65. The visualization ("when" and "where") of the EMIC-driven EPP as IPA provides key insights for understanding the dynamical space weather effects on the atmospheres not only for Earth, but also for the Jovian planets.
We used version Level 2A data of ozone from a 9.6 µm wavelength measurement by SABER35 on board the TIMED satellite orbiting at an altitude of 625 km with an orbital inclination of 74.1°. The Ozone measurement provides 35% precision in the mesosphere and the upper stratosphere66. The advantages of limb sounding by SABER are a precise vertical resolution and a stable background, unlike the variable surface of the Earth. The effect of the low horizontal resolution is not a determinant to this study, because the IPA is seen over a wide (> 1000 km) longitudinal width in a narrow (hundreds of kilometers) latitude.
To visualize the ionospheric footprint of REP from the Earth's radiation belt, we used IPA events observed by the SSUSI34 aboard the DMSP satellites in an 840-km polar orbit. The Lyman α (121.6 nm) emission by SSUSI is enhanced in the region of the IPA. The signature of the Lyman α emission using a space-based imager that is unaffected by lower atmosphere emissions contributes to identifying the localization of REP with a wide field of view. The auroral image is projected to an altitude of 110 km. The advantage of space-borne observations of auroral images is that the spatial coverage is wide, and observations are unaffected by weather conditions. The ultraviolet wavelength used is also visible when lighting conditions for optical emissions are not favorable. The duration of the IPA and REP is determined from the duration of the related EMIC/Pc1 waves.
To confirm REP, we used data from the RBM of the MAXI36 aboard the ISS orbiting at an altitude of ~ 400 km with an orbital inclination of 51.6°. The MAXI/RBM consists of two sets of PIN diode detectors with horizontal (H) and zenith (Z) directions. The RBM-H and RBM-Z detectors are sensitive to relativistic electrons above 0.3 MeV. In this study, we used RBM-Z data as a proxy for REP into the atmosphere.
The MEPED instrument on board the POES satellite constellation has several energy channels for electrons, but we used the E4 channel data (300 keV–2.5 MeV) from the 0° detector as a typical proxy for the REP42.
We focused on EMIC/Pc1 geomagnetic pulsations detected by ground-based magnetometers for identifying EMIC waves in the magnetosphere. The stations are located at Athabasca (ATH, 54.7°N, 246.7°E, L = 4.5)40 and Fort Churchill (FCHU, 58.8°N, 265.9°E, L = 7.2)43 in Canada. Sampling frequencies of induction magnetometers at ATH and FCHU are 64 Hz and 20 Hz, respectively. The Pc1 wave intensity was calculated from \(\sqrt {H^{2} + D^{2} }\), where H is the magnetic north–south component and D is the east–west component. The short time Fourier transform was used for a 60-s waveform with an overlap of 95%.
When a Gaussian wave spectral density is assumed for an EMIC wave, the pitch angle diffusion rate of relativistic electrons by the EMIC wave is approximately written as25
where
\({\Omega }_{e}\) is the electron cyclotron frequency, \(E = \frac{{E_{k} }}{{m_{e} c^{2} }}\) is the dimensionless electron kinetic energy of the electron energy \(E_{k}\), \(m_{e}\) is the electron rest mass, \(c\) is the speed of light, \(R\) is the ratio of the magnetic field wave power density to the background magnetic field power density, \(\nu = \sqrt \pi\cdot {\text{erf}}\left( 1 \right) \approx 1.49\), \({\text{erf}}\) is the error function, \(\alpha^{*} = \frac{{{\Omega }_{e}^{2} }}{{\omega_{pe}^{2} }}\), \(\omega_{pe}\) is the electron plasma frequency, and \(\varepsilon = \frac{{m_{e} }}{{m_{p} }}\) is the mass ratio between the electron rest mass \(m_{e}\) and the proton rest mass \(m_{p}\). The exponential function in Eq. (1) is derived from a Gaussian wave spectral density with
where \({\Omega }_{H + }\) is the proton cyclotron frequency, \(\omega_{m}\) is the center frequency of an EMIC wave, and \(\delta \omega\) is the bandwidth of an EMIC wave. We define the density ratio of each ion as \(\eta_{H + } = \frac{{n_{H + } }}{{n_{e} }}\), \(\eta_{He + } = \frac{{n_{He + } }}{{n_{e} }}\), and \(\eta_{O + } = \frac{{n_{O + } }}{{n_{e} }}\), where \(n_{e}\), \(n_{H + }\), \(n_{He + }\), and \(n_{O + }\) are the number densities of the cold electrons, protons, He + ions, and O + ions, respectively. We used the magnetic field wave power density of EMIC waves of 10 nT for a Gaussian spectrum with \(\omega_{m} = 0.1875\) Hz and \(\delta \omega = 0.0625\) Hz for the June 22, 2015, event, and the magnetic field wave power density of 1 nT with \(\omega_{m} = 0.1275\) Hz and \(\delta \omega = 0.0575\) Hz for the August 12, 2014, event in Fig. 4. The wave normal angle was simply assumed to be that of parallel plane wave propagation with respect to the background geomagnetic field line.
TIMED SABER data are available from the website http://saber.gats-inc.com/data.php. DMSP SSUSI data used in this study are publicly available from https://ssusi.jhuapl.edu/data_products. MAXI RBM data were obtained from https://data.darts.isas.jaxa.jp/pub/maxi/rbm/. POES MEPED data were obtained from https://www.ngdc.noaa.gov/stp/satellite/poes/dataaccess.html. EMIC/Pc1 wave data at ATH and FCHU were obtained from https://stdb2.isee.nagoya-u.ac.jp/magne/ and https://www.carisma.ca/, respectively. The Dst index used in this study was provided by the WDC for Geomagnetism, Kyoto (http://wdc.kugi.kyoto-u.ac.jp/wdc/Sec3.html, https://isds-datadoi.nict.go.jp/wds/10.17593__14515-74000.html).
The geomagnetic map was drawn using the mapping package M_Map for MATLAB, available online at www.eoas.ubc.ca/~rich/map.html.
Rodger, C. J. et al. Radiation belt electron precipitation into the atmosphere: recovery from a geomagnetic storm. J. Geophys. Res. 112, A11307. https://doi.org/10.1029/2007JA012383 (2007).
Article ADS CAS Google Scholar
Rozanov, E., Calisto, M., Egorova, T., Peter, T. & Schmutz, W. Influence of the precipitating energetic particles on atmospheric chemistry and climate. Surv. Geophys. 33, 483–501. https://doi.org/10.1007/s10712-012-9192-0 (2012).
Article ADS Google Scholar
Seppälä, A., Matthes, K., Randall, C. E. & Mironova, I. A. What is the solar influence on climate? Overview of activities during CAWSES-II. Prog. Earth Planet. Sci. 1, 24. https://doi.org/10.1186/s40645-014-0024-3 (2014).
Article ADS Google Scholar
Andersson, M. et al. Missing driver in the Sun-Earth connection from energetic electron precipitation impacts mesospheric ozone. Nat Commun. 5, 5197. https://doi.org/10.1038/ncomms6197 (2014).
Article ADS CAS PubMed Google Scholar
Heath, D. F., Kruger, A. J. & Crutzen, P. J. Solar proton events: influence on stratospheric ozone. Science 453(197), 886–889 (1977).
Article ADS Google Scholar
Jackman, C. H. et al. Northern hemisphere atmospheric effects due to the July 2000 solar proton events. Geophys. Res. Lett. 28, 2883–2886 (2001).
Article ADS CAS Google Scholar
Funke, B. et al. Composition changes after the "Halloween" solar proton event: the high-energy particle precipitation in the atmosphere (HEPPA) model versus MIPAS data intercomparison study. Atmos. Chem. Phys. 11, 9089–9139. https://doi.org/10.5194/acp-11-9089-2011 (2011).
Article ADS CAS Google Scholar
Solomon, S., Rusch, D. W., Gerard, J.-C., Reid, G. C. & Crutzen, P. J. The effect of particle precipitation events on the neutral and ion chemistry of the middle atmosphere. II—odd hydrogen. Planet. Space Sci. 29, 885–893. https://doi.org/10.1016/0032-0633(81)90078-7 (1981).
Article ADS CAS Google Scholar
Turunen, et al. Impact of different energies of precipitating particles on NOx generation in the middle and upper atmosphere during geomagnetic storms. J. Atmos. Solar Terr. Phys. 71(10–11), 1176–1189 (2009).
Article ADS CAS Google Scholar
Randall, C. E. et al. Enhanced NOx in 2006 linked to strong upper stratospheric Arctic vortex. Geophys. Res. Lett. 33, L18811. https://doi.org/10.1029/2006GL027160 (2006).
Article ADS CAS Google Scholar
Andersson, M. E. et al. Polar ozone response to energetic particle precipitation over decadal time scales: the role of medium-energy electrons. J. Geophys. Res.: Atmos. 123, 607–622. https://doi.org/10.1002/2017JD027605 (2018).
Article ADS CAS Google Scholar
Jia, J., Kero, A., Kalakoski, N., Szeląg, M. E. & Verronen, P. T. Is there a direct solar proton impact on lower-stratospheric ozone?. Atmos. Chem. Phys. 20, 14969–14982. https://doi.org/10.5194/acp-20-14969-2020 (2020).
Article ADS CAS Google Scholar
Baumgaertner, A. J. G., Seppälä, A., Jöckel, P. & Clilverd, M. A. Geomagnetic activity related NOx enhancements and polar surface air temperature variability in a chemistry climate model: modulation of the NAM index. Atmos. Chem. Phys. 11, 4521–4531. https://doi.org/10.5194/acp-11-4521-2011 (2011).
Article ADS CAS Google Scholar
Seppälä, A., Lu, H., Clilverd, M. A. & Rodger, C. J. Geomagnetic activity signatures in wintertime stratosphere wind, temperature, and wave response. J. Geophys. Res. Atmos. 118, 2169–2183. https://doi.org/10.1002/jgrd.50236 (2013).
Article ADS Google Scholar
Bornebusch, J. P., Wissing, J. M. & Kallenrode, M.-B. Solar particle precipitation into the polar atmosphere and their dependence on hemisphere and local time. Adv. Space Res. 45, 5. https://doi.org/10.1016/j.asr.2009.11.008 (2010).
Article CAS Google Scholar
Turunen, E. et al. Mesospheric ozone destruction by high-energy electron precipitation associated with pulsating aurora. J. Geophys. Res. Atmos. 121, 11852–11861. https://doi.org/10.1002/2016JD025015 (2016).
Article ADS CAS Google Scholar
Seppälä, A. et al. Relativistic electron microburst events: modeling the atmospheric impact. Geophys. Res. Lett. 45, 1141–1147. https://doi.org/10.1002/2017GL075949 (2018).
Article ADS Google Scholar
Verronen, P. T. et al. Diurnal variation of ozone depletion during the October–November 2003 solar proton events. J. Geophys. Res. 110, A09S32 (2005). https://doi.org/10.1029/2004JA010932
Article CAS Google Scholar
Yahnin, A. G., Yahnina, T. A. & Frey, H. U. Subauroral proton spots visualize the Pc1 source. J. Geophys. Res. 112, A10223. https://doi.org/10.1029/2007JA012501 (2007).
Article ADS CAS Google Scholar
Miyoshi, Y. et al. Precipitation of radiation belt electrons by EMIC waves, observed from ground and space. Geophys. Res. Lett. 35, L23101. https://doi.org/10.1029/2008GL035727 (2008).
Article ADS CAS Google Scholar
Sakaguchi, K. et al. Visualization of ion cyclotron wave and particle interactions in the inner magnetosphere via THEMIS-ASI observations. J. Geophys. Res. 117, A10204. https://doi.org/10.1029/2012JA018180 (2012).
Article ADS CAS Google Scholar
Van Allen, J. A. & Frank, L. A. Radiation around the earth to a radial distance of 107,400 km. Nature 183, 430–434 (1959).
Article ADS Google Scholar
Usanova, M. E. et al. Conjugate ground and multisatellite observations of compression-related EMIC Pc1 waves and associated proton precipitation. J. Geophys. Res. 115, A07208. https://doi.org/10.1029/2009JA014935 (2010).
Article ADS Google Scholar
Jordanova, V. K., Spasojevic, M. & Thomsen, M. F. Modeling the electromagnetic ion cyclotron wave-induced formation of detached subauroral proton arcs. J. Geophys. Res. 112, A08209. https://doi.org/10.1029/2006JA012215 (2007).
Article ADS CAS Google Scholar
Summers, D. & Thorne, R. M. Relativistic electron pitch-angle scattering by electromagnetic ion cyclotron waves during geomagnetic storms. J. Geophys. Res. 108, 1143. https://doi.org/10.1029/2002JA009489,A4 (2003).
Article Google Scholar
Omura, Y. & Zhao, Q. Relativistic electron microbursts due to nonlinear pitch angle scattering by EMIC triggered emissions. J. Geophys. Res. Space Phys. 118, 5008–5020. https://doi.org/10.1002/jgra.50477 (2013).
Article ADS CAS Google Scholar
Cornwall, J. M. Cyclotron instabilities and electromagnetic emission in the ultra low frequency and very low frequency ranges. J. Geophys. Res. 70(1), 61–69. https://doi.org/10.1029/JZ070i001p00061 (1965).
Article ADS Google Scholar
Anderson, B. J., Erlandson, R. E., & Zanetti, L. J. A statistical study of Pc 1–2 magnetic pulsations in the equatorial magnetosphere: 1. Equatorial occurrence distributions. J. Geophys. Res. 97(A3), 3075–3088. https://doi.org/10.1029/91JA02706 (1992).
Pickett, J. S. et al. Cluster observations of EMIC triggered emissions in association with Pc1 waves near Earth's plasmapause. Geophys. Res. Lett. 37, L09104. https://doi.org/10.1029/2010GL042648 (2010).
Article ADS CAS Google Scholar
Nomura, R. et al. Pulsating proton aurora caused by rising tone Pc1 waves. J. Geophys. Res. Space Phys. 121, 1608–1618. https://doi.org/10.1002/2015JA021681 (2016).
Article ADS Google Scholar
Ozaki, M. et al. Fast modulations of pulsating proton aurora related to subpacket structures of Pc1 geomagnetic pulsations at subauroral latitudes. Geophys. Res. Lett. 43, 7859–7866. https://doi.org/10.1002/2016GL070008 (2016).
Article ADS Google Scholar
Loto’aniu, T. M., Fraser, B. J. & Waters, C. L. Propagation of electromagnetic ion cyclotron wave energy in the magnetosphere. J. Geophys. Res. 110, A07214. https://doi.org/10.1029/2004JA010816 (2005).
Article ADS Google Scholar
Ozaki, M., Shiokawa, K., Horne, R. B., Engebretson, M. J., Lessard, M., Ogawa, Y., et al. Magnetic conjugacy of Pc1 waves and isolated proton precipitation at subauroral latitudes: Importance of ionosphere as intensity modulation region. Geophys. Res. Lett.. 48, e2020GL091384. (2021). https://doi.org/10.1029/2020GL091384
Paxton, L. J., Morrison, D., Zhang, Y., Kil, H., Wolven, B., Ogorzalek, B. S., et al. Validation of remote sensing products produced by the Special Sensor Ultraviolet Scanning Imager (SSUSI): A far UV-imaging spectrograph on DMSP F-16. Proc. SPIE 4485, Optical Spectroscopic Techniques, Remote Sensing, and Instrumentation for Atmospheric and Space Research IV, 4485, 338–348. (2002). https://doi.org/10.1117/12.454268
Russell, J. M. III, M. G. Mlynczak, L. L. Gordley, J. Tansock, and R. Esplin. An overview of the SABER experiment and preliminary calibration results. Proc. SPIE 3756, 277–288. https://doi.org/10.1117/12.366382 (1999).
Matsuoka, M. et al. The MAXI Mission on the ISS: science and instruments for monitoring all-sky X-ray images. Publ. Astron. Soc. Jpn. 61(5), 999–1010 (2009).
Article ADS CAS Google Scholar
Alken, P. et al. International geomagnetic reference field: the thirteenth generation. Earth Planets Space. 73, 49. https://doi.org/10.1186/s40623-020-01288-x (2021).
Article ADS Google Scholar
Tsyganenko, N. A. A model of the magnetosphere with a dawn-dusk asymmetry, 1 Mathematical structure. J. Geophys. Res. 107, A8. https://doi.org/10.1029/2001JA000219 (2002).
Article Google Scholar
Oyama, S.-I. et al. Lower thermospheric wind variations in auroral patches during the substorm recovery phase. J. Geophys. Res. Space Physics 121, 3564–3577. https://doi.org/10.1002/2015JA022129 (2016).
Article ADS Google Scholar
Shiokawa, K. et al. Ground-based instruments of the PWING project to investigate dynamics of the inner magnetosphere at subauroral latitudes as a part of the ERG-ground coordinated observation network. Earth Planets Space. 69, 160. https://doi.org/10.1186/s40623-017-0745-9 (2017).
Article ADS Google Scholar
Tsyganenko, N. A. & Sitnov, M. I. Modeling the dynamics of the inner magnetosphere during strong geomagnetic storms. J. Geophys. Res. 110, A03208. https://doi.org/10.1029/2004JA010798 (2005).
Article ADS Google Scholar
Peck, E. D., Randall, C. E., Green, J. C., Rodriguez, J. V. & Rodger, C. J. POES MEPED differential flux retrievals and electron channel contamination correction. J. Geophys. Res. Space Physics. 120, 4596–4612. https://doi.org/10.1002/2014JA020817 (2015).
Article ADS Google Scholar
Mann, I. R. et al. The upgraded CARISMA magnetometer array in the THEMIS era. Space Sci Rev. 141, 413–451. https://doi.org/10.1007/s11214-008-9457-6 (2008).
Article ADS Google Scholar
Gallagher, D. L., Craven, P. D. & Comfort, R. H. Global core plasma model. J. Geophys. Res. 105(A8), 18819–18833. https://doi.org/10.1029/1999JA000241 (2000).
Article ADS CAS Google Scholar
Ni, B. et al. Resonant scattering of outer zone relativistic electrons by multiband EMIC waves and resultant electron loss time scales. J. Geophys. Res. Space Phys. 120, 7357–7373. https://doi.org/10.1002/2015JA021466 (2015).
Article ADS Google Scholar
Albert, J. M. Evaluation of quasi-linear diffusion coefficients for EMIC waves in a multispecies plasma. J. Geophys. Res. 108, 1249. https://doi.org/10.1029/2002JA009792,A6 (2003).
Article Google Scholar
Ni, B. et al. Hot plasma effects on the cyclotron-resonant pitch-angle scattering rates of radiation belt electrons due to EMIC waves. Geophys. Res. Lett. 45, 21–30. https://doi.org/10.1002/2017GL076028 (2018).
Article ADS Google Scholar
Cao, X. et al. Scattering of ultra-relativistic electrons in the Van Allen radiation belts accounting for hot plasma effects. Sci Rep. 7, 17719. https://doi.org/10.1038/s41598-017-17739-7 (2017).
Article ADS CAS PubMed PubMed Central Google Scholar
Usanova, M. E. et al. Effect of EMIC waves on relativistic and ultrarelativistic electron populations: ground-based and Van Allen Probes observations. Geophys. Res. Lett. 41, 1375–1381. https://doi.org/10.1002/2013GL059024 (2014).
Article ADS CAS Google Scholar
Shprits, Y. et al. Wave-induced loss of ultra-relativistic electrons in the Van Allen radiation belts. Nat Commun 7, 12883. https://doi.org/10.1038/ncomms12883 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Chen, L., Thorne, R. M., Bortnik, J. & Zhang, X.-J. Nonresonant interactions of electromagnetic ion cyclotron waves with relativistic electrons. J. Geophys. Res. Space Phys. 121, 9913–9925. https://doi.org/10.1002/2016JA022813 (2016).
Article ADS Google Scholar
Kubota, Y., Omura, Y. & Summers, D. Relativistic electron precipitation induced by EMIC-triggered emissions in a dipole magnetosphere. J. Geophys. Res. Space Phys. 120, 4384–4399. https://doi.org/10.1002/2015JA021017 (2015).
Article ADS Google Scholar
Kubota, Y. & Omura, Y. Rapid precipitation of radiation belt electrons induced by EMIC rising tone emissions localized in longitude inside and outside the plasmapause. J. Geophys. Res. Space Phys. 122, 293–309. https://doi.org/10.1002/2016JA023267 (2017).
Article ADS Google Scholar
Shumko, M. et al. Proton aurora and relativistic electron microbursts scattered by electromagnetic ion cyclotron waves. Front. Astron. Space Sci. 9, 975123. https://doi.org/10.3389/fspas.2022.975123 (2022).
Article ADS Google Scholar
Ozaki, M. et al. Discovery of 1 Hz range modulation of isolated proton aurora at subauroral latitudes. Geophys. Res. Lett. 45, 1209–1217. https://doi.org/10.1002/2017GL076486 (2018).
Article ADS Google Scholar
Sato, K. et al. On the origins of mesospheric gravity waves. Geophys. Res. Lett. 36, L19801. https://doi.org/10.1029/2009GL039908 (2009).
Article ADS Google Scholar
Summers, D., Ni, B., & Meredith, N. P. Timescales for radiation belt electron acceleration and loss due to resonant wave-particle interactions: 2. Evaluation for VLF chorus, ELF hiss, and electromagnetic ion cyclotron waves. J. Geophys. Res., 112, A04207. https://doi.org/10.1029/2006JA011993 (2007).
Article ADS Google Scholar
Hsieh, Y.-K., Omura, Y., & Kubota, Y. Energetic electron precipitation induced by oblique whistler mode chorus emissions. J. Geophys. Res.: Space Phys. 127, e2021JA029583, (2022). https://doi.org/10.1029/2021JA029583
Bortnik, J., Albert, J. M., Artemyev, A., Li, W., Jun, C.-W., Grach, V. S., & Demekhov, A. G. Amplitude dependence of nonlinear precipitation blocking of relativistic electrons by large amplitude EMIC waves. Geophys. Res. Lett. 49, e2022GL098365 (2022). https://doi.org/10.1029/2022GL098365
Hendry, A. T., Seppälä, A., Rodger, C. J., & Clilverd, M. A. Impact of EMIC-wave driven electron precipitation on the radiation belts and the atmosphere. J. Geophys. Res.: Space Phys. 126, e2020JA028671 (2021). https://doi.org/10.1029/2020JA028671
Semenova, N. V., Yahnin, A. G., Yahnina, T. A. & Demekhov, A. G. Properties of localized precipitation of energetic protons equatorward of the isotropic boundary. Geophys. Res. Lett. 46, 10932–10940. https://doi.org/10.1029/2019GL085373 (2019).
Article ADS CAS Google Scholar
Morley, S. K., Sullivan, J. P., Henderson, M. G., Blake, J. B. & Baker, D. N. The global positioning system constellation as a space weather monitor: comparison of electron measurements with Van Allen Probes data. Space Weather 14, 76–92. https://doi.org/10.1002/2015SW001339 (2016).
Article ADS Google Scholar
Verronen, P. T., Andersson, M. E., Marsh, D. R., Kovács, T. & Plane, J. M. C. WACCM-D: whole atmosphere community climate model with D-region ion chemistry. J. Adv. Model. Earth Syst. 8, 954–975. https://doi.org/10.1002/2015MS000592 (2016).
Article ADS Google Scholar
Danabasoglu, G., Lamarque, J.-F., Bacmeister, J., Bailey, D. A., DuVivier, A. K., Edwards, J., et al. The community earth system model version 2 (CESM2). J. Adv. Model. Earth Syst. 12, e2019MS001916, (2020). https://doi.org/10.1029/2019MS001916
Yao, Z. et al. Revealing the source of Jupiter's x-ray auroral flares. Sci. Adv. 7, 28. https://doi.org/10.1126/sciadv.abf0851 (2021).
Article Google Scholar
Rong, P. P., Russell, J. M., Mlynczak, M. G., Remsberg, E. E., Marshall, B. T., Gordley, L. L., & López-Puertas, M. Validation of thermosphere ionosphere mesosphere energetics and dynamics/sounding of the atmosphere using broadband emission radiometry (TIMED/SABER) v1.07 ozone at 9.6 μm in altitude range 15–70 km. J. Geophys. Res. 114, D04306. https://doi.org/10.1029/2008JD010073 (2009).
Download references
The present study was supported by the Japan Society for the Promotion of Science, KAKENHI grants JP16H06286 and JP20H02162. We are grateful to the TIMED/SABER, DMSP/SSUSI, and POES/MEPED teams for providing the data. This study has made use of MAXI data provided by RIKEN, JAXA, and the MAXI team. The authors thank the CARISMA team for data. CARISMA is operated by the University of Alberta, funded by the Canadian Space Agency.
Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa, Japan
Mitsunori Ozaki, Satoshi Yagitani & Shion Hashimoto
Institute for Space-Earth Environmental Research, Nagoya University, Nagoya, Japan
Kazuo Shiokawa & Yuichi Otsuka
National Institute of Polar Research, Tachikawa, Japan
Ryuho Kataoka
Department of Polar Science, The Graduate University for Advanced Studies, SOKENDAI, Tachikawa, Japan
Ryuho Kataoka
NASA Langley Research Center, Hampton, VA, USA
Martin Mlynczak
The Johns Hopkins University Applied Physics Laboratory, Laurel, MD, USA
Larry Paxton
Athabasca University Observatories, Athabasca, AB, Canada
Martin Connors
Department of Physics and Astronomy, University of Calgary, Calgary, AB, Canada
Martin Connors
Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, Sagamihara, Japan
Satoshi Nakahira
Department of Physics, University of Alberta, Edmonton, AB, Canada
Ian Mann
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
M.O. conceived of the present study using the IPA, analyzed satellite and ground data, performed numerical computations of EMIC wave-particle interactions, and wrote the manuscript. K.S. is the project leader for the induction magnetometer at ATH and contributed to data interpretation. R.K. contributed to the model calculation of the direct ionization by REP and interpretation. M.M. is the principal investigator of the TIMED/SABER mission and contributed to the SABER data analysis. L.P. is the principal investigator of the SSUSI mission and contributed to the SSUSI data analysis for Lyman \({\alpha }\) emissions from the IPA. M.C. contributed to the induction magnetometer observation at ATH and the revision of Introduction and Discussion sections. S.Y. contributed to the revision of the manuscript. S.H. performed the figure plotting. Y.O. contributed to the data handling of the induction magnetometer at ATH. S.N. gives the MAXI/RBM data. I.M. supplied the induction magnetometer data from FCHU and contributed to the interpretation of REP. All authors provided feedback to improve the manuscript, and helped share the analysis and the manuscript.
Correspondence to Mitsunori Ozaki.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Ozaki, M., Shiokawa, K., Kataoka, R. et al. Localized mesospheric ozone destruction corresponding to isolated proton aurora coming from Earth's radiation belt. Sci Rep 12, 16300 (2022). https://doi.org/10.1038/s41598-022-20548-2
Download citation
Received: 26 May 2022
Accepted: 14 September 2022
Published: 11 October 2022
DOI: https://doi.org/10.1038/s41598-022-20548-2
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
Earth, Planets and Space (2023)
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.