JGR:Space physics

Van Allen Probe observations indicate that whistler-mode hiss waves below 1 kHz are absorbed at low altitudes near magnetic equator. The lowest cutoff frequency of equatorial hiss is close to the gyrofrequency of hydrogen ions. The lowest cutoff altitude of global hiss is extracted when its occurrence rate is equal to 0.005 on the plane of altitude (L in RE) and magnetic local time (MLT). By fitting the lowest cutoff altitude of global hiss, we constructed the empirical model of the lowest cutoff altitude of equatorial hiss under geomagnetically quiet (AE < 200 nT) and active (AE ≥ 200 nT) conditions. The enhanced substorm activities reduce the lowest cutoff altitude of hiss waves on the dawnside (MLT ∼ 1–5 hr), whereas the lowest cutoff altitude of the dayside hiss is nearly fixed at ∼1.1 RE (MLT ∼ 6–20 hr). From the dayside to the nightside (MLT ∼ 0–6 hr and 20–24 hr), the lowest cutoff altitude of equatorial hiss raises gradually from 1.1 RE to 1.4 RE.

This paper investigates the impact of a powerful gamma ray burst (GRB) that occurred on 9 October 2022, on the Earth's environment using a very low frequency receiver (VLF) to probe the lower ionospheric region (the D region). In addition to the VLF data analysis, we employ numerical simulation through the Long Wavelength Propagation Capability code (LWPC) to derive the increase in the D− region electron density. Our results revealed discernible perturbations in amplitude and phase across all transmitter paths (NAA, DHO, ICV, and NSC) to the Algiers receiver persisting for 40 min. At the maximum of the signal perturbation, the LWPC simulation results showed a decrease in the mean new reference height h′ from 74 to 65.71 km, along with an increase in the sharpness factor β from 0.3 to 0.4875 km−1. Under these new conditions, the electron density increased from its ambient value (216.10 cm−3) to 33.7 103 cm−3.

The ability to understand and model ionospheric plasma flow on all spatial scales has important implications for operational space weather models. This study exploits a recently developed method to statistically separate large-scale and meso-scale contributions to probability density functions (PDFs) of ionospheric flow vorticity measured by the Super Dual Auroral Radar Network (SuperDARN). The SuperDARN vorticity data are first sub-divided depending on the Interplanetary Magnetic Field (IMF) direction, and the separation method is applied to PDFs of vorticity compiled in spatial regions of size 1° of geomagnetic latitude by 1 hr of magnetic local time, covering much of the high-latitude ionosphere in the northern hemisphere. The resulting PDFs are fit by model functions using maximum likelihood estimation (MLE) and the spatial variations of the MLE estimators for both the large-scale and meso-scale components are presented. The spatial variations of the large-scale vorticity estimators are ordered by the average ionospheric convection flow, which is highly dependent on the IMF direction. The spatial variations of the meso-scale vorticity estimators appear independent of the senses of vorticity and IMF direction, but have a different character in the polar cap, the cusp, the auroral region, and the sub-auroral region. The paper concludes by discussing the sources of the vorticity components in the different regions, and the consequences for the fidelity of ionospheric plasma flow models.

Updated summaries of the August 1972 and March 1989 space weather events have been constructed. The features of these two events are compared to the Carrington 1859 event and a few other major space weather events. It is concluded that solar active regions release energy in a variety of forms (X-rays, EUV photons, visible light, coronal mass ejection (CME) plasmas and fields) and they in turn can produce other energetic effects (solar energetic particles (SEPs), magnetic storms) in a variety of ways. It is clear that there is no strong one-to-one relationship between these various energy sinks. The energy is often distributed differently from one space weather event to the next. Concerning SEPs accelerated at interplanetary CME (ICME) shocks, it is concluded that the Fermi mechanism associated with quasi-parallel shocks is relatively weak and that the gradient drift mechanism (electric fields) at quasi-perpendicular shocks will produce harder spectra and higher fluxes. If the 4 August 1972 intrinsic magnetic cloud condition (southward interplanetary magnetic field instead of northward) and the interplanetary Sun to 1 au conditions were different, a 4 August 1972 magnetic storm and magnetospheric dawn-to-dusk electric fields substantially larger than the Carrington event would have occurred. Under these special interplanetary conditions, a Miyake et al. (2012), https://doi.org/10.1038/nature11123-like extreme SEP event may have been formed. The long duration complex 1989 storm was probably greater than the Carrington storm in the sense that the total ring current particle energy was larger.

Electron density enhancements in the ionospheric D-region due to the precipitation of high-energy electrons (>30 keV) have been measured as increases in cosmic radio noise absorption (CNA) using ground-based riometers. CNA has been studied since the 1960s. However, there have been few studies of the spatiotemporal development of CNA at multi-point ground stations distributed in longitude at subauroral latitudes, where plasma particles with a wide energy range are intermingled. In this study, we analyzed the longitudinal development of CNA steep increases using simultaneous riometer observations at six stations at subauroral latitudes in Canada, Alaska, Russia, and Iceland over 3 years from 2017 to 2020. The results revealed that the occurrence rate of steep increases in CNA was highest at midnight at 22-08 magnetic local time (MLT), and lowest near dusk at 17–21 MLT. We also showed statistically that the CNA steep increases expanded eastward on the dawn side and westward on the dusk side. The CNA expansion velocity was slightly faster than the results of previous studies in the auroral zone. Correlation and superposed epoch analyses of CNA with solar wind and geomagnetic parameters revealed that CNA intensity was dependent on the Interplanetary Magnetic Field Bz, Interplanetary Electric Field Ey, SYM-H index, and SME index. These results indicate that the CNA at subauroral latitudes is closely related to solar wind and geomagnetic activities, and its propagation characteristics correspond to the dynamics of high energy electrons in the inner magnetosphere.

Electromagnetic Ion Cyclotron (EMIC) wave scattering has been proved to be responsible for the fast loss of both radiation belt (RB) electrons and ring current (RC) protons. However, its role in the concurrent dropout of these two co-located populations remains to be quantified. In this work, we study the effect of EMIC wave scattering on both populations during the 27 February 2014 storm by employing the global physics-based RAM-SCB model. Throughout this storm event, MeV RB electrons and 100s keV RC protons experienced simultaneous dropout following the occurrence of intense EMIC waves. By implementing data-driven initial and boundary conditions, we perform simulations for both populations through the interplay with EMIC waves and compare them against Van Allen Probes observations. The results indicate that by including EMIC wave scattering loss, especially by the He-band EMIC waves, the model aligns closely with data for both populations. Additionally, we investigate the simulated pitch angle distributions (PADs) for both populations. Including EMIC wave scattering in our model predicts a 90° peaked PAD for electrons with stronger losses at lower pitch angles, while protons exhibit an isotropic PAD with enhanced losses at pitch angles above 40°. Furthermore, our model predicts considerable precipitation of both particle populations, predominantly confined to the afternoon to midnight sector (12 hr < MLT < 24 hr) during the storm's main phase, corresponding closely with the presence of EMIC waves.

A new full model of the atmospheric transport of cosmogenic 10Be is presented based on the specialized SOCOL-AERv2-BE chemistry-climate model coupled with the CRAC:10Be isotope production model. The model includes all the relevant atmospheric processes and allows computing the isotope concentration at any given location and time. The full model is directly compared with 10Be isotope measurements in five Antarctic and Greenland ice cores for the period 1980–2007. The model reasonably well reproduces the average concentration and solar-cycle dependency or the lack of it for most observational sites but does not perfectly catch the interannual variability at sites with complex orography likely due to the coarse model grid. This implies that the model correctly reproduces the large-scale atmospheric dynamics but effectively averages out synoptic-scale variability. It is found that the dominant source of 10Be is located in the middle stratosphere (25–40 km), in the tropical (<30° latitudes) and polar (>60°) regions, as produced by galactic cosmic rays and solar energetic particles, respectively. It is shown that >60% (90%) of 10Be produced in the atmosphere reaches the Earth's surface within one (two) years, respectively. For practical purposes, a simple parameterization of the full-model results is presented which agrees with the full model within 20% in polar regions. This parameterization allows one to make a quick estimate of near-ground 10Be concentrations based only on production rates without heavy calculations. This practical approach can be applied to studies of solar and geomagnetic variability using cosmogenic isotopes.

Pulsating Aurora (PsA) is one of the major classes of diffuse aurora associated with precipitation of a few to a few tens of keV electrons from the magnetosphere. Recent studies suggested that, during PsA, more energetic (i.e., sub-relativistic/relativistic) electrons precipitate into the ionosphere at the same time. Those electrons are considered to be scattered at the higher latitude part of the magnetosphere by whistler-mode chorus waves propagating away from the magnetic equator. However, there have been no actual cases of simultaneous observations of precipitating electrons causing PsA (PsA electrons) and chorus waves propagating toward higher latitudes; thus, we still do not quite well understand under what conditions PsA electrons become harder and precipitate to lower altitudes. To address this question, we have investigated an extended interval of PsA on 12 January 2021, during which simultaneous observations with the Arase satellite, ground-based all-sky imagers and the European Incoherent SCATter (EISCAT) radar were conducted. We found that, when the PsA shape became patchy, the PsA electron energy increased and Arase detected intense chorus waves at magnetic latitudes above 20°, indicating the propagation of chorus waves up to higher latitudes along the field line. A direct comparison between the irregularities of the magnetospheric electron density and the emission intensity of PsA patches at the footprint of the satellite suggests that the PsA morphology and the energy of PsA electrons are determined by the presence of “magnetospheric density ducts,” which allow chorus waves to travel to higher latitudes and thereby precipitate more energetic electrons.

