🔭 astrophysics

Hour-scale response of low-Earth-orbit satellite altitude to extreme geomagnetic storms revealed by the NinjaSat CubeSat during the 2024 solar maximum

Using data from the NinjaSat CubeSat during the 2024 solar maximum, this study reveals that hour-scale LEO satellite altitude decay during extreme geomagnetic storms is primarily driven by polar Joule heating indexed by the PC index, which leads orbital response by approximately two hours, thereby offering a critical window for drag forecasting.

Original authors: Hayato Tanaka, Takao Nakagawa, Satoshi Nozawa, Toru Tamagawa, Teruaki Enoto, Takao Kitaguchi, Wataru Iwakiri, Tatehiro Mihara, Tomoshi Takeda, Naoyuki Ota, Amira Aoyama, Satoko Iwata, Takuya Takahashi

Published 2026-06-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Hayato Tanaka, Takao Nakagawa, Satoshi Nozawa, Toru Tamagawa, Teruaki Enoto, Takao Kitaguchi, Wataru Iwakiri, Tatehiro Mihara, Tomoshi Takeda, Naoyuki Ota, Amira Aoyama, Satoko Iwata, Takuya Takahashi, Kaede Yamasaki

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the Earth is wrapped in a giant, invisible blanket of air called the thermosphere. This blanket isn't static; it breathes. When the Sun is calm, the blanket is thin and steady. But when the Sun gets angry and sends out a massive burst of energy (a geomagnetic storm), the blanket suddenly puffs up, gets thicker, and becomes "heavier."

For satellites flying in Low Earth Orbit (LEO), this is a problem. A thicker, heavier blanket creates more friction (drag), causing the satellites to slow down and lose altitude faster than expected. If they drop too low, they could crash or burn up.

This paper is like a high-speed detective story. The researchers used a tiny satellite called NinjaSat (a "CubeSat" about the size of a shoebox) to watch how fast this "blanket" puffs up during the two biggest solar storms of 2024. They wanted to solve a mystery: What exactly triggers the blanket to puff up, and how long does it take?

Here is the breakdown of their findings using simple analogies:

1. The Detective Work: Separating the Noise

Satellites naturally lose altitude over time, just like a car slowly losing air from its tires. This is the "background noise." The researchers had to filter out this slow, steady loss to see the sudden, sharp drops caused by the storms.

The Method: They compared the satellite's actual GPS position (what it really did) against a computer prediction of where it should be if only the slow, steady drag existed. The difference between the two told them exactly how much extra "drag" the storm added, hour by hour.

2. The Two Suspects: The "Polar Cap" vs. The "Ring Current"

Scientists have two main "gauges" (indices) they use to measure solar storms. The researchers wanted to know which one actually causes the satellite to drop.

Suspect A: The PC Index (The Polar Cap). Think of this as a thermometer at the North and South Poles. It measures how much solar energy is being dumped directly into the polar atmosphere, heating it up like a blowtorch.
Suspect B: The Dst Index. Think of this as a global mood ring for the Earth's magnetic field. It measures a giant electric current (the ring current) that circles the Earth far above the atmosphere, but it doesn't measure the heat directly.

3. The Big Reveal: Who is the Culprit?

The researchers looked at the timing of the events like a race.

