-
Auroral, Ionospheric and Ground Magnetic Signatures of Magnetopause Surface Modes
Authors:
M. O. Archer,
M. D. Hartinger,
L. Rastaetter,
D. J. Southwood,
M. Heyns,
J. W. B. Eggington,
A. N. Wright,
F. Plaschke,
X. Shi
Abstract:
Surface waves on Earth's magnetopause have a controlling effect upon global magnetospheric dynamics. Since spacecraft provide sparse in situ observation points, remote sensing these modes using ground-based instruments in the polar regions is desirable. However, many open conceptual questions on the expected signatures remain. Therefore, we provide predictions of key qualitative features expected…
▽ More
Surface waves on Earth's magnetopause have a controlling effect upon global magnetospheric dynamics. Since spacecraft provide sparse in situ observation points, remote sensing these modes using ground-based instruments in the polar regions is desirable. However, many open conceptual questions on the expected signatures remain. Therefore, we provide predictions of key qualitative features expected in auroral, ionospheric, and ground magnetic observations through both magnetohydrodynamic theory and a global coupled magnetosphere-ionosphere simulation of a magnetopause surface eigenmode. These show monochromatic oscillatory field-aligned currents, due to both the surface mode and its non-resonant Alfvén coupling, are present throughout the magnetosphere. The currents peak in amplitude at the equatorward edge of the magnetopause boundary layer, not the open-closed boundary as previously thought. They also exhibit slow poleward phase motion rather than being purely evanescent. We suggest the upward field-aligned current perturbations may result in periodic auroral brightenings. In the ionosphere, convection vortices circulate the poleward moving field-aligned current structures. Finally, surface mode signals are predicted in the ground magnetic field, with ionospheric Hall currents rotating perturbations by approximately (but not exactly) 90° compared to the magnetosphere. Thus typical dayside magnetopause surface modes should be strongest in the East-West ground magnetic field component. Overall, all ground-based signatures of the magnetopause surface mode are predicted to have the same frequency across L-shells, amplitudes that maximise near the magnetopause's equatorward edge, and larger latitudinal scales than for field line resonance. Implications in terms of ionospheric Joule heating and geomagnetically induced currents are discussed.
△ Less
Submitted 2 March, 2023;
originally announced March 2023.
-
How a realistic magnetosphere alters the polarizations of surface, fast magnetosonic, and Alfvén waves
Authors:
M. O. Archer,
D. J. Southwood,
M. D. Hartinger,
L. Rastaetter,
A. N. Wright
Abstract:
System-scale magnetohydrodynamic (MHD) waves within Earth's magnetosphere are often understood theoretically using box models. While these have been highly instructive in understanding many fundamental features of the various wave modes present, they neglect the complexities of geospace such as the inhomogeneities and curvilinear geometries present. Here we show global MHD simulations of resonant…
▽ More
System-scale magnetohydrodynamic (MHD) waves within Earth's magnetosphere are often understood theoretically using box models. While these have been highly instructive in understanding many fundamental features of the various wave modes present, they neglect the complexities of geospace such as the inhomogeneities and curvilinear geometries present. Here we show global MHD simulations of resonant waves impulsively-excited by a solar wind pressure pulse. Although many aspects of the surface, fast magnetosonic (cavity/waveguide), and Alfvén modes present agree with the box and axially symmetric dipole models, we find some predictions for large-scale waves are significantly altered in a realistic magnetosphere. The radial ordering of fast mode turning points and Alfvén resonant locations may be reversed even with monotonic wave speeds. Additional nodes along field lines that are not present in the displacement/velocity occur in both the perpendicular and compressional components of the magnetic field. Close to the magnetopause the perpendicular oscillations of the magnetic field have the opposite handedness to the velocity. Finally, widely-used detection techniques for standing waves, both across and along the field, can fail to identify their presence. We explain how all these features arise from the MHD equations when accounting for a non-uniform background field and propose modified methods which might be applied to spacecraft observations.
△ Less
Submitted 31 January, 2022;
originally announced January 2022.
-
Magnetopause ripples going against the flow form azimuthally stationary surface waves
Authors:
M. O. Archer,
M. D. Hartinger,
F. Plaschke,
D. J. Southwood,
L. Rastaetter
Abstract:
Surface waves process the turbulent disturbances which drive dynamics in many space, astrophysical and laboratory plasma systems, with the outer boundary of Earth's magnetosphere, the magnetopause, providing an accessible environment to study them. Like waves on water, magnetopause surface waves are thought to travel in the direction of the driving solar wind, hence a paradigm in global magnetosph…
▽ More
Surface waves process the turbulent disturbances which drive dynamics in many space, astrophysical and laboratory plasma systems, with the outer boundary of Earth's magnetosphere, the magnetopause, providing an accessible environment to study them. Like waves on water, magnetopause surface waves are thought to travel in the direction of the driving solar wind, hence a paradigm in global magnetospheric dynamics of tailward propagation has been well-established. Here we show through multi-spacecraft observations, global simulations, and analytic theory that the lowest-frequency impulsively-excited magnetopause surface waves, with standing structure along the terrestrial magnetic field, propagate against the flow outside the boundary. Across a wide local time range (09-15h) the waves' Poynting flux exactly balances the flow's advective effect, leading to no net energy flux and thus stationary structure across the field also. Further down the equatorial flanks, however, advection dominates hence the waves travel downtail, seeding fluctuations at the resonant frequency which subsequently grow in amplitude via the Kelvin-Helmholtz instability and couple to magnetospheric body waves. This global response, contrary to the accepted paradigm, has implications on radiation belt, ionospheric, and auroral dynamics and potential applications to other dynamical systems.
△ Less
Submitted 6 October, 2021;
originally announced October 2021.
-
The STONE curve: A ROC-derived model performance assessment tool
Authors:
Michael W. Liemohn,
Abigail R. Azari,
Natalia Ganushkina,
Lutz Rastaetter
Abstract:
A new model validation and performance assessment tool is introduced, the sliding threshold of observation for numeric evaluation (STONE) curve. It is based on the relative operating characteristic (ROC) curve technique, but instead of sorting all observations in a categorical classification, the STONE tool uses the continuous nature of the observations. Rather than defining events in the observat…
▽ More
A new model validation and performance assessment tool is introduced, the sliding threshold of observation for numeric evaluation (STONE) curve. It is based on the relative operating characteristic (ROC) curve technique, but instead of sorting all observations in a categorical classification, the STONE tool uses the continuous nature of the observations. Rather than defining events in the observations and then sliding the threshold only in the classifier (model) data set, the threshold is changed simultaneously for both the observational and model values, with the same threshold value for both data and model. This is only possible if the observations are continuous and the model output is in the same units and scale as the observations, that is, the model is trying to exactly reproduce the data. The STONE curve has several similarities with the ROC curve, plotting probability of detection against probability of false detection, ranging from the (1,1) corner for low thresholds to the (0,0) corner for high thresholds, and values above the zero-intercept unity-slope line indicating better than random predictive ability. The main difference is that the STONE curve can be nonmonotonic, doubling back in both the x and y directions. These ripples reveal asymmetries in the data-model value pairs. This new technique is applied to modeling output of a common geomagnetic activity index as well as energetic electron fluxes in the Earth's inner magnetosphere. It is not limited to space physics applications but can be used for any scientific or engineering field where numerical models are used to reproduce observations.
△ Less
Submitted 22 April, 2020;
originally announced May 2020.
-
New magnetohydrodynamic model available at NASA Community Coordinated Modeling Center
Authors:
Ilja Honkonen,
Lutz Rastätter,
Alex Glocer
Abstract:
The Community Coordinated Modeling Center (CCMC) at NASA Goddard Space Flight Center is a multi-agency partnership to enable, support and perform research and development for next-generation space science and space weather models. CCMC currently hosts nearly 100 numerical models and a cornerstone of this activity is the Runs on Request (RoR) system which allows anyone to request a model run and an…
▽ More
The Community Coordinated Modeling Center (CCMC) at NASA Goddard Space Flight Center is a multi-agency partnership to enable, support and perform research and development for next-generation space science and space weather models. CCMC currently hosts nearly 100 numerical models and a cornerstone of this activity is the Runs on Request (RoR) system which allows anyone to request a model run and analyze/visualize the results via a web browser. CCMC is also active in the education community by organizing student research contests, heliophysics summer schools, and space weather forecaster training for students, government and industry representatives.
We present a generic magnetohydrodynamic (MHD) model - PAMHD - that has been added to the CCMC RoR system which allows the study of a variety of fluid and plasma phenomena in one, two and three dimensions using a dynamic point-and-click web interface. Flexible initial and boundary conditions allow experimentation with a variety of plasma physics problems such as shocks, instabilities, planetary magnetospheres and astrophysical systems. Experimentation with numerical effects, e.g. resolution, solution method and boundary conditions, is also possible and can provide valuable context for space weather forecasters when interpreting observations or modeling results.
We present an overview of the C++ implementation and show example results obtained through the CCMC RoR system, including the first to our knowledge MHD simulation of the interaction of the magnetospheres of Jupiter and Saturn in two dimensions.
△ Less
Submitted 21 March, 2017;
originally announced March 2017.
-
Ensemble modeling of CMEs using the WSA-ENLIL+Cone model
Authors:
M. L. Mays,
A. Taktakishvili,
A. A. Pulkkinen,
P. J. MacNeice,
L. Rastaetter,
D. Odstrcil,
L. K. Jian,
I. G. Richardson,
J. A. LaSota,
Y. Zheng,
M. M. Kuznetsova
Abstract:
Ensemble modeling of CMEs provides a probabilistic forecast of CME arrival time which includes an estimation of arrival time uncertainty from the spread and distribution of predictions and forecast confidence in the likelihood of CME arrival. The real-time ensemble modeling of CME propagation uses the WSA-ENLIL+Cone model installed at the CCMC and executed in real-time. The current implementation…
▽ More
Ensemble modeling of CMEs provides a probabilistic forecast of CME arrival time which includes an estimation of arrival time uncertainty from the spread and distribution of predictions and forecast confidence in the likelihood of CME arrival. The real-time ensemble modeling of CME propagation uses the WSA-ENLIL+Cone model installed at the CCMC and executed in real-time. The current implementation evaluates the sensitivity of WSA-ENLIL+Cone model simulations of CME propagation to initial CME parameters. We discuss the results of real-time ensemble simulations for a total of 35 CME events between January 2013 - July 2014. For the 17 events where the CME was predicted to arrive at Earth, the mean absolute arrival time prediction error was 12.3 hours, which is comparable to the errors reported in other studies. For predictions of CME arrival at Earth the correct rejection rate is 62% and the false-alarm rate is 38%. The arrival time was within the range of the ensemble arrival predictions for 8 out of 17 events. The Brier Score for CME arrival predictions is 0.15 (where 1 is a perfect forecast), indicating that on average, the predicted likelihood of CME arrival is fairly accurate. The reliability of ensemble CME arrival predictions is heavily dependent on the initial distribution of CME input parameters, particularly the median and spread. Preliminary analysis of the probabilistic forecasts suggests undervariability, indicating that these ensembles do not sample a wide enough spread in CME input parameters. Prediction errors can also arise from ambient model parameters, the accuracy of the solar wind background derived from coronal maps, or other model limitations. Finally, predictions of the Kp geomagnetic index differ from observed values by less than one for 11 out of 17 of the ensembles and Kp prediction errors computed from the mean predicted Kp show a mean absolute error of 1.3.
△ Less
Submitted 11 May, 2015; v1 submitted 16 April, 2015;
originally announced April 2015.