This research examines the plasma processes under penetration of the plasma clouds (plasmoids) across the magnetopause which is modeled as a tangential discontinuity (TD). Cases with the parallel magnetic field in both sides out of the TD are under investigation. Plasma parameters and magnetic field were chosen from the MMS mission and other spacecraft observations. The results are important for understanding the following basic space plasma physics problems: (a) plasma cloud deformation and strong phase mixing with magnetospheric plasma; (b) the transfer of mass, momentum and energy of magnetosheath and magnetic cloud plasma into magnetospheric plasmas; (c) necessary conditions for plasma cloud penetration via the magnetopause; (d) wave generation by plasma clouds inside the magnetopause.

We use MHD simulations to study the time sequence of magnetospheric responses to a synthetic event with a southward interplanetary magnetic field (IMF) turning. The onset of dayside magnetopause reconnection launches a weak rarefaction wave and sunward flow in the equatorial magnetosphere simultaneously with a tailward flow through the polar cap. This convection results in the accumulation of magnetic flux in the tail lobes and thinning of the tail current layer which provides favorable conditions for the onset of nightside reconnection. The onset of nightside reconnection about 40 min later closes the Dungey convection cycle, resulting in a second increase in the sunward flow in the equatorial plane. Variations of the magnetopause standoff distance as well as the size of the polar cap (PC) may indicate the onsets of the dayside and nightside reconnections. We compare the results of two MHD models and discuss their differences.

Utilizing magnetic field measurements made by the Iridium satellites and by ground magnetometers in North America we calculate the full ionospheric current system and investigate the substorm current wedge. The current estimates are independent of ionospheric conductance, and are based on estimates of the divergence-free (DF) ionospheric current from ground magnetometers and curl-free (CF) ionospheric currents from Iridium. The DF and CF currents are represented using spherical elementary current systems (SECS), derived using a new inversion scheme that ensures the current systems' spatial scales are consistent. We present 18 substorm events and find a typical substorm current wedge (SCW) in 12 events. Our investigation of these substorms shows that during substorm expansion, equivalent field-aligned currents (EFACs) derived with ground magnetometers are a poor proxy of the actual FAC. We also find that the intensification of the westward electrojet can occur without an intensification of the FACs. We present theoretical investigations that show that the observed deviation between FACs estimated with satellite measurements and ground-based EFACs are consistent with the presence of a strong local enhancement of the ionospheric conductance, similar to the substorm bulge. Such enhancements of the auroral conductance can also change the ionospheric closure of pre-existing FACs such that the ground magnetic field, and in particular the westward electrojet, changes significantly. These results demonstrate that attributing intensification of the westward electrojet to SCW current closure can yield false understanding of the ionospheric and magnetospheric state.

Using Time History of Events and Macroscale Interactions during Substorms (THEMIS) data, we studied the stepwise development in high-latitude geomagnetic perturbations and Pi1 and Pi2 pulsations during substorm onsets and their association with stepwise auroral onset arc development by analyzing four substorm events. We found that the geomagnetic perturbations and pulsations which are magnetic signatures of the substorm on the ground show stepwise changes and excitation similar to the development of the auroral onset arc which is the visual manifestation of the substorm. We observed minor to small changes in magnetic perturbations and excitation of Pi2 pulsations before initial brightening (IB), and the subsequent excitation of Pi1 and the second Pi2 at or around the further enhancement of onset arc (FE). Then, a steep fall in the magnetic northward component, and the largest-amplitude and highest-frequency Pi1 and Pi2 pulsations appeared at or after poleward expansion (PE). The appearance of FE in all four events and its association with magnetic perturbations and pulsations suggest that FE is an important step in addition to IB and PE. The detailed analysis of the FE step using ground- and space-based data may provide information on the substorm triggering mechanism, the sequence of mechanisms behind the substorm, as well as the mechanisms responsible for the excitation of Pi1 and Pi2 pulsations.

