Magnetopause Location Research Articles

Magnetopause location modeling using machine learning: inaccuracy due to solar wind parameter propagation

An intrinsic limitation of empirical models of the magnetopause location is a predefined magnetopause shape and assumed functional dependences on relevant parameters. We overcome this limitation using a machine learning approach (artificial neural networks), allowing us to incorporate general, purely data-driven dependences. For the training and testing of the developed neural network model, a data set of about 15,000 magnetopause crossings identified in the THEMIS A-E, Magion 4, Geotail, and Interball-1 satellite data in the subsolar region is used. A cylindrical symmetry around the direction of the impinging solar wind is assumed, and solar wind dynamic pressure, interplanetary magnetic field magnitude, cone angle, clock angle, tilt angle, and corrected Dst index are considered as parameters. The effect of these parameters on the magnetopause location is revealed. The performance of the developed model is compared with other empirical magnetopause models. Finally, we demonstrate and discuss the inaccuracy of magnetopause models due to the inaccurate information about the impinging solar wind parameters based on measurements near the L1 point. This inaccuracy imposes a theoretical limit on the precision of magnetopause predictions, a limit that our model closely approaches.

Finding the magnetopause location using soft X-ray observations and a statistical inverse method

Finding the magnetopause location using soft X-ray observations and a statistical inverse method

Extremely Distant Magnetopause Locations Caused by Magnetosheath Jets

AbstractMagnetopause position is controlled mainly by the solar wind dynamic pressure and north‐south interplanetary magnetic field component and these quantities are included in different empirical magnetopause models. We have collected about 50,000 of dayside magnetopause crossings observed by THEMIS in course of 2007–2019 and compared the observed magnetopause position with model prediction. The difference between observed and predicted magnetopause radial distance, Robs − Rmod is used for quantifying the model‐observation agreement. Its median values are well predicted for cases up to Robs ≈ 12 RE for all models but higher positive deviations are found for larger magnetopause distances, mainly under a nearly radial field and low dynamic pressure. The analysis reveals their connection with transient magnetopause displacements caused by strong sunward flows in the magnetosheath. We discuss the possible origin of the observed magnetosheath flow switching in terms of the interaction of magnetosheath jets with the magnetopause.

Properties of Magnetohydrodynamic Normal Modes in the Earth's Magnetosphere

AbstractThe Earth's magnetosphere supports a variety of Magnetohydrodynamic (MHD) normal modes with Ultra Low Frequencies (ULF) including standing Alfvén waves and cavity/waveguide modes. Their amplitudes and frequencies depend in part on the properties of the magnetosphere (size of cavity, wave speed distribution). In this work, we use ∼13 years of Time History of Events and Macroscale Interactions during Substorms satellite magnetic field observations, combined with linearized MHD numerical simulations, to examine the properties of MHD normal modes in the region L > 5 and for frequencies <80 mHz. We identify persistent normal mode structure in observed dawn sector power spectra with frequency‐dependent wave power peaks like those obtained from simulation ensemble averages, where the simulations assume different radial Alfvén speed profiles and magnetopause locations. We further show with both observations and simulations how frequency‐dependent wave power peaks at L > 5 depend on both the magnetopause location and the location of peaks in the radial Alfvén speed profile. Finally, we discuss how these results might be used to better model radiation belt electron dynamics related to ULF waves.

Tracking Magnetopause Motion Using Cold Plasmaspheric Ions

AbstractWe demonstrate that any plasmaspheric/cold ions accelerated in the vicinity of the magnetopause boundary, can proxy the local magnetopause motion over many minutes. The timeseries of this motion capture local structures such as waves on the boundary. We determine cold ion velocities normal to full magnetopause boundary crossings for three events with varying distances to the predicted reconnection X‐line, thus, providing a proof‐of‐concept study demonstrating the potential for using cold ion velocities to track magnetopause motion over a long period of time. Obtaining the time history of the (local) motion of the magnetopause relative to the spacecraft is determined by integrating the bulk (<100 eV for H+) ion velocities normal to the boundary. Timeseries of these tracked cold ion accelerations may be used to investigate boundary layer thicknesses, potential wave structures on the magnetopause, and their evolution beyond the boundary crossing. This method generally tracks magnetopause motion out to distances of ∼1–2RE away from the spacecraft during quasi‐steady space weather conditions.

Comparing magnetopause predictions from two MHD models during a geomagnetic storm and a quiet period

