Turbulence in the Local Interstellar Medium and the IBEX Ribbon
Pith reviewed 2026-05-24 19:54 UTC · model grok-4.3
The pith
Turbulence in the VLISM at scales of 100 au or larger produces an IBEX ribbon whose large-scale structure mismatches observations, while scales of 10 au or smaller yield a smoother ribbon consistent with data.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The inclusion of turbulent fluctuations at scales ≳100 au with power consistent with Voyager 1 observations produces a ribbon whose large-scale structure is inconsistent with IBEX observations, whereas restricting the fluctuations to ∼10 au or smaller produces a smoother ribbon structure similar to IBEX observations. Different turbulence realizations produce different small-scale features ≲10° in the ribbon, but its large-scale structure is robust if the maximum fluctuation size is ∼10 au. The ribbon thickness is considerably larger if the large-scale mean field is draped around the heliosphere, and the magnetic mirror force still plays an important role in trapping particles even though the
What carries the argument
Motion of charged particles under the magnetic mirror force in a turbulent magnetic field superposed on a large-scale mean field (uniform or MHD-derived), with turbulence modeled as a homogeneous random field whose power spectrum is set by Voyager 1 measurements.
If this is right
- The observed IBEX ribbon constrains the dominant turbulent fluctuation scale in the VLISM to be no larger than about 10 au.
- The heliosphere-VLISM interaction, rather than pure interstellar turbulence, sets the magnetic field geometry that shapes the ribbon at scales below 10 au.
- Particle trapping by magnetic mirroring remains effective for ribbon formation even when the ribbon lacks a double-peaked profile.
- MHD-derived draped mean fields around the heliosphere increase ribbon thickness compared with uniform-field assumptions.
Where Pith is reading between the lines
- Ribbon maps could serve as an indirect probe of the spatial transition between heliospherically influenced fields and truly interstellar turbulence.
- Future particle-tracing models may need to incorporate spatially varying turbulence spectra that decay away from the heliopause rather than assuming homogeneity throughout the VLISM.
- High-resolution IBEX or future mission data on ribbon small-scale features could distinguish between different turbulence realizations once the maximum scale is fixed at 10 au.
Load-bearing premise
The turbulent magnetic field component can be treated as a homogeneous random field whose power spectrum and amplitude at the tested scales are correctly given by Voyager 1 measurements.
What would settle it
A direct measurement in the VLISM of magnetic fluctuation power at scales between 10 and 100 au that is substantially lower than the Voyager 1 spectrum used in the simulations, combined with a ribbon map that still shows the large-scale inconsistencies predicted for the high-power case.
Figures
read the original abstract
The effects of turbulence in the very local interstellar medium (VLISM) have been proposed by Giacalone & Jokipii (2015) to be important in determining the structure of the Interstellar Boundary Explorer (IBEX) ribbon via particle trapping by magnetic mirroring. We further explore this effect by simulating the motion of charged particles in a turbulent magnetic field superposed on a large-scale mean field, which we have considered to be either spatially-uniform or a mean field derived from a 3D MHD simulation. We find that the ribbon is not double-peaked, in contrast to Giacalone & Jokipii (2015). However, the magnetic mirror force still plays an important role in trapping particles. Furthermore, the ribbon$'$s thickness is considerably larger if the large-scale mean field is draped around the heliosphere. Voyager 1 observations in the VLISM show a turbulent field component that is stronger than previously thought, which we test in our simulation. We find that the inclusion of turbulent fluctuations at scales ${\gtrsim}$100 au and power consistent with Voyager 1 observations produces a ribbon whose large-scale structure is inconsistent with IBEX observations. However, restricting the fluctuations to ${\sim}$10 au or smaller produces a smoother ribbon structure similar to IBEX observations. Different turbulence realizations produce different small-scale features ${\lesssim}10{\deg}$ in the ribbon, but its large-scale structure is robust if the maximum fluctuation size is ${\sim}$10 au. This suggests that the magnetic field structure at scales ${\lesssim}$10 au is determined by the heliosphere$'$s interaction with the VLISM and cannot entirely be represented by homogeneous interstellar turbulence.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses particle-tracing simulations to examine how turbulence in the VLISM affects the IBEX ribbon. Turbulent magnetic fluctuations with power spectra and amplitudes taken from Voyager 1 are superposed on either a uniform mean field or an MHD-derived draped field. The authors report that the ribbon is not double-peaked (in contrast to Giacalone & Jokipii 2015), that magnetic mirroring remains important for particle trapping, and that ribbon thickness increases when the mean field is draped. The central result is that including fluctuations at scales ≳100 au produces large-scale ribbon morphology inconsistent with IBEX data, while restricting the maximum fluctuation scale to ∼10 au yields smoother ribbons whose large-scale structure matches observations; different realizations affect only small-scale (<10°) features when the cutoff is ∼10 au. The authors conclude that field structure at ≲10 au scales is set by heliosphere-VLISM interaction and cannot be fully captured by homogeneous interstellar turbulence.
Significance. If the scale-dependent results are robust, the work supplies a concrete constraint on the turbulence scales that can be present in the VLISM while remaining consistent with the observed IBEX ribbon. It also supplies evidence that Voyager 1 spectra at the largest scales may contain heliospheric interaction signatures rather than pure interstellar turbulence. The use of both uniform and MHD mean fields plus multiple turbulence realizations provides a useful check on the robustness of the large-scale morphology.
major comments (2)
- [Abstract and simulation setup] Abstract and simulation-setup paragraph: the claim that power at scales ≳100 au produces ribbon structure inconsistent with IBEX rests on the assumption that the Voyager 1 spectrum at those scales represents a homogeneous random field of interstellar origin. The manuscript provides no test or discussion of whether those large-scale fluctuations could instead contain draping, compression, or other heliosphere-induced signatures, which would invalidate the homogeneous-superposition runs used to demonstrate inconsistency.
- [Results on ribbon morphology] Results section on ribbon morphology: the statement that restricting fluctuations to ∼10 au produces a ribbon “similar to IBEX observations” is presented without quantitative metrics (e.g., angular width, intensity contrast, or goodness-of-fit measures) or convergence tests with respect to the number of particles or grid resolution, making it difficult to assess how strongly the data support the scale cutoff.
minor comments (2)
- [Abstract] The abstract contains a typographical artifact (“ribbon$'$s”).
- Notation for the turbulence cutoff scale is given as both “∼10 au or smaller” and “maximum fluctuation size is ∼10 au”; a single, explicit definition would improve clarity.
Simulated Author's Rebuttal
We thank the referee for their constructive and insightful comments. We respond to each major comment below and will incorporate revisions to address the points raised.
read point-by-point responses
-
Referee: [Abstract and simulation setup] Abstract and simulation-setup paragraph: the claim that power at scales ≳100 au produces ribbon structure inconsistent with IBEX rests on the assumption that the Voyager 1 spectrum at those scales represents a homogeneous random field of interstellar origin. The manuscript provides no test or discussion of whether those large-scale fluctuations could instead contain draping, compression, or other heliosphere-induced signatures, which would invalidate the homogeneous-superposition runs used to demonstrate inconsistency.
Authors: We agree that the assumption of homogeneous turbulence is central and merits explicit discussion. Our results show that superposing homogeneous fluctuations at ≳100 au (with Voyager amplitudes) produces large-scale ribbon morphology inconsistent with IBEX, while smaller scales do not. This outcome itself indicates that the VLISM field at those scales cannot be purely homogeneous interstellar turbulence and must incorporate heliospheric interaction effects. We will revise the abstract and add a paragraph in the discussion section to clarify the modeling assumptions, note the possible heliospheric contributions to Voyager spectra at large scales, and discuss how this supports our conclusion that ≲10 au structure is set by the heliosphere-VLISM interaction. revision: yes
-
Referee: [Results on ribbon morphology] Results section on ribbon morphology: the statement that restricting fluctuations to ∼10 au produces a ribbon “similar to IBEX observations” is presented without quantitative metrics (e.g., angular width, intensity contrast, or goodness-of-fit measures) or convergence tests with respect to the number of particles or grid resolution, making it difficult to assess how strongly the data support the scale cutoff.
Authors: We acknowledge that the current manuscript relies on qualitative visual comparison. In the revised version we will add quantitative metrics, including the measured angular width (FWHM) of the ribbon and the peak-to-background intensity contrast, for the uniform-field and draped-field cases with different turbulence cutoffs. We will also report results from a convergence test varying the number of traced particles (e.g., 10^5 to 10^6) to confirm that the reported large-scale morphology is insensitive to particle count above our nominal value. revision: yes
Circularity Check
No significant circularity; forward simulations from external inputs
full rationale
The paper performs particle-tracing simulations that take the turbulent power spectrum and amplitude as direct inputs from Voyager 1 VLISM observations, superposed on either uniform or independently-derived MHD mean fields. Ribbon morphology for different fluctuation scale cutoffs is then compared to IBEX data. No step reduces a claimed prediction to a fitted parameter by construction, nor does any load-bearing premise collapse to a self-citation chain; the Giacalone & Jokipii (2015) reference is cited only for motivation and is explicitly contrasted with the new results. The derivation remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Charged-particle trajectories obey the Lorentz force including magnetic mirror term in a superposition of mean plus turbulent field.
- domain assumption Turbulent fluctuations at the tested scales can be represented as a statistically homogeneous random field whose power matches Voyager 1 observations.
Reference graph
Works this paper leans on
-
[1]
Batchelor, G. K. 1960, The Theory of Homogeneous Turbulence (2 nd ed.; Cambridge: Cambridge Univ. Press)
work page 1960
-
[2]
Burlaga, L. F., Florinski, V., & Ness, N. F. 2015, ApJL, 804, L31
work page 2015
-
[3]
Burlaga, L. F., Florinski, V., & Ness, N. F. 2018, ApJ, 854, 20
work page 2018
-
[4]
A., Czechowski, A., & Grygorczuk, J
Bzowski, M., Kubiak, M. A., Czechowski, A., & Grygorczuk, J. 2017, ApJ, 845, 15
work page 2017
-
[5]
Bzowski, M., Möbius, E., Tarnopolski, S., Izmodenov, V., & Gloeckler, G. 2009, 143, 177
work page 2009
-
[6]
Bzowski, M., Sokół, J. M., Tokumaru, M., et al. 2013, in Cross-Calibration of Far UV Spectra of Solar System Objects and the Heliosphere, ed. E. Quémerais, M. Snow, & R.-M. Bonnet (New York: Springer), 67
work page 2013
- [7]
- [8]
-
[9]
Florinski, V., Heerikhuisen, J., Niemiec, J., & Ernst, A. 2016, ApJ, 826, 197
work page 2016
-
[10]
P., Heerikhuisen, J., Hu, Q., & Khazanov, I
Florinski, V., Zank, G. P., Heerikhuisen, J., Hu, Q., & Khazanov, I. 2010, ApJ, 719, 1097
work page 2010
-
[11]
Fraternale, F., Pogorelov, N. V., Richardson, J. D., & Tordella, D. 2019, ApJ, 872, 40
work page 2019
-
[12]
O., Allegrini, F., Bochsler, P., et al
Funsten, H. O., Allegrini, F., Bochsler, P., et al. 2009, SSRv, 146, 75
work page 2009
-
[13]
A., Bochsler, P., Chornay, D., et al
Fuselier, S. A., Bochsler, P., Chornay, D., et al. 2009, SSRv, 146, 117
work page 2009
-
[14]
V., Heerikhuisen, J., & Rassoul, H
Gamayunov, K. V., Heerikhuisen, J., & Rassoul, H. 2017, ApJ, 845, 63
work page 2017
-
[15]
V., Heerikhuisen, J., & Rassoul, H
Gamayunov, K. V., Heerikhuisen, J., & Rassoul, H. 2019, ApJL, 876, L21
work page 2019
- [16]
-
[17]
Giacalone, J., & Jokipii, J. R. 1999, ApJ, 520, 204 18
work page 1999
-
[18]
Giacalone, J., & Jokipii, J. R. 2015, ApJL, 812, L9
work page 2015
- [19]
-
[20]
Grygorczuk, J., Ratkiewicz, R., Strumik, M., & Grzedzielski, S. 2011, ApJL, 727, L48
work page 2011
-
[21]
Gurnett, D. A., Kurth, W. S., Burlaga, L. F., Ness, N. F. 2013, Sci, 341, 1489
work page 2013
-
[22]
Heerikhuisen, J., & Pogorelov, N. V. 2010, in ASP Conf. Ser. 429, Numerical Modeling of Space Plasma Flows, ASTRONUM-2009, ed. N. V. Pogorelov, E. Audit, & G. P. Zank (San Francisco, CA: ASP), 227
work page 2010
-
[23]
Heerikhuisen, J., Pogorelov, N. V., & Zank, G. P. 2013, in ASP Conf. Ser. 474, Numerical Modeling of Space Plasma Flows, ASTRONUM-2012, ed. N. V. Pogorelov, E. Audit, & G. P. Zank (San Francisco, CA: ASP), 195
work page 2013
-
[24]
Heerikhuisen, J., & Pogorelov, N. V. 2011, ApJ, 738, 29
work page 2011
-
[25]
Heerikhuisen, J., Pogorelov, N. V., Zank, G. P., et al. 2010, ApJL, 708, L126
work page 2010
-
[26]
Heerikhuisen, J., Zirnstein, E. J., Funsten, H. O., Pogorelov, N. V., & Zank, G. P. 2014, ApJ, 784, 73
work page 2014
-
[27]
Isenberg, P. A. 2014, ApJ, 787, 76
work page 2014
-
[28]
Isenberg, P. A. 2015, JPCS, 577, 012014
work page 2015
- [29]
-
[30]
A., Swaczyna, P., Bzowski, M., et al
Kubiak, M. A., Swaczyna, P., Bzowski, M., et al. 2016, ApJS, 223, 25
work page 2016
-
[31]
Lallement, R., Quémerais, E., Koutroumpa, D., et al. 2010, in AIP Conf. Ser. 1216, Twelfth International Solar Wind Conference, ed. M. Maksimovic, K. Issautier, N. Meyer-Vernet, M. Moncuquet, & F. Pantellini (Melville, NY: AIP), 555
work page 2010
- [32]
-
[33]
J., Allegrini, F., Bzowski, M., et al
McComas, D. J., Allegrini, F., Bzowski, M., et al. 2014, ApJS, 213, 20
work page 2014
-
[34]
McComas, D. J., Bzowski, M., Fuselier, S. A., et al. 2015, ApJS, 220, 22
work page 2015
-
[35]
McComas, D. J., Zirnstein, E. J., Bzowski, M., et al. 2017, ApJS, 229, 41 Möbius, E., Bochsler, P., Bzowski, M., et al. 2009, Sci, 326, 969 Möbius, E., Liu, K., Funsten, H., Gary, S. P., & Winske, D. 2013, ApJ, 766, 129
work page 2017
-
[36]
Niemiec, J., Florinski, V., Heerikhuisen, J., Nishikawa, K.-I. 2016, ApJ, 826, 198
work page 2016
- [37]
-
[38]
Pogorelov, N. V., Heerikhuisen, J., Zank, G. P., et al. 2011, ApJ, 742, 104
work page 2011
-
[39]
Pogorelov, N. V., Zank, G. P., Borovikov, S. N., et al. 2008, in ASP Conf. Ser. 385, Numerical Modeling of Space Plasma Flows, ed. N. V. Pogorelov, E. Audit, & G. P. Zank (San
work page 2008
-
[40]
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 2002, Numerical Recipes in C: The Art of Scientific Computing (2nd ed.; Cambridge: Cambridge Univ. Press)
work page 2002
- [41]
-
[42]
Schwadron, N. A., Bzowski, M., Crew, G. B., et al. 2009, Sci, 326, 966
work page 2009
- [43]
-
[44]
Schwadron, N. A., Möbius, E., McComas, D. J., et al. 2016, ApJ, 828, 81 19 Sokół, J. M., Bzowski, M., & Tokumaru, M. 2019, ApJ, in press, arXiv:1809.09823v2 Sokół, J. M., Swaczyna, P., Bzowski, M., & Tokumaru, M. 2015, SoPh, 290, 2589
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[45]
Stone, E. C., Cummings, A. C., McDonald, F. B., et al. 2013, Sci, 341, 150
work page 2013
-
[46]
Strumik, M., Ben-Jaffel, L., Ratkiewicz, R., & Grygorczuk, J. 2011, ApJL, 741, L6
work page 2011
-
[47]
Summerlin, E. J., Viñas, A. F., Moore, T. E., Christian, E. R., & Cooper, J. F. 2014, ApJ, 793, 93
work page 2014
-
[48]
Zank, G. P. 1999, SSRv, 89, 413
work page 1999
-
[49]
Zank, G. P. 2015, ARAA, 53, 449
work page 2015
- [50]
-
[51]
J., Heerikhuisen, J., & Dayeh, M
Zirnstein, E. J., Heerikhuisen, J., & Dayeh, M. A. 2018, ApJ, 855, 30
work page 2018
-
[52]
J., Heerikhuisen, J., Funsten, H
Zirnstein, E. J., Heerikhuisen, J., Funsten, H. O., et al. 2016, ApJL, 818, L18
work page 2016
-
[53]
Zirnstein, E. J., Heerikhuisen, J., Zank, G. P., et al. 2017, ApJ, 836, 238
work page 2017
-
[54]
Zirnstein, E. J., McComas, D. J., Schwadron, N. A., et al. 2019, ApJ, 876, 92
work page 2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.