Tracking Magnetic Topological Change in a Time-Dependent Coronal Model

Cooper Downs; Emily I. Mason; Jon A. Linker; Roberto Lionello; Viacheslav Titov

arxiv: 2604.21639 · v1 · submitted 2026-04-23 · 🌌 astro-ph.SR

Tracking Magnetic Topological Change in a Time-Dependent Coronal Model

Emily I. Mason , Cooper Downs , Roberto Lionello , Jon A. Linker , Viacheslav Titov This is my paper

Pith reviewed 2026-05-08 14:09 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords magnetic topologyinterchange reconnectionopen fluxsolar coronaMHD simulationtime-dependent modelsquashing factorsolar wind

0 comments

The pith

A time-dependent coronal model indicates that interchange reconnection processes 3.5% of the open magnetic flux each day.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper applies a mapping technique to follow how magnetic field lines in the Sun's corona change over time in a simulation. By tracking individual magnetic elements instead of static connections, it measures both the open flux and the portion that undergoes interchange reconnection. The simulation, representing solar minimum conditions, finds that this reconnection handles a median of 3.5% of the total open flux in any 24-hour period. Such a finding matters because it quantifies a process thought to contribute to the slow solar wind and the overall evolution of the heliospheric magnetic field. The approach also allows visualization of these changes in a simple 2D plot.

Core claim

The slip-back mapping applied to the thermodynamic MHD simulation tracks topological changes by following advecting magnetic elements across the map. This recovers the flux processed through interchange reconnection, with involved field lines following high squashing factor lines. In the model scaled to solar minimum-like activity, a median value of 3.5% of the total open flux is processed through interchange reconnection in any given 24-hour interval, representing a high proportion of the total open flux changes.

What carries the argument

The slip-back mapping method, which follows individual magnetic elements as they move to determine changes in open and closed field regions over time.

If this is right

The method effectively tracks open flux and recovers interchange reconnection flux in time-dependent models.
Field lines in these processes align with high squashing factor lines, consistent with theories of slow solar wind acceleration.
A median 3.5% of open flux is processed daily, indicating significant topological evolution over time.
Temporal changes can be communicated via an intuitive 2D plot that reduces visual complexity.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

This daily processing rate suggests interchange reconnection could be a major contributor to maintaining the open flux balance in the heliosphere.
Applying the method to simulations with varying activity levels could reveal how the rate changes with the solar cycle.
Correlating these modeled events with spacecraft measurements of solar wind composition might test the link to slow wind streams.

Load-bearing premise

The slip-back mapping applied to the thermodynamic MHD simulation accurately captures physical interchange reconnection without significant contamination from numerical diffusion or grid-scale artifacts in the time-dependent run.

What would settle it

A measurement from solar magnetograms showing a substantially different daily percentage of open flux undergoing reconnection, or a higher-resolution version of the same simulation producing a different median value, would settle the claim.

Figures

Figures reproduced from arXiv: 2604.21639 by Cooper Downs, Emily I. Mason, Jon A. Linker, Roberto Lionello, Viacheslav Titov.

**Figure 1.** Figure 1: Diagram showing the algorithm for both the primary (a) and dual (b) mapping methods detailed here. Each step of the mapping process is outlined in the figure, with steps 1, 3, and 5 being field line traces and 2 and 4 being advections back or forward in time between t0 and t1. The footpoints whose identity determine the SBM code are marked with P (primary) or D (dual) and a number, and are shown in the sam… view at source ↗

**Figure 2.** Figure 2: (a) Example plot of an interchange reconnection code, with labeled footpoints. P2 and D2 are located in the past (t0), while P1, P3, and D3 are in the present (t1), but all are co-plotted for context. Please note that some footpoints are very close to each other, making it difficult to distinguish between them. (b) and (c) illustrate interchange reconnection in an MHD simulation, where the open red and clo… view at source ↗

**Figure 3.** Figure 3: a shows the evolution of these categories of SBM fluxes over time, using a 6-hr cadence for the mapping on the full 721-hour simulation. The FO and FC categories show much higher overall fluxes and more variability than the DC category, as expected from a SBM applied to the inner slip surface. Overall, the FO slightly exceeds the FC category, as can be seen more clearly in Figure 3b view at source ↗

