JGR:Space physics

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.

Flux transfer events (FTEs) are a type of magnetospheric phenomena that exhibit distinctive observational signatures from the in situ spacecraft measurements. They are generally believed to possess a magnetic field configuration of a magnetic flux rope and formed through magnetic reconnection at the dayside magnetopause, sometimes accompanied with enhanced plasma convection in the ionosphere. We examine two FTE intervals under the condition of southward interplanetary magnetic field (IMF) with a dawn-dusk component. We apply the Grad-Shafranov (GS) reconstruction method to the in situ measurements by the Magnetospheric Multiscale (MMS) spacecraft to derive the magnetic flux contents associated with the FTE flux ropes. In particular, given a cylindrical magnetic flux rope configuration derived from the GS reconstruction, the magnetic flux content can be characterized by both the toroidal (axial) and poloidal fluxes. We then estimate the amount of magnetic flux (i.e., the reconnection flux) encompassed by the area “opened” in the ionosphere, based on the ground-based Super Dual Auroral Radar Network (SuperDARN) observations. We find that for event 1, the FTE flux rope is oriented in the approximate dawn-dusk direction, and the amount of its total poloidal magnetic flux falls within the range of the corresponding reconnection flux. For event 2, the FTE flux rope is oriented in the north-south direction. Both the FTE flux and the reconnection flux have greater uncertainty. We provide a detailed description about a formation scenario of sequential magnetic reconnection between adjacent field lines based on the FTE flux rope configurations from our results.

Energetic electron dynamics in the Earth's radiation belts and near-Earth plasma sheet are controlled by multiple processes operating on very different time scales: from storm-time magnetic field reconfiguration on a timescale of hours to individual resonant wave-particle interactions on a timescale of milliseconds. The most advanced models for such dynamics either include test particle simulations in electromagnetic fields from global magnetospheric models, or those that solve the Fokker-Plank equation for long-term effects of wave-particle resonant interactions. The most prospective method, however, would be to combine these two classes of models, to allow the inclusion of resonant electron scattering into simulations of electron motion in global magnetospheric fields. However, there are still significant outstanding challenges that remain regarding how to incorporate the long term effects of wave-particle interactions in test-particle simulations. In this paper, we describe in details two approaches that incorporate electron scattering in test particle simulations: stochastic differential equation (SDE) approach and the mapping technique. Both approaches assume that wave-particle interactions can be described as a probabilistic process that changes electron energy, pitch-angle, and thus modifies the test particle dynamics. To compare these approaches, we model electron resonant interactions with field-aligned whistler-mode waves in dipole magnetic fields. This comparison shows advantages of the mapping technique in simulating the nonlinear resonant effects, but also underlines that more significant computational resources are needed for this technique in comparison with the SDE approach. We further discuss applications of both approaches in improving existing models of energetic electron dynamics.

We present a case study of plasma and magnetic field observations in the Martian magnetotail using data from the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission during an orbit when the spacecraft was in the optical shadow, past the dusk terminator and downstream of the strongest crustal magnetic fields. In this region, we observed multiple magnetic field rotations (a signature of currents) closely associated with energized (up to 100 eV) electron populations. Several transitions between closed and draped magnetic topologies also occur in this region, which are likely to be caused by magnetic reconnection between the interplanetary magnetic field (IMF) and crustal magnetic fields. We also observe two regions of energized, counter-streaming electrons, which are rare in the magnetotail, but twice as likely to occur downstream of strong crustal magnetic fields when they are near the evening terminator. Together, the multiple magnetic field rotations, topological changes, and counter streaming electrons suggest the presence of an electric potential structure similar to those observed above the auroral arc regions at Earth.