Magnetopause location is an important prediction of numerical simulations of the magnetosphere, yet the models can err, either under-predicting or over-predicting the motion of the boundary. This study compares results from two of the most widely used magnetohydrodynamic (MHD) models, the Lyon–Fedder–Mobarry (LFM) model and the Space Weather Modeling Framework (SWMF), to data from the GOES 13 and 15 satellites during the geomagnetic storm on 22 June 2015, and to THEMIS A, D, and E during a quiet period on 31 January 2013. The models not only reproduce the magnetopause crossings of the spacecraft during the storm, but they also predict spurious magnetopause motion after the crossings seen in the GOES data. We investigate the possible causes of the over-predictions during the storm and find the following. First, using different ionospheric conductance models does not significantly alter predictions of the magnetopause location. Second, coupling the Rice Convection Model (RCM) to the MHD codes improves the SWMF magnetopause predictions more than it does for the LFM predictions. Third, the SWMF produces a stronger ring current than LFM, both with and without the RCM and regardless of the LFM spatial resolution. During the non-storm event, LFM predicts the THEMIS magnetopause crossings due to the southward interplanetary magnetic field better than the SWMF. Additionally, increasing the LFM spatial grid resolution improves the THEMIS predictions, while increasing the SWMF grid resolutions does not.

Study of Extreme Magnetopause Distortions Under Varying Solar Wind Conditions

AbstractTo first order, the magnetopause (MP) is defined by a pressure balance between the solar wind and the magnetosphere. The boundary moves under the influence of varying solar wind conditions and transient foreshock phenomena, reaching unusually large and small distances from the Earth. We investigate under which solar wind conditions such extreme MP distortions occur. Therefore, we construct a database of magnetopause crossings (MPCs) observed by the THEMIS spacecraft in the years 2007 to mid‐2022 using a simple Random Forest Classifier. Roughly 7% of the found crossing events deviate beyond reported errors in the stand‐off distance from the Shue et al. (1998, https://doi.org/10.1029/98JA01103) MP model and thus are termed extreme distortions. We find the occurrence of these extreme events in terms of expansion or compression of the MP to be linked to different solar wind parameters, most notably to the IMF magnitude, cone angle, velocity, Alfvén Mach number and temperature. Foreshock transients like hot‐flow anomalies and foreshock bubbles could be responsible for extreme magnetospheric expansions. The results should be incorporated into future magnetopause models and may be helpful for the reconstruction of the MP locations out of soft x‐ray images, relevant for the upcoming SMILE mission.

Variation of the Jovian Magnetopause Under Constant Solar Wind Conditions: Significance of Magnetodisc Dynamics

AbstractIt is generally believed that variations in the upstream solar wind (SW) and interplanetary magnetic field (IMF) conditions are the main cause of changes of Jupiter's magnetopause (JM) location. However, most previous pressure balance models for the JM are axisymmetric and do not consider internal drivers, for example, the dynamics of the magnetodisc. We use three‐dimensional global magnetosphere simulations to investigate the variation of the JM under constant SW/IMF conditions. These simulations show that even without variations in the upstream driving conditions, the JM can exhibit dynamic variations, suggesting a range as large as 50 Jupiter radii in the subsolar location. Our study shows that the interchange structures in the Jovian magnetodisc will introduce significant radial dynamic pressure, which can drive significant variation in the JM location. The results provide important new context for interpreting the JM location and dynamics, with key implications for other internally mass‐loaded and/or rapidly rotating systems.

Investigating Potential Causes for the Prediction of Spurious Magnetopause Crossings at Geosynchronous Orbit in MHD Simulations

AbstractDuring intense geomagnetic storms, the magnetopause can move in as far as geosynchronous orbit, leaving the satellites in that orbit out in the magnetosheath. Spacecraft operators turn to numerical models to predict the response of the magnetopause to solar wind conditions, but the predictions of the models are not always accurate. This study investigates four storms with a magnetopause crossing by at least one GOES satellite, using four magnetohydrodynamic models at NASA's Community Coordinated Modeling Center to simulate the events, and analyzes the results to investigate the reasons for errors in the predictions. Two main reasons can explain most of the erroneous predictions. First, the solar wind input to the simulations often contains features measured near the L1 point that did not eventually arrive at Earth; incorrect predictions during such periods are due to the solar wind input rather than to the models themselves. Second, while the models do well when the primary driver of magnetopause motion is a variation in the solar wind density, they tend to overpredict or underpredict the integrated Birkeland currents and their effects during times of strong negative interplanetary magnetic field (IMF) Bz, leading to poorer prediction capability. Coupling the MHD codes to a ring current model, when such a coupling is available, generally will improve the predictions but will not always entirely correct them. More work is needed to fully characterize the response of each code under strong southward IMF conditions as it relates to prediction of magnetopause location.

Magnetopause Standoff Position Changes and Geosynchronous Orbit Crossings: Models and Observations

AbstractThis research examines the ability of current physics‐based models to predict the magnetopause location. We use the Run‐On‐Request capabilities at the Community Coordinated Modeling Center at NASA GSFC with 4 magnetohydrodynamic (MHD) models. The magnetopause position prediction and response time to the solar wind changes is then compared to extreme solar wind conditions where magnetopause crossing of geosynchronous orbit have been observed by the Geostationary Operational Environmental Satellite (GOES) satellites. Rigorous analysis/comparison of observations and models is critical in determining magnetosphere dynamics for model validation. This paper is a preliminary effort defining the metrics necessary to understand the current magnetosphere model capabilities and challenges. Results show that there are discrepancies between the MHD models' standoff positions of the dayside magnetopause for the same solar wind conditions on events that included (a) an increase in solar wind dynamic pressure and a step function in the Interplanetary Magnetic Field Bz component; (b) nominal solar wind conditions (values of approximated 400 km/s for solar wind speed and 5 nT for the magnetic field magnitude) with a northward IMF; and (c) compression caused by several coronal mass ejections impacting the near Earth environment. Overall, the models predicted different magnetopause subsolar locations sometimes in the order of 3 RE. Contingency tables were calculated to show model performance in comparison with the data observed with the GOES 13/15 geosynchronous orbit for extreme events and skill scores were calculated.

Interplanetary Magnetic Field By Controls the Magnetopause Location

AbstractWe present a statistical study of interplanetary magnetic field (IMF) effects on the magnetopause location based on nearly 15,000 magnetopause crossings identified in the THEMIS A‐E, Magion 4, Geotail, and Interball‐1 satellite data. In order to exclude the dominant controlling role of the solar wind dynamic pressure, differences between observed and empirical model magnetopause distances are analyzed. We find that the magnetopause distance is significantly controlled by the IMF clock angle, not considered in the used empirical models. Furthermore, during very low Alfvén Mach numbers (MA < 4), IMF By component can induce a considerable dawn‐dusk asymmetry of the magnetopause shape. The magnetopause distance is found to be larger in the direction of the IMF vector, that is, not only the magnitude but also the orientation of the IMF By component is eventually important. The obtained results are qualitatively consistent with a magnetohydrodynamic model run at the Community Coordinated Modeling Center.

The Relationship between Solar Wind Charge Exchange Soft X-ray Emission and the Tangent Direction of Magnetopause in an XMM–Newton Event

With the advent of soft X-ray imaging enabling global magnetopause detection, it is critical to use reconstruction techniques to derive the 3-dimensional magnetopause location from 2-dimensional X-ray images. One of the important assumptions adopted by most techniques is that the direction with maximum soft X-ray emission is the tangent direction of the magnetopause, which has not been validated in observation so far. This paper analyzes a magnetospheric solar wind charge exchange (SWCX) soft X-ray event detected by XMM–Newton during relatively stable solar wind and geomagnetic conditions. The tangent direction of the magnetopause is determined by an empirical magnetopause model. Observation results show that the maximum SWCX soft X-ray intensity gradient tends to be the tangent of the magnetopause’s inner boundary, while the maximum SWCX soft X-ray intensity tends to be the tangent of the magnetopause’s outer boundary. Therefore, it is credible to use the assumption that the tangent direction of the magnetopause is the maximum SWCX soft X-ray intensity or its gradient when reconstructing the 3-dimensional magnetopause location. In addition, since these two maxima tend to be the inner and outer boundaries of the magnetopause, the thickness of magnetopause can also be revealed by soft X-ray imaging.

Quantifying the global solar wind-magnetosphere interaction with the Solar-Terrestrial Observer for the Response of the Magnetosphere (STORM) mission concept

Much of what we know about the solar wind’s interaction with the Earth’s magnetosphere has been gained from isolated in-situ measurements by single or multiple spacecraft. Based on their observations, we know that reconnection, whether on the dayside magnetopause or deep within the Earth’s magnetotail, controls the bulk flow of solar wind energy into and through the global system and that nightside activity provides the energized particles that power geomagnetic storms. But by their very nature these isolated in-situ measurements cannot provide an instantaneous global view of the entire system or its cross-scale dynamics. To fully quantify the dynamics of the coupled solar wind-magnetosphere requires comprehensive end-to-end global imaging of the key plasma structures that comprise the magnetosphere which have spatial resolutions that exceeds anything possible with multi-point or constellation situ measurements. Global, end-to-end, imaging provides the pathway to understanding the system as a whole, its constituent parts, and its cross-scale processes on a continuous basis, as needed to quantify the flow of solar wind energy through the global magnetospheric system. This paper describes how a comprehensively-instrumented single spacecraft in a high-altitude, high-inclination orbit coupled with ground-based instruments provides the essential observations needed to track and quantify the flow of solar wind energy through the magnetosphere. This includes observations of the solar wind plasma and magnetic field input, the magnetopause location in soft X-rays, the auroral oval in far ultraviolet, the ring current in energetic neutrals, the plasmasphere in extreme ultraviolet, the exosphere in Lyman-α, and the microstructure of the nightside auroral oval from ground-based all sky cameras.

Juno Magnetometer Observations at Ganymede: Comparisons With a Global Hybrid Simulation and Indications of Magnetopause Reconnection

AbstractJuno's flyby of Ganymede on 7 June 2021, provides a unique opportunity to explore the moon's magnetosphere. By means of Magnetometer (MAG) observations and a hybrid numerical simulation, we provide a global description of this environment and analyze the upstream magnetopause in detail. In particular, LATMOS Hybrid Simulation results suggest Juno observed open (one foot‐point connected to Ganymede, the other to Jupiter) and closed magnetic field lines along its trajectory. Additionally, we determine that the upstream magnetopause location and orientation seen by Juno MAG are consistent with previous Galileo observations. We observe a non‐zero normal component of the magnetic field across the magnetopause, along with flux rope structures embedded within the boundary's current sheet. Both signatures are strong indicators that magnetic reconnection was occurring along Ganymede's magnetopause. Based on this crossing, we calculate a shear angle of ∼78° across the magnetopause, with a dimensionless reconnection rate of ∼0.12.

Statistics of geomagnetic storms: Global simulations perspective

We present results of 131 geomagnetic storm simulations using the University of Michigan Space Weather Modeling Framework Geospace configuration. We compare the geomagnetic indices derived from the simulation with those observed, and use 2D cuts in the noon-midnight planes to compare the magnetopause locations with empirical models. We identify the location of the current sheet center and look at the plasma parameters to deduce tail dynamics. We show that the simulation produces geomagnetic index distributions similar to those observed, and that their relationship to the solar wind driver is similar to that observed. While the magnitudes of the Dst and polar cap potentials are close to those observed, the simulated AL index is consistently underestimated. Analysis of the magnetopause position reveals that the subsolar position agrees well with an empirical model, but that the tail flaring in the simulation is much smaller than that in the empirical model. The magnetotail and ring currents are closely correlated with the Dst index, and reveal a strong contribution of the tail current beyond 8REto the Dst index during the storm main phase.

Laser-driven, ion-scale magnetospheres in laboratory plasmas. II. Particle-in-cell simulations

Ion-scale magnetospheres have been observed around comets, weakly magnetized asteroids, and localized regions on the Moon and provide a unique environment to study kinetic-scale plasma physics, in particular in the collision-less regime. In this work, we present the results of particle-in-cell simulations that replicate recent experiments on the large plasma device at the University of California, Los Angeles. Using high-repetition rate lasers, ion-scale magnetospheres were created to drive a plasma flow into a dipolar magnetic field embedded in a uniform background magnetic field. The simulations are employed to evolve idealized 2D configurations of the experiments, study highly resolved, volumetric datasets, and determine the magnetospheric structure, magnetopause location, and kinetic-scale structures of the plasma current distribution. We show the formation of a magnetic cavity and a magnetic compression in the magnetospheric region, and two main current structures in the dayside of the magnetic obstacle: the diamagnetic current, supported by the driver plasma flow, and the current associated with the magnetopause, supported by both the background and driver plasmas with some time-dependence. From multiple parameter scans, we show a reflection of the magnetic compression, bounded by the length of the driver plasma, and a higher separation of the main current structures for lower dipolar magnetic moments.

Massive Multi‐Mission Statistical Study and Analytical Modeling of the Earth's Magnetopause: 4. On the Near‐Cusp Magnetopause Indentation

AbstractThe shape and location of the magnetopause current sheet in the near‐cusp region is still a debated question. Over time, several observations led to contradictory conclusions regarding the presence of an indentation of the magnetopause in that region. As a result several empirical models consider the surface is indented in that region, while some others do not. To tackle this issue, we fit a total of 17,230 magnetopause crossings to various indented and non‐indented analytical models. The results show that while all models describe the magnetopause position and shape equivalently far from the cusp region, the non‐indented version over‐estimate the radial position of the near‐cusp magnetopause. Among indented models, we show that the one designed from non‐linearmagneto hydrodynamic simulations fits well the near‐cusp magnetopause location, while the other underestimate its position probably because their design was possibly based on magnetopause crossing catalogs that contain cusp inner boundary crossings.

Massive Multi‐Mission Statistical Study and Analytical Modeling of the Earth's Magnetopause: 2. Shape and Location

AbstractThe Earth magnetopause is the boundary between the magnetosphere and the shocked solar wind. Its location and shape are primarily determined by the properties of the solar wind and interplanetary magnetic field (IMF) but the nature of the control parameters and to what extent they impact the stand‐off distance, the flaring, and the symmetries, on the dayside and night side, is still not well known. We present a large statistical study of the magnetopause location and shape based an extensive multi‐mission magnetopause database, cumulating 17,230 crossings on 17 different spacecraft, from the dayside to lunar nightside distances. The results confirm the power‐law dependency of the stand‐off position on the solar wind pressure. The IMF clock angle itself (all amplitudes combined) does not impact the stand‐off distance, nor does the cone angle. However, the magnetopause is found to move Earthward as the IMF gets stronger and more southward. All upstream conditions combined, it is found that the function used at the root of several analytical models still holds at lunar distances. We find that the equatorial flaring is larger than the meridional one. However, the meridional flaring is found to depend on the seasonal tilt conditions, being larger in the summer hemisphere. The flaring is also found to depend on the IMF clock angle. Meridional flaring increases as the IMF turns south and is then larger than the equatorial flaring. The equatorial flaring barely changes or weakly increases as the IMF turns northward, and is larger than the meridional flaring for northward conditions. The results of the study pave the way for the elaboration of a new analytical empirical expression of the magnetopause location and shape.

Massive Multi‐Mission Statistical Study and Analytical Modeling of the Earth's Magnetopause: 3. An Asymmetric Non Indented Magnetopause Analytical Model

AbstractIn a companion statistical study, we showed that the expression of the magnetopause surface as a power law of an elliptic function of the zenith angle θ holds at lunar distances, that the flaring of the magnetopause surface is influenced by the Interplanetary Magnetic Field (IMF) By component and that the IMF Bx component had no influence on the stand‐off distance. As a follow‐up to these statistical results, this paper presents a new empirical analytical asymmetric and non‐indented model of the magnetopause location and shape. This model is obtained from fitting of 15,349 magnetopause crossings using 17 different spacecraft and is parametrized by the upstream solar wind dynamic and magnetic pressures, the IMF clock angle and the Earth dipole tilt angle. The constructed model provides a more accurate prediction of the magnetopause surface location than current Magnetopause surface models, especially on the night side of the magnetosphere.

Validation of the use of THEMIS-B and THEMIS-C as a near-Earth solar wind monitor

Two THEMIS (Time History of Events and Macroscale Interactions during Substorms) spacecraft, B and C, began orbiting the Moon in 2011 and have since provided routine measurements of the plasma conditions in the lunar orbit. In this study, we systematically compare these measurements in near-Earth space with solar wind measurements obtained from the Lagrangian 1 (L1) point and propagated to the Earth, including measurements in the OMNI database and from the Wind spacecraft. A statistical comparison between THEMIS data and data from the OMNI database from September 2011 to December 2017 showed that the <em>Y</em> and <em>Z</em> components of the magnetic field and the flow speed are generally consistent between the two data sets. The ion number density and the dynamic pressure measured by THEMIS in near-Earth space are lower than those in the OMNI database, suggesting possible variation in the solar wind environment while propagating from the L1 point to near-Earth space. We further show two examples in which near-Earth solar wind measurements are applied in calculating the magnetopause location and in quantifying the magnetic field response to interplanetary shocks. Both examples suggest that using solar wind data from near-Earth space achieves better results than using solar wind data from the L1 point. These results provide validation of THEMIS-B and THEMIS-C as an alternative monitor of the near-Earth solar wind environment.