**Figure 4.** Figure 4: Top: plot showing the hourly change in open flux evolution for the three categories of SBM events (red), as well as the standard open flux calculation for the simulation (black). The unsigned difference between these two quantities is plotted as the unsigned interchange reconnection in blue. The orange line corresponds to the ratio of the hourly unsigned IR to the total open flux, with associated units of … view at source ↗

**Figure 5.** Figure 5: Overlay of crossing points of every field line involved in slip-back mapping associated with flux opening, closing, or disconnecting. The bottom layer of the image is a signed log(Q) mapping of the magnetic field at the height listed in the title of the panel. These images show how the field lines that undergo topological change are entirely isolated to high-Q regions, as predicted by theory. An animation … view at source ↗

read the original abstract

We apply the slip-back mapping method of Titov et al. 2009 and Lionello et al. 2020 to a thermodynamic MHD simulation to track topological changes in the magnetic field at a range of temporal cadences. The method constitutes the logical successor to a simple open-field map for a steady-state model, as it tracks changes in the open and closed fields for a time-dependent model by tracking individual magnetic elements as they advect across the map, rather than simply tracing field line connectivity from each cell. Through careful categorization of the slip-back mapping values and analysis of the flux changes, we not only effectively track the open flux but can recover the flux processed through interchange reconnection as well. The field lines involved in these processes are shown to follow lines of high squashing factor, as proposed by interchange reconnection-driven slow solar wind theory. The time-dependent model, which is scaled to solar minimum-like activity, projects that a median value of 3.5% of the total open flux in any given 24-hour interval has been processed through interchange reconnection. This corresponds to a relatively high proportion of the total open flux changes over time in the heliosphere. Our results show that not only is this method a useful tool for accurately tracking topological change in time-dependent simulations, but that its inherent complexity can be visually reduced into an intuitive 2D plot that simply and effectively communicates temporal changes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper extends slip-back mapping to a time-dependent coronal MHD run and reports a median 3.5% of open flux processed by interchange reconnection, but the number may include numerical diffusion effects.

read the letter

The core result here is a concrete 3.5% median for open flux that has undergone interchange reconnection over 24-hour intervals in their solar-minimum thermodynamic MHD simulation. They achieve this by applying the existing slip-back mapping technique to track how magnetic elements move across the open-closed boundary as the field evolves, then categorizing the changes to isolate reconnection-processed flux and displaying the results in simplified 2D plots. This is a straightforward extension of the method from steady-state models, and the plots do make the temporal evolution easier to follow than raw connectivity maps would. The link they draw to high-squashing-factor field lines also lines up with existing interchange-reconnection ideas for slow solar wind. That part is useful for people who already run or analyze coronal MHD simulations and want a practical way to quantify topology changes over time. The main soft spot is the lack of any reported checks on whether numerical diffusion is inflating the 3.5% figure. In a finite-resolution MHD run, field lines can slip across separatrices through grid-scale resistivity or diffusion, and the mapping will register those changes exactly as it does physical reconnection. Without resolution-convergence tests, an ideal-MHD comparison run, or a scan of explicit diffusivity, it is difficult to bound how much of the reported flux is artifact. The abstract workflow is logical, but that missing step directly affects in the headline percentage. This work is aimed at modelers who study coronal magnetic topology and solar-wind source regions. Readers already familiar with the Titov-Lionello slip-back papers will see the incremental advance clearly and may pick up the visualization trick for their own runs. It is solid enough on its own terms to merit a serious referee, mainly because it supplies a specific number from a time-dependent case rather than another steady-state snapshot. I would recommend sending it to peer review with a request that the authors add at least a brief resolution or diffusivity test to support the flux accounting.

Referee Report

1 major / 2 minor

Summary. The manuscript applies the slip-back mapping technique (Titov et al. 2009; Lionello et al. 2020) to a thermodynamic MHD simulation of the solar corona scaled to solar-minimum conditions. By tracking individual magnetic elements across time-dependent open/closed maps and categorizing the resulting slip-back values, the authors quantify both open-flux evolution and the subset of flux that has undergone interchange reconnection. They report a median value of 3.5% of the total open flux processed through interchange reconnection in any given 24-hour interval, show that the relevant field lines align with high-squashing-factor regions, and present a reduced 2D visualization of the temporal topology changes.

