REVIEW 2 major objections 5 minor 31 references
Opposite polarity magnetic field and convective downflows in a simulated sunspot penumbra
T0 review · 2 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Simulated sunspot penumbra contains about twice as much opposite-polarity magnetic field as current telescopes can detect.
desk verdict Useful quantitative catalog of how resolution and noise hide opposite-polarity penumbral fields, but the reference fractions rest on an artificially inclined simulated field that may inflate them. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is a synthetic observation pipeline built on forward modeling. The SPINOR code with STOPRO computes LTE Stokes profiles of the Fe I 6301.5 Å and Fe I 6302.5 Å lines at the full resolution of the MHD simulations; those profiles are then convolved with telescope point-spread functions (including the Hinode spider obscuration), smeared in wavelength, and given $10^{-3}$ Ic polarimetric noise, before being interpreted with the standard observational recipes used for real sunspot data: line-bisector velocities at the 80% level, center-of-gravity magnetic field, far-wing magnetograms, and the three-lobe Stokes V selection. The undegraded reference is the physical data extracted on τ = 1 and τ = 0.1 surfaces in the simulation. This pipeline lets the same trade-offs of spatial resolution, spectral sampling, and noise thresholds act on simulated data, so the gap between native and degraded filling factors can be attributed to observational degradation rather than to the physics of the model.
What would settle it
Degrade the same synthetic profiles to a very high resolution, low-noise telescope (around 0.03 arcsec, comparable to DKIST) and measure the opposite-polarity fraction: if it stays near 11 percent instead of rising toward the native-resolution 21 percent, the claim that current observations hide most of the reversed field would be falsified. Alternatively, rerun the native-resolution analysis with the artificial top-boundary inclination replaced by a potential-field-matched boundary and check whether the opposite-polarity fraction falls below the degraded observational value.
Extended reading notes
Core claim
The paper establishes that observed penumbral magnetograms substantially underestimate the true opposite-polarity magnetic flux and convective downflow coverage, because both quantities are carried by small patches with steep gradients that typical spatial smearing and polarimetric noise erase. At the τ=1 surface (where the continuum forms), the simulation at native resolution puts 21 percent of the penumbral area in opposite polarity, while Hinode-class (0.5 m) degraded maps show only 11 percent; far-wing magnetograms drop from 17 percent native to 8 percent at 0.5 m. The downflow filling factor falls from 43 percent (τ-surface analysis) or 42 percent (bisector velocity) at native resolution to 32 percent or 14 percent at 0.5 m. Meanwhile, the apparent association of opposite polarity with downflows rises from roughly 70 percent native to 88 percent at 0.5 m in the τ-surface analysis and from 67 percent to 72 percent in the bisector/magnetogram analysis, because spatial smearing preferentially deletes small isolated reversed-field patches not tied to downflows. These quantities change little between simulations with 32 km and 12 km grid spacing, so the loss is attributed to observational degradation, not to numerical resolution.
Load-bearing premise
The load-bearing premise is that the simulated sunspot, whose top boundary condition artificially steepens the magnetic field inclination compared with a potential-field extrapolation, produces a realistic spatial distribution of opposite-polarity field and downflows; if that artificial inclination inflates the amount or location of reversed field, the reference fractions against which observations are judged would be biased.
Editorial extensions
If this is right
- Hinode-class measurements of penumbral opposite-polarity flux should be read as lower limits; the true small-scale reversed flux may be roughly twice as large.
- Spatial smearing systematically biases the measured opposite-polarity/downflow association upward, so reported strong associations partially reflect a selection effect and not just physical co-location.
- Detecting the elongated reversed-field patches along penumbral filament sides requires roughly 1 m class telescope resolution, or advanced inversion of 0.5 m data.
- The penumbral opposite-polarity and downflow fractions are robust to simulation grid spacing (32 km vs 12 km), so remaining disagreement between models and observations likely lies in telescope resolution and noise rather than numerical resolution.
- Future DKIST-class resolution and polarimetric sensitivity should reveal substantially more opposite-polarity field and convective downflows in the inner and middle penumbra than current magnetograms show.
Reading between the lines
- If instrumental degradation hides about half of the reversed flux, statistical studies of small-scale flux cancellation in sunspot penumbrae based on Hinode data may systematically underestimate the rate of flux submergence and the role of convection in recycling penumbral magnetic flux.
- The same forward-modeling degradation test could be applied to infrared spectral lines or to inversion codes; calibrating inversions on these synthetic degraded cubes could provide a way to recover the hidden fractions from real observations.
- The paper's comparison suggests that the especially high opposite-polarity/downflow association reported in some observational studies is partly method-induced; applying the three-lobe method and a spatially coupled 2D inversion to the same synthetic data would quantify that method bias.
- Because the simulation's top boundary artificially increases the inclination of the magnetic field, the native-resolution 21 percent value might itself be an overestimate for real sunspots; a high-resolution observation that recovers a substantially lower opposite-polarity fraction would indicate the absolute level is model-dependent, not just obscured.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a forward-modeling study of synthetic Stokes profiles from 3D MHD simulations of a sunspot penumbra (Rempel 2012) at two horizontal resolutions (32 km and 12 km). The synthetic profiles are degraded to mimic Hinode SP, 1 m, and 1.5 m telescope observations, and the authors use four diagnostic methods (tau-surface masks, far-wing Stokes V magnetograms, line-bisector velocities, and COG/COG+3-lobe profiles) to measure the area fraction of opposite-polarity magnetic field, the downflow filling factor, and their spatial association. The central result is that these quantities are strongly affected by spatial smearing and noise: for example, the opposite-polarity fraction at tau=1 drops from 21% at native resolution to 11% at Hinode/1 m resolution (Table 1), and the downflow filling factor in bisector maps drops from 42% to 14% at Hinode resolution (Table 2). The authors conclude that a significant fraction of opposite-polarity field and downflows is hidden in typical observations, and that these quantities are robust between the 32 km and 12 km simulation grids. The paper also discusses the performance of the 3-lobe method and the effect of noise.
Significance. If its quantitative conclusions hold, the paper provides a useful calibration of the observational visibility of small-scale opposite-polarity fields and convective downflows in sunspot penumbrae, and it supports the need for high-resolution, high-sensitivity instruments such as DKIST. The study combines full Stokes synthesis, realistic instrument degradation, and multiple retrieval methods in a single framework, and the use of two simulation resolutions is a positive internal check. The authors are transparent about the non-potential top boundary condition in Section 2, which is a known limitation of the Rempel (2012) models. However, the central quantitative claim (that observations hide a large fraction of the opposite-polarity/downflow signal) depends on the reference simulation's magnetic geometry at tau=1, and the paper does not quantify how much the artificial top boundary inflates the reported fractions. Because of this, the paper's contribution is a valuable but currently provisional estimate rather than a definitive number.
major comments (2)
- [Section 2 and Section 3.1] The top boundary condition of the reference simulations, acknowledged in Section 2 as artificially increasing the inclination angle of the magnetic field compared to a potential field extrapolation, may directly bias the central diagnostic: the area fraction of pixels with vertical (line-of-sight) field opposite to the spot polarity at tau=1. An artificially more horizontal background field reduces the initial vertical-field magnitude, so convective downflows need to bend field lines through a smaller angle to produce a sign reversal. Consequently, the reference fractions (21% native, 11% Hinode at tau=1; Table 1) and the associated conclusion that a significant fraction of opposite-polarity field is hidden in observations could be inflated relative to a sunspot with a more potential-like upper boundary. The 32 km vs 12 km comparison in Tables 1 and 2 changes only numerical resolution while keeping the same boundary condition, so it does not bound this systematic uncertainty. Please provide a sensitivity test (e.g., a run with a different top boundary or a comparison with a potential-field extrapolation) or explicitly quantify the possible bias and soften the conclusions accordingly.
- [Section 3.1, Tables 1-3] All quoted percentages are computed from a single snapshot for each resolution (one 32 km run and one 12 km run). The abstract and discussion claim that these quantities are 'robust within the simulations' based solely on the comparison of these two snapshots. This comparison conflates grid resolution with temporal evolution, since the two snapshots are taken from runs evolved for different durations (26 and 15 minutes, Section 2). To support the robustness claim, the authors should either analyze multiple snapshots (reporting mean and standard deviation or the full temporal spread) or explicitly state that the results are single-snapshot values and refrain from making a general robustness claim. This is load-bearing because the robustness statement is one of the four main results listed in the abstract.
minor comments (5)
- [Throughout] Please correct typographical errors: 'spatialy' to 'spatially', 'methodes' to 'methods', '1,5' to '1.5', 'uplows' to 'upflows', 'Feanz' to 'Franz', and 'GREGORE' to 'GREGOR'.
- [Abstract] The sentence 'forward modeling of synthetic Stokes profiles of the Fe I 6301.5 Å and Fe I 6302.5 Å lines)' contains an unmatched closing parenthesis; remove it.
- [Section 3.2] In the paragraph reporting downflow fractions, 'Degraded Hinode (0.5 m), 1 m and 1.5 m occupy 14%, 35%and 38%, respectively' has a missing space before 'and' and inconsistent spacing around percent signs; please standardize.
- [Section 3.3] The sentence 'The fraction of opposite polarity in downflows is lower in our analysis' is ambiguous: lower than what? Please specify the comparison (e.g., lower than in the tau=1 reference or lower than in Franz & Schlichenmaier 2013) so the reader can assess the claim.
- [Table 2] The table caption uses the phrase 'WITH AS WELL AS WITH OUT NOISE'; this should read 'with and without noise' for clarity.
Circularity Check
No significant circularity: the detectability claim is a forward-model output, not a fitted or self-referential quantity.
full rationale
The paper's derivation chain is self-contained as a forward-modeling study: it takes pre-existing MHD sunspot simulations (Rempel 2012), computes synthetic Stokes profiles, degrades them with instrument PSFs, spectral resolution, and noise, applies standard retrieval methods (bisector, COG, wing magnetograms, 3-lobe), and then measures opposite-polarity fractions and downflow associations. No parameter is fitted to observational data, and no 'prediction' is equivalent to an input by construction. The central conclusion that smearing and noise hide a significant fraction of opposite-polarity field and downflows follows directly from comparing degraded synthetic maps with the native-resolution simulation maps; it is not a renamed fit. The self-citations involved (Rempel 2012 for the simulations, Bharti et al. 2011 for the forward-modeling approach) are load-bearing only insofar as any use of a prior simulation as ground truth is load-bearing, but the cited simulations are external published products with their own assumptions and are not invoked as a uniqueness theorem or ansatz that makes the present results true by definition. The disclosed artificial top-boundary inclination (Section 2) is a legitimate physical-correctness concern about whether the reference simulations represent real sunspots, but it does not make the detectability analysis circular: the paper does not define opposite-polarity fraction in terms of the degraded observations, nor does it fit the simulation to the observations. Accordingly, no circular step meets the evidentiary standard of the review rules.
Assumptions & free parameters
free parameters (6)
- Far wing offset for Stokes V magnetograms =
20-30 mÅ
- Stokes V detection threshold =
3 sigma
- Bisector intensity level =
80% of continuum
- Velocity exclusion threshold =
100 m/s (absolute)
- Instrumental polarimetric noise =
10^-3 Ic
- Spectral and spatial smearing parameters =
Gaussian FWHM 21.54 mÅ for Hinode; 45 mÅ for 1 m/1.5 m; PSF mix chosen to match 12%/8% RMS contrast
assumptions (6)
- domain assumption LTE line formation with SPINOR/STOPRO is adequate for Fe I 6301.5 and 6302.5 lines
- domain assumption Vertical line-of-sight and disk-center geometry
- domain assumption Constant tau surfaces from a vertical ray approximate formation heights
- ad hoc to paper Top boundary condition artificially increases field inclination but does not distort the relevant penumbral structure
- domain assumption The simulation snapshots represent statistically typical penumbral states
- domain assumption Mask selection isolates penumbra without bias
Cite this review
Pith. "Pith review of Opposite polarity magnetic field and convective downflows in a simulated sunspot penumbra." pith.science (2026). https://pith.science/paper/HREMMDI4
@misc{pith2026190806439,
author = {Pith},
title = {Pith review of: Opposite polarity magnetic field and convective downflows in a simulated sunspot penumbra},
year = {2026},
howpublished = {\url{https://pith.science/paper/HREMMDI4}},
note = {Machine review of arXiv:1908.06439}
}
read the original abstract
Recent numerical simulations and observations of sunspots show a significant amount of opposite polarity magnetic field within the sunspot penumbra. Most of the opposite polarity field is associated with convective downflows. We present an analysis of 3D MHD simulations through forward modeling of synthetic Stokes profiles of the Fe\sci 6301.5 \AA~ and Fe\sci 6302.5 \AA~ lines). The synthetic Stokes profiles are spatially and spectrally degraded considering typical instrument properties. Line bisector shifts of the Fe\sci 6301.5 \AA~ line are used to determine line-of-sight velocities. Far wing magnetograms are constructed from the Stokes V profiles of the Fe\sci 6302.5 \AA~ line. While we find an overall good agreement between observations and simulations, the fraction of opposite polarity magnetic field, the downflow filling factor and the opposite polarity-downflow association are strongly affected by spatial smearing and presence of strong gradients in the line-of-sight magnetic field and velocity. A significant fraction of opposite polarity magnetic field and downflows are hidden in the observations due to typical instrumental noise. Comparing simulations that differ by more than a factor of two in grid spacing we find that these quantities are robust within the simulations.
Reference graph
Works this paper leans on
-
[1]
Berdyugina, S. V ., Solanki, S. K., Frutiger, C., 2003, A&A, 4 12, 513
work page 2003
-
[2]
H., Rempel, M., Hirzberger, J., Sola nki, S
Bharti, L., Cameron, R. H., Rempel, M., Hirzberger, J., Sola nki, S. K. 2012, ApJ, 752, 128
work page 2012
- [3]
- [4]
-
[5]
Borrero, J. M., & Ichimoto, K. 2011, Living Rev. Sol. Phys., 8 , 4 del Toro Iniesta, J. C., Bellot Rubio, L. R., & Collados, M. 20 01, ApJ, 549, L139
work page 2011
-
[6]
E., Bellot Rubio, L., de la Cruz Rodríguez
Esteban, P . E., Bellot Rubio, L., de la Cruz Rodríguez. 2015, ApJ, 803, 93
work page 2015
- [7]
- [8]
Show all 31 references
-
[9]
Ichimoto, K. eta al. 2007, PASJ, 59, 593
2007
-
[10]
K ., Aznar Cuadrado, R., Merenda, L
Joshi, Jayant, Pietarila, A., Hirzberger, J., Solanki, S. K ., Aznar Cuadrado, R., Merenda, L. 2011, ApJ, 734, L18
2011
-
[11]
K., Tiwari, S
Joshi, J., Lagg, A., Hirzberger, J., Solanki, S. K., Tiwari, S. K. 2017, A&A, 599, A35
2017
-
[12]
B., Kiselman, D., Löfdahl, M
Langhans, K., Scharmer, G. B., Kiselman, D., Löfdahl, M. G., & Berger, T. E. 2005, A&A, 436, 1087
2005
-
[13]
2011, IAU symposium series, 273, 8-14, ed
Rempel, M. 2011, IAU symposium series, 273, 8-14, ed. Prasad Choudhary, D. & Strassmeier, K. G
2011
-
[14]
2011, Living Rev
Rempel, M., & Schlichenmaier, R. 2011, Living Rev. Sol. Phys ., 8, 3
2011
-
[15]
H., & Knölker, M
Rempel, M., Schüssler, M., Cameron, R. H., & Knölker, M. 2009 b, Science, 325, 171 Ruiz Cobo, B. & Asensio Ramos, A. 2013, A&A, 549, L4
2009
-
[16]
Scharmer, G. B. 2014, A&A, 561, 31
2014
-
[17]
B., de la Cruz Rodriguez, J., Sütterlin, P ., Hen riques, V
Scharmer, G. B., de la Cruz Rodriguez, J., Sütterlin, P ., Hen riques, V . M. J. 2013, A&A, 553, 63
2013
-
[18]
Scharmer, G. B. & Henriques, V . M. J. 2012, A&A, 540, A19
2012
-
[19]
2012, Astronomische Nachrichten, 333, 796
Schmidt, W., von der Lühe, O., V olkmer, R., et al. 2012, Astronomische Nachrichten, 333, 796
2012
-
[20]
B., Henriques, V
Scharmer, G. B., Henriques, V . M. J., Kiselman, D., & de la Cru z Rodríguez, J. 2011, Science, 333, 316
2011
-
[21]
B., Bjelksjö, K., Korhonen, T
Scharmer, G. B., Bjelksjö, K., Korhonen, T. K., Lindberg, B. , & Petterson, B. 2003, in Society of Photo-Optical Instrumentation Engin eers (SPIE) Conference Series, V ol. 4853, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. S. L. Keil & S. V ...
2003
-
[22]
K., Berdyugina, S
Shelyag, S., Schüssler, M., Solanki, S. K., Berdyugina, S. V ., Vögler, A. 2004, A&A, 427, 335
2004
-
[23]
Solanki, S. K. 2003, A&ARv, 11, 15
2003
-
[24]
H., & Weiss, N
Thomas, J. H., & Weiss, N. O. 2004, ARAA, 42, 517
2004
-
[25]
H., & Weiss, N
Thomas, J. H., & Weiss, N. O. 2008, Sunspots and Starspots (Ca mbridge University Press)
2008
-
[26]
K., van Noort, M., Lagg, A., Solanki, S
Tiwari, S. K., van Noort, M., Lagg, A., Solanki, S. K. 2013, A& A, 557, 25
2013
-
[27]
K., van Noort, M., Solanki, S
Tiwari, S. K., van Noort, M., Solanki, S. K., Lagg, A. 2015, A& A, 583, A119
2015
-
[28]
Tiwari, S. K. 2017, 2017arXiv171207174T
2017
-
[29]
2003, ApJ, 592, 1225 van Noort, M., Lagg, A., Tiwari, S
Uitenbroek, H. 2003, ApJ, 592, 1225 van Noort, M., Lagg, A., Tiwari, S. K., Solanki, S. K. 2013, A& A, 557, 24 van Noort, M. 2012, A&A, 548, 5 Westendorp Plaza, C., del Toro Iniesta, J. C., Ruiz Cobo, B., et al. 1997, Nature, 389, 47 Westendorp Plaza, C., del Toro Iniesta, J. ...
2003
-
[30]
2001, ApJ, 547, 1148
Pillet, V . 2001, ApJ, 547, 1148
2001
-
[31]
L., Solanki,S
Zakharov, V ., Hirzberger, J., Riethmüller, T. L., Solanki,S. K., Kobel, P . 2008, A&A, 488, L17
2008
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.