{"id":"168563a6-9597-472c-a6a8-3fe0155c12f2","arxiv_id":"2507.16283","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The magnetic field and plasma-beta are mapped along isolated sunspot fan loops from the photosphere to the corona for the first time.","lead":"Using the Sun's own 3-minute oscillations as tracers, researchers followed individual magnetic loops from a sunspot up into the corona and measured how the magnetic field weakens along the way. The field drops from about 2,000 Gauss at the surface to roughly 200 Gauss at the loop's base in the corona, giving new inputs for sunspot and wave models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The derived B(h) and beta(h) rest on correlation-contour areas at 94±2% of maximum correlation, calibrated at only the coronal footpoint; if 3-min oscillation patches include background umbral coherence, the area expansion and hence the headline field strengths are systematically biased.","rationale":"The reader's weakest assumption points to the correlation-contour areas, and I agree that this is the load-bearing step. The area enters directly in Eq. 5, so a systematic error in A(h) changes B(h) and beta(h) in a way that no amount of careful DEM analysis or model scaling can fix. The paper does have some supporting evidence: the coronal-footpoint contour is calibrated to an independent FWHM measurement (within errors), the derived coronal fields (158-236 G) are consistent with independent estimates by Brooks et al. (2021) and Gupta & Nayak (2022), and the overall beta(h) pattern resembles the Gary (2001) model. These give the result plausibility, but they do not validate the lower-atmosphere areas, because the calibration is at one height only and the Gary model is not a per-loop measurement. I did not elevate the density-scaling factor (dividing the Fontenla model by 2) to primary status because it affects only beta, not the B-field claim, and beta is already very small in the corona; however, it does add uncertainty to the beta=1 height, which is quoted without error bars. Also, the photospheric area Ap is sub-pixel for HMI, so Bp is an average over a larger region; this is a secondary concern. The proposed forward-model test would directly determine whether the wave-correlation patches measure flux-tube cross-sections or the coherence length of umbral oscillations, and would settle the concern. Given the currently unresolved nature of this assumption, the reader's CONDITIONAL verdict is appropriate; I would not change it.","tokens_in":14074,"tokens_out":10655,"duration_ms":121042,"concrete_test":"Perform a controlled numerical experiment: take a model umbral atmosphere with a known expanding flux tube (prescribed A(h), B(h), density, temperature), inject 3-min slow magnetoacoustic waves guided by the tube together with a realistic spatially coherent background umbral oscillation, synthesize IRIS/AIA-like data cubes at the observed passbands, cadence, resolution and noise, and apply the exact correlation-contour pipeline (Section 3.3, Appendix C). Compare the recovered A(h) and B(h) with the known truth; if the recovered coronal-footpoint B deviates by more than the quoted error bars or beta(h) crosses 1 incorrectly, the correlation-patch identification is biased and the central claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 3.3 and Appendix C, lower-atmosphere loop cross-sections are obtained by summing pixels inside contours at 94±2% of the maximum correlation coefficient in 3-min filtered correlation images. The contour level is calibrated to the FWHM area at the coronal footpoint only (92% for that height), and the same level is assumed to represent the true flux-tube cross-section at every lower height. The 3-min band is dominated by omnipresent umbral oscillations; a correlation patch can be broadened by the spatial coherence of this background oscillation field or by the limited spatial resolution of IRIS/AIA (0.332 arcsec/pixel) rather than by the true magnetic flux-tube width. Since B(h) is derived from B(h)=Bp Ap/A(h) (Eq. 5), any systematic error in A(h) propagates linearly into B(h) and quadratically into plasma-beta (Eq. 2). A further red flag is that Ap ≈ 0.22 arcsec^2 at the photosphere is smaller than one HMI 0.5-arcsec pixel, so the measured Bp=1928±13 G is itself a coarse average. The paper gives no independent verification of the lower-atmosphere areas (e.g., from an extrapolation or a synthetic test) and no estimate of the sensitivity of B(h) to the chosen contour level. If the correlation patches are contaminated by background oscillations, the headline decrease from 1596-2269 G to 158-236 G, and the beta<1 conclusion, could be systematically wrong.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper derives magnetic field strength and plasma-β along individual umbral fan loops from the photosphere to the corona. The authors use IRIS and SDO observations, trace loops in the lower atmosphere via 3-min slow-wave correlation analysis, obtain cross-sectional areas at multiple heights, and then apply flux conservation, B(h) = Bp Ap / A(h), with Bp from the HMI magnetogram. For the representative Loop 2 in AR 12470 they report a photospheric field of 1928 ± 13 G, a coronal footpoint field of 158 ± 50 G, β ≈ 0.3 at the photospheric footpoint and β ≈ 0.0006 at the coronal footpoint, and β < 1 along the whole loop with the β ≈ 1 layer at sub-photospheric heights. Similar results are summarized for six loops in two active regions, and the trends are compared with the Gary (2001) umbral model.","tokens_in":14406,"tokens_out":4147,"duration_ms":47492,"significance":"If the area estimates are reliable, this is a genuinely novel combination of wave seismology and flux-tube theory: it gives the first height-resolved magnetic-field and plasma-β profiles along individual umbral fan loops spanning photosphere to corona. The method is empirically grounded (Bp is measured, not fitted) and the paper reports error bars and a transparent propagation of uncertainties. The results are potentially important for wave-heating models and for constraining umbral atmosphere models. The main weakness is that the decisive lower-atmosphere cross-sectional areas rest on a correlation-contour calibration performed at a single height, so the headline B(h) and β(h) profiles require additional validation before the claims can be considered robust.","major_comments":[{"comment":"The lower-atmosphere cross-sectional areas, which enter Eq. (5) linearly, are defined by contours at 94 ± 2% of the maximum correlation coefficient, but the contour level is calibrated to the FWHM-based area only at the coronal footpoint (92%). The manuscript provides no sensitivity test of the chosen level and no independent verification of these areas. Because the 3-min band is dominated by umbral oscillations, the correlation patches could be broadened by background spatial coherence or by the 0.332 arcsec/pixel resolution rather than by the true flux-tube cross-section. A systematic error in A(h) propagates linearly into B(h) via Eq. (5) and quadratically into β via Eq. (2). I request a quantitative sensitivity study (e.g., repeating the analysis for contour levels from 90% to 98%) and, if possible, a comparison with an independent estimate such as magnetic flux conservation from HMI at the footpoint or a modeled flux-tube expansion.","section":"Section 3.3, Appendix C, Eq. (5)"},{"comment":"The density in the lower atmosphere is taken from the Fontenla sunspot model and then arbitrarily divided by a factor of 2 to match the DEM density at the coronal footpoint. This scaling factor has no quoted uncertainty and is not propagated into the β estimates. Since β ∝ N, the reported βp values in Table 1 and the qualitative conclusion β < 1 along the whole loop depend directly on this ad hoc factor. The paper should either justify the scaling factor with a fitting procedure over the full height range or explicitly include it in the systematic-error budget and show how βp changes if the factor is varied over a plausible range.","section":"Section 3.1, Fig. 2, Eq. (2)"},{"comment":"The photospheric magnetic field Bp is obtained from HMI at a formation height of approximately 269 km, but the photospheric loop area Ap ≈ 0.22 arcsec² is smaller than one HMI pixel (0.5 arcsec ≈ 0.25 arcsec²). The manuscript does not state whether Bp is the unsigned line-of-sight field or a vector-field magnitude, nor how the coarse HMI spatial resolution and the reported µ ≈ 0.90 projection geometry affect the value inserted into Eq. (5). This is important because Bp sets the absolute normalization of B(h) and therefore of β(h). Please clarify the precise quantity used and assess the resulting systematic uncertainty.","section":"Section 3.4, Table 1"},{"comment":"The statement that the β ≈ 1 layer lies at sub-photospheric heights is an extrapolation, not a measurement: the authors fit a straight line to β values between the temperature-minimum region and the photosphere, and extend the line to β = 1 at negative heights. This extrapolation is model-dependent and no uncertainty is given for the quoted heights in Table 1 (e.g., −617 km for Loop 2). Since β < 1 at all measured heights, the sub-photospheric location of the β = 1 layer should be presented as a tentative inference, with an error estimate or at least a clear caveat.","section":"Section 3.4, Fig. 4"}],"minor_comments":[{"comment":"The final sentence of the abstract, 'we provide estimates on the magnetic field and plasma-β variations along the various fan loops traced from the photosphere to the corona using 3-min waves and along the corona', contains a redundant 'and along the corona' and should be rewritten for clarity.","section":"Abstract"},{"comment":"The paper states that the temperature sensitivity of each passband is 'well known' and then assigns a single formation height to each filter. Since IRIS/AIA passbands have broad response functions, a brief discussion of the systematic uncertainty in the assigned formation heights would strengthen the density and β profiles.","section":"Section 3.1"},{"comment":"The assumed 15% error in density and temperature is stated without justification or a reference for this specific value; given that these errors dominate the β uncertainty, a short justification would be helpful.","section":"Section 3.1"},{"comment":"The black solid lines representing the Gary (2001) umbral model are useful, but the paper does not quantify the discrepancy between the model β values and the loop β values; adding a brief quantitative statement would make the comparison more informative.","section":"Figure 4"}],"recommendation":"major_revision","confidential_remarks":"This is a promising paper with a novel method, but the central claim rests on correlation-contour areas calibrated at a single height and on an ad hoc density scaling factor. Both issues are fixable with additional analysis, so I recommend major revision rather than rejection. I would not reject based on disagreement with the Gary (2001) model; the internal validation is the main concern."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe one thing to know: this is the first time anyone has produced B(h) and beta(h) profiles along individual sunspot fan loops from the photosphere to the corona. The underlying wave-tracing method was published before by the same group, and flux conservation is textbook, but the combination — using 3-min slow waves to define loop cross-sections at multiple heights, then applying B(h) = Bp Ap / A(h) — is a real step forward. The coronal values (158 ± 50 G at the footpoint) agree with independent spectroscopic and extrapolation estimates, which gives the result some external anchor.\n\nThe paper does a good job of comparing with Gary's (2001) active-region model and shows a similar trend with different absolute beta values, which is physically plausible for isolated umbral loops as opposed to a region average. The data handling is careful: co-alignment, background subtraction, DEM analysis, and error propagation are all described in enough detail to reproduce.\n\nThe soft spot is the lower-atmosphere area measurement. The 94 ± 2% correlation contour is calibrated at the coronal footpoint only, where it matches the FWHM area, and then used at every height below. If the 3-min correlation patches at those heights are broadened by the background umbral oscillation field or by spatial resolution rather than by the actual flux-tube width, the area expansion is systematically wrong, and B(h) scales linearly with that error. The photospheric footpoint area (0.22 arcsec^2) is smaller than one HMI pixel, so the 'photospheric Bp' is really an average over a larger region. There is also an arbitrary factor of 2 in the density scaling to match DEM at the coronal footpoint, and the beta=1 layer heights are extrapolated without uncertainties. None of these is fatal on its own — the overall result, beta < 1 along the whole loop, is robust to plausible changes in the density scale — but the absolute values of B and beta in the lower atmosphere could be off by tens of percent.\n\nWho is this for? Solar physicists working on sunspot waves, loop seismology, and coronal heating. It is a solid within-subfield contribution, not a paradigm shift.\n\nRecommendation: send it to peer review. The method deserves scrutiny, and a good referee will ask for a sensitivity analysis of the contour level and a synthetic test of the area recovery. That should be a major revision, not a rejection.","headline":"First full-height B and beta profiles along individual sunspot fan loops, with a credible but under-tested lower-atmosphere area calibration.","tokens_in":14937,"tokens_out":2949,"would_cite":true,"duration_ms":28688,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Umbral fan loops traced by 3-min slow waves show magnetic field falling from about 2000 G at the photosphere to about 200 G in the corona, with plasma-beta below 1 along the whole loop.","keywords":["Solar atmosphere","Sunspots","Umbral fan loops","3-min slow waves","Plasma beta","Magnetic field strength","Magnetohydrodynamics","Coronal seismology"],"falsifier":"Compare the predicted $B(h)$ with an independent magnetic field measurement at a chromospheric or transition-region height along the same loops, e.g., spectropolarimetry of a chromospheric line or microwave imaging; if the measured field deviates from $B_p A_p / A(h)$ by more than the error bars at any height, the area-tracing assumption fails.","tokens_in":13882,"feed_emoji":"☀️","tokens_out":6131,"duration_ms":59212,"temperature":0.7,"pith_summary":"This paper claims that sunspot fan loops can be traced continuously from the photosphere to the corona by cross-correlating 3-minute slow oscillations across AIA, IRIS, and HMI passbands, and that this tracing makes it possible to map the magnetic field strength and plasma-beta along individual loops for the first time. Using flux conservation, the photospheric field measured by HMI (1596-2269 G across the sampled loops) is scaled by the measured loop cross-sectional area to obtain the field at every height, reaching 158-236 G at the coronal footpoints. The resulting plasma-beta is below unity all along the loops, with the beta=1 layer located 200-1000 km below the photosphere. If correct, the method turns the ubiquitous 3-min umbral oscillations into a seismological tool for probing the magnetic structure of sunspot atmospheres and for testing MHD wave models in expanding waveguides.","feed_headline":"3-min sunspot waves map magnetic field from photosphere to corona","feed_subtitle":"Loop tracing shows the field falling from ~2000 G to ~200 G, with plasma-beta below 1 throughout.","key_machinery":"The load-bearing identity is magnetic flux conservation, $B(h)A(h)=B_p A_p$, applied along a loop whose cross-sectional area $A(h)$ is measured at each atmospheric height. In the lower atmosphere the area comes from closed contours at $94\\pm2\\%$ of the maximum correlation-coefficient value in correlation images of 3-min filtered light curves, calibrated at the coronal footpoint against the FWHM-derived diameter; in the corona, $A(h)$ comes from Gaussian fits to cross-loop intensity profiles. The photospheric anchor $B_p$ is the RMS HMI magnetogram value, and the temperature and density needed for plasma-$\\beta$ come from a sunspot model at lower heights and from differential emission measure analysis in the corona, with the loop inclination fixed by comparing observed and hydrostatic density scale heights.","core_discovery":"The authors establish that the magnetic field and plasma-beta vary along individual umbral fan loops from the photosphere to the corona, using the loops' own 3-min slow oscillations as tracers. For the representative Loop 2, the field drops from 1928 +/- 13 G at the photospheric footpoint to 158 +/- 50 G at the coronal footpoint; across all studied loops the photospheric range is 1596-2269 G and the coronal range 158-236 G. Plasma-beta at the footpoints is 0.2-0.5 at the photosphere and 0.0001-0.001 in the corona, and the loops remain below beta=1 throughout, implying the magnetized plasma is magnetic-pressure-dominated in the whole visible atmosphere. The obtained patterns resemble the previously established active-region model, but with lower chromospheric beta values, which the authors attribute to tracing isolated loops rather than the integrated umbra.","pith_inferences":["If the $\\beta=1$ depth extrapolation is sensitive to the chosen linear fit, the true transition could be several hundred kilometres higher or lower; testing this requires measuring the plasma-beta profile independently in the temperature-minimum region.","The assumption that 3-min correlation patches represent the flux-tube cross-section at every height could be checked by comparing correlation-contour areas at chromospheric heights with radio or spectropolarimetric maps of the same umbral loops; disagreement would indicate contamination by the background umbral oscillation field.","Applied to loops without visible coronal emission, the technique might extend $B(h)$ mapping to faint or newly forming loops, where coronal FWHM calibration would need to be replaced by a different anchor."],"forward_implications":["A height-resolved $B(h)$ and $\\beta(h)$ profile is now available for individual umbral loops, not just for the global corona or for isolated loop segments, which is the input needed to model wave propagation in density-stratified, expanding flux tubes.","The sub-photospheric $\\beta=1$ layer found here identifies the region where mode conversion of slow waves can occur, an important constraint for umbral helioseismology.","Coronal fields of 158-236 G along these loops imply Alfv\\'en speeds and energy fluxes much larger than typical global coronal values, strengthening the case for Alfv\\'en-wave heating models powered by umbral oscillations.","The same correlation-tracing recipe can be applied to other multi-passband sunspot observations to build a statistical sample of loop magnetic parameters."],"supporting_citations":[{"why":"Supplies the technique of tracing fan loops from the corona to the photosphere using 3-min oscillations and obtaining cross-sectional areas at lower heights.","marker":"A. Rawat & G. Gupta 2023"},{"why":"Provides the sunspot atmosphere model that gives temperature and total number density as functions of height in the lower atmosphere.","marker":"J. Fontenla et al. 1999"},{"why":"Provides the reference active-region plasma-beta variation model used for comparison and the definition of plasma-beta.","marker":"G. A. Gary 2001"},{"why":"Supplies the differential emission measure tool used to obtain coronal loop temperature and density.","marker":"I. G. Hannah & E. P. Kontar 2012"},{"why":"Provides spectroscopic coronal loop magnetic field and plasma-beta values against which the derived coronal field is compared.","marker":"D. H. Brooks et al. 2021"},{"why":"Supplies the density scale-height method used to derive the loop inclination from the observed versus expected scale heights.","marker":"M. J. Aschwanden et al. 1999"},{"why":"Describes the HMI instrument whose magnetogram data give the photospheric magnetic field strength $B_p$.","marker":"P. H. Scherrer et al. 2012"},{"why":"Describes the IRIS instrument whose slit-jaw images provide the lower-atmosphere passbands used for loop tracing.","marker":"B. De Pontieu et al. 2014"}],"fun_headline_variants":["Sunspot 3-min waves trace field drop from 2000 G to 200 G","Umbral loops mapped: plasma-beta below 1 from photosphere to corona","3-min slow waves reveal magnetic field and beta along sunspot loops","Seismology of umbral loops: field falls 10x, beta < 1 throughout","Fan loop profiling: 3-min waves gauge magnetic field and plasma-beta"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The cross-sectional area of each loop at lower atmospheric heights is taken to be the 94 +/- 2% maximum-correlation contour of the 3-min oscillation signal, a level calibrated at only the coronal footpoint and assumed to hold at all heights.","fun_headline_variants_meta":{"raw":{"variants":["Sunspot 3-min waves trace field drop from 2000 G to 200 G","Umbral loops mapped: plasma-beta below 1 from photosphere to corona","3-min slow waves reveal magnetic field and beta along sunspot loops","Seismology of umbral loops: field falls 10x, beta < 1 throughout","Fan loop profiling: 3-min waves gauge magnetic field and plasma-beta"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000214,"raw_usage":{"total_tokens":1488,"prompt_tokens":1073,"completion_tokens":415,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":689,"completion_tokens_details":{"reasoning_tokens":309}},"tokens_in":689,"tokens_out":415,"duration_ms":4808,"temperature":1.0,"reasoning_tokens":309,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:12:54.695531+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the predicted $B(h)$ with an independent magnetic field measurement at a chromospheric or transition-region height along the same loops, e.g., spectropolarimetry of a chromospheric line or microwave imaging; if the measured field deviates from $B_p A_p / A(h)$ by more than the error bars at any height, the area-tracing assumption fails.","supporting_citations":[{"cited_title":"H., Warren , H","cited_arxiv_id":null,"evidence_quote":"Provides spectroscopic coronal loop magnetic field and plasma-beta values against which the derived coronal field is compared."}],"review_version":1}