JGR:Space physics

The proton radiation belt contains high fluxes of adiabatically trapped protons varying in energy from ∼one to hundreds of megaelectron volts (MeV). At large radial distances, magnetospheric field lines become stretched on the nightside of Earth and exhibit a small radius of curvature R C near the equator. This leads protons to undergo field line curvature (FLC) scattering, whereby changes to the first adiabatic invariant accumulate as field strength becomes nonuniform across a gyroorbit. The outer boundary of the proton belt at a given energy corresponds to the range of magnetic L shell over which this transition to nonadiabatic motion takes place, and is sensitive to the occurrence of geomagnetic storms. In this work, we first find expressions for nightside equatorial R C and field strength B e as functions of Dst and L* to fit the TS04 field model. We then apply the Tu et al. (2014, https://doi.org/10.1002/2014ja019864) condition for nonadiabatic onset to solve the outer boundary L*, and refine our expression for R C to achieve agreement with Van Allen Probes observations of 1–50 MeV proton flux over the 2014–2018 era. Finally, we implement this nonadiabatic onset condition into the British Antarctic Survey proton belt model (BAS-PRO) to solve the temporal evolution of proton fluxes at L ≤ 4. Compared with observations, BAS-PRO reproduces storm losses due to FLC scattering, but there is a discrepancy in mid-2017 that suggests a ∼5 MeV proton source not accounted for. Our work sheds light on outer zone proton belt variability at 1–10 MeV and demonstrates a useful tool for real-time forecasting.

Dust storms and the atmospheric waves can play a significant role in the dynamics in the upper neutral atmosphere of Mars. Recent observations found that a periodic variation of neutral H and O exists in the upper atmosphere, which is likely associated with atmospheric waves that occurred during the regional dust storm on 4 September 2016. However, such periodic variations accompanying the dust storm are far from understood in terms of the Martian ionized particles (i.e., ionosphere). Here we investigated the periodic variation of the ion density in the Martian ionosphere during the regional dust storm in September 2016 based on the MAVEN/STATIC observations. Assuming a simple Chapman layer model, we implemented numerical fitting for the ion density altitude profile to retrieve the peak ion densities in the Martian ionosphere. We then applied periodogram analysis to these peak ion densities in order to identify the peak periodicities together with their confidence levels. We identified several distinct peak periodicities with a scale from ∼1 day up to ∼20 days. The peak periodicities around 6–9 days are comparable to those seen in the Martian neutral atmosphere. In addition, the other peak periodicities likely correspond to the periodic crossings of the local crustal fields and periodic variation of the upstream solar wind, indicating a strong regional coupling between the lower atmosphere and the upstream solar wind.

The coupling between the ionosphere and plasmasphere is not known well during the quiet and storm times. The topside ionosphere ion composition reflects the exchange between the ionosphere and plasmasphere. We used the DMSP satellite observations to show the changes of the 840 km altitude H+ and O+ ions at ∼3–6 LT in the mid- and low-latitudes on 20–31 August 2018 where a storm occurred on 25 August. The H+ was the main ion at 3–4 LT, and its concentration presented two peaks around 40°N/S magnetic latitudes on 20–24 August. The H+ concentration strongly lessened for more than five days around 40°N/S magnetic latitudes in the recovery phase of the storm on 26–30 August. Further, the enhancement and reduction of the H+ concentration persisted mainly in the Pacific and Asian sectors. The O+ density increased in Northern Hemisphere at 6:36 LT, and the H+ concentration showed the maximum around the magnetic equator during the quiet times, while the H+ concentration still mainly reduced around 40°N/S and stronger at 40°S in the recovery phase. The plasmaspheric total electron content showed a peak around the magnetic equator on 20–24 August, and the peak was weakened or disappeared on 26–30 August. The plasmasphere might be the principal sources for the mid-latitude H+ enhancement along the geomagnetic field lines during the quite times and the storm-time suppressed plasmasphere plasma density would weaken the supply of the H+ at 840 km altitude.

Energetic ions outflowing from Jupiter's atmosphere was observed during Juno's 12th perijove crossing (PJ12) in the vicinity of Io's auroral footprint and reported by prior studies. It was hypothesized that Wave-Particle Interactions (WPI) with ion cyclotron waves observed coincident with the ion outflow may be responsible for the heating and subsequent outflow. This study uses numerical simulation and data model comparison to test whether ion cyclotron resonant heating is indeed a plausible mechanism to explain the intense ion outflow observed. Our simulations assume that the wave heating is of limited duration due to Io's footprint motion. The simulations are moreover compared to the previously published Jupiter Energetic Particle Detector Instruments (JEDI) observations at high energies, and the lower energy Jovian Auroral Distributions Experiment (JADE) observations that were not previously reported. We find that the ion cyclotron resonant heating mechanism can indeed lead to ion conic formation and strong vertical transport under certain assumptions about the distribution of wave power with altitude. We also find that the ion outflow is energized quickly with very rapid formation of the ion conic distribution. The implications of the intense ion outflow are also examined and it is found that such strong wave heating can lead to a depletion of the topside ionosphere.

Whistler mode waves are a common type of electromagnetic waves in the Martian induced magnetosphere. Using high-resolution magnetic field data from the Magnetometer (MAG) instrument onboard Mars Atmosphere and Volatile Evolution (MAVEN) from October 2014 to November 2022, we perform a detailed analysis of the statistical distribution of the occurrence rate, averaged amplitude, peak frequency, wave normal angle and ellipticity of left-hand and right-hand polarized whistler mode waves in the Martian induced magnetosphere. Our results show that whistler mode waves are mainly observed in the subsolar and magnetic pileup region, with the occurrence rate of right-hand mode waves higher than that of left-hand mode waves. The averaged wave amplitude ranges from 0.02 to 0.13 nT and peak wave frequency ranges from 2 to 9 Hz. We also find that the wave normal angles for both left-hand and right-hand polarized whistler waves are relatively larger in the subsolar region and magnetic pileup region where the corresponding wave ellipticity is closer to the linear polarization. Our results are valuable to in-depth understanding of the generation mechanism of whistler mode waves as well as their contributions to the electron dynamics in the Martian induced magnetosphere.

We use the three-dimensional (3-D) global hybrid code ANGIE3D to simulate the interaction of four solar wind tangential discontinuities (TDs) observed by ARTEMIS P1 from 0740 UT to 0800 UT on 28 December 2019 with the bow shock, magnetosheath, and magnetosphere. We demonstrate how the four discontinuities produce foreshock transients, a magnetosheath cavity-like structure, and a brief magnetopause crossing observed by THEMIS and MMS spacecraft from 0800 UT to 0830 UT. THEMIS D observed entries into foreshock transients exhibiting low density, low magnetic field strength, and high temperature cores bounded by compressional regions with high densities and high magnetic field strengths. The MMS spacecraft observed cavities with strongly depressed magnetic field strengths and highly deflected velocity in the magnetosheath downstream from the foreshock. Dawnside THEMIS A magnetosheath observations indicate a brief magnetosphere entry exhibiting enhanced magnetic field strength, low density, and decreased and deflected velocity (sunward flow). The solar wind inputs into the 3-D hybrid simulations resemble those seen by ARTEMIS. We simulate the interaction of four oblique TDs with properties similar to those in the observation. We place virtual spacecraft at the locations where observations were made. The hybrid simulations predict similar characteristics of the foreshock transients, a magnetosheath cavity, and a magnetopause crossing with characteristics similar to those observed by the multi-spacecraft observations. The detailed and successful comparison of the interaction involving multiple TDs will be presented.

During geomagnetic storms relativistic outer radiation belt electron flux exhibits large variations on rapid time scales of minutes to days. Many competing acceleration and loss processes contribute to the dynamic variability of the radiation belts; however, distinguishing the relative contribution of each mechanism remains a major challenge as they often occur simultaneously and over a wide range of spatiotemporal scales. In this study, we develop a new comprehensive model for storm-time radiation belt dynamics by incorporating electron wave-particle interactions with parallel propagating whistler mode waves into our global test-particle model of the outer belt. Electron trajectories are evolved through the electromagnetic fields generated from the Multiscale Atmosphere-Geospace Environment (MAGE) global geospace model. Pitch angle scattering and energization of the test particles are derived from analytical expressions for quasi-linear diffusion coefficients that depend directly on the magnetic field and density from the magnetosphere simulation. Using a study of the 17 March 2013 geomagnetic storm, we demonstrate that resonance with lower band chorus waves can produce rapid relativistic flux enhancements during the main phase of the storm. While electron loss from the outer radiation belt is dominated by loss through the magnetopause, wave-particle interactions drive significant atmospheric precipitation. We also show that the storm-time magnetic field and cold plasma density evolution produces strong, local variations of the magnitude and energy of the wave-particle interactions and is critical to fully capturing the dynamic variability of the radiation belts caused by wave-particle interactions.

The aurora often appears as an approximately oval shape surrounding the magnetic poles, and is a visible manifestation of the intricate coupling between the Earth's upper atmosphere and the near-Earth space environment. While the average size of the auroral oval increases with geomagnetic activity, the instantaneous shape and size of the aurora is highly dynamic. The identification of auroral boundaries holds significant value in space physics, as it serves to define and differentiate regions within the magnetosphere connected to the aurora by magnetic field lines. In this work, we demonstrate a new method to estimate the spatiotemporal variations of the poleward and equatorward boundaries in global UV images. We apply our method, which is robust against outliers and occasional bad data, to 2.5 years of UV imagery from the Imager for Magnetopause-to-Aurora Global Exploration satellite. The resulting data set is compared to recently published boundaries based on the same images (Chisham et al., 2022, https://doi.org/10.1029/2022JA030622), and shown to give consistent results on average. Our data set reveals a root mean square boundary normal velocity of 149 m/s for the poleward boundary and 96 m/s for the equatorward boundary and the velocities are shown to be stronger on the nightside than on the dayside. Interestingly, our findings demonstrate an absence of correlation between the amount of open magnetic flux and the amount of flux enclosed within the auroral oval.

Solar wind protons can charge exchange with the extensive hydrogen corona of Mars, resulting in a significant flux of energetic neutral atoms (ENAs). As these solar wind hydrogen ENAs precipitate into the upper atmosphere, they can experience electron attachment or detachment, resulting in populations of H− and H+, respectively, with upstream velocity. We seek to characterize the behavior of H− in the ionosphere of Mars through a combination of in situ data analysis and mathematical models. Observations indicate that measurable H− precipitation in the ionosphere of Mars is rare, occurring during only 1.8% of available observations. These events occur primarily during high energy solar wind conditions near perihelion. We also compare H− fluxes to those of H+ and find that H− fluxes are ∼4.5 times less than H+, indicating preferential conversion of hydrogen ENAs to H+. We develop a simple model describing the evolution of the charged and neutral fraction of ENAs and H− ions versus altitude. We find that 0.29%–0.78% of ENAs are converted to H− for solar wind energies 1–3 keV. We also predict that the effects of photodetachment on the H-H− system are non-negligible.

The sporadic E (Es) layer virtual height h′Es ${h}^{\prime }Es$ and the ordinary wave critical frequency foEs $foEs$ are routinely measured ionogram parameters, used to characterize the Es altitude occurrence and intensity. It has become common practice to take the real height hEs $hEs$ to be about equal to h′Es ${h}^{\prime }Es$ by assuming that signal propagation delays in the E region plasma below the layer are small and can be neglected. Although this applies for nighttime, during daytime it may overestimate hEs $hEs$ significantly. The present paper relies on true height analysis theory to devise a simplified method and propose an algorithm that can estimate Es real heights reasonably well. The method relies on h′Es ${h}^{\prime }Es$ and foEs $foEs$ ionosonde measurements and E region electron density profiles obtained from the International Reference Ionosphere model. The algorithm is applied to a typical set of Digisonde observations to compute hEs $hEs$ and examine real height variations and functional dependencies. Whereas hEs $hEs$ ≃ h′Es ${h}^{\prime }Es$ at nighttime, during daytime there are notable h′Es ${h}^{\prime }Es$ − hEs $hEs$ differences taking values less than 10 km for most of the observed layers. During the early morning and early afternoon hours, however, when weak layers appear at upper heights, the virtual to real height differences become larger reaching 20–25 km. The method proposed here for the estimation of hEs $hEs$ can be easily applied to improve the accuracy of the results of sporadic E layer studies.

The potential of deep learning for the investigation of medium scale traveling ionospheric disturbances (MSTIDs) has been exploited through the Sodankylä rapid-run ionosonde in this statistical study. The complementing observations of the Sodankylä ionosonde with those of the Sodankylä meteor radar reveals the diurnal and seasonal occurrence rate of high-latitude MSTIDs in the recent low solar activity period, 2018–2020. In our results, the daytime, nighttime and dusk MSTIDs are predominantly identified during winter, summer, and equinoctial months, respectively. The winter daytime higher (lower) occurrence rate is well correlated with the lower (higher) altitude of the height of the F2-layer peak (hmF2), and the low occurrence rate of the summer daytime is well correlated with the mesosphere-lower-thermosphere wind shear and higher gradient of temperature. Relatively high occurrence rate (>0.4) of summer nighttime MSTIDs has a general—but not one-to-one agreement—with post-noon to evening IU (eastward auroral current index) inferred ionospheric conductivity. Rather, we see a one-to-one relationship between the summer nighttime MSTIDs and zonal wind shear suggesting that the wind shear-induced electrodynamic processes could play significant roles for higher occurrence rate of MSTIDs. Furthermore, significant MSTIDs with ∼0.4 occurrence rate are so far revealed during spring and autumn transition periods. The enhanced nighttime MSTID amplitudes during the equinox are observed to be well correlated with IL index (westward auroral current indicator) suggesting that the particle precipitation during substorms could be the primary cause.

The relative motion of an artificial plasma stream in an ambient plasma background transverse to the geomagnetic field produces a polarization current and an E⇀×B⇀ $\overset{\rightharpoonup }{\mathrm{E}}\times \overset{\rightharpoonup }{\mathrm{B}}$ drift in the direction of the injected photoionized neutral cloud. The polarization current couples the ion cloud momentum to the ambient plasma causing the stream to brake or “skid” before coming to rest. A multi-fluid five-moment resistive magnetohydrodynamic (MHD) model has been previously developed and is used in this study to simulate an artificial barium cloud released at 8 km/s in the lower ionosphere transverse to the magnetic field. The MHD model’s governing equations are used to derive 2D analytical solutions to the equations of motion that capture oscillatory behavior relating the cloud ion cyclotron frequency, artificial and ambient plasma density, initial release velocity, and photoionization rate. The analytical solutions are compared to the MHD results for time t < 1 s and previous 1D momentum coupling models. It is shown that motion in the cloud release direction is consistent with previous models while motion in the polarization direction is consistent with numerical results.

The solar wind sputtering in the magnetospheric polar cusp is an important source of heavy atoms in Mercury’s exosphere and magnetosphere. However, the majority of ejected atoms are neutral, undergoing an extended period before photoionization occurs. In this study, we employ an ab initio simulation to investigate the behavior of sodium (Na) atoms prior to their photoionization. Our results reveal that overall only approximately 2.7% of the sputtered atoms contribute to magnetospheric ions, while the vast majority of these ions (∼82.9%) escape into interplanetary space. The remaining fraction (14.4%) eventually returns to the planetary surface. For Na atoms ionized inside the magnetosphere, a larger proportion of Na+ (53.5%) is supplied to the magnetotail compared to the polar cusp (39.4%), which is due to the tailward acceleration caused by solar radiation. Additionally, the remaining Na+ (7.1%) contributes to the dayside ring current region, as demonstrated by the observation of the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft. Our research introduces a perspective on Na+ transport in the magnetosphere that complements and coexists with traditional mechanisms.

The influence of atmospheric planetary waves on the occurrence of irregularities in the low latitude ionosphere is investigated using Whole Atmosphere Community Climate Model with thermosphere-ionosphere eXtension (WACCM-X) simulations and Global Observations of the Limb and Disk (GOLD) observations. GOLD observations of equatorial plasma bubbles (EPBs) exhibit a ∼6–8 day periodicity during January–February 2021. Analysis of WACCM-X simulations, which are constrained to reproduce realistic weather variability in the lower atmosphere, reveals that this coincides with an amplification of the westward propagating wavenumber-1 quasi-six day wave (Q6DW) in the mesosphere and lower thermosphere (MLT). The WACCM-X simulated Rayleigh-Taylor (R-T) instability growth rate, considered as a proxy of EPB occurrence, is found to exhibit a ∼6-day periodicity that is coincident with the enhanced Q6DW in the MLT. Additional WACCM-X simulations performed with fixed solar and geomagnetic activity demonstrate that the ∼6-day periodicity in the R-T instability growth rate is related to the forcing from the lower atmosphere. The simulations suggest that the Q6DW influences the day-to-day formation of EPBs through interaction with the migrating semidiurnal tide. This leads to periodic oscillations in the zonal winds, resulting in periodic variability in the strength of the prereversal enhancement, which influences the R-T instability growth rate and EPBs. The results demonstrate that atmospheric planetary waves, and their interaction with atmospheric tides, can have a significant impact on the day-to-day variability of EPBs.

Realistic modeling of the dynamics and variability in the upper mesosphere and lower thermosphere (UMLT) is critical to understand the coupling between different layers of the whole atmosphere system. Here we present simulations of the UMLT temperatures at ∼100 km altitude for one year during 2014 by the Whole Atmosphere Community Climate Model with thermosphere-ionosphere extension (WACCM-X) constrained below ∼90 km using meteorological analysis products of the high-altitude version of Navy Global Environmental Model (NAVGEM-HA). The model results are sampled at the same times and latitudes and longitudes as the satellite observations from Thermosphere Ionosphere and Mesosphere Electric Dynamics/Sounding of the Atmosphere using Broadband Emission Radiometry (TIMED/SABER). Comparisons show that the observed and modeled daily zonal mean temperatures are correlated (r ∼0.5–0.7) at most latitudes between ±50°. Both the observations and simulations show an annual variation at mid-latitudes in two hemispheres with the temperature maximum in summer and the minimum in winter, and at lower latitudes the semiannual variation becomes stronger having the temperature maximums at equinoxes and the minimums during solstices. However, the temperatures observed are on average ∼5–10 K (3%–5%) smaller than the model and the observations show a larger variability. Moreover, migrating tidal amplitudes are mostly overestimated by the model. Though differences are noticed, the WACCM-X simulations with NAVGEM-HA meteorological analyses are overall consistent with the SABER observations. These results support that whole atmosphere models informed by high altitude observations would help to simulate the UMLT variability and the atmosphere and ionosphere coupling.

Magnetosonic (MS) waves are electromagnetic emissions from a few to 100 Hz primarily confined near the magnetic equator both inside and outside the plasmasphere. Previous studies proved that MS waves can transport equatorially mirroring electrons from an equatorial pitch angle of 90° down to lower values by bounce resonance. However, the dependence of the bounce resonance effect on wave or background plasma parameters is still unclear. Here, we applied a test particle simulation to investigate electron transport coefficients, including diffusion and advection coefficients in energy and pitch angle, due to bounce resonance with MS waves. We investigate five wave field parameters, including wave frequency width, wave center frequency, latitudinal distribution width, wave normal angle root-mean-square of wave magnetic amplitude, and two background parameters, L-shell value and plasma density. We find different transport coefficient peaks resulting from different bounce resonance harmonics. Higher-order harmonic resonances exist, but the effect of fundamental resonance is much stronger. As the wave center frequency increases, higher-order harmonics start to dominate. With wave frequency width increasing, the energy range of effective bounce resonance broadens, but the effect itself weakens. The bounce resonance effect will increase when we decrease the wave normal angle, or increase the wave amplitude, latitudinal distribution width, L-shell value, and plasma density. The parametric study will advance our understanding of the favorable conditions of bounce resonance.

Measurements of CO2, Ar, N2, and O densities between 150 and 200 km from the Mars Atmosphere and Volatile Evolution Neutral Gas and Ion Mass Spectrometer during February 2015 to February 2022 are analyzed to provide a comprehensive analysis of their longitudinal wavenumber k = 1, 2, and 3 components. Variations in density amplitudes (A k ) with solar flux are marginally detectable during this period. The A k binned and averaged in latitude, local solar time and Ls are referenced to diurnal- and zonal-mean backgrounds in accord with how tides and stationary planetary waves (SPWs) are defined in theory and modeling. The resulting global A k distributions are the interference patterns formed by superposition of diurnal tides, SPWs and/or semidiurnal tides; consequently, a simple dependence on species mass consistent with thermal expansion (diffusive equilibrium) that might exist for some individual wave components is obscured. Additionally, vertical winds likely contribute to deviations from diffusive equilibrium. Complementary analyses of the Mars Climate Database indicate that the major contributors to the A k are DE2, SE1, DE1, and SPW1, 2, 3; support the absence of significant variability due to solar flux; and indicate a more well-defined sensitivity to species mass. The A k and their phases (longitudes of maxima) for the whole data set are available as part of Supporting Information S1.

Gravity waves (GWs) play an important role in the dynamics and energetics of the mesosphere. Geomagnetic activity is a known source of GWs in the upper atmosphere. However, how deep the effects of geomagnetic activity induced GWs penetrate into the mesosphere remains an open question. We use temperature measurements from the SABER/TIMED instrument between 2002 and 2018 to study the variations of mesospheric GW activity following intense geomagnetic disturbances identified by AE and Dst indices. By considering several case studies, we show for the first time that the GWs forced by geomagnetic activity can propagate down to about 80 km in the high latitude mesosphere. Only regions above 55° latitudes show a clear response. The fraction of cases in which there is an unambiguous enhancement in GW activity following the onset of geomagnetic disturbance is smaller during summer than other seasons. Only about half of the events show an unambiguous increase in GW activity during non-summer periods and about one quarter of the events in summer show an enhancement in GWs. In addition, we also find that the high latitude mesopause is often seen to descend in altitude following onset of geomagnetic activity in the non-summer high latitude region.

Most ionospheric models cannot sufficiently reproduce the observed electron density profiles in the E-region ionosphere, since they usually underestimate electron densities and do not match the profile shape. Mitigation of these issues is often addressed by increasing the solar soft X-ray flux which is ineffective for resolving data-model discrepancies. We show that low-resolution cross sections and solar spectral irradiances fail to preserve structure within the data, which considerably impacts radiative processes in the E-region, and are largely responsible for the discrepancies between observations and simulations. To resolve data-model inconsistencies, we utilize new high-resolution (0.001 nm) atomic oxygen (O) and molecular nitrogen (N2) cross sections and solar spectral irradiances, which contain autoionization and narrow rotational lines that allow solar photons to reach lower altitudes and increase the photoelectron flux. This work improves upon Meier et al. (2007, https://doi.org/10.1029/2006gl028484) by additionally incorporating high-resolution N2 photoionization and photoabsorption cross sections in model calculations. Model results with the new inputs show increased O+ production rates of over 500%, larger than those of Meier et al. (2007, https://doi.org/10.1029/2006gl028484) and total ion production rates of over 125%, while N2+ ${\mathrm{N}}_{2}^{+}$ production rates decrease by ∼15% in the E-region in comparison to the results obtained using the cross section compilation from Conway (1988, https://apps.dtic.mil/sti/pdfs/ADA193866.pdf). Low-resolution molecular oxygen (O2) cross sections from the Conway compilation are utilized for all input cases and indicate that O2+ ${\mathrm{O}}_{2}^{+}$ is a dominant contributor to the total ion production rate in the E-region. Specifically, the photoionization contributed from longer wavelengths is a main contributor at ∼120 km.

This Commentary article deals with the important role of large-amplitude, short-duration (<1 hr) and southwardly directed magnetic field incursions within the Sheath region of an interplanetary coronal mass ejection, which were recently shown to have led to extreme auroral activity during the early part of the Halloween storm (Ohtani, 2022, https://doi.org/10.1029/2022ja030596). For such largely geoeffective magnetic field structures some suggestions are given toward their possible interplanetary causes, which could also be associated with the origin of similar Sheath-structures observed during other events with very intense geomagnetic activity. A particular attention is given to a Sheath-incorporated and largely geoeffective flux rope-hypothesis. At the end we add some comments on further related magnetospheric and space weather issues.