Energy transfer and transport in the terrestrial magnetotails are primarily driven by dipolarization fronts (DFs) embedded inside plasma jets. The DF-driven energy transfer has hitherto been believed to occur locally at the fronts. Different from the traditional knowledge, here we present the first observation of persistent energy conversion extended far behind a DF. The persistent energy conversion, which was dominated by energy loads and mainly contributed by electron currents, developed inside a turbulent, decaying flux pileup region (FPR), nearly 10 d DF (DF’s thickness) behind the DF. The energy transfer chain may be initiated by interaction between the ion flow and ambient plasmas and closed by electron dynamics, leading to electron acceleration perpendicular to magnetic field. These results highlight that electron physics in turbulent FPRs plays a crucial role in the energy transport in the planetary magnetospheres.

Observational data sets for the high latitude middle atmosphere are key to understand the dynamics over those latitudes and the coupling between the lower and middle atmosphere. Utilizing long-term data sets from an all-sky imager at Tromsø, Norway (69.6°N, 19.2°E), the characteristics of 18 mesospheric frontal events in the Arctic winter mesosphere from 2011 to 2015 were studied. These frontal events exhibit horizontal extensions exceeding 500 km and were characterized by a sharp leading front, sometimes followed by a quasi-monochromatic wave train or a turbulent region. A subset of these frontal gravity wave events has been identified in the past as “bores.” While there have been numerous previous reports from low- and mid-latitude sites, and also from southern high latitudes, there have been a few from northern high latitudes. This study focuses on the frontal events in the northern high latitudes and provides new insights into the characteristics of these events. Their horizontal wavelengths primarily ranged from 20 to 40 km, and they exhibited phase speeds in the range 30–80 m/s. Most events were observed before local midnight. No clear link between these events and auroral activity was found. The majority of fronts were found propagating in the north-west direction, which might be due to the wind filtering effects.

Utilizing the multi-point observations by Van Allen Probe A, GOES 13 and 15, we analyzed the competing influences of earthward convection and azimuthal drift loss on the pitch angle distributions of energetic electrons during the simultaneous increases in solar wind flow velocity and pressure. The increase in solar wind speed amplifies the dawn-dusk convection electric field and causes the earthward transport of energetic electrons, and meanwhile the enhancement of solar wind dynamic pressure causes the inward displacement of dayside magnetopause and triggers the azimuthal drift loss of energetic electrons. The earthward convection of low-energy electrons (<60 keV) is much faster than their azimuthal drift loss at most pitch angles, and the fast earthward convections make the butterfly-like electron pitch angle distributions formed early become pancake-like distributions. The 60–530 keV electrons maintain the butterfly-like pitch angle distributions during the earthward convections, whereas the high-energy electrons above 530 keV are not transported to the low-L shells because of fast drift loss in the high-L source region. The competition between the earthward convection and the azimuthal drift loss finally determines the pitch angle distributions of energetic electrons near the trapping boundary during the increases in solar wind flow speed and pressure.

High energy resolution DEMETER satellite observations from the Instrument for the Detection of Particle (IDP) are analyzed during an electromagnetic ion cyclotron (EMIC)-induced electron precipitation event. Analysis of an Interval Pulsation with Diminishing Periods (IPDP)-type EMIC wave event, using combined satellite observations to correct for incident proton contamination, detected an energy precipitation spectrum ranging from ∼150 keV to ∼1.5 MeV. While inconsistent with many theoretical predictions of >1 MeV EMIC-induced electron precipitation, the finding is consistent with an increasing number of experimentally observed events detected using lower resolution integral channel measurements on the POES, FIREBIRD, and ELFIN satellites. Revised and improved DEMETER differential energy fluxes, after correction for incident proton contamination shows that they agree to within 40% in peak flux magnitude, and 85 keV (within 40%) for the energy at which the peak occurred as calculated from POES integral channel electron precipitation measurements. This work shows that a subset of EMIC waves found close to the plasmapause, that is, IPDP-type rising tone events, can produce electron precipitation with peak energies substantially below 1 MeV. The rising tone features of IPDP EMIC waves, along with the association with the high cold plasma density regime, and the rapidly varying electron density gradients of the plasmapause may be an important factor in the generation of such low energy precipitation, co-incident with a high energy tail. Our work highlights the importance of undertaking proton contamination correction when using the high-resolution DEMETER particle measurements to investigate EMIC-driven electron precipitation.