The exchange of kinetic and electromagnetic energy by precipitation and/or outflow and through field-aligned currents are two aspects of the ionosphere-magnetosphere coupling. A thorough investigation of these processes is required to better understand magnetospheric dynamics. Building on our previous study using the Defense Meteorological Satellite Program spectrometer data, here we use Swarm vector field magnetometer data to describe the auroral zone morphology in terms of east-west magnetic field perturbations. We define a threshold for detecting magnetic fluctuations based on the power spectral density of ΔB EW, and derive the disturbed magnetic field occurrence probability (dBOP) at low [0.1–1 Hz] and high [2.5–5 Hz] frequencies. High-frequency distributions of dBOP reveal a dayside-nightside asymmetry, whereas low-frequency dBOP exhibits a persistent morphological asymmetry between the dawn-to-noon and the dusk-to-midnight sectors, peaking at dawn. Notably, weak solar wind conditions are associated with an increase in the dBOP asymmetric patterns. At low frequency in particular, while the dBOP seems to be primarily constant at dawn, the dusk dBOP decreases during quiet times, inducing a relatively larger dawn-dusk asymmetry in such conditions. Furthermore, based on a comparison analysis between low-frequency dBOP and auroral electron precipitation occurrence probability, we suggest that both types of distribution offer a more comprehensive understanding of the large-scale auroral zone when considered concurrently. Our interpretation is that the dBOP at low frequencies reflects a quasi-steady state circulation of energy, while the high-frequency dBOP reflects the regions of rapid changes in the magnetosphere. The dBOP is therefore a crucial source of information regarding the magnetosphere-ionosphere coupling.

