Shocks, compressible perturbations, and intermittency in the very local interstellar medium: Voyager 1 and 2 observations and numerical modeling

This paper combines Voyager 1 and 2 observations with MHD modeling to demonstrate that solar-cycle-driven compressions propagating through the heliopause explain recent magnetic field anomalies in the very local interstellar medium, while turbulence analysis reveals intermittent structures and enables predictions of future magnetic activity through 2030.

The Gist

Imagine the Sun as a giant, glowing campfire in the middle of a vast, cold forest. This fire blows a steady stream of hot air and sparks outward in all directions, creating a protective bubble around the campfire. This bubble is called the Heliosphere.

For decades, the Voyager 1 and 2 spacecraft have been like brave explorers walking away from this campfire, past the edge of the bubble, and into the cold, quiet forest known as the Interstellar Medium.

This paper is a report from a team of scientists who used powerful computer simulations to understand what these explorers are seeing as they walk through that forest. Here is the story, told in simple terms:

1. The Mystery of the "Magnetic Hump"

When Voyager 1 crossed the edge of the Sun's bubble in 2012, it expected to find calm, quiet space. But recently, something strange happened. In 2020, the spacecraft hit a sudden "wall" of pressure (a shockwave). But instead of the magnetic field dropping back to normal afterward, it did something weird:

• It jumped up.
• It stayed high for a long time.
• Then, it formed a giant, slow "hump" that peaked in 2021 and hasn't gone down yet.

It was like walking past a storm, expecting the wind to die down, but instead, the wind kept getting stronger and forming a giant, lingering hill of air that just wouldn't flatten out. Scientists were confused: Why is the magnetic field still so strong? Is Voyager 1 in a completely different part of the universe?

2. The Solution: The Sun's "Breathing"

The scientists built a massive, 3D computer model of the Sun and the space around it. They fed it real data about the Sun's activity over the last 30 years.

They discovered that the Sun doesn't just blow wind; it breathes. Every 11 years, the Sun goes through a "solar cycle," getting very active (lots of sunspots and flares) and then quiet.

• The Analogy: Imagine the Sun is a person blowing bubbles. When they blow hard (solar maximum), they send out a big, fast gust of wind. When they blow gently (solar minimum), the wind is slow.
• The Collision: When the fast wind from the "active" Sun catches up to the slow wind from the "quiet" Sun, they crash into each other. This creates a giant compression wave.

The paper explains that the weird "hump" Voyager 1 saw wasn't a new part of the universe. It was simply the result of the Sun's 11-year breathing cycle (specifically Cycle 24) sending a massive wave of pressure that hit the edge of the Sun's bubble, bounced off, and traveled all the way to Voyager 1.

3. Why Voyager 1 and 2 See Different Things

You might wonder: If it's the same wave, why did Voyager 2 (which is in a different direction) see something different?

• The Terrain Analogy: Imagine the edge of the Sun's bubble (the Heliopause) isn't a smooth balloon. It's more like a crinkled, wavy sheet of fabric.
• Voyager 1 walked into a "valley" where the wave hit hard and piled up, creating that strong magnetic hump.
• Voyager 2 walked into a different spot where the wave was smoother and less intense. It also crossed the edge of the bubble at a slightly different time, so it missed the main "hump" event entirely.

The computer model showed that the Sun's waves hit the bubble at different angles, creating different effects depending on where you are standing.

4. The "Intermittency" Puzzle

Scientists also looked at how "turbulent" the space is. Think of turbulence like the difference between a smooth river and a river full of rapids and rocks.

• The Mystery: For a while, Voyager 1 saw lots of "rapids" (turbulence). Then, around 2022, the water seemed to become perfectly smooth and calm. Some scientists thought this meant Voyager 1 had entered a totally new, pristine region of space.
• The Reality: The paper says, "Not so fast!" The turbulence didn't disappear because Voyager 1 moved to a new place. It disappeared because the Sun's "breathing" is currently in a quiet phase. The waves are getting weaker as they travel further out. The "smoothness" is just the calm after the storm, not a new world.

