{"id":"89cb8d04-79f4-4c1a-a213-3cdc2fb72b0d","arxiv_id":"1908.06439","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Synthetic observations of simulated sunspot penumbrae show that spatial smearing and instrumental noise hide most opposite-polarity magnetic field patches and downflows, while overestimating their association.","lead":"This paper takes a computer simulation of a sunspot and converts it into synthetic telescope images to test how well real instruments can see small regions of reversed magnetic field and downward flows in the penumbra. It shows that typical smearing and noise hide a significant fraction of these features and that correctly measuring them needs higher resolution or advanced inversion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The artificial top-boundary inclination disclosed in Section 2 may inflate the reference opposite-polarity fraction, so the 'hidden in observations' conclusion is not yet tied to real sunspot geometry.","rationale":"The paper is a clean forward-modeling exercise: take a known MHD snapshot, synthesize Stokes profiles, degrade them to instrument properties, and compare detection diagnostics against a simulation truth. The central quantitative claim is that smearing and noise hide a significant fraction of opposite-polarity field and downflows. I looked first for internal inconsistency or circularity. The analysis pipeline is self-consistent, the tables support a strong resolution dependence, and the 32 km versus 12 km comparison is a genuine, if limited, convergence check. The non-monotonic behavior at 1.5 m resolution in Table 1 is odd, but it is presented as a smearing effect and does not overturn the qualitative conclusion. The one premise on which everything else depends is the simulated reference. The authors themselves flag in Section 2 that the top boundary artificially increases the inclination angle relative to a potential-field extrapolation. Since opposite polarity is defined by the sign of the vertical component of an already strongly inclined field, changing the mean inclination can directly change how much of the penumbra appears reversed. The grid-spacing robustness check does not address this, because both runs share the same boundary treatment. Thus the absolute fractions (21%, 17%, and the degraded values) and the inferred 'hidden fraction' are conditionally valid at best. The reader's CONDITIONAL verdict already captures this boundary-condition uncertainty, so no verdict change is needed; the condition should be made explicit. I see no basis for REJECT: the forward-modeling logic is sound, and the qualitative claim that resolution hides fine-scale opposite polarity is likely to survive even if the absolute numbers shift.","tokens_in":11562,"tokens_out":8219,"duration_ms":92663,"concrete_test":"Re-run the 32 km simulation with the top boundary replaced by a potential-field-compatible condition, or extend the vertical domain so the field can relax before the boundary, keeping line synthesis, masking, and degradation identical; then recompute the native and Hinode-degraded opposite-polarity fractions and downflow associations in Tables 1 and 2. If the native tau = 1 fraction moves by more than about 5 percentage points from 21%, or the Hinode-degraded value moves outside the 8-11% range, the central 'hidden fraction' conclusion is boundary-condition-dependent and the paper should be revised to state that dependence explicitly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing concern is that the reference simulation is not a neutral test bed for the opposite-polarity fraction. Section 2 states that the computational box uses a top boundary about 700 km above the photosphere that 'artificially increases the inclination angle of the magnetic field compared to a potential field extrapolation.' The core diagnostic in Section 3.1 is the area fraction of pixels where the vertical (line-of-sight) field at tau = 1 has sign opposite to the spot polarity. In a background field that is artificially more horizontal, a convective downflow must bend field lines through a smaller angle before B_z changes sign, so the measured 21% native and 11% Hinode-degraded opposite-polarity fractions, and the associated downflow associations, could be inflated relative to a sunspot whose upper boundary is more potential-like. The 32 km versus 12 km resolution check changes only numerical resolution while keeping the same boundary condition, so it does not bound this systematic error. Without a sensitivity study of the top boundary, the paper's statement that a significant fraction of opposite-polarity field and downflows are hidden in typical observations rests on an unverified reference geometry.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":11769,"tokens_out":7538,"duration_ms":70328,"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":[{"comment":"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":"Section 2 and Section 3.1"},{"comment":"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.","section":"Section 3.1, Tables 1-3"}],"minor_comments":[{"comment":"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'.","section":"Throughout"},{"comment":"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":"Abstract"},{"comment":"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":"Section 3.2"},{"comment":"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.","section":"Section 3.3"},{"comment":"The table caption uses the phrase 'WITH AS WELL AS WITH OUT NOISE'; this should read 'with and without noise' for clarity.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed forward-modeling study, but the boundary-condition issue is the main risk to the central quantitative claim. I recommend major revision and encourage the authors to add a sensitivity test or restrict their conclusions to the Rempel (2012) model class. The paper is within scope for the journal and would be a useful contribution once this systematic uncertainty is addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely useful: it takes a known simulation, runs it through a realistic forward-modeling pipeline, and produces a quantitative catalog of how opposite-polarity magnetic field and downflow fractions change at native, Hinode, 1 m, and 1.5 m resolution using four different retrieval methods. The multi-resolution and multi-method comparison is new, and the 32 km versus 12 km robustness check is a legitimate addition. The tables are consistent with the stated trends, and the authors are transparent about several limitations. The central qualitative conclusion, that ordinary observations hide a significant fraction of the small-scale opposite-polarity flux and downflows, holds up in the sense that degrading resolution and adding noise consistently reduces detected fractions across methods and both simulation resolutions.\n\nThe main soft spot is exactly what the stress-test note flags. Section 2 states that the top boundary sits about 700 km above the photosphere and 'artificially increases the inclination angle of the magnetic field compared to a potential field extrapolation.' The paper's reference fractions at tau=1 are defined by sign changes in the vertical field. If the background field is artificially more horizontal, downflows have to bend field lines through a smaller angle to produce opposite polarity, so the 21% native and 11% Hinode-degraded fractions are plausibly inflated relative to a more potential-like upper boundary. The 12 km run uses the same boundary condition and therefore does not bound this systematic error. I do not think this kills the paper, because the direction of the resolution effect is robust, but the absolute calibration numbers should not be quoted as real-sunspot quantities without a sensitivity study or much stronger caveats.\n\nTwo smaller issues: the percentages come from single snapshots with no uncertainty estimates, and the detection thresholds are chosen from the same data, so the exact numbers carry more weight than they should. The absence of public code and data is a reproducibility cost, though not a fatal one since the pipeline is described in reasonable detail.\n\nThis paper is for observers who want to calibrate what Hinode and 1 m class instruments are missing in penumbral magnetograms, and for simulators comparing model predictions to those observations. It deserves a serious referee. I would send it out with a request for a sensitivity test or an explicit discussion of how the top-boundary inclination affects the reference fractions, and I would not let the current version pass as is.","headline":"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.","tokens_in":12321,"tokens_out":1437,"would_cite":true,"duration_ms":18214,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulated sunspot penumbra contains about twice as much opposite-polarity magnetic field as current telescopes can detect.","keywords":["sunspot penumbra","opposite polarity magnetic field","convective downflows","Stokes polarimetry","forward modeling","spatial smearing","MHD simulations","filling factors"],"falsifier":"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.","tokens_in":11310,"feed_emoji":"🌞","tokens_out":5675,"duration_ms":55397,"temperature":0.7,"pith_summary":"This paper asks how much of the small-scale opposite-polarity magnetic field and convective downflows predicted by sunspot simulations would actually be visible with present-day telescopes. It forward models synthetic Stokes profiles from two MHD sunspot simulations, degrades them to Hinode (0.5 m), 1 m, and 1.5 m telescope resolutions with realistic noise, and measures how the recovered fractions change. The central finding is that spatial smearing and noise hide a substantial part of the opposite-polarity field and downflows: at the τ=1 surface, the fraction of penumbral area covered by opposite polarity drops from 21 percent at native simulation resolution to about 11 percent at 0.5 m resolution. Degradation also inflates the apparent association between opposite polarity and downflows, because only the strongest, largest patches survive detection. If correct, published penumbral magnetograms systematically underestimate small-scale reversed flux and its role in sunspot convection.","feed_headline":"Telescope blur hides half the reversed field in sunspots","feed_subtitle":"Simulated penumbra has 21 percent opposite-polarity coverage; Hinode-class views catch only about 11 percent.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the MHD sunspot simulations (32/16 km and 12/8 km grid spacings) whose penumbral opposite-polarity field and downflows are analyzed.","marker":"Rempel (2012)"},{"why":"Establishes the forward-modeling and degradation approach for synthetic Fe I profiles compared with different telescope resolutions.","marker":"Bharti et al. (2011)"},{"why":"Provides the far-wing magnetogram and three-lobe Stokes V method used to retrieve opposite-polarity fields, along with the observed baseline to compare against.","marker":"Franz & Schlichenmaier (2013)"},{"why":"Supplies the wing magnetogram construction technique and the earlier observed association between opposite polarity and downflows.","marker":"Ichimoto et al. (2007)"},{"why":"Provides the Hinode point-spread function with central obscuration and spider used for the 0.5 m spatial convolution.","marker":"van Noort (2012)"},{"why":"Contributes the SPINOR/STOPRO radiative transfer code that computes the synthetic LTE Stokes profiles.","marker":"Berdyugina et al. (2003)"},{"why":"Supplies high-resolution 1 m observations showing elongated opposite-polarity patches, used as the comparison for 1 m degraded synthetic maps.","marker":"Scharmer et al. (2013)"},{"why":"Offers independent simulation and Hinode results at different optical depths that the paper finds consistent with its own filling-factor findings.","marker":"Joshi et al. (2017)"}],"fun_headline_variants":["Telescope smearing hides half of sunspot's reversed magnetic field","Observational blur conceals most opposite-polarity flux in sunspot penumbra","Sunspot simulations: true opposite-polarity field double what telescopes see","Half of sunspot reversed-field patches vanish in real telescope views"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Telescope smearing hides half of sunspot's reversed magnetic field","Observational blur conceals most opposite-polarity flux in sunspot penumbra","Sunspot simulations: true opposite-polarity field double what telescopes see","Half of sunspot reversed-field patches vanish in real telescope views"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001437,"raw_usage":{"total_tokens":5838,"prompt_tokens":1035,"completion_tokens":4803,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":651,"completion_tokens_details":{"reasoning_tokens":4725}},"tokens_in":651,"tokens_out":4803,"duration_ms":34399,"temperature":1.0,"reasoning_tokens":4725,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:45:52.754196+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"2011, ApJ, 739, 35","cited_arxiv_id":null,"evidence_quote":"Establishes the forward-modeling and degradation approach for synthetic Fe I profiles compared with different telescope resolutions."},{"cited_title":"& Schlichenmaier, R","cited_arxiv_id":null,"evidence_quote":"Provides the far-wing magnetogram and three-lobe Stokes V method used to retrieve opposite-polarity fields, along with the observed baseline to compare against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the wing magnetogram construction technique and the earlier observed association between opposite polarity and downflows."},{"cited_title":"V ., Solanki, S","cited_arxiv_id":null,"evidence_quote":"Contributes the SPINOR/STOPRO radiative transfer code that computes the synthetic LTE Stokes profiles."},{"cited_title":"B., de la Cruz Rodriguez, J., Sütterlin, P ., Hen riques, V","cited_arxiv_id":null,"evidence_quote":"Supplies high-resolution 1 m observations showing elongated opposite-polarity patches, used as the comparison for 1 m degraded synthetic maps."},{"cited_title":"K., Tiwari, S","cited_arxiv_id":null,"evidence_quote":"Offers independent simulation and Hinode results at different optical depths that the paper finds consistent with its own filling-factor findings."}],"review_version":1}