pith. sign in

arxiv: 2604.16098 · v1 · submitted 2026-04-17 · 🌌 astro-ph.SR

Towards a Robust Estimate of the Solar Photospheric Poynting Flux and Helicity Flux

Jiayi Liu , Xudong Sun , Peter W. Schuck , Lars K. S. Daldorff This is my paper

Pith reviewed 2026-05-10 07:33 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords solar photospherePoynting fluxmagnetic helicityDoppler velocityelectric fieldactive regionflux estimationHelmholtz decomposition
0
0 comments X

The pith

Differences in Doppler velocity handling cause major discrepancies in estimates of solar Poynting and helicity fluxes.

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

The paper examines three common methods for calculating the Poynting flux and helicity flux from solar photospheric observations of magnetic fields and velocities. These methods produce substantially different results for the total energy and helicity injected into an active region, sometimes disagreeing even on the sign. By applying the Helmholtz-Hodge decomposition, the analysis reveals that the Doppler velocity contributes to the fluxes primarily through its effect on the non-inductive electric field. The discrepancies are traced to the distinct, ad hoc ways each method incorporates the Doppler and transverse velocity components. Accurate flux estimates matter for understanding the buildup of magnetic energy that powers solar eruptions.

Core claim

On NOAA active region 12673, the PDFI, DAVE4VM, and DAVE4VMwDV methods yield inconsistent accumulated Poynting and helicity fluxes. The Helmholtz-Hodge decomposition demonstrates that Doppler velocities significantly influence the non-inductive electric field contribution to these fluxes, and the ad hoc treatments of velocities in the methods account for the observed differences in values and signs.

What carries the argument

Helmholtz-Hodge decomposition of the observed velocity field to separate contributions to the electric field, highlighting the role of the non-inductive component from Doppler velocity.

If this is right

  • Improved velocity observations are required to better constrain the electric field estimates.
  • Standardizing the treatment of Doppler and transverse velocities across methods would reduce inconsistencies in flux calculations.
  • The choice of method affects conclusions about energy and helicity injection in active regions.
  • Future work should focus on observations that minimize uncertainties in velocity components.

Where Pith is reading between the lines

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

  • These discrepancies could impact models of solar flare and coronal mass ejection triggering that rely on photospheric flux inputs.
  • Applying the decomposition to other active regions might reveal if the Doppler contribution is generally significant.
  • Simulations with known ground-truth velocities could test which method recovers the true fluxes most accurately.

Load-bearing premise

The Helmholtz-Hodge decomposition cleanly separates the inductive and non-inductive electric field contributions from the observed velocities without substantial interference from measurement errors or the particular properties of the active region.

What would settle it

High-precision, simultaneous measurements of all velocity components in a well-observed active region where independent verification of the total energy input is possible, or controlled numerical simulations where the true Poynting flux is known.

Figures

Figures reproduced from arXiv: 2604.16098 by Jiayi Liu, Lars K. S. Daldorff, Peter W. Schuck, Xudong Sun.

