Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geomagnetic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Solar Energetic Proton Access to the Inner Magnetosphere During the September 7–8, 2017 Event

Here, the access of solar energetic protons into the inner magnetosphere on September 7–8, 2017 is investigated by following reversed proton trajectories to compute the proton cutoff energy using the Dartmouth geomagnetic cutoff code (Kress et al., 2010, https://doi.org/10.1029/2009sw000488). The cutoff energies for protons coming from the west and east direction, the minimum and maximum cutoff energy respectively, are calculated every 5 min along the orbit of Van Allen Probes using TS07 and the Lyon-Fedder-Mobarry (LFM) MHD magnetic field model. The result shows that the cutoff energy increases significantly as the radial distance decreases, and that the cutoff energy decreases with the building up of the ring current during magnetic storms. Solar wind dynamic pressure also affects cutoff suppression (Kress et al., 2004, https://doi.org/10.1029/2003gl018599). The LFM-RCM model shows stronger suppression of cutoff energy than TS07 during strong solar wind driving conditions. The simulation result is compared with proton flux measurements, showing consistent variation of the cutoff location during the September 7–8, 2017 geomagnetic storm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Examination of Radiation Belt Dynamics During Substorm Clusters: Activity Drivers and Dependencies of Trapped Flux Enhancements

Here, dynamical variations of radiation belt trapped electron fluxes are examined to better understand the variability of enhancements linked to substorm clusters. Analysis is undertaken using the Substorm Onsets and Phases from Indices of the Electrojet substorm cluster algorithm for event detection. Observations from low earth orbit are complemented by additional measurements from medium earth orbit to allow a major expansion in the energy range considered, from medium energy energetic electrons up to ultra-relativistic electrons. The number of substorms identified inside a cluster does not depend strongly on solar wind drivers or geomagnetic indices either before, during, or after the cluster start time. Clusters of substorms linked to moderate (100 nT < AE ≤ 300 nT) or strong AE (AE ≥ 300 nT) disturbances are associated with radiation belt flux enhancements, including up to ultra-relativistic energies by the strongest substorms (as measured by strong southward Bz and high AE). These clusters reliably occur during times of high speed solar winds streams with associated increased magnetospheric convection. However, substorm clusters associated with quiet AE disturbances (AE ≤ 100 nT) lead to no significant chorus whistler mode intensity enhancements, or increases in energetic, relativistic, or ultra-relativistic electron flux in the outer radiation belts. In these cases the solar wind speed is low, and the geomagnetic Kp index indicates a lack of magnetospheric convection. Our study clearly indicates that clusters of substorms occurring outside of high speed wind streams are not by themselves sufficient to drive acceleration, which may be due to the lack of pre-cluster convection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Inner Belt Electron Decay Timescales: A Comparison of Van Allen Probes and DREAM3D Losses Following the June 2015 Storm

Abstract NASA's Van Allen Probes observed significant, long‐lived fluxes of inner belt electrons up to ∼1 MeV after geomagnetic storms in March and June 2015. Reanalysis of Magnetic Electron Ion Spectrometer (MagEIS) data with improved background correction showed a clearer picture of the relativistic electron population that persisted through 2016 and into 2017 above the Fennell et al. (2015, https://doi.org/10.1002/2014gl062874 ) limit. The intensity and duration of these enhancements allow estimation of decay timescales for comparison with simulated decay rates and theoretical lifetimes. We compare decay timescales from these data and DREAM3D simulations based on them using geomagnetic activity‐dependent pitch angle diffusion coefficients derived from plasmapause‐indexed wave data (Malaspina et al., 2016, https://doi.org/10.1002/2016gl069982 , 2018, https://doi.org/10.1029/2018gl078564 ) and phase space densities derived from MagEIS observations. Simulated decay rates match observed decay rates more closely than the theoretical lifetime due to significantly nonequilibrium pitch angle distributions in simulation and data. We conclude that nonequilibrium effects, rather than a missing diffusion or loss process, account for observed short decay rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multitaper Magnitude‐Squared Coherence for Time Series With Missing Data: Understanding Oscillatory Processes Traced by Multiple Observables

To explore the hypothesis of a common source of variability in two time series, observers may estimate the magnitude-squared coherence (MSC), which is a frequency-domain view of the cross correlation. For time series that do not have uniform observing cadence, MSC can be estimated using Welch's overlapping segment averaging. However, multitaper has superior statistical properties to Welch's method in terms of the tradeoff between bias, variance, and bandwidth. The classical multitaper technique has recently been extended to accommodate time series with underlying uniform observing cadence from which some observations are missing. This situation is common for solar and geomagnetic data sets, which may have gaps due to breaks in satellite coverage, instrument downtime, or poor observing conditions. We demonstrate the scientific use of missing-data multitaper magnitude-squared coherence by detecting known solar mid-term oscillations in simultaneous, missing-data time series of solar Lyman α flux and geomagnetic Disturbance Storm Time index. Due to their superior statistical properties, we recommend that multitaper methods be used for all heliospheric time series with underlying uniform observing cadence.

Astro-statistics techniques (1886)↗
DOE OSTI2016

Energy Estimation of Cosmic Rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy density is determined from the radio pulses at each observer position and is interpolated using a two dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge excess emission components. We found that the spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy – corrected for geometrical effects – is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. Finally we find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Constraining the p¯/p ratio in TeV cosmic rays with observations of the Moon shadow by HAWC

An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed toward the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a function of the energy and charge of the cosmic rays. The displacement provides a natural tool for momentum/charge discrimination that can be used to study the composition of cosmic rays. Using 33 months of data comprising more than 80 billion cosmic rays measured by the High Altitude Water Cherenkov observatory, we have analyzed the Moon shadow to search for TeV antiprotons in cosmic rays. Here, we present our first upper limits on the p¯/p fraction, which in the absence of any direct measurements provide the tightest available constraints of ~1 % on the antiproton fraction for energies between 1 and 10 TeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Trial Implementation of a Multihazard Risk Assessment Framework for High-Impact Low-Frequency Power Grid Events

The Pacific Northwest National Laboratory developed a risk framework for modeling high-impact, low-frequency power grid events to support risk-informed decisions. In this research article, we briefly recap the framework and demonstrate its implementation for seismic and geomagnetic hazards using a benchmark reliability test system. We describe integration of a collection of models implemented to perform hazard analysis, fragility evaluation, consequence estimation and post-event restoration. We demonstrate the value of the framework as a multi-hazard power grid risk assessment and management tool. The research will benefit transmission planners and emergency planners by improving their ability to maintain a resilient grid infrastructure against impacts from major events.

power grid, risk assessment, seismic hazard, dynam↗
DOE OSTI2019

Volcanic stratigraphy and age model of the Kimama deep borehole (Project Hotspot): Evidence for 5.8 million years of continuous basalt volcanism, central Snake River Plain, Idaho

The Snake River Plain of central Idaho represents the world’s best example of a mantle hotspot track impinging upon continental crust and provides a record of bimodal volcanism extending over 12 Ma to the present. Project Hotspot recovered almost 2 km of continuous drill core from the Kimama borehole, located in central Idaho on the axial volcanic zone. The Kimama drill core represents the most complete record of mafic volcanism along the Yellowstone–Snake River Plain hotspot track. A total of 432 basalt flow units, representing 183 basalt flows, 78 basalt flow groups, and 34 super groups, along with 42 sediment interbeds are recognized using volcanic facies observations, stratigraphic relationships, borehole geophysical logs, and paleosecular variation in magnetostratigraphy. Rhyolite and other non-basaltic volcanic materials were not encountered in the drill core. Ages for six basalt lava flows were determined by 40 Ar/ 39 Ar using incremental heating experiments. Paleomagnetic inclination was measured on over 1200 samples collected at roughly 2-m-depth intervals, yielding mean values of paleosecular variation between ±50° to ±70° in Kimama flow groups, close to the expected 61° axial dipole average for the Kimama borehole location. Twenty-three magnetic reversals were identified and correlated to dated geomagnetic chrons and subchrons and compared with the 40 Ar/ 39 Ar radiometric ages. A linear fit to 40 Ar/ 39 Ar dates, geomagnetic chron and subchron boundaries, and volcanogenic zircon U-Pb ages defines a mean accumulation rate of ~320 m/m.y. and extrapolates to a bottom hole age of 6.3 Ma. Average thicknesses of lithologic units increase from 2.7 m (sediment), 4 m (flow units), 10 m (flows), 23 m (flow groups), to 53 m (super groups). On average, one lava flow inundated the Kimama borehole location every 33 k.y. Intercalated sediments, ranging from 0.06 to 24.5 m thick, make up roughly 6% of the drill core and indicate lulls in local volcanic activity that may have lasted up to 77 k.y. Neutron and gamma-ray logs supplement observations from the drill cores: neutron logs document individual flow units through the contrast between massive flow interiors and more porous flow surfaces, and gamma-ray logs document the depth and thickness of sedimentary interbeds and high–K-Fe basalts. The 5.8 m.y. duration of basaltic volcanism in the Kimama drill core implies a steady rate of volcanism, indicating a relatively stable rate of mantle upflow along the lithosphere-mantle boundary in the wake of Yellowstone–Snake River Plain plume volcanism.

58 GEOSCIENCES↗
DOE OSTI2015

Modeling of different scenarios of thin current sheet equilibria in the Earth’s magnetotail

The Earth’s magnetosphere is an open dynamic system permanently interacting with the solar wind, i.e., the plasma flow from the Sun. Some plasma processes in the magnetosphere are of spontaneous explosive character, while others develop rather slowly as compared to the characteristic times of plasma particle motion in it. The large-scale current sheet in the magnetotail can be in an almost equilibrium state both in quiet periods and during geomagnetic perturbations, and its variations can be considered quasistatic. Thus, under some conditions, the magnetotail current sheet can be described as an equilibrium plasma system. Its state depends on various parameters, in particular, on those determining the dynamics of charged particles. Knowing the main governing parameters, one can study the structure and properties of the current sheet equilibrium. This work is devoted to the self-consistent modeling of the equilibrium thin current sheet (TCS) of the Earth’s magnetotail, the thickness of which is comparable with the ion gyroradius. The main objective of this work is to examine how the TCS structure depends on the parameters characterizing the particle dynamics and magnetic field geometry. A numerical hybrid self-consistent TCS model in which the tension of magnetic field lines is counterbalanced by the inertia of ions moving through the sheet is constructed. The ion dynamics is considered in the quasi-adiabatic approximation, while the electron motion, in the conductive fluid approximation. Depending on the values of the adiabaticity parameter κ (which determines the character of plasma particle motion) and the dimensionless normal component of the magnetic field b{sub z}, the following two scenarios are considered: (A) the adiabaticity parameter is proportional to the particle energy and b{sub z} = const and (B) the particle energy is fixed and the adiabaticity parameter is proportional to b{sub z}. The structure of the current sheet and particle dynamics in it are studied as functions of the parameters κ and b{sub z}. It is shown that, in scenario A, the current sheet thickness decreases with increasing adiabaticity parameter due to a decrease in the ion gyroradius. Accordingly, the radius of curvature of magnetic field lines decreases, which leads to an increase in the contribution of electron drift currents near the neutral plane z = 0. Numerical simulations demonstrate that current equilibria can exist if the adiabaticity parameter lies in the range 0.05 ≤ κ ≤ 0.7. At κ ∼ 0.7, the contribution of electron drift currents to the total current density is much larger than the contribution of ions and the ion motion becomes chaotic. At larger values of the adiabaticity parameter, no equilibrium solutions were found in the framework of the given one-dimensional model. Therefore, the value κ = 0.7 corresponds to the upper applicability limit of the quasi-adiabatic model of the current sheet. In scenario B, an increase in the parameter κ leads to the appearance of a large number of quasi-trapped ions in the current sheet, due to which the current sheet thickens and the amplitude of the current density decreases. As a result, equilibrium solutions exist in a much narrower range of the adiabaticity parameter, 0.05≤ κ ≤ 0.25. Consequences of the existence of parametric boundaries of equilibrium solutions for the TCS under actual geomagnetic conditions are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Modeling of the propagation of low-frequency electromagnetic radiation in the Earth’s magnetosphere

A numerical algorithm for solving the set of differential equations describing the propagation of low-frequency electromagnetic radiation in the magnetospheric plasma, including in the presence of geomagnetic waveguides in the form of large-scale plasma density inhomogeneities stretched along the Earth’s magnetic field, has been developed. Calculations of three-dimensional ray trajectories in the magnetosphere and geomagnetic waveguide with allowance for radiation polarization have revealed characteristic tendencies in the behavior of electromagnetic parameters along the ray trajectory. The results of calculations can be used for magnetospheric plasma diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Predicting cutoff L-shells of solar protons using the GPPSn particle dataset

Solar energetic protons (SEPs) arriving at the Earth trigger severe radiation storms in the near-Earth space, directly impacting space missions operating at various altitudes. Therefore, monitoring SEP events and predicting the penetration depths of solar protons are critical for aerospace sectors. Building on previous efforts, here we demonstrate the feasibility of using proton measurements from the Global Prompt Proton Sensor network (GPPSn), enabled by Los Alamos National Laboratory developed combined X-ray dosimeters aboard GPS satellites, to characterize and predict the penetration of solar protons into the geomagnetic field. The inclined medium-Earth-orbits (MEOs) of the global GPS constellation offer a unique advantage of allowing simultaneous measurements of penetrating solar protons inside both open- and closed-field line regions. Therefore, the L-profiles of ∼10s–100 MeV solar protons and their associated cutoff L-shells can be determined from the GPPSn dataset, using predefined threshold proton flux values rather than traditional flux ratios. After examining a list of SEP event intervals across solar cycles 23, 24 and 25—including the 2024 Mother’s Day superstorm, we showcase how the latest GPPSn proton dataset (release v1.10), reprocessed and calibrated, can not only be used to monitor solar proton distributions inside the dynamic geomagnetic field for individual events, but also to derive a new empirical model linking cutoff L-shells with several key space weather parameters. This newly developed SEPCL-MEO model demonstrates high predictive performance; for example, predictions for > 30 MeV solar protons yield a correlation coefficient of 0.85 and performance efficiency of 0.67 when validated against GPPSn observations. Results from this pilot study underscores the scientific and operational value of the GPPSn dataset, and this dataset—when paired with machine-learning techniques—can play a critical role in observing and predicting the effects of future incoming SEP events, including extreme ones.

58 GEOSCIENCES↗
DOE OSTI2016

MONITORING SHORT-TERM COSMIC-RAY SPECTRAL VARIATIONS USING NEUTRON MONITOR TIME-DELAY MEASUREMENTS

Neutron monitors (NMs) are ground-based detectors of cosmic-ray showers that are widely used for high-precision monitoring of changes in the Galactic cosmic-ray (GCR) flux due to solar storms and solar wind variations. In the present work, we show that a single neutron monitor station can also monitor short-term changes in the GCR spectrum, avoiding the systematic uncertainties in comparing data from different stations, by means of NM time-delay histograms. Using data for 2007–2014 from the Princess Sirindhorn Neutron Monitor, a station at Doi Inthanon, Thailand, with the world’s highest vertical geomagnetic cutoff rigidity of 16.8 GV, we have developed an analysis of time-delay histograms that removes the chance coincidences that can dominate conventional measures of multiplicity. We infer the “leader fraction” L of neutron counts that do not follow a previous neutron count in the same counter from the same atmospheric secondary, which is inversely related to the actual multiplicity and increases for increasing GCR spectral index. After correction for atmospheric pressure and water vapor, we find that L indicates substantial short-term GCR spectral hardening during some but not all Forbush decreases in GCR flux due to solar storms. Such spectral data from Doi Inthanon provide information about cosmic-ray energies beyond the Earth’s maximum geomagnetic cutoff, extending the reach of the worldwide NM network and opening a new avenue in the study of short-term GCR decreases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Sheath-accumulating Propagation of Interplanetary Coronal Mass Ejection

Fast interplanetary coronal mass ejections (ICMEs) are the drivers of strong space weather storms such as solar energetic particle events and geomagnetic storms. The connection between the space-weather-impacting solar wind disturbances associated with fast ICMEs at Earth and the characteristics of causative energetic CMEs observed near the Sun is a key question in the study of space weather storms, as well as in the development of practical space weather prediction. Such shock-driving fast ICMEs usually expand at supersonic speeds during the propagation, resulting in the continuous accumulation of shocked sheath plasma ahead. In this paper, we propose a “sheath-accumulating propagation” (SAP) model that describes the coevolution of the interplanetary sheath and decelerating ICME ejecta by taking into account the process of upstream solar wind plasma accumulation within the sheath region. Based on the SAP model, we discuss (1) ICME deceleration characteristics; (2) the fundamental condition for fast ICMEs at Earth; (3) the thickness of interplanetary sheaths; (4) arrival time prediction; and (5) the super-intense geomagnetic storms associated with huge solar flares. We quantitatively show that not only the speed but also the mass of the CME are crucial for discussing the above five points. The similarities and differences between the SAP model, the drag-based model, and the“snow-plow” model proposed by Tappin are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Magnetosphere dynamics during the 14 November 2012 storm inferred from TWINS, AMPERE, Van Allen Probes, and BATS-R-US–CRCM

During the 14 November 2012 geomagnetic storm, the Van Allen Probes spacecraft observed a number of sharp decreases ("dropouts") in particle fluxes for ions and electrons of different energies. In this paper, we investigate the global magnetosphere dynamics and magnetosphere–ionosphere (M–I) coupling during the dropout events using multipoint measurements by Van Allen Probes, TWINS, and AMPERE together with the output of the two-way coupled global BATS-R-US–CRCM model. We find different behavior for two pairs of dropouts. For one pair, the same pattern was repeated: (1) weak nightside Region 1 and 2 Birkeland currents before and during the dropout; (2) intensification of Region 2 currents after the dropout; and (3) a particle injection detected by TWINS after the dropout. The model predicted similar behavior of Birkeland currents. TWINS low-altitude emissions demonstrated high variability during these intervals, indicating high geomagnetic activity in the near-Earth tail region. For the second pair of dropouts, the structure of both Birkeland currents and ENA emissions was relatively stable. The model also showed quasi-stationary behavior of Birkeland currents and simulated ENA emissions with gradual ring current buildup. We confirm that the first pair of dropouts was caused by large-scale motions of the OCB (open–closed boundary) during substorm activity. We show the new result that this OCB motion was associated with global changes in Birkeland (M–I coupling) currents and strong modulation of low-altitude ion precipitation. The second pair of dropouts is the result of smaller OCB disturbances not related to magnetospheric substorms. The local observations of the first pair of dropouts result from a global magnetospheric reconfiguration, which is manifested by ion injections and enhanced ion precipitation detected by TWINS and changes in the structure of Birkeland currents detected by AMPERE. This study demonstrates that multipoint measurements along with the global model results enable the reconstruction of a more complete system-level picture of the dropout events and provides insight into M–I coupling aspects that have not previously been investigated.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Statistical and temporal characteristics of sawtooth events

Magnetospheric sawtooth events are characterized by periodic particle injections and magnetic dipolarizations spread quasi-simultaneously across a wide range of magnetic local times. We present a comprehensive statistical study of magnetospheric sawtooth events (STEs) during solar cycle 24 (2008–2016), extending previous catalogs and enabling solar cycle comparisons. Our results confirm that STEs predominantly occur during the rising and declining phases of the solar cycle, and are strongly associated with geomagnetic storms. Superposed epoch analysis reveals near-simultaneous particle injections across all magnetic local time sectors, but magnetic field dipolarization confined to the midnight region. These results support a scenario in which nightside tail reconnection and enhanced convection are the primary drivers of sawtooth oscillations. The localization of magnetic dipolarizations during STEs challenges global instability interpretations and suggests that STEs represent a stormtime substorm mode triggered under specific solar wind and magnetotail conditions. Superposed epoch analyses also show enhanced oxygen content in the magnetosphere during sawtooth events, but do not show a significant difference from geomagnetic storms that do not exhibit periodic behavior.

58 GEOSCIENCES↗
NASA NTRS1960

Observations of the Van Allen Radiation Regions During August and September 1959

An integrating ionization chamber and a single Geiger counter were flown on United States satellite Explorer VI in an elliptical orbit extending to 48,000 km. In addition to the Van Allen inner zone and the great outer zone, a stable and distinct intermediate zone was detected throughout August and September 1959. The outer-zone intensity showed a large decrease following the sudden commencement of a geomagnetic storm. Later in the storm the outer zone increased to much in excess of its prestorm level. During stable periods the outer zone was fairly constant and less intense than it had been observed to be with Pioneer III or Pioneer IV or the first Soviet cosmic rocket. Cosmic-ray background counting rates were reached on most passes in August and Septmeber near apogee of the satellite. The radiation ?dumped? from the outer zone during the geomagnetic storm fits very well with the intensity and latitude distribution required to account for balloon observations of auroral X rays made during the IGY period. This paper is based on preliminary analysis of Explorer VI data.

Arnoldy, R. L.↗