5. What Happens Next? (The Forecast)

Using their model, the scientists made some predictions for the future:

• Voyager 1: It will likely stay in this "strong magnetic field" zone until about 2030. After that, things will get quieter again.
• Voyager 2: It is in a more active spot. It should hit several more "waves" and pressure fronts before 2026, and then a big one around 2030 (from the next solar cycle).
• New Horizons: This other spacecraft is still inside the Sun's bubble but is getting close to the edge. The model predicts it will cross the boundary (the Termination Shock) around 2031, roughly 80 times the distance from the Earth to the Sun.

The Big Takeaway

The universe isn't as chaotic or random as it seems. The strange things Voyager 1 saw—the strong magnetic fields, the "hump," and the sudden calm—are all connected to the Sun's own rhythm. The Sun is the conductor, and the space around it is the orchestra. Even though Voyager is far away, it is still listening to the music the Sun is playing.

The paper confirms that we don't need to invent new physics to explain these observations; we just need to understand how the Sun's 11-year cycle pushes and pulls on the fabric of space.

Technical Summary

1. Problem Statement

The Voyager 1 (V1) and Voyager 2 (V2) spacecraft have provided unique in situ measurements of the Very Local Interstellar Medium (VLISM) since crossing the heliopause (HP). However, recent observations present significant challenges to existing theoretical models:

• The V1-pf2 Anomaly: In 2020.4, V1 observed a pressure front (pf2) followed by a distinct "magnetic hump" and a persistently strong magnetic field that has not recovered to pre-event levels. Previous models struggled to reproduce the non-decreasing magnetic field strength ($|B|$) immediately behind the shock and the specific "hump" structure.
• Alternative Hypotheses: The inability of standard models to fit these data led to alternative hypotheses, such as V1 penetrating a new, pristine region of the LISM or the event not being of solar origin.
• Intermittency Disappearance: Observations since 2022 suggested a "vanishing" of magnetic turbulence intermittency, which some interpreted as a fundamental change in VLISM properties.
• V2 Discrepancies: V2 observed different magnetic field profiles (smoother, wave-like) and did not detect a strong analog to the V1-pf2 event, despite being in the same general region.

2. Methodology

The authors employed a comprehensive approach combining advanced numerical modeling with detailed data analysis:

Numerical Modeling:

• Framework: Used the Multi-Scale Fluid–Kinetic Simulation Suite (MS-FLUKSS), a 3D, global, data-driven Magnetohydrodynamic (MHD) model.
• Physics: A five-fluid model treating the plasma mixture, pickup ions (PUIs) as a separate comoving fluid, and four populations of neutral hydrogen (including those born in the outer heliosheath).
• Boundary Conditions (BCs): Inner BCs at 1 au were derived from a combination of OMNI, Interplanetary Scintillation (IPS), Ulysses, and Wilcox Solar Observatory (WSO) data. This allowed for high-cadence, time-dependent inputs covering multiple solar cycles.
• Simulation Scope: The study simulated four solar cycles (approx. 1980–2025) to capture self-consistent, long-term global responses. Two primary cases were run with different Interstellar Magnetic Field (ISMF) strengths ($2.93 \, \mu\text{G}$ and $3.50 \, \mu\text{G}$).
• Trajectory Analysis: Virtual spacecraft trajectories were shifted to account for systematic underestimations of the heliosheath thickness in the models, ensuring accurate comparison with Voyager data.

Data Analysis:

• Turbulence Analysis: Applied wavelet transforms (Morlet) and structure-function analysis to Voyager MAG data (48s resolution).
• Intermittency Quantification: Calculated running kurtosis of magnetic field increments to identify non-Gaussian statistics (intermittency) across various time scales (from 288s to 60 days).
• Foreshock Geometry: Analyzed the geometric connection between the spacecraft and shock surfaces to identify potential foreshock regions.

3. Key Contributions

