JGR:Space physics

In this work, IMF B y effects on field-aligned currents (FACs) are examined in different local time sectors, seasons, and hemispheres. At dusk and 09–14 MLT, when the eastward polar electrojet (PEJ) prevails, the northern FACp (poleward side FACs) are stronger when IMF B y  < 0 than when IMF B y  > 0. Conversely, at dawn, 21–02 MLT, and 09–14 MLT with westward PEJ, the northern FACp are stronger with IMF B y  > 0 compared to IMF B y  < 0. The southern FACp shows a reversed relationship with IMF B y direction. The dependence of FACe (equatorward side FACs) on IMF B y is weaker, except for the midday FACe, which shows opposite variations with respect to IMF B y when compared to FACp. Stronger IMF B y effect is observed in local summer in most of local times. The northern FACs are located at higher latitude for IMF B y  > 0 than for IMF B y  < 0 in local times with eastward PEJ, while the opposite trend is observed in other local times and in the Southern Hemisphere. The hemispheric difference in the peak latitude of FACs demonstrates an inverse relationship with its intensity, with stronger FACs located at lower latitudes. Overall, the local time and hemispheric differences in FACs strength and latitude are discussed in the context of interhemispheric field-aligned currents linked to IMF B y .

This paper conducts a multi-instrument analysis of a latitudinal plasma density peak at the middle latitudes during the early recovery phase of the April 2023 geomagnetic storm. The total electron content (TEC), peak density of the F layer, and the in situ plasma density from Swarm and Defense Meteorological Satellite Program (DMSP) satellites all capture this peak feature. This narrow latitudinal peak structure appeared around 50°N and extended from 40°E to 150°E in longitude with a prolonged duration of about 8 hr from sunset to midnight. This mid-latitude peak reveals a noticeable equatorward motion and a slight westward shift. According to the plasma composition observations from DMSP satellites, this peak structure shows an O+ ions dominance, which means that this peak is more likely to be formed by an internal rather than an external source from the plasmasphere. Meanwhile, the middle latitude Fabry–Perot interferometer (FPI) observed strong equatorward thermospheric winds, and the peak height of the F layer presented a visible elevation, which suggests that the equatorward wind lifting caused the plasma density enhancement. Besides, the O/N2 ratio significantly decreased at lower and middle latitudes, and ion drift observations showed a distinct subauroral westward channel. Based on these simultaneous measurements, this structure's sharp equatorward and poleward boundaries might be related to the O/N2 ratio change and the subauroral polarization stream (SAPS) flow separately.

The ion foreshock, filled with backstreaming foreshock ions, is very dynamic with many transient structures that disturb the bow shock and the magnetosphere-ionosphere system. It has been shown that foreshock ions can be generated through either solar wind reflection at the bow shock or leakage from the magnetosheath. While solar wind reflection is widely believed to be the dominant generation process, our investigation using Time History of Events and Macroscale Interactions during Substorms mission observations reveals that the relative importance of magnetosheath leakage has been underestimated. We show from case studies that when the magnetosheath ions exhibit field-aligned anisotropy, a large fraction of them attains sufficient field-aligned speed to escape upstream, resulting in very high foreshock ion density. The observed foreshock ion density, velocity, phase space density, and distribution function shape are consistent with such an escape or leakage process. Our results suggest that magnetosheath leakage could be a significant contributor to the formation of the ion foreshock. Further characterization of the magnetosheath leakage process is a critical step toward building predictive models of the ion foreshock, a necessary step to better forecast foreshock-driven space weather effects.

We conducted an all-sky imaging transient search with the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) data collected during the Perseid meteor shower in 2018. The data collection during the meteor shower was motivated to conduct a search for intrinsic radio emission from meteors below 60 MHz known as the meteor radio afterglows (MRAs). The data collected were calibrated and imaged using the core array to obtain lower angular resolution images of the sky. These images were input to a pre-existing LWA transient search pipeline to search for MRAs as well as cosmic radio transients. This search detected 5 MRAs and did not find any cosmic transients. We further conducted peeling of bright sources, near-field correction, visibility differencing and higher angular resolution imaging using the full array for these 5 MRAs. These higher angular resolution images were used to study their plasma emission structures and monitor their evolution as a function of frequency and time. With higher angular resolution imaging, we resolved the radio emission size scales to less than 1 km physical size at 100 km heights. The spectral index mapping of one of the long duration event showed signs of diffusion of plasma within the meteor trails. The unpolarized emission from the resolved radio components suggest resonant transition radiation as the possible radiation mechanism of MRAs.