The present study provides an evidence for the generation of harmonics of magnetosonic waves in the Martian magnetosheath region. The wave signatures are manifested in the magnetic field measurements recorded by the fluxgate magnetometer instrument onboard the Mars Atmosphere and Volatile Evolution missioN (MAVEN) spacecraft in the dawn sector around 5–10 LT at an altitude of 4,000–6,000 kms. The wave that is observed continuously from 19.1 to 20.7 UT below the proton cyclotron frequency (f ci  ≈ 46 mHz) is identified as fundamental mode of the magnetosonic wave. Whereas harmonics of the magnetosonic wave are observed during 19.7–20.3 UT at frequencies that are multiple of f ci . The ambient solar wind proton density and plasma flow velocity are found to vary with a fundamental mode frequency of 46 mHz. It is noticed that the fundamental mode is mainly associated with the left-hand (LH), and higher frequency harmonics are associated with the right-hand (RH) circular polarizations. A clear difference in the polarization and ellipticity is noticed during the time of occurrence of harmonics. The magnetosonic wave harmonics are found to propagate in the quasi-perpendicular directions to the ambient magnetic field. The results of linear theory and Particle-In-Cell simulation performed here are in agreement with the observations. The present study provides a conclusive evidence for the occurrence of harmonics of magnetosonic wave in the close vicinity of the magnetosheath region of the unmagnetized planet Mars.

Solar eclipse traveled across South China in the afternoon on 21 June 2020. Five ionosondes located from mid-to low-latitudes and on both north and south of the eclipse path were applied to investigate the ionospheric responses. Both the zonal and meridional ranges of the observation region have exceeded 1,000 km. All the five ionosondes had observed the Intermediate Descending Layers (IDLs) simultaneously just after the eclipse maximum and this is a very small probability event. During the solar eclipse, the multi-hop echoes above the Es, the rising Es to 150 km altitude, the plasma flux from above F2-layer were also observed and analyzed. The descending trend of the IDLs and the peak height of F2-layer (h m F2) shows great consistency, indicating the close relationship between the eclipse induced plasma flux and the IDLs. The traces of gravity waves were also found in the IDLs and F-layer. The plasma flux may carry the ions to valley region and the eclipse produced gravity waves were responsible for the formation of the IDLs.

Ultra-low frequency (ULF) waves transfer energy and momentum into the ionosphere-thermosphere system. To quantify this energy, this paper first presents a new method to quantitatively detect ULF waves in Incoherent Scatter Radar (ISR) data based on 2D fast-Fourier transforms and subsequent reconstruction of the wave. In parallel with other data sets, including optical, magnetometer, satellite, and models, we present the first full ionospheric energy dissipation rates for a ULF wave, split into electromagnetic (EM) and kinetic fluxes. The EM energy deposition is calculated from the use of the Poynting theorem, looking at Joule and frictional heating rates, where both rates show the same order of magnitude (1.24 × 1013 and 7.3 × 1012 J) respectively when integrated over the wave lifetime of 2 hr 15 min and an area of 4° magnetic latitude × 74° magnetic longitude. However, contrary to the common assumption that the EM flux is dominant, we determined the kinetic flux, to be almost equal in magnitude (8.7 × 1012 J). This indicates that previous papers might have underestimated the total energy dissipation by ULF waves. Compared to the substorm energy budget, we find that locally, the ULF wave event studied here makes up approximately 10% of a typical substorm cycle budget.

For the first time, characteristics of the geographical and seasonal distribution of the quasi-diurnal lunar O1 tide were derived from a time series of ionospheric total electron content (TEC) maps provided by International Global Navigation Satellite System Service (IGS). The data analysis is focused on solar minimum in 2008 and 2009 where disturbing influences of geomagnetic and solar activity were minimal. We found that the magnitude of the O1 tide is as strong as the “dominant” semidiurnal lunar M2 tide. Relative amplitudes of 10% and larger are observed in some regions for the O1 component in TEC. The O1 component is particularly strong in northern hemispheric winter over the west coast of South America. There, two maxima occur which are northward and southward of the magnetic equator in the Equatorial Ionization Anomaly (EIA) crest regions. Following Yamazaki et al. (2017, https://doi.org/10.1002/2017ja024601), it might be assumed that a longitudinal anomaly of ionospheric conductivities in the Peruvian sector leads to a stronger modulation of the equatorial electrojet by the lunar tides. Electrodynamic lifting of plasma and transport to the EIA crests may explain the variations of the O1 component in TEC. Contrary to many studies, we find the O1 component (period 25.82 hr) more important than the M1 component (period 24.84 hr, a lunar day). We show that the geographical distribution of the O1 component is totally different from that of the M1 component which is smaller. The seasonal variation of O1 shows maximal amplitudes in northern hemispheric winter and minimal amplitudes in southern hemispheric winter.