Dusk Flank Research Articles

Lower hybrid and solitary waves in dusk flank region of the Earth’s magnetosphere

A variety of plasma waves have been detected in the vicinity of the magnetopause by various spacecraft missions. In this paper, utilizing high-resolution data from the Magnetospheric Multiscale (MMS) mission, we present new observations of simultaneous lower hybrid and ion solitary waves in the dusk flank region of Earth’s magnetosphere. All four MMS spacecraft consistently observed this wave activity during their traversal from Earth’s magnetosphere to the magnetosheath. Our analysis suggests that the lower hybrid drift waves, driven by lower-hybrid drift instability, were observed in correlation with density gradients. Furthermore, the analysis indicates that the entire ion bulk population drifts, which drives ion solitary waves. It is found that these waves play a crucial role in particle heating in the dusk flank region of Earth’s magnetosphere.

Plasmapause Surface Waves Triggered by Substorms

AbstractIn this study, we present in‐situ measurement on board Van Allen Probes showing that surface wave can cause oscillations of the plasmapause and modulate hiss and electron cyclotron harmonic (ECH) waves. Plasmapause surface wave (PSW) in Ps6 band was observed during an intense substorm on 16 July 2017, accompanied with quasi‐periodic emissions of hiss (positively correlated with plasma density) and ECH waves (anticorrelated with density). Phase relation between magnetic field and velocity perturbations indicate that the PSW was an eigenmode between southern and northern ionosphere. Measurements from ground‐based magnetometers suggest that such PSW propagates sunward at the dusk flank and were excited around the expansion phase of an intense substorm. Additional cases of PSWs are also presented to demonstrate that such waves are commonly observed near plasmapause and are likely to be related to the substorm activity.

Kinetic simulations of solar wind plasma irregularities crossing the Hermean magnetopause

Context. The physical mechanisms that favor the access of solar wind plasma into the magnetosphere have not been entirely elucidated to date. Studying the transport of finite-sized magnetosheath plasma irregularities across the magnetopause is fundamentally important for characterizing the Hermean environment (of Mercury) as well as for other planetary magnetic and plasma environments. Aims. We investigate the kinetic effects and their role on the penetration and transport of localized solar wind or magnetosheath plasma irregularities within the Hermean magnetosphere under the northward orientation of the interplanetary magnetic field. Methods. We used three-dimensional (3D) particle-in-cell (PIC) simulations adapted to the interaction between plasma elements (irregularities or jets) of a finite spatial extent and the typical magnetic field of Mercury’s magnetosphere. Results. Our simulations reveal the transport of solar wind plasma across the Hermean magnetopause and entry inside the magnetosphere. The 3D plasma elements are braked and deflected in the equatorial plane. The entry process is controlled by the magnetic field gradient at the magnetopause. For reduced jumps of the magnetic field (i.e., for larger values of the interplanetary magnetic field), the magnetospheric penetration is enhanced. The equatorial dynamics of the plasma element is characterized by a dawn-dusk asymmetry generated by first-order guiding center drift effects. More plasma penetrates into the dusk flank and advances deeper inside the magnetosphere than in the dawn flank. Conclusions. The simulated solar wind or magnetosheath plasma jets can cross the Hermean magnetopause and enter into the magnetosphere, as described by the impulsive penetration mechanism.

Simulated Trapping of Solar Energetic Protons for the 8–10 March 2012 Geomagnetic Storm: Impact on Inner Zone Protons as Measured by Van Allen Probes