Mercury is the only planetary magnetosphere which does not possess a significantly conducting ionosphere, yet magnetic field observations by the Mariner 10 and MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) missions revealed a stable presence of large-scale field-aligned currents (FACs). Several empirical magnetic field models have been developed to describe Mercury’s magnetospheric currents and their associated magnetic fields, but none attempted to include an explicit description of large-scale FACs. Here, we describe a dynamic FAC magnetic field model for Mercury based on the analytical solution to conical currents, which are shielded by a dynamic magnetopause boundary. The model contains free parameters setting the FACs' latitudinal extent and amplitude as functions of magnetospheric activity. The parameters are fit by minimizing the root-mean-square (RMS) differences between the model and MESSENGER Magnetometer data. During magnetically quiet conditions, the modeled FACs have an intensity of ∼10 nA/m2 and extend in latitude from ∼83°N to ∼40°N. They intensify to ∼20 nA/m2 and expand equatorward to ∼28°N during the most active times. The inclusion of the FAC model reduces low-altitude RMS residuals by 8.1% when compared to prior models. The model effectively captures the azimuthal component of the magnetic field present in the MESSENGER low-altitude data, but largely misses the radial and co-latitudinal components, indicating other large-scale physics remain missing in the mathematical descriptions of Mercury's magnetosphere.

The Plasma Environment Modeling in the Earth's Magnetosphere (PEMEM) is a European Space Agency activity supporting the development of a new specification model for the spacecraft surface charging risk assessment. This paper presents a description of the basic model version: the PEMEM percentile model. The model is intended to be used for space missions with near-equatorial orbits. The model is based on the Van Allen Probes particle measurements inside the geostationary orbit. The model's primary input is a planned spacecraft trajectory. It outputs statistical characteristics of the plasma environment which are expected to be encountered during a mission lifetime. These characteristics include differential electron and proton flux percentiles for a set of energies (percentile spectra), and percentiles of the integrated electron flux. The model covers the energy range of 1–100 keV for electrons and 40 eV–51 keV for protons. Since extreme spacecraft charging usually occurs in the eclipse, the same characteristics can be separately output for the periods when the spacecraft is shadowed by the Earth.

Chorus waves are a kind of intense electromagnetic emission wave in magnetized planets and can play important roles in the kinetic electron dynamics in planetary magnetospheres. Rapid changes of the ring electron current belt in Mercury’s magnetosphere and the contribution of chorus waves have remained long-standing scientific issues from the first Mercury flyby observations by Mariner 10 in 1970s because of the small size of the magnetosphere. Based on theoretical analyses and simulations successfully reconstructing Earth’s chorus wave properties, we report on possible generation regions of chorus waves in Mercury’s magnetosphere. The theoretical analysis for low-temperature-anisotropy electrons shows a clear asymmetric day–night spatial distribution of the possible chorus generation region because of the difference in the nonlinear convective wave growth along the magnetic field lines. Simulation results show a rapid enhancement of the ring electron current belt by resonant interactions with repetitive chorus waves. Our study suggests that energetic electrons in Mercury’s magnetosphere can be enhanced locally by nonlinear chorus wave–particle interactions.

We show that Traveling Ionospheric Disturbances (TIDs) may affect the Farley Buneman Instability (FBI) and Gradient Drift Instability (GDI) echoes referred to as the Near Range Echoes (NREs) in the SuperDARN radar backscatter from the lower part of the E-region. TIDs and NREs are observed concomitantly by the Zhongshan SuperDARN radar (69.38°S, 76.38°E) in the far and near ranges, respectively. At the moment, there is no study about the effects of TIDs on the NREs caused by the FBI using the SuperDARN radars. The GDI are more likely to occur at a lower altitude while FBI occurs at a slightly higher altitude in the lower part of the ionospheric E-region. We use the Spearman Correlation Coefficient (SCC) to show that a part of the NREs backscatter power could be statistically explained by the MSTIDs backscatter power received by the same radar. We also investigate the simultaneous occurrence rate of the NREs and MSTIDs during the 24th solar cycle. Seasonal variability shows that MSTIDs-NREs events over Zhongshan mostly occur in summer and equinoxes during local night and morning. The majority of these events lasted between ∼4 and ∼8 hr. Most events disappeared early in the morning. Statistics of the Spearman correlation coefficient values show that ∼9% of NRE amplitude modulation could be due to the MSTIDs. There are almost equal numbers of negative and positive Spearman correlation coefficient values. The relative velocity between the E-region NREs and the F-region MSTIDs switching the electric field polarities between the crests and troughs could be the cause of those equal number of the Spearman correlation coefficient values. The orientation of the ionospheric current relative to the MSTID polarization electric field may also play a significant role in the reported Spearman correlation coefficient values. We argue that in some cases, the TIDs might have been close enough to the NREs altitude to modulate them directly by transporting the plasma up and down through shear or compression.

