skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Predicting Geostationary (GOES) 4.1–30 keV Electron Flux Over All MLT Using LEEMYR Regression Models
Abstract Regression models (LEEMYR: Low Energy Electron MLT geosYnchronous orbit Regression) predict hourly 4.1–30 keV electron flux at geostationary orbit (GOES‐16) using solar wind, IMF, and geomagnetic index parameters. Multiplicative interaction and polynomial terms describe synergistic and nonlinear effects. We reduce predictors to an optimal set using stepwise regression, resulting in models with validation comparable to a neural network. Models predict 1, 3, 6, 12, and 24 hr into the future. Validation correlations are as high as 0.78 (4.1 and 11 keV, 1 hr prediction) and Heidke Skill scores (HSS) up to 0.66. A 3 hr ahead prediction is more practical, with slightly lower validation correlation (0.75) and HSS (0.61). The addition of location (MLT: magnetic local time) as a covariate, including multiplicative interaction terms, accounts for location‐dependent flux differences and variation of parameter influence, and allows prediction over the full orbit. Adding a substorm index (SME) provides minimal increase in validation correlation (0.81) showing that other parameters are good proxies for an unavailable real time substorm index. Prediction intervals on individual values provide more accurate assessments of model quality than confidence intervals on the mean values. An inverse N‐weighted least squares approach is impractical as it increases false positive warnings. Physical interpretations are not possible as spurious correlations due to common cycles are not removed. However,SME, Bz, Kp, and Dst are the highest correlates of electron flux, with solar wind velocity, density, and pressure, and IMF magnitude being less well correlated.  more » « less
Award ID(s):
2246912
PAR ID:
10589823
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Space Weather
Volume:
22
Issue:
8
ISSN:
1542-7390
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Electron fluxes (20 eV–2 MeV, RBSP‐A satellite) show reasonable simple correlation with a variety of parameters (solar wind, IMF, substorms, ultralow frequency (ULF) waves, geomagnetic indices) over L‐shells 2–6. Removing correlation‐inflating common cycles and trends (using autoregressive and moving average terms in an ARMAX analysis) results in a 10 times reduction in apparent association between drivers and electron flux, although many are still statistically significant (p < 0.05). Corrected influences are highest in the 20 eV–1 keV and 1–2 MeV electrons, more modest in the midrange (2–40 keV). Solar wind velocity and pressure (but not number density), IMF magnitude (with lower influence ofBz), SME (a substorm measure), a ULF wave index, and geomagnetic indices Kp and SymH all show statistically significant associations with electron flux in the corrected individual ARMAX analyses. We postulate that only pressure, ULF waves, and substorms are direct drivers of electron flux and compare their influences in a combined analysis. SME is the strongest influence of these three, mainly in the eV and MeV electrons. ULF is most influential on the MeV electrons. Pressure shows a smaller positive influence and some indication of either magnetopause shadowing or simply compression on the eV electrons. While strictly predictive models may improve forecasting ability by including indirect driver and proxy parameters, and while these models may be made more parsimonious by choosing not to explicitly model time series behavior, our present analyses include time series variables in order to draw valid conclusions about the physical influences of exogenous parameters. 
    more » « less
  2. Abstract The magnetotail lobe region at Mercury is characterized by low plasma density and low magnetic field variability compared to the nightside magnetosheath and central plasma sheet. At Mercury, as well as other planets, lobe magnetic fields play a crucial role in storing and releasing magnetic flux in response to changing upstream solar wind conditions such as interplanetary magnetic field (IMF) orientation and solar wind dynamic pressure (Pdyn). This makes the region significant for studying the magnetospheric interaction with the intense solar wind conditions at Mercury's orbit. Here, we identify and analyze magnetotail lobe observations made by the Mercury Surface, Space Environment, Geochemistry and Ranging (MESSENGER) spacecraft during its 4 years orbital phase. We empirically determined a set of criteria using magnetometer (MAG) and the Fast Imaging Plasma Spectrometer instruments onboard MESSENGER to identify lobe magnetic field intervals. From 3,332 MESSENGER orbits, we identify 1,242 lobe field intervals. We derive an expression for the average lobe magnetic field strength in nanotesla with respect to radial distance downtail:Blobe(r) = (135 ± 8) * r(−2.1±0.3) + (31 ± 8). The lobe magnetic field exhibits both small‐scale (∼3 min) and orbit‐to‐orbit (∼8–12 hr) variability in magnetic field strength compared to this averaged field strength expression. The orbit‐to‐orbit variability in lobe field strength is not significantly correlated with estimated IMF orientation, but is directly correlated withPdyn. Thus, our findings provide evidence for the pressure balance between the inward facingPdynon the nightside magnetopause and the outward facing magnetic pressure supplied by the lobes. 
    more » « less
  3. Pulsating aurora are common diffuse-like aurora. Studies have suggested that they contain higher energy particles than other types and are possibly linked to substorm activity. There has yet to be a quantitative statistical study of the variation in pulsating aurora energy content related to substorms. We analyzed the inverted energy content from 53 events using the Poker Flat Incoherent Scatter Radar. To reduce the uncertainty, we split the differential energy flux into low and high energy using the limit of 30 keV. We also analyzed the lower altitude boundary of the electron density profile, characterized by a number density of > 1 0 10 m −3 , and used this as a proxy for high energy. We compared both of these to magnetic local time (MLT), AE index, and temporal proximity to substorm onset. There was a slight trend in MLT, but a much stronger one in relation to both substorm onset and AE index. For higher AE and closer to onset the total energy flux and flux above 30 keV increased. In addition, this higher energy remained enhanced for an hour after substorm onset. Our results confirm the high energy nature of pulsating aurora, demonstrate the connection to substorms, and imply their importance to coupling between the magnetosphere and atmosphere. 
    more » « less
  4. Abstract During 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. 
    more » « less
  5. Abstract Radiation Belt Storm Probes (RBSP) data show that seed electrons generated by sub‐storm injections play a role in amplifying chorus waves in the magnetosphere. The wave‐particle interaction leads to rapid heating and acceleration of electrons from 10's of keV to 10's of MeV energies. In this work, we examined the changes in the radiation belt during geomagnetic storm events by studying the RBSP REPT, solar wind, AL, SML, and Dst data in conjunction with the WINDMI model of the magnetosphere. The field‐aligned current output from the model is integrated to generate a proxy E index for various energy bands. These E indices track electron energization from 40 KeV to 20 MeV in the radiation belts. The indices are compared to RBSP data and GOES data. Our proxy indices correspond well to the energization data for electron energy bands between 1.8 and 7.7 MeV. Each E index has a unique empirical loss rate term (τL), an empirical time delay term (τD), and a gain value, that are fit to the observations. These empirical parameters were adjusted to examine the delay and charging rates associated with different energy bands. We observed that theτLandτDvalues are clustered for each energy band.τLandτDconsistently increase going from 1.8 to 7.7 MeV in electron energy fluxEeand the dropout interval increases with increasing energy level. The average trend of ΔτD/ΔEewas 4.1 hr/MeV and the average trend of ΔτL/ΔEewas 2.82 hr/MeV. 
    more » « less