Significance. If the mapping procedure cleanly isolates physical interchange reconnection, the 3.5% median supplies a concrete, falsifiable benchmark for the contribution of interchange reconnection to open-flux variability and slow-wind source regions. The work extends steady-state open-field mapping to fully time-dependent models without introducing new free parameters and demonstrates a practical reduction of complex 3-D topology diagnostics to intuitive 2-D plots. These strengths would make the paper a useful reference for both observational and modeling studies of coronal magnetic evolution.

major comments (1)

[analysis of flux changes and categorization of slip-back mapping values] The central 3.5% median result is obtained by categorizing slip-back mapping values to isolate flux that has undergone interchange reconnection. The manuscript does not report resolution-convergence tests, ideal-MHD comparison runs, or explicit resistivity scans. Because the underlying thermodynamic MHD simulation admits numerical diffusion at grid scales that produces connectivity changes indistinguishable from physical reconnection in the mapping, any unquantified numerical contribution directly scales the reported percentage. This issue is load-bearing for the headline quantitative claim.

minor comments (2)

[Abstract] The abstract states that the 3.5% value 'corresponds to a relatively high proportion of the total open flux changes'; a brief quantitative comparison to the total open-flux variation rate or to steady-state expectations would make this statement precise.
[section describing the intuitive 2D plot] The description of the 2-D visualization would benefit from an explicit statement of which mapping quantities are plotted on each axis and how the color scale encodes the categorized reconnection flux.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive evaluation of the work's significance and for identifying the need to clarify the role of numerical effects in the reported interchange reconnection rate. We respond to the major comment below and will make corresponding revisions to the manuscript.

read point-by-point responses

Referee: The central 3.5% median result is obtained by categorizing slip-back mapping values to isolate flux that has undergone interchange reconnection. The manuscript does not report resolution-convergence tests, ideal-MHD comparison runs, or explicit resistivity scans. Because the underlying thermodynamic MHD simulation admits numerical diffusion at grid scales that produces connectivity changes indistinguishable from physical reconnection in the mapping, any unquantified numerical contribution directly scales the reported percentage. This issue is load-bearing for the headline quantitative claim.

Authors: We agree that the absence of explicit resolution-convergence tests, ideal-MHD comparisons, and resistivity scans is a limitation that affects interpretation of the 3.5% median value. The thermodynamic MHD simulation employs a finite-volume scheme with inherent numerical diffusion at grid scales, and the slip-back mapping cannot distinguish these from physical reconnection. Our categorization of slip-back values was chosen to select only those changes aligned with high-squashing-factor regions, consistent with theoretical expectations for interchange reconnection, but this does not fully eliminate the numerical contribution. In revision we will add a new paragraph to the Methods section describing the numerical resistivity and effective magnetic Reynolds number of the model, and we will insert a caveats subsection in the Results noting that the reported percentage represents an upper bound that includes any numerical effects. These additions will constitute a partial revision; we do not plan to perform new simulations but will provide this contextual information to allow readers to assess the claim. revision: partial

Circularity Check

0 steps flagged

No circularity: 3.5% flux result computed from applying prior mapping method to independent MHD simulation

full rationale

The paper derives its headline median value of 3.5% open flux processed through interchange reconnection by applying the slip-back mapping procedure (cited from Titov et al. 2009 and Lionello et al. 2020) to output from a new thermodynamic MHD simulation run. This produces a post-processed statistic on tracked magnetic-element advection and flux-change categorization, which is not equivalent to any input parameter, fitted quantity, or self-referential definition within the present work. Although author overlap exists with the cited method papers, the method functions as an external tool whose validity is presupposed rather than re-derived here; the simulation itself supplies the independent data on which the percentage is evaluated. No self-definitional reductions, fitted-input predictions, ansatz smuggling, or renaming of known results occur in the derivation chain. The result remains falsifiable against the simulation's time-dependent field evolution and external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard MHD equations and the validity of the slip-back mapping algorithm from earlier papers; no new free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (2)

standard math The thermodynamic MHD equations govern the evolution of the coronal magnetic field and plasma.
Invoked implicitly as the basis for the simulation to which the mapping is applied.
domain assumption Slip-back mapping correctly identifies interchange reconnection events when applied to time-dependent output.
This is the central methodological premise carried over from the cited 2009 and 2020 papers.