In planetary radiation belts, the Kennel-Petschek flux limit is expected to set an upper limit on trapped electron fluxes at 80–600 keV in the presence of efficient electron loss through pitch-angle diffusion by whistler-mode chorus waves generated around the magnetic equator by the same 80–600 keV electron population. Comparisons with maximum measured fluxes have been relatively successful, but several key assumptions of the Kennel-Petschek model have not been experimentally tested. The Kennel-Petschek model notably assumes an exponential growth of chorus waves as the trapped electron flux increases, and a fixed maximum wave power gain of about 3. Here, we describe a method for inferring the near-equatorial wave power gain using only measurements of trapped, precipitating, and backscattered electron fluxes at low altitude. Next, we make use of Electron Losses and Fields Investigation (ELFIN) CubeSats measurements of such electron fluxes during two moderate geomagnetic storms with sustained electron injections to infer the corresponding chorus wave power gains as a function of time, energy, and equatorial trapped electron flux. We show that wave power increases exponentially with trapped flux, with a wave power gain roughly proportional to the theoretical linear convective gain, and that the maximum inferred gain near the upper flux limit is roughly 10, with a factor of 2 uncertainty. Therefore, two key theoretical underpinnings of the Kennel-Petschek model are borne out by the present results, although the strong inferred gains should correspond to higher flux limits than in traditional estimates.

Precipitation of relativistic electrons into the Earth's atmosphere regulates the outer radiation belt fluxes and contributes to magnetosphere-atmosphere coupling. One of the main drivers of such precipitation is electron scattering by whistler-mode waves. Such waves typically originate at the equator, where they can resonate with and scatter sub-relativistic (tens to a few hundred keV) electrons. However, they can occasionally propagate far away from the equator along field lines, reaching middle latitudes, where they can resonate with and scatter relativistic (>500 keV) electrons. Such a propagation is typical for the dayside, but statistically has not been found on the nightside where the waves are quickly damped along their propagation due to Landau damping. Here we explore two events of relativistic electron precipitation from low-altitude observations on the nightside. Combining measurements of whistler-mode waves from ground observatories, relativistic electron precipitation from low-altitude satellites, total electron content maps from GPS receivers, and magnetic field and electron flux from equatorial satellites, we show wave ducting by plasma density gradients is the possible channel that allows the waves to reach middle latitudes and scatter relativistic electrons. We suggest that both whistler-mode wave generation and ducting can be driven by equatorial mesoscale (with spatial scales of about one Earth radius) transient structures during nightside injections. We also compare these nightside events with observations of ducted waves and relativistic electron precipitation at the dayside, where wave generation and ducting are driven by ultra-low-frequency waves. This study demonstrates the potential importance of mesoscale transients in relativistic electron precipitation, but does not however unequivocally establish that ducted whistler-mode waves are the primary cause of the observed electron precipitation.

As activity in Earth orbit continues to grow, it is important to characterize the environment of near-Earth space. One means of remotely sensing lower thermospheric neutrals is by measurement of O and N 2 density through the observation of far-ultraviolet (FUV) airglow of atomic oxygen at 135.6 nm and the N 2 Lyman-Birge-Hopfield (LBH) bands (~130–180 nm), as has been done on the Ionospheric Connection Explorer (ICON), Global-scale Observations of the Limb and Disk (GOLD), and Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) missions. This technique is not without limitations, however, as the FUV measurements suffer from contamination by ionospheric emissions at low latitudes and auroral emissions excited by precipitating energetic electrons and protons at high latitudes. Previous work has shown the potential for making measurements of O and N 2 density in the lower-middle thermosphere using observations of extreme-ultraviolet (EUV) airglow. This measurement approach has a potential advantage in that it does not have an inherent ionospheric emission that must be accounted for. Additionally, these emissions are primarily excited directly by solar UV rather than electron impact and thus have the potential to enable expansion of neutral density observations into the auroral zone and polar cap where the FUV measurement cannot be applied. This article demonstrates a new approach and algorithm designed to retrieve thermospheric O and N 2 density from 150 to 400 km using measurements from the ICON EUV instrument. The retrieval results throughout 2020 are summarized and compared to measurements from ICON FUV, GOLD, and SWARM.