• Unified Solar-Cycle Explanation: The paper provides a self-consistent explanation linking V1-pf2, the magnetic hump, and the sustained high magnetic field to Solar Cycle 24 (SC-24) global compressions, rather than a change in the pristine LISM.
• Reproduction of the "Hump": For the first time, the model successfully reproduces the non-decreasing magnetic field behind a shock and the subsequent "hump" feature. This is attributed to the "pileup" (merging) of multiple compressible wave fronts generated at the heliopause between 2016–2018.
• Clarification of Intermittency: The study demonstrates that the apparent disappearance of intermittency after 2022 is due to turbulence weakening (lower signal-to-noise ratio and declining solar activity) rather than a fundamental change in the medium's nature.
• V1 vs. V2 Differences: The model explains the differences between V1 and V2 observations as a result of latitudinal effects (V2 is in a more compressed southern hemisphere region) and timing (V2 crossed the HP after the SC-24 perturbation was already generated, missing the initial shock front).

4. Key Results

A. The V1-pf2 Event and Magnetic Hump

• Origin: V1-pf2 is identified as a global compression wave resulting from the merger of shocks driven by the increase in solar wind dynamic pressure post-solar maximum (SC-24).
• Magnetic Field Behavior: The model reproduces the observed "hump" and sustained high $|B|$ as a result of wave pileup. Multiple waves generated at the heliopause merged, creating a shallower ramp and preventing the magnetic field from relaxing to pre-shock levels.
• Prediction: The magnetic field at V1 is predicted to remain elevated until $\sim$2030. Models with an ISMF strength of $3.5 \, \mu\text{G}$ provide the best fit to observations.

B. Voyager 2 Observations

• Smooth Profiles: V2's smoother, wave-like profiles are attributed to its location in the southern hemisphere where the outer heliosheath is more compressed, and to the fact that it encountered the tail of the SC-24 perturbation rather than the leading shock front.
• Future Events: V2 is predicted to encounter multiple solar-driven compressions before 2026 and a major SC-25 event around 2030.

C. Turbulence and Intermittency

• Time Dependence: Intermittency is highly time-dependent. Significant fine-scale intermittency (scales $<1$ hour) was observed in two specific intervals: $\sim$2014.5–2015.5 (associated with V1-sh1) and 2018–2021.5.
• Post-2022 Behavior: The near-Gaussian statistics observed after 2022 reflect a weakening of turbulence intensity, not a change in the medium's fundamental properties.
• Foreshock Connection: The study identifies a broad foreshock region for V1-pf2, suggesting V1 may have been magnetically connected to the shock as early as 2017.2. This extended foreshock likely explains the enhanced fine-scale turbulence and energetic particle fluxes observed prior to the shock arrival.

D. New Horizons (NH) Predictions

• The model predicts NH will cross the Termination Shock (TS) around 2031 at a heliocentric distance of $80 \pm 2$ au.
• The simulations show good agreement with Voyager TS crossing times (within 3% in distance).

5. Significance

• Validation of Global Models: The study validates the capability of global, data-driven MHD models to reproduce complex, transient interstellar phenomena driven by solar cycles, resolving long-standing discrepancies between theory and Voyager data.
• Refutation of "New Regime" Hypotheses: It strongly refutes the hypothesis that V1 has entered a pristine, undisturbed LISM region, confirming instead that the VLISM is dynamically coupled to the Sun via global compressions.
• Future Mission Planning: The predictions for V1, V2, and New Horizons provide critical guidance for interpreting future data and planning for the IMAP mission.
• Methodological Advancement: The work highlights the necessity of multi-cycle simulations and uncertainty quantification (ensemble runs) to account for the intrinsic stochasticity of the heliosphere and the sensitivity of shock arrival times to heliopause corrugation.

In conclusion, the paper establishes that the complex magnetic and plasma structures observed in the VLISM are primarily driven by the Sun's 11-year cycle, mediated by global compressions and shock merging, rather than by intrinsic changes in the interstellar medium itself.