Figure 1
Figure 1. Figure 1: Overview of the evolution of the magnetic field and Doppler velocity NOAA AR 12673 from HMI observation on four different times. Left column: Vector magnetic field maps. The background gray map shows the vertical magnetic field Bz. The red (cyan) arrows show the horizontal field vectors in positive- (negative-) Bz regions with magnetic field strength |B| > 250 G. Right column: Doppler velocity maps. Blue (… view at source ↗
Figure 2
Figure 2. Figure 2: The time evolution of the location, magnetic flux and Doppler velocity for NOAA AR 12673. Top: The time evolution of µ (cosine of the heliocentric angle) for the centroid of the AR. Middle: The time evolution of magnetic flux Φ. Red, blue, and black curves refer to positive, absolute negative, and net magnetic fluxes. The shadow region marks the time when the center of the active region is 45◦ away from th… view at source ↗
Figure 3
Figure 3. Figure 3: The time evolution of estimated Poynting flux and helicity flux for NOAA AR 12673 from PDFI (black), DAVE4VM (blue), and DAVE4VMwDV (red). Top: The time evolution of Poynting flux (left) and accumulated Poynting flux (right). Here and after, data gaps are assigned a zero value for integration. Bottom: The time evolution of helicity flux (left) and accumulated helicity flux (right). Uncertainties due to the… view at source ↗
Figure 4
Figure 4. Figure 4: The time evolution of Poynting flux contributed from inductive and non-inductive electric field for NOAA AR 12673. Top: The time evolution of inductive Poynting flux (left) and accumulated inductive Poynting flux (right). Bottom: The time evolution of non-inductive Poynting flux (left) and accumulated non-inductive Poynting flux (right). 1.00 0.75 0.50 0.25 0.00 0.25 0.50 0.75 1.00 H elicit y Flu x [× 1 0 … view at source ↗
Figure 5
Figure 5. Figure 5: Similar to [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The difference of inductive electric field E I between PDFI, DAVE4VM, and DAVE4VMwDV with the background of Doppler velocity vl map. Left: The difference in E I between PDFI and DAVE4VM. Right: The difference in E I between PDFI and DAVE4VMwDV. The black arrows show the difference of horizontal inductive electric field ∆E I h. The background map shows the time derivative of the vertical magnetic field ∂Bz/… view at source ↗
Figure 7
Figure 7. Figure 7: A comparison between observed ∂Bz/∂t and the derived (∇ ×E)·zˆ from DAVE4VM (left) and DAVE4VMwDV (right). From top to bottom, we show the comparison at 2017-09-03 23:54:00 UT and 2017-09-05 23:54:00 UT. The Spearman coefficient (ρ), Pearson coefficient (r), and slope (S) are also shown on the plots. horizontal coordinates z = Z(x, y) and thus the induction equation for the observations on the τ = 1 surfac… view at source ↗
Figure 8
Figure 8. Figure 8: The distributions of r, normalized residual of ∂Bz/∂t (Equation (28)) at September 3, 23:54UT (top) and September 5, 23:54UT (bottom). Left: the distributions of r based on ∂Bz/∂t calculated with five-point stencil. Right: The distributions for r based on the least-square fit of ∂Bz/∂t. A Gaussian distribution N ∼ (0, 1) is plotted as black curve. The best fit t-distribution for DAVE4VM are plotted as red … view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of S I z and S NI z after adopting the internal fitting. Top: The time evolution of S I z (left) and accumulated S I z (right). Bottom: The time evolution of S NI z (left) and accumulated S NI z (right). Black, blue, and red curves refer to the results from PDFI, DAVE4VM internal fitting, and DAVE4VMwDV internal fitting. when the Doppler velocity becomes larger. The worse fit of the induction eq… view at source ↗
Figure 10
Figure 10. Figure 10: PSD for selected variables as a function of spatial scale d at September 3, 23:54 UT (left) and September 6, 23:54 UT (right). The black, blue, and red lines represent the time derivative of the vertical magnetic field from observation, internal fitting of DAVE4VM and DAVE4VMwDV. The green line represents the time derivative of the vertical magnetic field from DAVE4VMwDV calculated from velocity. in Secti… view at source ↗
Figure 11
Figure 11. Figure 11: Comparison of non-inductive Poynting flux distributions. From left to right, we show the S NI z maps from DAVE4VM, DAVE4VMwDV, and PDFI, respectively. The top, middle, and bottom rows are for September 3 23:54 UT, September 4 23:54 UT, and September 7 12:54 UT, respectively. The red contours mark the regions with Doppler velocity |vl| > 0.5 km s−1 . One way to remove the effect of first term is to use the… view at source ↗
Figure 12
Figure 12. Figure 12: The two-dimension (2D) histograms of non-inductive Poynting flux S NI z from different methods. From left to right, we show the 2D histogram of S NI z between DAVE4VMwDV and PDFI, and DAVE4VMwDV and PDFI+DAVE4VMwDV. The top panel shows the comparison on September 3, 23:54UT. The bottom panel shows the comparison on September 4, 23:54UT. The Pearson coefficient (r), Spearman coefficient (ρ), and the slope … view at source ↗
Figure 13
Figure 13. Figure 13: Similar to [PITH_FULL_IMAGE:figures/full_fig_p019_13.png] view at source ↗
read the original abstract

The observed solar photospheric magnetic fields and Doppler velocities are frequently used to quantify the Poynting flux and helicity flux. Multiple methods have been developed for this purpose, but their estimates of the Poynting flux and helicity flux often differ from one another. Here we study the performance of three widely used methods on NOAA active region 12673: "PTD-Doppler-FLCT Ideal" (PDFI), "Differential Affine Velocity Estimator for Vector Magnetograms" (DAVE4VM), and an extension of the latter with Doppler velocity constraint (DAVE4VMwDV). We find that the values of the accumulated energy and helicity differ significantly between the three methods, even in signs. Using the Helmholtz-Hodge decomposition, we show that Doppler velocity can contribute significantly to the Poynting flux and helicity flux through the non-inductive (curl-free) electric field. The different, ad hoc treatments of the Doppler and transverse velocities in three methods are directly responsible for the discrepancies. We discuss the desired future observations that can better constrain these methods.

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.

Referee Report

2 major / 2 minor

Summary. The manuscript compares three methods (PDFI, DAVE4VM, and DAVE4VMwDV) for computing photospheric Poynting and helicity fluxes from vector magnetograms and Doppler velocities observed in NOAA AR 12673. It reports large differences in accumulated energy and helicity (including sign reversals) across the methods and applies the Helmholtz-Hodge decomposition to the velocity fields to show that Doppler velocities contribute substantially to the fluxes through the curl-free (non-inductive) electric-field component. The central claim is that the discrepancies originate directly from the differing, ad-hoc treatments of Doppler versus transverse velocities in the three algorithms.

Significance. If the attribution holds, the work is significant because Poynting and helicity fluxes are central inputs to models of coronal heating and eruptive activity; clarifying why standard methods disagree improves the reliability of these quantities. The explicit use of Helmholtz-Hodge decomposition on real data provides a useful diagnostic that can guide future algorithm development. The analysis is grounded in a well-observed active region and employs existing, widely used velocity estimators rather than introducing new free parameters.

major comments (2)
  1. [Helmholtz-Hodge decomposition analysis (results section)] The attribution of all sign and magnitude discrepancies to velocity-handling differences rests on the Helmholtz-Hodge decomposition isolating a substantial non-inductive contribution from the Doppler velocities. However, the manuscript provides no error propagation through the decomposition, no synthetic-data tests with controlled noise levels (typical photospheric Doppler uncertainties are several hundred m/s), and no assessment of leakage from measurement noise or AR-specific flows into the curl-free component. This leaves the isolation unverified and weakens the claim that velocity-treatment differences are the sole or primary cause.
  2. [Method comparison and accumulated-flux results] The quantitative comparison of accumulated energy and helicity across the three methods does not include a controlled isolation of the effect of each velocity constraint (e.g., the Doppler constraint added in DAVE4VMwDV versus the unconstrained DAVE4VM). Without such a breakdown or sensitivity runs that vary only the Doppler/transverse handling while holding other algorithmic assumptions fixed, it remains unclear how much of the reported discrepancy is truly due to the ad-hoc treatments versus other differences in electric-field reconstruction.
minor comments (2)
  1. [Abstract] The abstract states that the methods 'differ significantly' but does not report the actual numerical values or ratios of the accumulated fluxes, making it harder for readers to gauge the practical size of the discrepancies.
  2. [Figures and methods section] Figure captions and the data-processing description lack explicit statements of the spatial and temporal resolution, the exact time interval used for accumulation, and whether any smoothing or interpolation was applied to the velocity fields before decomposition.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed review of our manuscript. The comments highlight important aspects of the analysis that we will address in the revision. We respond to each major comment below.

read point-by-point responses

  1. Referee: [Helmholtz-Hodge decomposition analysis (results section)] The attribution of all sign and magnitude discrepancies to velocity-handling differences rests on the Helmholtz-Hodge decomposition isolating a substantial non-inductive contribution from the Doppler velocities. However, the manuscript provides no error propagation through the decomposition, no synthetic-data tests with controlled noise levels (typical photospheric Doppler uncertainties are several hundred m/s), and no assessment of leakage from measurement noise or AR-specific flows into the curl-free component. This leaves the isolation unverified and weakens the claim that velocity-treatment differences are the sole or primary cause.

    Authors: We acknowledge that the manuscript does not include formal error propagation through the Helmholtz-Hodge decomposition or synthetic tests with controlled noise levels. The decomposition is applied directly to the velocity fields derived from the three methods on the observed data of NOAA AR 12673 to illustrate that the curl-free (non-inductive) component associated with Doppler velocities contributes substantially to both the Poynting and helicity fluxes. This provides a diagnostic explanation for the sign and magnitude differences seen in the accumulated quantities. We agree, however, that the robustness of this isolation against typical Doppler uncertainties and possible leakage has not been quantified. In the revised manuscript we will add a dedicated discussion of these limitations together with a set of synthetic tests that inject controlled noise levels into the velocity fields before decomposition. revision: yes

  2. Referee: [Method comparison and accumulated-flux results] The quantitative comparison of accumulated energy and helicity across the three methods does not include a controlled isolation of the effect of each velocity constraint (e.g., the Doppler constraint added in DAVE4VMwDV versus the unconstrained DAVE4VM). Without such a breakdown or sensitivity runs that vary only the Doppler/transverse handling while holding other algorithmic assumptions fixed, it remains unclear how much of the reported discrepancy is truly due to the ad-hoc treatments versus other differences in electric-field reconstruction.

    Authors: We note that DAVE4VM and DAVE4VMwDV share the same core differential affine velocity estimator and differ only in the explicit inclusion of the Doppler velocity constraint; their direct comparison therefore isolates the effect of that constraint. PDFI, by contrast, employs an independent Poloidal-Toroidal decomposition that incorporates both Doppler and FLCT transverse velocities under different assumptions. Nevertheless, we agree that an explicit breakdown or additional sensitivity experiments that vary only the Doppler/transverse handling while freezing other algorithmic choices would make the attribution clearer. We will add such a controlled comparison and sensitivity analysis to the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity; analysis is empirical comparison on observational data

full rationale

The paper applies three established velocity-inversion methods (PDFI, DAVE4VM, DAVE4VMwDV) and the standard Helmholtz-Hodge decomposition directly to vector magnetogram and Doppler data from NOAA 12673. Discrepancies in accumulated Poynting and helicity fluxes are shown by explicit numerical evaluation on the same input observations; no parameter is fitted to a subset and then relabeled as a prediction, no equation reduces by construction to its own input, and no load-bearing premise rests on a self-citation chain. The derivation chain is therefore self-contained against external benchmarks and receives the default non-circularity finding.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard solar-physics assumptions about the applicability of the ideal MHD electric field and the validity of the Helmholtz-Hodge decomposition for separating curl-free and divergence-free components from photospheric vector fields.

axioms (1)
  • domain assumption The Helmholtz-Hodge decomposition separates the electric field into inductive (curl) and non-inductive (curl-free) parts that can be computed from observed velocities and magnetic fields.
    Invoked to show that Doppler velocity contributes significantly to the Poynting flux through the non-inductive component.

pith-pipeline@v0.9.0 · 5506 in / 1281 out tokens · 36170 ms · 2026-05-10T07:33:20.621267+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

    work page

  3. [3]

    - [1] #1 = = ^ ^ ^ .\!\!^ d .\!\!^ h .\!\!^ m .\!\!^ s .\!\!^ @mss

    thebibliography [1] 20pt to REFERENCES 6pt =0pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key on reference command E...

    work page doi:10.5281/zenodo.831784 2021