The PC Index (The Heat): When the solar wind hit the poles, the "Polar Cap" gauge spiked. About 2 to 3 hours later, the satellite started dropping rapidly.
- The Analogy: Imagine turning on a heater in a room (the poles). It takes a couple of hours for the warm air to rise and fill the whole room (the satellite's altitude). The heater (PC Index) turned on first, and the room got hot second. This proves the heat from the poles is what pushes the satellite down.
The Dst Index (The Mood Ring): The "Mood Ring" gauge also spiked, but it did so after the satellite had already started dropping.
- The Analogy: The satellite dropped before the mood ring changed color. This means the mood ring isn't the cause of the drop; it's just a side effect that happens later. The researchers concluded that the Dst index is a "lagging indicator"—it tells you a storm is happening, but it's too slow to predict the drag.

4. The "Where" Matters: The May vs. October Storms

The researchers also compared NinjaSat (which flies over the poles) with another satellite, XRISM (which flies near the equator).

May Storm: Both satellites dropped. This means the "heater" was so strong it warmed the air all the way from the poles down to the equator.
October Storm: Only NinjaSat (the polar flyer) dropped. XRISM (the equatorial flyer) didn't budge. This means the "heater" was localized; it only warmed the air near the poles, leaving the equator cool.

5. Why This Matters (According to the Paper)

The paper claims that for the next few years, as the Sun gets more active, we need to stop relying on the slow "Mood Ring" (Dst) to predict satellite crashes. Instead, we should watch the "Polar Cap Thermometer" (PC Index).

Because the PC Index spikes 2–3 hours before the satellite starts to fall, it gives operators a "warning window." It's like seeing the smoke before the fire alarm goes off. If you know the heater is on, you can move the satellite to a safer spot before the blanket gets too thick.

In summary: The paper proves that solar storms heat the air at the poles first, and that heat takes a few hours to travel up and push satellites down. The "Polar Cap" gauge is the best early warning system we have, while the "Ring Current" gauge is just a late notification.

Technical Summary: Hour-scale Response of Low-Earth-Orbit Satellite Altitude to Extreme Geomagnetic Storms

Problem Statement
The rapid growth of the Low-Earth Orbit (LEO) satellite population necessitates accurate, short-latency predictions of thermospheric density to ensure collision avoidance and orbital sustainability. While secular orbital decay is well-understood, the impulsive density enhancements driven by geomagnetic storms remain poorly resolved on hour-scale timescales. Conventional operational models (e.g., NRLMSISE-00, DTM-2013, JB2008) rely on geomagnetic indices with coarse temporal resolution (3-hour Kp/Ap or hourly Dst) or fail to capture the specific causal mechanisms of storm-time heating. A critical unresolved question is which geomagnetic observable serves as the most direct causal proxy for the impulsive heating that drives rapid orbital decay, and what the specific transport timescale is from energy input to atmospheric expansion.

Methodology
The study utilizes high-cadence Global Navigation Satellite System (GNSS) telemetry from two X-ray astronomy satellites operating during the 2024 solar maximum:

NinjaSat: A 6U CubeSat in a sun-synchronous orbit ( $97^\circ$ inclination, $\sim$ 500 km altitude) providing 10-minute cadence position data.
XRISM: An X-ray Imaging and Spectroscopy Mission satellite in a low-inclination orbit ( $31^\circ$ , $\sim$ 560 km altitude) providing 30-second cadence data.

To isolate storm-driven perturbations from the background, the authors constructed an altitude residual ( $\Delta a$ ) defined as the difference between the GNSS-derived altitude ( $a_{GNSS}$ ) and the altitude propagated from NORAD Two-Line Elements (TLEs) using the SGP4 model ( $a_{SGP4}$ ). Since SGP4 captures secular decay but not impulsive density spikes, a positive $\Delta a$ indicates enhanced drag.

The analysis involved:

Baseline Removal: A linear regression of the smoothed altitude change rate ( $\dot{g}\Delta a$ ) against GOES-18 EUV irradiance (140.5 nm) during geomagnetically quiet periods (July–September 2024) to model and subtract the EUV-driven background.
Cross-Correlation: The residual storm-time altitude change rate was cross-correlated with the Polar Cap (PC) index (a proxy for solar-wind energy input to the high-latitude ionosphere) and the Disturbance Storm Time (Dst) index (a proxy for the ring current).
Latitudinal Comparison: By comparing the response of NinjaSat (polar coverage) and XRISM (low-latitude coverage), the study assessed the latitudinal extent of thermospheric heating during two major 2024 storms: the May "Gannon" event and the October event.

Key Results

PC Index as the Primary Driver: The storm-time altitude change rate correlates strongly with the PC index ( $r \gtrsim 0.8$ ). Crucially, the PC index leads the altitude response by approximately 2 to 3 hours (combined inverse-variance weighted mean lag $\tau = -2.2 \pm 2.6$ h; October event specifically $\tau = -3.17 \pm 3.17$ h). This lead time aligns with the physical timescale required for Joule heating deposited at 110–130 km to propagate vertically and meridionally to the $\sim$ 500 km thermosphere.
Dst Index as a Consequence, Not a Cause: While the altitude change rate also correlates strongly with the magnitude of the Dst index ($|Dst|$, $r \approx 0.82–0.84$ ), the temporal ordering is reversed. The altitude drop precedes the Dst minimum by $\sim$ 3 hours. This indicates that the Dst-altitude correlation arises from a common solar-wind driver rather than direct thermal coupling from the ring current, which resides hundreds of thousands of kilometers above the thermosphere.
Latitudinal Dependence:
- May 2024 (Gannon Storm): Both NinjaSat and XRISM exhibited significant orbital decay, indicating that the heating extended to low latitudes ( $<30^\circ$ ).
- October 2024: Only NinjaSat responded significantly; XRISM showed no deviation from its secular trend. This suggests the heating was confined to the polar cap and auroral zone, demonstrating that storm impacts are not uniform and can be polar-localized even during events with comparable Dst minima.
Energetics: The PC-derived hemispheric Joule-heating power ( $Q_J$ ) during the storms exceeded the concurrent EUV heating baseline by factors of 2.5 to 3, establishing geomagnetic Joule heating as the dominant incremental energy source during these events.

Significance and Claims
The paper claims to provide the first quantitative, in-orbit demonstration that high-latitude ionospheric Joule heating is the decisive coupling channel for storm-time LEO orbital decay. The findings challenge the utility of the Dst index as a timely predictor for drag, noting that operational models using Dst or 3-hour Kp/Ap indices systematically underpredict density enhancements because they average out impulsive heating and lag the actual atmospheric response.

The study proposes a concrete pathway for improving operational drag forecasting: replacing Kp, Ap, or Dst proxies with the real-time, 1-minute cadence PC index. Because the PC index leads the orbital response by 2–3 hours, it offers a natural operational window for satellite operators to recompute conjunction assessments and plan corrective maneuvers before the peak density surge arrives. The authors emphasize that while the results are robust for the two analyzed events, a multi-storm analysis across the descending phase of Solar Cycle 25 is required to fully establish the universality of the lag times and amplitude dependencies.