pith-pipeline@v0.9.0 · 5564 in / 1419 out tokens · 34965 ms · 2026-05-08T14:09:14.618395+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

10 extracted references · 10 canonical work pages

[1]

Otherwise, the connectivity is recorded as closed (C, if both footpoints of the field line are atR 0) or disconnected (D, if both footpoints are atR 1)

Att=t 1, we trace a magnetic field line from the initial point P1 located on one of the slip surfaces and determine the connectivity: if the field line reaches the other slip surface, the connectivity is open (O). Otherwise, the connectivity is recorded as closed (C, if both footpoints of the field line are atR 0) or disconnected (D, if both footpoints are atR 1)

work page
[2]

We label the final point P2

The end-point of the field line is then advected back in time fromt 1 tot 0 with the˜ v(t) flow. We label the final point P2

work page
[3]

Att=t 0, we trace a magnetic field line from P2 and record the connectivity (O, C, or D)

work page
[4]

We label the final point P3

We advect the end-point of the field line forward in time fromt 0 tot 1 using˜ v(t). We label the final point P3

work page
[5]

Notice that P2 and P3 may be a little above or below but still close to the associated slip surface

Att=t 1, we trace a magnetic field line from P3 and record the connectivity (O, C, or D). Notice that P2 and P3 may be a little above or below but still close to the associated slip surface. Analogously, we obtain the dual slip-back mapping (Fig. 1b) through the following steps 14

work page
[6]

Att=t 1, we trace a magnetic field line from the initial point P1 and record the connectivity (O, C, or D) as in the primary mapping

work page
[7]

We label the final point D2

P1 is then advected back in time fromt 1 tot 0 with the˜ v(t) flow. We label the final point D2

work page
[8]

Att=t 0, we trace a magnetic field line from D2 and record the connectivity (O, C, or D)

work page
[9]

We label the final D3

We advect the end-point of the field line forward in time fromt 0 tot 1 using˜ v(t). We label the final D3

work page
[10]

Again, D2 and D3 are generally close but removed from the associated slip surface

Att=t 1, we trace a magnetic field line from D3 and record the connectivity (O, C, or D). Again, D2 and D3 are generally close but removed from the associated slip surface. We now have all the elements to obtain the codes described in§3

work page

[1] [1]

Otherwise, the connectivity is recorded as closed (C, if both footpoints of the field line are atR 0) or disconnected (D, if both footpoints are atR 1)

Att=t 1, we trace a magnetic field line from the initial point P1 located on one of the slip surfaces and determine the connectivity: if the field line reaches the other slip surface, the connectivity is open (O). Otherwise, the connectivity is recorded as closed (C, if both footpoints of the field line are atR 0) or disconnected (D, if both footpoints are atR 1)

work page

[2] [2]

We label the final point P2

The end-point of the field line is then advected back in time fromt 1 tot 0 with the˜ v(t) flow. We label the final point P2

work page

[3] [3]

Att=t 0, we trace a magnetic field line from P2 and record the connectivity (O, C, or D)

work page

[4] [4]

We label the final point P3

We advect the end-point of the field line forward in time fromt 0 tot 1 using˜ v(t). We label the final point P3

work page

[5] [5]

Notice that P2 and P3 may be a little above or below but still close to the associated slip surface

Att=t 1, we trace a magnetic field line from P3 and record the connectivity (O, C, or D). Notice that P2 and P3 may be a little above or below but still close to the associated slip surface. Analogously, we obtain the dual slip-back mapping (Fig. 1b) through the following steps 14

work page

[6] [6]

Att=t 1, we trace a magnetic field line from the initial point P1 and record the connectivity (O, C, or D) as in the primary mapping

work page

[7] [7]

We label the final point D2

P1 is then advected back in time fromt 1 tot 0 with the˜ v(t) flow. We label the final point D2

work page

[8] [8]

Att=t 0, we trace a magnetic field line from D2 and record the connectivity (O, C, or D)

work page

[9] [9]

We label the final D3

We advect the end-point of the field line forward in time fromt 0 tot 1 using˜ v(t). We label the final D3

work page

[10] [10]

Again, D2 and D3 are generally close but removed from the associated slip surface

Att=t 1, we trace a magnetic field line from D3 and record the connectivity (O, C, or D). Again, D2 and D3 are generally close but removed from the associated slip surface. We now have all the elements to obtain the codes described in§3

work page