Energetic particle precipitation is the major source of electron production that controls the ionospheric Pedersen and Hall conductances at high latitudes. Many studies use empirical formulas to estimate conductances. The particle precipitation spectra measured by the Defense Meteorological Satellite Program (DMSP) Special Sensor J are often used as the input to the empirical formulas. In this study, we evaluate the empirical formulas of ionospheric conductances during four different types of auroral precipitation conditions based on 63 conjugate events observed by DMSP and EISCAT. The conductances calculated from the DMSP data with the empirical formulas are compared with those based on EISCAT measurements with the standard equations. The best correlation between these two is found when the empirical Robinson formulas (Robinson et al., 1987, https://doi.org/10.1029/ja092ia03p02565) are used in the presence of diffuse electron precipitation without ions. In the presence of ion precipitation, the correlation coefficients are smaller, but the correlation improves when the Galand formulas (Galand & Richmond, 2001, https://doi.org/10.1029/1999ja002001) are used to estimate the contribution of ion precipitation to the conductances. We also found that pure ion precipitation can cause the increase of conductances up to 2–7 S for Pedersen and 2.5–10 S for Hall conductances, which is positively correlated with the auroral electrojet index. Overall, the empirical formulas applied to the DMSP particle spectra underestimate the ionospheric conductances.

The turbulent foreshock region upstream of the quasi-parallel bow shock is dominated by waves and reflected particles that interact with each other and create a large number of different foreshock transients. The structures with the enhanced magnetic field, Short Large Amplitude Magnetic Structures, and density spikes named plasmoids are frequently observed. They are one of the suggested sources of transient flux enhancements or jets in the magnetosheath. Using measurements of the Magnetospheric Multiscale Spacecraft (MMS) and OMNI solar wind database between the 2015 and 2018 years, we have found that there is a category of events exhibiting both magnetic field and density enhancements simultaneously and we introduce the term “mixed structures” for them. Consequently, we divided our set of observations into three groups of events and present their comparative statistical analysis in the subsolar foreshock. Based on our results and previous research, we discuss the properties, possible origin and occurrence rates of these events under different upstream conditions and their possible relation to the jet and plasmoid formation in the magnetosheath.

We produced a database of over 41,000 ionospheric cusp locations using 40 years of energetic particle measurements from 14 Defense Meteorological Satellite Program (DMSP) satellites. We limited the database to periods when the Auroral Electrojet (AE) was <100 nT and the Interplanetary Magnetic Field (IMF) measurements were available, then calculated the magnetic latitude (MLAT), magnetic local time (MLT) and the dipole tilt (Λ) for each boundary, multiplying Λ by −1 in the southern hemisphere. We then binned the data in running one-hour bins from 8.5 to 15.5 MLT. To obtain the dependence of the MLAT of the cusp boundary on Λ and IMF By, we performed linear fits on the MLAT versus Λ distributions, separated by ±IMF By. The dependence (degrees MLAT/degrees Λ) is asymmetric in both MLT and hemisphere as well as IMF By. In the northern hemisphere, the dependence for By > 0 (By < 0) is greatest in the afternoon (morning) sector for both poleward and equatorward boundaries. The opposite is true in the southern hemisphere. However, in the noon sector, the dependence is nearly the same for all boundaries and By sign. In addition, separation by Λ sign shows dramatically greater dependence variance for Λ > 0 than for Λ < 0 in the afternoon and morning sectors (but not in the noon sector). The observed asymmetries revealed in this work are thought to be due to the effect of By on reconnection and the location of the reconnection X line.

The magnetotail lobe region at Mercury is characterized by low plasma density and low magnetic field variability compared to the nightside magnetosheath and central plasma sheet. At Mercury, as well as other planets, lobe magnetic fields play a crucial role in storing and releasing magnetic flux in response to changing upstream solar wind conditions such as interplanetary magnetic field (IMF) orientation and solar wind dynamic pressure (P dyn ). This makes the region significant for studying the magnetospheric interaction with the intense solar wind conditions at Mercury's orbit. Here, we identify and analyze magnetotail lobe observations made by the Mercury Surface, Space Environment, Geochemistry and Ranging (MESSENGER) spacecraft during its 4 years orbital phase. We empirically determined a set of criteria using magnetometer (MAG) and the Fast Imaging Plasma Spectrometer instruments onboard MESSENGER to identify lobe magnetic field intervals. From 3,332 MESSENGER orbits, we identify 1,242 lobe field intervals. We derive an expression for the average lobe magnetic field strength in nanotesla with respect to radial distance downtail: B lobe (r) = (135 ± 8) * r (−2.1±0.3) + (31 ± 8). The lobe magnetic field exhibits both small-scale (∼3 min) and orbit-to-orbit (∼8–12 hr) variability in magnetic field strength compared to this averaged field strength expression. The orbit-to-orbit variability in lobe field strength is not significantly correlated with estimated IMF orientation, but is directly correlated with P dyn . Thus, our findings provide evidence for the pressure balance between the inward facing P dyn on the nightside magnetopause and the outward facing magnetic pressure supplied by the lobes.

In the mesosphere and lower thermosphere, diurnal tides are responsible for the dynamics and structures at low latitudes since they have largest amplitudes there. Based on the 20-year (2002−2021) observations from Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics Doppler interferometer (TIDI), we investigate the seasonal variations of three diurnal tidal components (DW1, DE3, DW2), and their responses to stratospheric quasi-biennial oscillation (SQBO) and solar cycle globally. The results show that: (a) DW1, DE3 and DW2 show prominent semiannual (SAO) and annual oscillations (AO) in their peak regions where the ratio of their annual mean amplitude to their maximum annual mean amplitude is larger than 0.8. DW1 also exhibit strong terannual oscillations (TAO) especially in meridional wind. (b) The responses of the amplitudes of seasonal variations of DW1, DE3 and DW2 to SQBO and solar cycle are comparable in magnitude to those of their annual mean amplitudes (for response to solar cycle, even stronger in some cases), and thus cannot be neglected. (c) In their respective peak regions, the responses of annual mean amplitudes of these diurnal tidal components to SQBO are uniformly positive except DW2, and their responses to solar cycle are uniformly negative. For these diurnal tides, the amplitudes of the dominant seasonal variations exhibit consistent response patterns with those of annual means to the SQBO/solar cycle. (d) Empirical formulas are given, which well describe the seasonal and interannual variations of dominant diurnal tidal components in their peak regions.

Joule heating is the primary source of the high latitude thermospheric energy dissipation during geomagnetic storm period. The nitric oxide (NO) emission at 5.3 μm accounts for the majority of Joule heating energy due to its radiative nature. It is the dominating cooling agent above 100 km that effectively regulates the thermospheric temperature. We studied the relationship between the response time of NO cooling emission and Joule heating rate during geomagnetic storm periods by using the NO emission observations by Sounding of Atmosphere using Broadband Emission Radiometry (SABER) onboard NASA’s Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics (TIMED) satellite and the Joule heating rate measured by the European incoherent scatter (EISCAT) radar over Tromsø (geographic:69.59°N, 19.22°E), Norway, and the thermosphere-ionosphere-electrodynamics general circulation model (TIEGCM) simulation. We selected seven geomagnetic storms during which there were continuous measurements of the Pederson conductivity and electric fields when the TIMED/SABER satellite was in the northview mode. The TIEGCM gives a fair description of NO cooling. However, it was found that the Joule heating rates obtained with the TIEGCM often do not show good agreement with those from EISCAT observations. The Joule heating rate peaks with 3–7 hr of storm' onset. The NO cooling flux takes longer time to respond to the energy input during storm period. The fastest response of the NO cooling emission to the Joule heating rate is found during the strongest storm. The weakest storm does not have the longest response time. No correlation between the response time of NO cooling flux with respect to the Joule heating rate and storm’s intensity was observed.

While the recent Van Allen Probes mission has provided a wealth of trapped particle measurements, questions still arise from analysis of their data. In particular, 10–100s of keV electrons exhibit dynamics not well understood with current data. Injections of 33–80 keV electrons can occur during both storm and quiet times as shown by Van Allen Probes data. However, due to the Probes' orbit, they can not distinguish precipitation during these events, necessary to quantify injection rates. Analysis of a Low Earth Orbit (LEO) mission measuring these electron populations found enhancements not explained with current injection sources. Future measurements including the quasi-trapped and precipitating populations are necessary to resolve these dynamics. We present the Medium Energy Electron Telescope to resolve these uncertainties surrounding 10–100s of keV electrons in the inner belt. This solid-state particle telescope is optimized for these measurements in the high-flux inner belt environment, with flight heritage from prior instruments. However, novel instrument design is required to measure these populations at fine energy resolution and fit the instrument into a 1U volume, including that of the detectors, electronics, instrument housing, and collimator components. The design is guided by Geant4 analysis and consideration for expected fluxes using the AE9 and AP9 models for a LEO mission and associated tradeoffs are discussed. We develop 59 energy channels with fine nominal resolution (<20%) for 30–800 keV electrons and show proton measurement capabilities for 1.1–>60 MeV populations. Instrument saturation and proton contamination are quantified by analysis of the instrument's response.

We investigate the generation of low latitude daytime E region field-aligned irregularities (FAIs) and focus on whether the FAIs are exclusively linked with sporadic E (Es) and whether the FAIs are driven by wind-driven gradient-drift instability (GDI). For this investigation, we use observations of FAIs made by the 30 MHz Gadanki Ionospheric Radar Interferometer (GIRI) and Es made by the DPS-4D digital ionosonde, both collocated at Gadanki. We have also examined simultaneous GIRI and Ionospheric Connection Explorer wind observations to substantiate our results. Results show that daytime E region FAIs are closely linked with Es activity and FAIs are not formed when Es is absent. FAIs are formed when blanketing frequency of Es (fbEs) exceeds frequency of the E layer. Both signal-to-noise ratio (SNR) and spectral width of the FAIs echoes display close relationship with top frequency of Es (ftEs). Further, SNR and spectral width of the FAIs echoes display a near-linear relationship with minimum deviation when Es activity is moderate and show dispersion from near-linear relationship when Es activity is strong. Zonal drift of the FAIs are predominantly eastward displaying both positive and negative vertical shear. Results also show that Es activity in terms of fbEs and ftEs is found to be closely related to the vertical shear in the zonal drift of the FAIs. Considering the height region of FAIs being highly collisional, we argue that zonal drifts are dominated by zonal wind and vertical shear in the zonal neutral wind is responsible for the formation/maintenance of Es layer and the eastward neutral wind is primarily responsible for the generation of the FAIs through the GDI.

The magnetopause deformation due to the upstream magnetosheath pressure perturbations is important to understand the solar wind-magnetosphere coupling process, but how to identify such events from in situ spacecraft observations is still challenging. In this study, we investigate magnetopause crossing events with fast-moving cold ions in the magnetosphere from Magnetospheric Multiscale observations, and find when fast-moving cold ions are present at the magnetopause, they are closely associated with the magnetopause deformation, which is featured by fast magnetopause motion and significant magnetopause normal deflection from model predictions. Therefore, fast-moving cold ions can be a useful indicator to search for magnetopause deformation events. By integrating the cold ion speed, the inferred magnetopause deformation amplitude varies from 0.1 to 2.0 RE. Further statistics indicate that such magnetopause deformation events prefer to occur under quasi-radial interplanetary magnetic field and fast solar wind conditions, suggesting high-speed magnetosheath jets could be one direct cause of magnetopause deformations.