AbstractSolar energetic protons (SEPs) have been shown to contribute significantly to the inner zone trapped proton population for energies <100 MeV and L > 1.3 (Selesnick et al., 2007, https://doi.org/10.1029/2006sw000275). The Relativistic Electron Proton Telescope (REPT) on the Van Allen Probes launched 30 August 2012 observed a double‐peaked (in L) inner zone population throughout the 7‐year lifetime of the mission. It has been proposed that a strong SEP event accompanied by a CME‐shock in early March 2012 provided the SEP source for the higher L trapped proton population, which then diffused radially inward to be observed by REPT at L ∼ 2. Here, we follow trajectories of SEP protons launched isotropically from a sphere at 7 Re in 15 s cadence fields from an Lyon‐Fedder‐Mobarry coupled to Rice Convection Model global magnetohydrodynamic (MHD) simulation driven by measured upstream solar wind parameters. The timescale of the interplanetary shock arrival is captured, launching a magnetosonic impulse propagating azimuthally along the dawn and dusk flanks inside the magnetosphere, shown previously to produce SEP trapping. The MHD‐test particle simulation uses Geostationary Operational Environmental Satellite (GOES) proton energy spectra to weight the initial radial profile required for the radial diffusion calculation over the following 2 years. GOES proton measurements also provide a dynamic outer boundary condition for radial diffusion. A direct comparison with REPT measurements 20 months following the trapping event in March 2012 supports this novel combination of short‐term and long‐term evolution of the newly trapped protons.

Kelvin-Helmholtz Vortices as an Interplay of Magnetosphere-Ionosphere Coupling

The solar wind-magnetosphere interaction drives diverse physical processes on the flanks of Earth’s magnetopause, and in turn these processes couple to the ionosphere. We investigate simultaneous multipoint in-situ spacecraft and ground-based measurements to determine the role of Kelvin-Helmholtz waves at the Earth’s magnetopause and the low-latitude boundary layer in the magnetosphere-ionosphere coupling process. Nonlinear Kelvin-Helmholtz waves develop into flow vortices that twist and/or shear flux tube magnetic fields, thereby generating localized field-aligned currents. Kelvin-Helmholtz vortices on the dusk (dawn) flanks of the magnetosphere generate clockwise (counter-clockwise) rotations and upward (downward) field-aligned currents inside the flux tubes, consistent with the region-1 field-aligned current. We present in-situ MMS and Cluster spacecraft observations of Kelvin-Helmholtz vortices at the magnetopause that map to the poleward edge of the auroral regions. The FAST spacecraft and the ground-based magnetometers from which spherical elementary currents (acting as a proxy for vertical currents) can be calculated observe corresponding field-aligned current signatures. This study demonstrates the role played by the Kelvin-Helmholtz waves in linking magnetopause boundary fluctuations to ionospheric phenomena.

Multi-Spacecraft Observations of Fluctuations Occurring Along the Dusk Flank Magnetopause, and Testing the Connection to an Observed Ionospheric Bead

During 2018 November 06, 11:30—18:00 UT, the MMS constellation, the Cluster set of spacecraft, and the Geotail spacecraft were all situated near the dusk flank magnetopause. Large scale fluctuations were observed by the available and operating science instruments at these various spacecraft (i.e., magnetic field, plasma moment, and energy flux measurements). Similar fluctuations were not observed by upstream solar wind monitors, suggesting that the waves were initiated at the magnetopause. A localized emission ‘bead’ from the post-noon ionosphere was also observed from low Earth orbit. The nature and relation of the fluctuations observed at all of these spacecraft at the magnetosphere boundary and the connection to the post-noon high-latitude ionosphere are investigated in this study.

Polarization Properties of 3‐D Field Line Resonances

AbstractRecent work has shown that field line resonances (FLRs) can form with intermediate magnetic field polarizations, that is, between toroidal (azimuthal) and poloidal (radial). In this case, they are referred to as “3‐D”, given that the equilibrium and therefore associated wave solutions vary in all three dimensions. Such 3‐D variations in the medium can be expected in the magnetosphere in particular during the presence of a plasmaspheric drainage plume in the dusk flank. In this case, strong field‐aligned currents driven by FLRs can form along contours of Alfvén frequency which are very asymmetric azimuthally. In this study, we use magnetohydrodynamic (MHD) simulations performed in a background dipole magnetic field, to consider what the satellite signatures would be for a spacecraft observing such non‐toroidal FLRs under different magnetopause driving conditions. We consider how the 3‐D nature impacts the amplitude and phase relations traditionally expected for FLRs. We find that in the 3‐D FLR regions of intermediate polarization, the FLR has strong contributions from both the radial and azimuthal field components (velocity, magnetic or electric) and these are in or out of phase. We show how hodograms can be used to infer the FLR polarization in these cases as well as provide predictions for the variation of the FLR polarization with MLT across a plasmaspheric plume. We further show how the presence of a significant global fast mode can skew the resulting polarization estimates.

The Kelvin-Helmholtz Instability From the Perspective of Hybrid Simulations

The flow shear-driven Kelvin-Helmholtz (KH) instability is ubiquitous in planetary magnetospheres. At Earth these surface waves are important along the dawn and dusk flanks of the magnetopause boundary while at Jupiter and Saturn the entire dayside magnetopause boundary can exhibit KH activity due to corotational flows in the magnetosphere. Kelvin-Helmholtz waves can be a major ingredient in the so-called viscous-like interaction with the solar wind. In this paper, we review the KH instability from the perspective of hybrid (kinetic ions, fluid electrons) simulations. Many of the simulations are based on parameters typically found at Saturn’s magnetopause boundary, but the results can be generally applied to any KH-unstable situation. The focus of the discussion is on the ion kinetic scale and implications for mass, momentum, and energy transport at the magnetopause boundary.

Solar Wind—Magnetosphere Coupling During Radial Interplanetary Magnetic Field Conditions: Simultaneous Multi‐Point Observations

AbstractIn‐situ spacecraft missions are powerful assets to study processes that occur in space plasmas. One of their main limitations, however, is extrapolating such local measurements to the global scales of the system. To overcome this problem at least partially, multi‐point measurements can be used. There are several multi‐spacecraft missions currently operating in the Earth's magnetosphere, and the simultaneous use of the data collected by them provides new insights into the large‐scale properties and evolution of magnetospheric plasma processes. In this work, we focus on studying the Earth's magnetopause (MP) using a conjunction between the Magnetospheric Multiscale and Cluster fleets, when both missions skimmed the MP for several hours at distant locations during radial interplanetary magnetic field (IMF) conditions. The observed MP positions as a function of the evolving solar wind conditions are compared to model predictions of the MP. We observe an inflation of the magnetosphere (∼0.7 RE), consistent with magnetosheath pressure decrease during radial IMF conditions, which is less pronounced on the flank (0.2 RE). There is observational evidence of magnetic reconnection in the subsolar region for the whole encounter, and in the dusk flank for the last portion of the encounter, suggesting that reconnection was extending more than 15 RE. However, reconnection jets were not always observed, suggesting that reconnection was patchy, intermittent or both. Shear flows reduce the reconnection rate up to ∼30% in the dusk flank according to predictions, and the plasma β enhancement in the magnetosheath during radial IMF favors reconnection suppression by the diamagnetic drift.

Characteristics of Jupiter's X‐Ray Auroral Hot Spot Emissions Using Chandra

AbstractTo help understand and determine the driver of jovian auroral X‐rays, we present the first statistical study to focus on the morphology and dynamics of the jovian northern hot spot (NHS) using Chandra data. The catalog we explore dates from December 18, 2000 up to and including September 8, 2019. Using a numerical criterion, we characterize the typical and extreme behavior of the concentrated NHS emissions across the catalog. The mean power of the NHS is found to be 1.91 GW with a maximum brightness of 2.02 Rayleighs (R), representing by far the brightest parts of the jovian X‐ray spectrum. We report a statistically significant region of emissions at the NHS center which is always present, the averaged hot spot nucleus (AHSNuc), with mean power of 0.57 GW and inferred average brightness of 1.2 R. We use a flux equivalence mapping model to link this distinct region of X‐ray output to a likely source location and find that the majority of mappable NHS photons emanate from the pre‐dusk to pre‐midnight sector, coincident with the dusk flank boundary. A smaller cluster maps to the noon magnetopause boundary, dominated by the AHSNuc, suggesting that there may be multiple drivers of X‐ray emissions. On application of timing analysis techniques (Rayleigh, Monte Carlo, Jackknife), we identify several instances of statistically significant quasi‐periodic oscillations (QPOs) in the NHS photons ranging from 2.3 to 36.4 min, suggesting possible links with ultra‐low frequency activity on the magnetopause boundary (e.g., dayside reconnection, Kelvin‐Helmholtz instabilities).

Effects of the IMF Direction on Atmospheric Escape From a Mars‐like Planet Under Weak Intrinsic Magnetic Field Conditions

AbstractDirection of the upstream interplanetary magnetic field (IMF) significantly changes the magnetospheric configuration, influencing the atmospheric escape mechanism. This paper investigates effects of IMF on the ion escape mechanism from a Mars‐like planet that has a weak dipole magnetic field directing northward on the equatorial surface. The northward (parallel to the dipole at subsolar), southward (antiparallel), and Parker‐spiral IMFs under present solar wind conditions are compared based on multispecies magnetohydrodynamics simulations. In the northward IMF case, molecular ions escape from the high‐latitude lobe reconnection region, where ionospheric ions are transported upward along open field lines. Atomic oxygen ions originating either in the ionosphere or oxygen corona escape through a broader ring‐shaped region. In the southward IMF case, the escape flux of heavy ions increases significantly and has peaks around the equatorial dawn and dusk flanks. The draped IMF can penetrate into the subsolar ionosphere by erosion, and the IMF becomes mass‐loaded as it is transported through the dayside ionosphere. The mass‐loaded draped IMF is carried to the tail, contributing to ion escape. The escape channels in the northward and southward IMF cases are different from those in the Parker‐spiral IMF case. The escape rate is the lowest in the northward IMF case and comparable in the Parker‐spiral and southward IMF cases. In the northward IMF case, a weak intrinsic dipole forms a magnetosphere configuration similar to that of Earth, quenching the escape rate, while the Parker‐spiral and southward IMFs cause reconnection and erosion, promoting ion escape from the upper atmosphere.

Field‐Aligned Currents in Auroral Vortices

AbstractAuroral bright spots have been observed at Earth, Jupiter, and Saturn in regions that map to the magnetopause boundary layer. It has been suggested that the bright spots are associated with the Kelvin‐Helmholtz (KH) instability. We utilize a quasistatic magnetosphere‐ionosphere coupling model driven by a vortex in the boundary layer to determine how the field‐aligned current structure depends on ionospheric and boundary layer parameters. We compare vortex induced currents with shear‐flow induced currents. We find that the strength of the maximum currents are comparable, but the structure is significantly different. For a vortex, the current and electron precipitation maximize when the vortex size mapped to the ionosphere is approximately 1.5L, where is the auroral scale length, ΣP is the Pedersen conductivity, and κ is the Knight parameter. For a vortex, the current width provides a direct measure of the size, Δ, of the boundary layer structure, while the current width of shear‐flow aurora is generally determined by the larger of Δ or L. For comparison with observations, an event is considered where auroral bright spots in the ionosphere are detected by DMSP SSUSI FUV when KH structures are observed on the dusk flank by THEMIS.

Formation and Release of the Harang Reversal Relating With the Substorm Onset Process

AbstractWith the interplanetary magnetic field (IMF) that turns from northward to southward, the global simulation successively reproduces the growth phase, the onset, and the expansion phase of the substorm. The calculated ionospheric convection for the growth phase reproduces the development of the Harang reversal (HR), the upward field‐aligned current (FAC), and the upcoming onset point as observed. Magnetic field lines traced from the center of the nightside upward FAC are open, while the magnetic field line traced from the onset point is closed. These open magnetic field lines map to flow shear just outside the O/C boundary. Seen from the magnetic configuration, the growth phase proceeds as the replacement process of nulls. Two new nulls appear on the dayside under the southward IMF, while two old nulls under the northward IMF retreat tailward forming two bifurcation regions on the dawn and dusk flanks. Flow shear around the O/C boundary forms by magnetospheric convection that returns to the dayside via bifurcation regions. The expansion phase proceeds through a topological change by the near‐earth neutral line (NENL). The NENL occurs inside the thinned structure of the northward IMF remnant, on the low‐latitude side of the flow shear, and projects down to the low‐latitude edge of the upward FAC. Associated with the NENL, the convection return path changes to the center of the plasma sheet and reveals in the ionosphere as the release of the HR. By the shrinkage of projecting magnetic field line, the O/C boundary migrates poleward in the ionosphere.

Prediction and Understanding of Soft-proton Contamination in XMM-Newton: A Machine Learning Approach

Abstract One of the major and unfortunately unforeseen sources of background for the current generation of X-ray telescopes are few tens to hundreds of keV (soft) protons concentrated by the mirrors. One such telescope is the European Space Agency’s (ESA) X-ray Multi-Mirror Mission (XMM-Newton). Its observing time lost due to background contamination is about 40%. This loss of observing time affects all the major broad science goals of this observatory, ranging from cosmology to astrophysics of neutron stars and black holes. The soft-proton background could dramatically impact future large X-ray missions such as the ESA planned Athena mission (http://www.the-athena-x-ray-observatory.eu/). Physical processes that trigger this background are still poorly understood. We use a machine learning (ML) approach to delineate related important parameters and to develop a model to predict the background contamination using 12 yr of XMM-Newton observations. As predictors we use the location of the satellite and solar and geomagnetic activity parameters. We revealed that the contamination is most strongly related to the distance in the southern direction, Z (XMM-Newton observations were in the southern hemisphere), the solar wind radial velocity, and the location on the magnetospheric magnetic field lines. We derived simple empirical models for the first two individual predictors and an ML model that utilizes an ensemble of the predictors (Extra-Trees Regressor) and gives better performance. Based on our analysis, future missions should minimize observations during times associated with high solar wind speed and avoid closed magnetic field lines, especially at the dusk flank region in the southern hemisphere.

Latitudinal Dependence of the Kelvin‐Helmholtz Instability and Beta Dependence of Vortex‐Induced High‐Guide Field Magnetic Reconnection

AbstractWe investigate both large‐ and small‐scale properties of a Kelvin‐Helmholtz (KH) event at the dusk flank magnetopause using Magnetospheric Multiscale observations on 8 September 2015. We first use two types of 3‐D simulations (global and local) to demonstrate that Magnetospheric Multiscale is close to the most KH unstable region, and so the occurrence of vortex‐induced reconnection may be expected. Because they produce low‐shear current sheets, KH vortices constitute a perfect laboratory to investigate magnetic reconnection with large guide field and low asymmetry. Recent works suggest that magnetic reconnection may be suppressed when a current sheet combines large guide field and pressure gradient (which induces a diamagnetic drift). We thus perform a statistical analysis of high‐resolution data for the 69 KH‐induced low‐shear magnetic reconnection events observed on that day. We find that the suppression mechanism is not at work for most of the observed reconnecting current sheets, as predicted, but we also find that almost all nonreconnecting current sheets should be reconnecting according to this model. This confirms the fact that the model provides a necessary but not sufficient condition for reconnection to occur. Finally, based on the same data set, we study the latitudinal distribution of these magnetic reconnection events combined with global magnetospheric modeling. We find that reconnection associated with KH vortices occurs over a significant range of latitudes at the flank magnetopause. It is not confined to the plane where the growth rate is maximum, in agreement with recent 3‐D simulations.

Characteristics of the Flank Magnetopause: MMS Results

AbstractWe have used a large number of magnetopause crossings by the Magnetospheric Multiscale (MMS) mission to investigate macroscopic properties of this current sheet, with emphasis on the flanks of the magnetopause. Macroscopic features such as thickness, location, and motion of the magnetopause were calculated as a function of local time sector. The results show that the flanks of the magnetopause are significantly thicker than the dayside magnetopause. Thicknesses vary from about 650 km near noon to over 1,000 km near the terminator. Current densities vary in a similar manner, with average current densities around noon almost twice as high as near the terminator. We also find a dawn‐dusk asymmetry in many of the macroscopic parameters; the dawn magnetopause is thicker than at dusk, while the dusk flank is more dynamic, with a higher average normal velocity.

Generation of Turbulence in Kelvin‐Helmholtz Vortices at the Earth's Magnetopause: Magnetospheric Multiscale Observations

AbstractThe Kelvin‐Helmholtz instability (KHI) at Earth's magnetopause and associated turbulence are suggested to play a role in the transport of mass and momentum from the solar wind into Earth's magnetosphere. We investigate electromagnetic turbulence observed in Kelvin‐Helmholtz vortices encountered at the dusk flank magnetopause by the Magnetospheric Multiscale (MMS) spacecraft under northward interplanetary magnetic field (IMF) conditions in order to reveal its generation process, mode properties, and role. A comparison with another MMS event at the dayside magnetopause with reconnection but no KHI signatures under a similar IMF condition indicates that while high‐latitude magnetopause reconnection excites a modest level of turbulence in the dayside low‐latitude boundary layer, the KHI further amplifies the turbulence, leading to magnetic energy spectra with a power law index −5/3 at magnetohydrodynamic scales even in its early nonlinear phase. The mode of the electromagnetic turbulence is analyzed with a single‐spacecraft method based on Ampère's law, developed by Bellan (2016, https://doi.org/10.1002/2016JA022827), for estimating wave vectors as a function of spacecraft frame frequency. The results suggest that the turbulence does not consist of propagating normal‐mode waves but is due to interlaced magnetic flux tubes advected by plasma flows in the vortices. The turbulence at sub‐ion scales in the early nonlinear phase of the KHI may not be the cause of the plasma transport across the magnetopause but rather a consequence of three‐dimensional vortex‐induced reconnection, the process that can cause an efficient transport by producing tangled reconnected field lines.

Simulations of MHD Wave Propagation and Coupling in a 3‐D Magnetosphere

AbstractA novel simulation grid is devised that is optimized for studying magnetohydrodynamic (MHD) wave coupling and phase mixing in a dipole‐like magnetic field. The model also includes flaring on the dawn and dusk flanks. The location of the magnetopause is quite general. In particular, it does not have to coincide with a coordinate surface. Simulations indicate the central role of global fast waveguide modes. These switch from being azimuthally standing in nature at noon, to propagating antisunward on the flanks. The field line resonances (FLRs) seen in the simulation results are three dimensional and not strictly azimuthally polarized. When a plume is present, the FLRs cross a range of 2 in shell, and have a polarization that is midway between toroidal and poloidal.

Solar Wind Ion Entry Into the Magnetosphere During Northward IMF.

Extended periods of northward interplanetary magnetic field (IMF) lead to the formation of a cold, dense plasma sheet due to the entry of solar wind plasma into the magnetosphere. Identifying the paths that the solar wind takes to enter the magnetosphere, and their relative importance has remained elusive. Any theoretical model of entry must satisfy observational constraints, such as the overall entry rate and the dawn‐dusk asymmetry observed in the cold, dense plasma sheet. We model, using a combination of global magnetohydrodynamic and test particle simulations, solar wind ion entry into the magnetosphere during northward IMF and compare entry facilitated by the Kelvin‐Helmholtz instability to cusp reconnection. For Kelvin‐Helmholtz entry we reproduce transport rates inferred from observation and kinetic modeling and find that intravortex reconnection creates buoyant flux tubes, which provides, through interchange instability, a mechanism of filling the central plasma sheet with cold magnetosheath plasma. For cusp entry we show that an intrinsic dawn‐dusk asymmetry is created during entry that is the result of alignment of the westward ion drift with the dawnward electric field typically observed during northward IMF. We show that both entry mechanisms provide comparable mass but affect entering plasma differently. The flank‐entering plasma is cold and dawn‐dusk symmetric, whereas the cusp‐entering plasma is accelerated and preferentially deflected toward dawn. The combined effect of these entry mechanisms results in a plasma sheet population that exhibits dawn‐dusk asymmetry in the manner that is seen in nature: a two‐component (hot and cold) dusk flank and hotter, broadly peaked dawn population.

Prolonged Kelvin–Helmholtz Waves at Dawn and Dusk Flank Magnetopause: Simultaneous Observations by MMS and THEMIS

Abstract The Kelvin–Helmholtz (K-H) waves predominantly excited at the Earth’s low-latitude magnetopause were suggested to be dawn–dusk asymmetric. We report a prolonged simultaneous observations of the K-H waves on the dawn and dusk magnetopause by Magnetospheric Multiscale (MMS) and THEMIS-A (THA) spacecraft, respectively. The quasi-periodic K-H waves on both flanks have unambiguous low-density and high-speed patterns. The wave periods vary gradually on both flanks, with similar average periods (303 ± 107 s for MMS and 266 ± 102 s for THA). The lag time between the variations of the wave periods is close to the wave propagation time from THA to MMS, which suggests that the K-H waves generate and propagate quasi-symmetrically on both flanks. Larger local magnetic shear angles are observed on the trailing edges by MMS than by THA, which is probably due to the strong magnetic field distortion during the tailward propagation. The increased magnetic shear may excite magnetic reconnection, thus contributing to the formation of the low-latitude boundary layer.