Study of debris-plasma interaction in the Earth’s ionosphere.
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Very-Low-Frequency (VLF) transmitters operate worldwide mostly at frequencies of 10–30 kilohertz for submarine communications. While it has been of intense scientific interest and practical importance to understand whether VLF transmitters can affect the natural environment of charged energetic particles, for decades there remained little direct observational evidence that revealed the effects of these VLF transmitters in geospace. Here we report a radially bifurcated electron belt formation at energies of tens of kiloelectron volts (keV) at altitudes of ~0.8–1.5 Earth radii on timescales over 10 days. Using Fokker-Planck diffusion simulations, we provide quantitative evidence that VLF transmitter emissions that leak from the Earth-ionosphere waveguide are primarily responsible for bifurcating the energetic electron belt, which typically exhibits a single-peak radial structure in near-Earth space. Since energetic electrons pose a potential danger to satellite operations, our findings demonstrate the feasibility of mitigation of natural particle radiation environment.
An existing technique to use broadband lightning waveforms to probe the D-region ionosphere (60–90 km altitude) is shown to be extendable to a global scale using the Earth Networks Total Lightning Network (ENTLN). This paper demonstrates the technique in detail on a region of the Southeastern United States. This demonstration shows that diurnal D-region height variation and smaller time-scale variations on the order of tens of minutes to hours are evident in the measurement. The technique is then extended to three additional global regions on this same day: Northeastern U.S., India, and Japan. The diurnal behavior between these different regions is compared to a D-region model from the International Reference Ionosphere.
Within the fully integrated magnetosphere-ionosphere system, many electrodynamic processes interact with each other. We review recent advances in understanding three major meso-scale coupling processes within the system: the transient field-aligned currents (FACs), mid-latitude plasma convection, and auroral particle precipitation. (1) Transient FACs arise due to disturbances from either dayside or nightside magnetosphere. As the interplanetary shocks suddenly compress the dayside magnetosphere, short-lived FACs are induced at high latitudes with their polarity successively changing. Magnetotail dynamics, such as substorm injections, can also disturb the current structures, leading to the formation of substorm current wedges and ring current disruption. (2) The mid-latitude plasma convection is closely associated with electric fields in the system. Recent studies have unraveled some important features and mechanisms of subauroral fast flows. (3) Charged particles, while drifting around the Earth, often experience precipitating loss down to the upper atmosphere, enhancing the auroral conductivity. Recent studies have been devoted to developing more self-consistent geospace circulation models by including a better representation of the auroral conductance. It is expected that including these new advances in geospace circulation models could promisingly strengthen their forecasting capability in space weather applications. The remaining challenges especially in the global modeling of the circulation system are also discussed.
Density irregularities play a vital role in determining how radio frequency (RF) waves travel through plasmas. In the Earth’s ionosphere, these density irregularities also impact radio communication. In this study, we conduct a detailed numerical analysis of RF wave propagation in small-scale ionospheric density irregularities using the advanced Petra-M code. We focus specifically on high-frequency (HF) waves, ranging from 3 to 30 MHz, which are essential for military, amateur radio operators, and emergency communications. By introducing density structures, such as equatorial plasma bubbles derived from fluid simulations, we demonstrate that HF waves can scatter in multiple directions when they encounter these irregularities. Additionally, we observe significant mode conversion, where incoming electromagnetic waves transform into electrostatic modes within the density gradient layer. This shows that smaller density irregularities can greatly weaken signals or cause complete signal loss for receivers, emphasizing the need for increased awareness and innovative solutions in radio communication transmission.
Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.
TRAPPIST-1e, an Earth-sized exoplanet in the habitable zone of the nearby M dwarf TRAPPIST-1, may experience magnetospheric responses that vary with stellar space weather, which could potentially influence both its habitability and radio emissions. Our objective is to investigate how different Earth-like magnetospheric configurations of TRAPPIST-1e – specifically variations in dipolar magnetic field strength and axial tilt – respond to diverse stellar space weather conditions, including events analogous to coronal mass ejections (CMEs), and to assess their implications for potential habitability and expected radio emissions. We conducted 3D magnetohydrodynamic simulations of the TRAPPIST-1e system using the PLUTO code in spherical coordinates. The planetary magnetic field was modelled as dipolar, with equatorial strengths from Earth-like to several times stronger. The dipole axis spans a representative range of axial tilts. We investigate four stellar wind environments, from sub-Alfvénic flow to CME-like disturbances. Planetary shielding was quantified based on the magnetopause standoff distance, and radio powers were estimated via empirical scaling laws. Our simulations show that both shielding and radio power depend strongly on the magnetic configuration. Stronger fields increase protection, while larger tilts reduce it. Radio power increases with both field strength and tilt across all wind regimes. An Earth-like magnetic field can provide effective shielding even under intense CMEs, whereas high tilts require stronger fields. Predicted radio powers reach ~10 20 erg s −1 during CMEs, making bursts from close-in, magnetised planets more detectable. However, for TRAPPIST-1e, the maximum cyclotron frequency lies below the Earth’s ionospheric cutoff (~10 MHz), making ground-based detection currently infeasible.
Over the duration of June 22 nd to July 15 th , I interned at Los Alamos National Laboratory within ISR-2 on the Los Alamos Portable Pulser (LAPP) Program. This program concerns a radar dish used for characterizing signal propagation through earth’s ionosphere. The bulk of my internship consisted of characterizing the delay of a Marx Bank Generator and writing hardware-interface C/C++ code for the pulse-recording Digitizer System.
We propose the use of the Earth as a transducer for ultralight dark-matter detection. In particular, we point out a novel signal of kinetically mixed dark-photon dark matter: a monochromatic oscillating magnetic field generated at the surface of the Earth. Similar to the signal in a laboratory experiment in a shielded box (or cavity), this signal arises because the lower atmosphere is a low-conductivity air gap sandwiched between the highly conductive interior of the Earth below and ionosphere or interplanetary medium above. At low masses (frequencies) the signal in a laboratory detector is usually suppressed by the size of the detector multiplied by the dark-matter mass. Crucially, in our case the suppression is by the radius of the Earth, and not by the (much smaller) height of the atmosphere. We compute the size and global vectorial pattern of our magnetic field signal, which enables sensitive searches for this signal using unshielded magnetometers dispersed over the surface of the Earth. In principle, the signal we compute exists for any dark photon in the mass range 10 –21 eV ≲ m A' ≲ 3 × 10 –14 eV . We summarize the results of our companion paper [M. A. Fedderke et al., Search for dark-photon dark matter in the SuperMAG geomagnetic field dataset, arXiv:2108.08852], in which we detail such a search using a publicly available dataset from the SuperMAG Collaboration: we report no robust signal candidates and so place constraints in the (more limited) dark-photon dark-matter mass range 2 × 10 –18 eV ≲ m A' ≲ 7 × 10 –17 eV (corresponding to frequencies 6 × 10 –4 Hz ≲ f ≲ 2 × 10 –2 Hz). These constraints are complementary to existing astrophysical bounds. Future searches for this signal may improve the sensitivity over a wide range of ultralight dark-matter candidates and masses.
In weakly collisional plasma environments with sufficiently low electron beta, Alfvénic turbulence transforms into inertial Alfvénic turbulence at scales below the electron skin depth, k ⊥ d e ≳ 1. We argue that, in inertial Alfvénic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as cos α k ∝ $k$$^{–1}_{⊥}$. The kinetic and magnetic energy spectra scale as ∝ $k$$^{–5/3}_{⊥}$ and ∝ $k$$^{–11/3}_{⊥}$, respectively. As a result of the dual direct cascade, the generalized helicity spectrum scales as ∝ $k$$^{–5/3}_{⊥}$, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the k ⊥ d e >> 1 range. Turbulent eddies exhibit a phase-space anisotropy k ∥ ∝ $k$$^{5/3}_{⊥}$, consistent with critically balanced inertial Alfvén fluctuations. Furthermore, our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth’s magnetosheath and ionosphere, solar corona, and nonrelativistic pair plasmas, as well as to strongly rotating nonionized fluids.
Spacecraft operating in low-Earth orbit are subjected to a number of hazardous environmental constituents that can lead to decreased system performance and reduced operational lifetimes. Due to their thermal, optical, and mechanical properties, polymers are used extensively in space systems; however they are particularly susceptible to material erosion and degradation as a result of exposure to the LEO environment. The focus of this research is to examine the material erosion and mass loss experienced by a custom Kaptonlike polyimide due to exposure in a simulated low-earth orbit environment. The deployable membrane telescope design discussed in this research is named Peregrine and is the scientific payload for FalconSat-7, a 3U cubesat designed and developed at the United States Air Force Academy (USAFA) for the purpose of demonstrating the capability of deployable membrane telescope technology on a nanosatellite platform. In addition to the polymer samples, chrome, silver and gold specimens will be examined to measure the oxidation rate and act as a control specimen, respectively. A magnetically filtered atomic oxygen plasma source has previously been developed and characterized for the purpose of simulating the low-Earth orbit environment. The plasma source can be operated at a variety of discharge currents and gas flow rates, of which the plasma parameters downstream of the source are dependent. The characteristics of the generated plasma were examined as a function of these operating parameters to optimize the production of O + ions with energy relevant to LEO applications. The erosion yield of the Kapton-like polyimide was experimentally determined to be 2.84 x 10 -24 cm 3 per atom which is in close agreement when compared to the on-orbit measurement for reaction efficiency of Kapton HN. We report the experimentally determined reaction rate for the Kapton-like polyimide was used to estimate the operational lifetime of the photon sieve during the solar conditions expected beginning in May 2019. The effective lifetime is estimated between 126 and 153 days.
In weakly-collisional plasma environments with sufficiently low electron beta, Alfvenic turbulence transforms into inertial Alfvenic turbulence at scales below the electron skin-depth. We argue that, in inertial Alfvenic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale-dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as proportional to the perpendicular wavenumber to the power of -1. The kinetic and magnetic energy spectra scale as the perpendicular wavenumber to the power of -5/3 and -11/3, respectively. As a result of the dual direct cascade, the generalized-helicity spectrum scales as the perpendicular wavenumber to the power of -5/3, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the sub-electron-skin-depth range. Turbulent eddies exhibit a phase-space anisotropy, with the parallel wavenumber scaling as the perpendicular wavenumber to the power of 5/3, consistent with critically-balanced inertial Alfven fluctuations. Our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth's magnetosheath and ionosphere, solar corona, non-relativistic pair plasmas, as well as to strongly rotating non-ionized fluids.
An electromagnetic (EM) wave originating at or near the earth’s surface can propagate upwards, traverse the ionosphere, and then be detected at a satellite. While the signal’s magnitude will affect its detectabilty, its polarization will play an important role as well. Furthermore the signal polarization can reveal properties about its source, and the ionospheric conditions through which it traveled. Thus, it is desirable to predict the polarization characteristics of trans-ionospheric EM signals. This situation is complicated by the birefringence imposed on an EM wave propagating in the magnetized plasma of the ionosphere. Each frequency component of the signal can excite two waves that propagate at different phase velocities due to their separate refractive characteristics.For this work, the amplitude and phase of each EM wave field spatial component is followed as it crosses the lower ionospheric boundary, traverses the ionosphere and then crosses the upper boundary. Results are for a single frequency, and can be easily extended to broadband waveforms using fourier transform techniques.
The foundation stone of atmospheric electricity is the existence of a Global Electric Circuit (GEC) influencing all of the planet. In fact, global thunderstorm activity acts as a voltage source that imposes a potential difference between the ionosphere (positively charged) and the Earth’s surface (negatively charged) of about ~300 kV. Such potential difference is discharged through the poorly conducting atmosphere in fair-weather regions, with atmospheric electric fields of ~100 V/m. First evidence of the GEC was gathered by scientists involved on the Carnegie vessel expeditions (early 20th Century), in which they describe a similar daily variation of the atmospheric electric field in different parts of the Pacific Ocean. Such variation become known as the Carnegie Curve (CC).
Geomagnetic substorms are major energy transfer events where energy stored in the Earth's magnetotail is released into the ionosphere. Substorm phenomena, including auroral activities, earthward Poynting flux, magnetic field dipolarization, etc, have been extensively studied. However, the complex interplay among them is not fully understood. In a fortuitous event on June 7, 2013, the twin Van Allen Probes (separated by 0.4 h in local time) observed bursts of earthward Alfvenic Poynting flux in the vicinity of the plasma sheet boundary layer (PSBL). The Poynting flux bursts correlate with enhancements of auroral brightness around the footpoints of both spacecraft. This indicates a temporal and spatial correlation between the auroral brightening and Poynting flux bursts, and that the auroral motion is directly linked to the perpendicular expansion of the Alfven wave. These observations suggest that the Alfvenic Poynting flux is a primary driver for the auroral electron acceleration. Furthermore, around the time of auroral brightening, a dipolarization was seen to propagate more than 4 h in local time during a 20 min period. The azimuthal phase speed of this dipolarization (2 deg/min) is too small to explain the azimuthal motion of the aurora (13.6 deg/min), but the dipolarization could be related to the generation of the Alfvenic Poynting flux through phase mixing at strong density gradients like those in the PSBL.
One of the prominent effects of space weather is the formation of rapid geomagnetic field variations on Earth's surface driven by the magnetosphere-ionosphere system. These geomagnetic disturbances (GMDs) cause geomagnetically induced currents to run through ground conducting systems. In particular, localized GMDs (LGMDs) can be high amplitude and can have an effect on scale sizes less than 100 km, making them hazardous to power grids and difficult to predict. In this study, we examine the ability of the Space Weather Modeling Framework (SWMF) to reproduce LGMDs in the 7 September 2017 event using both existing and new metrics to quantify the success of the model against observation. We show that the high-resolution SWMF can reproduce LGMDs driven by ionospheric sources, but struggles to reproduce LGMDs driven by substorm effects. We calculate the global maxima of the magnetic fluctuations to show instances when the SWMF captures LGMDs at the correct times but not the correct locations. To remedy these shortcomings we suggest model developments that will directly impact the ability of the SWMF to reproduce LGMDs, most importantly updating the ionospheric conductance calculation from empirical to physics-based.
Abstract The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite‐based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low‐Earth orbit. As a result, physics‐based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high‐fidelity approach for physics‐based space weather modeling based on the solution of the Navier–Stokes equations using a high‐order discontinuous Galerkin method, combined with a matrix‐free strategy suitable for high‐performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in nonhydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a nonsegregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics‐based models of the ionosphere‐thermosphere system. Finally, a one‐dimensional radial derivation of the physics‐based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions.