An Enhanced "Flux-Corrected Transport"-Based Plasmasphere Refilling Model

The Gist

The Big Picture: Rebuilding a Damaged Ring

Imagine the Earth is surrounded by a giant, invisible, donut-shaped ring made of cold, dense gas (plasma). Scientists call this the plasmasphere. Think of it like a protective bubble of air hugging our planet.

When a massive solar storm hits Earth, it's like a hurricane blowing through this bubble. It strips away most of the gas, leaving the ring thin and empty. Once the storm passes, the Earth needs to "refill" this ring. The gas comes from the layer of the atmosphere just below the ring (the ionosphere) and flows upward along invisible magnetic "straws" (flux tubes) to fill the gap.

The Old Model vs. The New Model

For a long time, scientists used computer models to predict how fast this refilling happens.

• The Old Way: Imagine trying to fill a bathtub, but you assume the water temperature stays exactly the same everywhere, no matter how much you pour in. The old models did this with the gas: they assumed the temperature was constant and unchanging along the magnetic straws.
• The New Way (This Paper): The authors, Jaden Fitzpatrick and colleagues, realized that in reality, the gas gets hotter or cooler depending on where it is and when. They upgraded their model to let the temperature change naturally, just like real water heating up or cooling down as it flows.

The "Two-Stage" Refilling Process

The most exciting discovery from this new, smarter model is that the refilling happens in two distinct stages, like a two-speed car.

1. Stage 1 (The Slow Start): At first, the gas flows up slowly. It's a bit of a trickle.
2. Stage 2 (The Rush): Suddenly, the flow speeds up dramatically, filling the ring much faster.

Why did the old model miss this?
Because the old model assumed the temperature was flat and boring. The new model shows that as the gas moves, it creates temperature differences (gradients). These differences act like a hidden engine. They create an invisible "push" (called an ambipolar electric field) that accelerates the gas, triggering that sudden switch from the slow stage to the fast stage.

The Cast of Characters: Different Types of Gas

The plasmasphere isn't just one type of gas; it's a mix of three main "characters": Hydrogen (H+), Helium (He+), and Oxygen (O+). The new model shows how each one behaves differently:

• Oxygen (The Heavy Hauler): At the very beginning (Stage 1), the heavy Oxygen ions get a big boost from the temperature changes. They rush up early, but then they slow down and can't quite reach the top.
• Hydrogen (The Main Filler): Hydrogen is the lightest and most common. It takes a little longer to get going, but once the second stage starts, it becomes the main worker, filling up the majority of the ring.
• Helium (The Middleman): Helium is the tricky one. The model shows that right when the switch happens between Stage 1 and Stage 2, Helium's presence spikes. It's like a temporary bridge that helps keep the system balanced while the Hydrogen catches up.

Why This Matters

The authors tested their model with different scenarios, like changing the starting amount of gas or simulating different seasons (winter vs. summer). They found that:

• The temperature changes are the secret sauce that makes the "two-stage" behavior happen. Without them, the model just shows a boring, steady flow.
• The model works well whether you look at the ring close to Earth or a bit further out, suggesting it's a solid tool for the future.

The Bottom Line

By letting the temperature change naturally in their computer simulation, the authors created a much more realistic picture of how Earth repairs its protective plasma ring after a solar storm. They proved that heat isn't just a background detail; it's a driver that controls how fast and in what order the different gases fill up the ring. This helps scientists better understand the complex dance of particles that happens in space after a storm.

Technical Summary

Technical Summary: An Enhanced "Flux‐Corrected Transport"‐Based Plasmasphere Refilling Model

Problem Statement
The plasmasphere, a torus of cold, dense plasma surrounding Earth, is frequently eroded by geomagnetic storms. Recovery occurs via "plasmasphere refilling," where ionospheric plasma flows upward along magnetic flux tubes to replenish the depleted region. While previous hydrodynamic models successfully reproduced general refilling dynamics using the Flux-Corrected Transport (FCT) method, they relied on a significant simplification: assuming a spatially uniform and constant electron temperature along the flux tube. This assumption forced the pressure gradient and ambipolar electric field—key drivers of ion transport—to be determined solely by density gradients. Consequently, these models could not fully capture the complex, self-consistent interactions between thermal evolution and multi-ion transport, nor could they adequately explain the variability observed in two-stage refilling behaviors reported in recent satellite data.

Methodology
The authors extended an existing multi-ion, two-stream hydrodynamic FCT model (originally developed by Chatterjee & Schunk, 2019, 2020a) by incorporating a self-consistent electron energy equation. The core modification involves solving the one-dimensional heat conduction equation for electron temperature ($T_e$) as a function of both space and time, rather than holding it constant.

Key methodological components include:

• Electron Energy Equation: The model solves a heat conduction equation (based on Khazanov et al., 1992) that accounts for thermal conduction ($T_e^{5/2}$ dependence) and a constant electron heating rate ($Q_e$). Radiative and inelastic cooling are neglected at plasmaspheric altitudes.
• Modified Ambipolar Electric Field: With $T_e$ now varying spatially and temporally, the ambipolar electric field ($E_{\parallel}$) is recalculated to include both density and temperature gradients: $E_{\parallel} = -\frac{1}{e n_e} \frac{\partial}{\partial s}(n_e k T_e)$. This replaces the density-only formulation used in prior constant-temperature models.
• Coupled System: The system iteratively couples the electron energy equation with the multi-ion continuity and momentum equations for $H^+$, $He^+$, and $O^+$. The ion temperature is assumed equal to the local electron temperature.
• Simulation Configuration: Simulations were conducted for L-shells 3 and 4, extending from $\pm 56^\circ$ magnetic latitude to 1,500 km altitude. Initial conditions included symmetric and asymmetric ion density profiles to simulate standard and seasonal (winter solstice) scenarios.

Key Contributions

1. Self-Consistent Temperature Evolution: The primary contribution is the removal of the constant-temperature assumption, allowing the model to resolve field-aligned temperature gradients that develop dynamically during refilling.
2. Enhanced Physical Representation: By coupling temperature evolution with ion transport, the model provides a more physically complete description of the pressure gradients and ambipolar electric fields that accelerate ions.
3. Investigation of Two-Stage Refilling: The extended framework enables a detailed analysis of the two-stage refilling process (early and late stages), which was previously difficult to distinguish or explain solely through density-driven models.

Results

• Temperature Evolution: A strong field-aligned temperature gradient develops rapidly (within the first hour), with equatorial temperatures differing from high-latitude boundaries by approximately 1,500 K. The system reaches a quasi-steady state in temperature within roughly 60 minutes.
• Refilling Dynamics: The inclusion of variable temperature significantly alters refilling trajectories. Compared to the constant-temperature model, the variable-temperature model predicts a longer total refilling time (extending from ~26 to ~37 hours) and a reduced average refilling rate. These changes bring the model's predictions closer to observational data from LANL and Van Allen Probes.
• Two-Stage Behavior: The model clearly reproduces the two-stage refilling process:
  • Early Stage (0–7 h): Characterized by supersonic, counter-streaming flows from both hemispheres. The model shows enhanced early-time contributions from heavy ions ($O^+$ and $He^+$) due to the stronger ambipolar electric field driven by temperature gradients.
  • Late Stage (>7 h): As densities rise, Coulomb collisions thermalize the streams, and the refilling rate increases significantly until saturation.
• Ion Species Sensitivity:
  • $H^+$: Initial $H^+$ concentration primarily determines the final saturated concentration and the duration of the late stage but has little effect on early-stage dynamics.
  • $He^+$: Initial $He^+$ concentration has a negligible impact on the overall refilling trajectory. However, a peak in the fractional abundance of $He^+$ occurs at the transition between stages, driven by the ambipolar field compensating for $H^+$ loss.
  • $O^+$: Initial $O^+$ concentration significantly influences the late-stage refilling length and final equilibrium density. Asymmetric reductions in $O^+$ (simulating seasonal effects) produce distinct outcomes compared to symmetric reductions, highlighting the role of heavy ions in shaping the recovery trajectory.
• L-Shell Dependence: Comparisons between L=3 and L=4 confirm that refilling stages are compressed in duration and reach higher concentrations at lower L-values, consistent with the shorter flux tube lengths.

Significance and Claims
The paper claims that incorporating spatiotemporal temperature variability is essential for accurately modeling plasmasphere refilling. The authors state that this enhancement:

• Provides a more realistic representation of the acceleration mechanisms for different ion species, particularly the early dominance of heavy ions.
• Offers a physical explanation for the transition between early and late refilling stages, specifically how the ambipolar electric field couples $H^+$ and $He^+$ concentrations.
• Improves the agreement between model-predicted refilling rates and durations with observational data, although the model still predicts faster rates than some long-duration observations (5+ days), suggesting that future inclusion of time-dependent boundary conditions and heating rates is necessary.
• Establishes a robust framework for future extension to three-dimensional geometries and coupling with global ionosphere-thermosphere models to better interpret complex multi-ion transport processes during geomagnetic storm recovery.

The authors maintain a modest stance, acknowledging that the current model still relies on simplified boundary conditions (constant heating rates and fixed boundary temperatures) and a hydrodynamic approach that assumes Maxwellian velocity distributions, which may not hold for very low densities in the early refilling phase. They posit that future work addressing these limitations will further refine the model's predictive capabilities.