The latitudinal location of the Equatorial Ionization Anomaly (EIA) crest has seasonal variation, and there are disagreements on the interpretation of such seasonal characteristic in previous studies. Some studies suggested that this seasonal characteristic is determined by the seasonal characteristic of the equatorial electric field. Others suggested that this seasonal characteristic is determined by the seasonal changes of the thermospheric wind. The current paper uses Total Electron Content (TEC) data and the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM) to analyze the seasonal variation of the northern EIA crest in the eastern Asian sector under low solar activity. Our results show that the monthly averaged latitudinal location of the northern EIA crest has a good linear relationship (r = 0.74) with the monthly averaged Equatorial Electrojet (EEJ) intensity, which is a commonly used proxy of the eastward electric field. However, TIEGCM simulations with and without F-region wind indicate that such a relationship might be attributed to wind effects. Additionally, the linear relationship between the EEJ intensity and the northern-southern EIA crest distance is not significant (r = 0.47) in the eastern Asian sector. Our results suggest that a good correspondence between the eastward electric field and the latitudinal location of the EIA crest is not assured annually, as the seasonally varying F-region wind significantly influences EIA evolution.

In the mesosphere and lower thermosphere (MLT) region, semidiurnal and terdiurnal tides are dominant at middle and high latitudes and are important for the dynamics and structures. Using the global neutral horizontal wind data in the MLT region observed by TIMED Doppler interferometer from 2002 to 2021, the seasonal and interannual variations of six semidiurnal and terdiurnal tidal components (SW2, SW3, SW4, TW3, TW2 and TW4) are investigated. Particularly, the responses of these tidal components to the stratospheric quasi-biennial oscillation (SQBO) and the solar cycle are presented. The results indicate that: (a) seasonal and interannual variations of these tidal components are pronounced in their peak regions. The peak region is defined as the region where the ratio of the annual mean amplitude to the maximum annual mean amplitude is larger than 0.8. (b) In their peak regions, migrating tidal components exhibit strong seasonal variations at annual (annual oscillations (AO)), semiannual (semiannual oscillations (SAO)) and terannual (terannual oscillations (TAO)) periods, whereas nonmigrating tidal components only exhibit prominent AO and SAO. (c) In their peak regions, the responses of the annual mean amplitudes of these tidal components to SQBO are negative; the responses of migrating tidal components to solar cycle are negative whereas those of nonmigrating tidal components are positive. (d) For these tidal components, the responses of the amplitudes of seasonal variations to SQBO and solar cycle are weaker than their annual mean amplitudes, and the latitude × altitude response patterns are not always consistent with those of the annual mean amplitudes.

We compared the performance of DREAM3D simulations in reproducing the long-term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low-energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (D LL ), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics inside L* = 6, with a better model performance at lower μ than higher μ, where μ is the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data-driven outer boundary condition is critical to the model performance, while adding in the data-driven seed population doesn't further improve the performance. (b) The model shows comparable performance with D LL from Brautigam and Albert (2000, https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014, https://doi.org/10.1002/2013ja019204), and Liu et al. (2016, https://doi.org/10.1002/2015gl067398), while with D LL from Ali et al. (2016, https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data-based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higher μ using the Wang et al. (2019, https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under-reproduced in the current model.

This study uses over 2 years of 16 Hz density measurements, 50 Hz magnetic field data and ROTI data from the Swarm mission to perform long term statistics of plasma structuring in the polar ionosphere. The timeframe covers more than 2 years near the 24th solar cycle peak. We additionally use 3 years of data obtained from a timeframe close to solar minimum for discussion. We present power spectral densities (PSD) of electron density irregularities and magnetic field for 1-min intervals. These PSD have been characterized by the probability of a slope steepening, and by integrating the power deposited within frequency intervals corresponding to kilometer scales. For the electron density, we observe seasonal dependencies for both the integrated power and slope characteristics. While the dual slope probability, especially within the polar cap, varies with solar EUV-radiation, the integrated power is strongest around the equinoxes. Additionally, while we found similar results for the slope probability for both hemispheres, the integrated power exhibits strong hemispheric asymmetries with stronger enhancements within local summer in the southern hemisphere. The ROTI data shows a similar seasonal variability as the density PSD integrated power, in both seasonal dependency and interhemispheric variability. However, for the ROTI data the strongest fluctuations were found within the nightside auroral oval and the cusp. For the PSD of the magnetic field data, we obtain the strongest enhancements within the cusp for all seasons and all hemispheres. The fluctuations may indicate an increase in Alfvénic energy associated with a downward Poynting flux.

On 26–27 December 2022, Mars experienced an extremely low-density solar wind stream, which was encountered first by Earth because of the radial alignment of the two planets (i.e., Mars opposition). During this event, two important properties of the ionospheric and magnetospheric states changed significantly in response to the low solar wind ram pressure, as inferred from the superthermal electron observations from the Mars Atmospheric and Volatile EvolutioN (MAVEN) mission. The interface between the ionosphere and magnetosphere expanded to thousands of kilometers, outside of the nominal bow shock locations, coinciding with the expansion of the cold planetary ions. Meanwhile, the ambipolar electrostatic potential arising from the ionospheric electron pressure gradient increased from the nominal ∼ −0.7 to ∼ −2 V (relative to the lower ionosphere). This enhanced ambipolar potential likely facilitated the observed ionosphere expansion.

The moving solar terminator (ST) generates atmospheric disturbances, broadly termed solar terminator waves (STWs). Despite theoretically recurring daily, STWs remain poorly understood, partially due to measurement challenges near the ST. Analyzing Michelson Interferometer for Global High-resolution Thermospheric Imaging (MIGHTI) data from NASA's Ionospheric Connection Explorer (ICON) observatory, we present observations of STW signatures in thermospheric neutral winds, including the first reported meridional wind signatures. Seasonal analysis reveals STWs are most prominent during solstices, when they intersect the ST about ∼20° latitude from the equator in the winter hemisphere and have phase fronts inclined at a ∼40° angle to the ST. We also provide the first observed STW altitude profiles, revealing large vertical wavelengths above 200 km. Comparing these observations to four different models suggests the STWs likely originate directly or indirectly from waves from below 97 km. STWs may play an under-recognized role in the daily variability of the thermosphere-ionosphere system, warranting further study.

In a case study using Van Allen Probe B, we investigate chorus wave observations near the western edge of a plasmaspheric plume characterized by steep density gradients. Initially, wave vectors are oriented anti-Earthward, but they become very oblique and eastward as the probe approaches the plume boundary. Treating the plume boundary as an azimuthal density gradient, ray tracing can reproduce the observed wave vector directions. Ray tracing shows that the azimuthal density gradient strongly inclines the wave vectors eastward. Consequently, waves are reflected upon reaching the Gendrin angle and cannot enter the plume. We establish an analytical criterion for the azimuthal density enhancement, determining the condition for chorus waves to enter plumes near the equatorial region. Our results partly explain the oblique chorus near plumes observed by Hartley, Chen et al. (2022, https://doi.org/10.1029/2022GL098710), offering insight into wave-particle interactions by chorus waves with the influence of azimuthal density structures.

We examine a 6-day traversal of the magnetotail by the ARTEMIS satellites during an interval of prolonged northward IMF. The electrostatic analyzer (ESA) onboard the ARTEMIS spacecraft measures high ion and electron fluxes at approximately 60 R E downtail in regions of the magnetotail which would normally be the magnetotail lobe, containing open flux evacuated of plasma. We interpret these observations as trapped plasma on closed magnetic flux indicating that the magnetotail is closed or partially closed but extends at least as far as ∼60 R E downtail. We find that the occurrence of plasma in the magnetotail and the closure of the magnetosphere results in distinct changes to the magnetotail structure including a reduction in the magnetic field strength and pressure as well as a narrowing of the tail by approximately 20 R E .

Dayside magnetosphere interactions are essential for energy and momentum transport between the solar wind and the magnetosphere. In this study, we investigate a new phenomenon within this regime. Sudden enhancements of ion fluxes followed by repeating dropouts and recoveries were observed by Magnetospheric Multiscale on 5 November 2016, which is the very end of the recovery phase from a moderate geomagnetic storm. These repetitive flux variations display energy-dispersive characteristics with periods relevant to ion bounce motion, suggesting they are corresponding echoes. Alongside the flux variations, bipolar electric field impulses originating from external sources were detected. We traced the source region of the initial injection and found it is located near the spacecraft's position. To elucidate the underlying physics, a test-particle simulation is conducted. The results reveal that radial transport resulting from impulse-induced acceleration can give rise to these echoes. Observations demonstrate dayside magnetosphere interactions are more common than we previously considered, which warrants further research.