{"id":"c505e079-57a8-4e84-9f81-c31c908e8ddd","arxiv_id":"2608.05575","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Using 150 active regions, the paper finds helicity and energy are concentrated low in the corona and proposes an 81 Mm extrapolation cutoff for 97% retention.","lead":"Magnetic helicity and energy in the Sun's corona are concentrated near the surface, according to a study of 150 active regions. The authors recommend stopping magnetic field extrapolations at about 81 Mm to retain 97% of both quantities while saving about 38% of computation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 97% retention claim is built on slicing a single 133 Mm extrapolation rather than re-extrapolating at shorter heights; since NLFFF solutions and the reference potential field depend on box height, the recommended 81 Mm cutoff is not validated.","rationale":"We read the paper in good faith. The study is well-organized, uses a long SMFT dataset, and applies established NLFFF and finite-volume helicity methods. The quality metrics (σ_J, ⟨|f_i|⟩, E_div/E) are appropriate and reported with IQRs. The resolution tests in §3 (Fig. 6) show internal consistency of the slicing procedure. However, these do not address the core assumption that a sliced subvolume of a tall extrapolation is a proxy for a genuine shorter extrapolation. NLFFF solutions are known to depend on the computational box and its boundary conditions (Wiegelmann & Sakurai 2012, 2021); the optimization objective (Eq. 1) and the force-free residual are global, so moving the top boundary from 133 to 81 Mm changes the solution in the lower layers. Furthermore, the relative helicity calculation (§2.3) constructs a reference potential field from the normal component on each subvolume's boundary; for a sliced subvolume, that boundary is an artifact of the tall-box solution, not a solution of the physical boundary-value problem. Thus the 97% retention figure is not evidence about what a real 81 Mm extrapolation would deliver. The paper itself states the slicing procedure explicitly, so the limitation is not hidden. But because the abstract and conclusions present 'at least 81 Mm' as a validated cutoff with a specific accuracy and cost saving, the central claim is over-stated. A modest revision could reframe the result as a truncation study of tall boxes, or add the re-extrapolation validation proposed above. Given the reader's verdict of REJECT with moderate confidence, we find no reason to change it; the concern is the same and it lands.","tokens_in":11953,"tokens_out":6128,"duration_ms":60025,"concrete_test":"Select a subsample of ~10 active regions (3–4 from each flux group). For each, run the same NLFFF code to convergence in four boxes with identical lateral dimensions and resolution but heights of 49, 81, 100, and 133 Mm (or the grid equivalents). Compute H_R and E from each full short box exactly as in §2.3–2.4, and compare these with the values obtained by slicing the 133 Mm run at the corresponding heights. If the short-box values differ from the sliced values by more than the 3% tolerance used to define the 97% threshold (or if the cumulative 97% heights shift by more than ~10%), the slicing proxy is invalid and the recommended heights must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim (abstract; §4) – that an extrapolation height of 81 Mm retains 97% of magnetic helicity/energy and saves ~38% cost – rests entirely on the procedure of §2.3: 'we do not re-extrapolate the field; instead, we slice the extrapolated 3D field to define specific subvolumes.' This assumes that the field in the lower 81 Mm of a 133 Mm NLFFF box is identical to the field that would be obtained from an actual NLFFF run with a box capped at 81 Mm. That assumption is untested and likely false. The NLFFF optimization (Eq. 1) minimizes the force-free and divergence-free residuals over the whole box; lateral and top boundary conditions differ when the box height changes, and these boundaries affect the solution inside, particularly for fields anchored near the top. Moreover, the relative helicity H_R (Eq. 4) depends on the reference potential field P that matches the normal component on the subvolume boundary (Eq. 5); for a sliced subvolume, P is computed from the tall-box field's normal component on the artificial top boundary, which is not the same as the potential field that would match the photospheric boundary in a genuine short-box run. Thus, even if the internal field were identical, H_R would differ. Because the 97% threshold and the 81 Mm/49 Mm recommendations derive from these sliced subvolumes, they do not establish the behavior of actual shorter extrapolations. The 38% cost saving is just volume arithmetic (1 - 81/133 ≈ 39%), but the accuracy claim is unsupported. A related unit-conversion inconsistency (grid height 111 ≈ 74 Mm, not 81 Mm, at the stated 200 points = 133 Mm spacing) further weakens the headline number.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes 150 SMFT active-region vector magnetograms (1988–2019) grouped by absolute flux, performs NLFFF extrapolations to a height of about 133 Mm, and computes relative magnetic helicity with the Coulomb–Yang finite-volume method and magnetic energy from the extrapolated field. By slicing the single extrapolated volume into nested subvolumes, the authors define cumulative helicity and energy ratios and find that 97% of the total is reached at about 81 Mm for helicity and 49 Mm for energy; they recommend these as lower bounds for future extrapolation heights and claim a computational saving of about 38%. The paper also reports quality metrics for the force-free and divergence-free solutions and tests the effect of reduced spatial resolution.","tokens_in":12322,"tokens_out":5567,"duration_ms":45398,"significance":"If the central recommendation were validated, the paper would supply a practical, observationally based calibration of NLFFF extrapolation heights for large statistical helicity studies, which is a genuinely useful goal given the computational cost of such surveys. The dataset of 150 active regions spanning three solar cycles is a valuable resource, and the authors report standard quality metrics (σ_J, divergence fraction, E_div/E) with medians and interquartile ranges, which is good practice. The energy part of the analysis is on firmer ground because magnetic energy density is local and additive, so the energy profile is a direct property of the extrapolated field. The paper also honestly exposes its main methodological shortcut by stating that it does not re-extrapolate the field, and it explicitly acknowledges the non-additivity of relative helicity.","major_comments":[{"comment":"The central quantitative claim — that an extrapolation height of 81 Mm retains 97% of magnetic helicity and energy — rests on the assumption that slicing a single 133 Mm NLFFF solution into nested subvolumes reproduces what a genuine extrapolation with a shorter box would yield. As stated in §2.3, 'we do not re-extrapolate the field; instead, we slice the extrapolated 3D field.' This assumption is not tested and is likely violated because the NLFFF optimization of Eq. (1) minimizes residuals over the entire box, and lateral/top boundary conditions change when the box height is changed, affecting the solution in the lower layers as well. The paper should re-extrapolate a representative subset of active regions with box heights near 49, 81, and 133 Mm and compare the helicity and energy in the overlapping lower volume; without such a test, the retention percentages in the abstract and §4 are not established for actual truncated extrapolations.","section":"§2.3 and §4"},{"comment":"The ratio H(h)/H_total used in Figure 5 and in the 97% statements is not the fraction of the full volume's magnetic helicity located below height h. For each subvolume, H(h) is the relative helicity computed with its own potential-field reference P satisfying Eq. (5) on the artificial top boundary, so H(h) is a different gauge-invariant quantity, not the restriction of the full-volume helicity to the lower layers. The authors acknowledge the non-additivity of relative helicity, but they still interpret H(h)/H_total as a contribution rate. This is a load-bearing interpretational issue: the statement 'extrapolation height of at least 81 Mm retains 97% of the total magnetic helicity' is not implied by the calculation. The authors should either adopt a quantity with an additive layer decomposition (e.g., field-line helicity densities or a gauge-invariant flux-weighted measure) or explicitly rephrase the abstract and conclusions as describing the ratio of subvolume-relative helicities, without claiming that this ratio equals the retained fraction of the total helicity.","section":"§2.3, Eq. (4)–(6), Fig. 5"},{"comment":"The choice of the 97% threshold is justified by stating that discrepancies among FV methods for computing helicity are typically within 3% (Valori et al. 2016). However, that 3% refers to differences between independent finite-volume methods applied to the same volume, not to the error introduced by truncating the extrapolation height. The sentence 'When the extrapolation reaches or exceeds these heights, the error in the computed energy will be less than 3%' is therefore not supported by the cited reference. The threshold is arbitrary in the sense that any quantile could be chosen; the paper should present the full cumulative curves (as it does) and let readers judge, rather than claiming a quantitative error bound that the data do not provide.","section":"§3, paragraph after Fig. 5"},{"comment":"The claimed ~38% computational cost saving is essentially volume arithmetic: 1 - 81/133 ≈ 0.39. The actual computational cost of the NLFFF optimization with multigrid and boundary weighting ω(x, y, z) in Eq. (1) is not necessarily proportional to the number of grid points; it may scale nonlinearly with box height and resolution, and the 38% figure should be supported by measured runtimes from actual shorter extrapolations, or explicitly labeled as a first-order volume-based estimate.","section":"§4 and abstract"}],"minor_comments":[{"comment":"The weight function ω(x, y, z) is described as 'set to 1 in the core region and gradually decreases to zero near the boundary,' but the exact functional form and the thickness of the boundary layer (10–20 grid points) are not given; please provide the precise prescription used, since it affects how subvolume boundaries are weighted in the sliced analysis.","section":"§2.2, Eq. (1)"},{"comment":"The 97% heights (60.9, 90.1, 111.2 grid units etc.) are reported as group averages without scatter; adding the interquartile range or standard deviation across the 50 active regions per group would make the recommendation more robust and allow readers to assess case-to-case variability.","section":"§3, Fig. 5 and §3 resolution test"},{"comment":"The description 'data scaled to 1/2 (pixel resolution is 2 arcseconds)' is ambiguous: please clarify whether the number of pixels is halved, the field of view is kept constant, and what the resulting physical box height is for the 200-point vertical grid in each resolution case.","section":"§3, resolution test and Fig. 6"},{"comment":"The sign convention in Eq. (4) uses (A + A_P)·(B − P), which is one common form, but the paper should cite the precise form from Yang et al. (2013, 2018) or state the corresponding boundary-term cancellation, so readers can verify the implementation.","section":"§2.3, Eq. (4)"},{"comment":"There are minor typographical issues: 'the heighth helicity' in §3, 'where the red, green, and blue curves' in the caption of Fig. 5 (the bottom panel appears to show fractional contributions), and 'magnetic energy, more than 97% accuracy in magnetic energy calculations' in §4 is redundant. These do not affect the science but should be cleaned up.","section":"§3 and §4"}],"recommendation":"major_revision","confidential_remarks":"The paper has a valuable dataset and a clearly stated methodology, but the central recommendation is not supported without re-extrapolation tests. I would suggest requiring a revision that either (a) re-extrapolates a subsample of active regions at shorter box heights and compares the lower-layer fields and helicity/energy values, or (b) substantially softens the claims in the abstract and conclusions to describe the results as properties of subvolume-relative helicity ratios rather than as retained fractions of total helicity. The energy recommendation is more robust because energy is additive, but the cost and 97%-accuracy statements need qualification. The manuscript is otherwise well organized and the quality metrics are reported transparently."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The 81 Mm/97% cutoff you see in the abstract is not supported by the evidence. The paper computes helicity and energy in nested subvolumes of one 133 Mm NLFFF extrapolation per region, but never re-extrapolates with a shorter box. That matters: NLFFF solutions and the potential field reference depend on the box height, so the lower 81 Mm of a tall run is not necessarily what you'd get from a genuine 81 Mm run. The authors are honest about the slicing in §2.3, but the abstract and conclusions use the sliced results to recommend a cutoff for future extrapolations, which is an overreach.\n\nWhat's genuinely new and useful is the statistical profile: 150 active regions, uniformly processed SMFT data, decent quality metrics, and a clear demonstration that both helicity and energy are concentrated at low coronal heights. The 97% retention heights (around 111 grid layers for helicity, 67 for energy) are a new quantitative result, and the resolution-downgrade test is a good sanity check. The paper also correctly notes the non-additivity of relative helicity when it computes layer contributions.\n\nSoft spots: (1) the re-extrapolation issue is load-bearing; a handful of shorter-box runs would test it, and they haven't done any. (2) There's a unit inconsistency: 200 grid points at 1 arcsec is ~145 Mm, not 133 Mm, and 111 grid layers at 0.665 Mm/layer is 74 Mm, not 81 Mm. That needs to be reconciled. (3) The 97% threshold is a choice, not derived, but it's openly stated.\n\nThis paper is a reasonable proof-of-concept for truncation of a given tall-box solution, not a validated prescription for box sizing. It deserves a serious referee, though—the dataset is large and the question is practical. I'd send it to review and ask the authors to re-extrapolate a subset of regions with shorter boxes and fix the unit conversion before accepting.","headline":"The 81 Mm/97% cutoff is not supported by the evidence—it's built on slicing one tall NLFFF box, not re-extrapolating—but the 150-region statistical profile is useful and the paper deserves a serious referee.","tokens_in":12833,"tokens_out":3584,"would_cite":false,"duration_ms":29634,"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":"Based on 150 solar active regions, this paper argues that extrapolating the coronal magnetic field to a height of 81 Mm retains 97% of the total magnetic helicity and energy, allowing computational costs to be cut by about 38%.","keywords":["magnetic helicity","magnetic energy","nonlinear force-free field","NLFFF extrapolation","solar active regions","coronal magnetic field","finite-volume helicity","SMFT magnetograms"],"falsifier":"Re-run the extrapolations with actual box heights of 81 Mm and 49 Mm using the same preprocessing, code, and boundary conditions, and compare the resulting relative helicity and energy with the values from the 133 Mm box. If the shorter boxes systematically return less than 97% of the full-box quantities, or if the lower layers of the field change substantially when the ceiling is lowered, the proposed calibration does not apply to real extrapolations. A complementary check is to extrapolate a few cases to a taller box (e.g., 160 Mm) to confirm that the 133 Mm reference is itself converged.","tokens_in":11756,"feed_emoji":"☀️","tokens_out":6517,"duration_ms":49516,"temperature":0.7,"pith_summary":"This paper asks how high a nonlinear force-free field (NLFFF) extrapolation must go to capture nearly all of an active region's magnetic helicity and energy. Using 150 SMFT vector magnetograms from 1988 to 2019, grouped by magnetic flux, the authors extrapolate each field once and then slice the volume into nested subvolumes to build vertical profiles. They find that 97% of the total relative magnetic helicity is reached at about 81 Mm, 97% of the magnetic energy at about 49 Mm, and that stopping at these heights reduces computational cost by roughly 38% relative to a 133 Mm box. The practical payoff is a physically motivated, rather than purely empirical, choice of extrapolation height for long-term statistical studies of helicity in active regions.","feed_headline":"81 Mm captures 97% of coronal helicity and energy","feed_subtitle":"150 active regions show a 38% computational saving with negligible accuracy loss.","key_machinery":"The method hinges on a nested-subvolume slicing of a single NLFFF extrapolation. For each magnetogram, one optimization-based NLFFF solution is computed in a box about 133 Mm tall, then subvolumes are defined from the photosphere upward. Relative magnetic helicity in each subvolume is computed with a Coulomb-gauge finite-volume method: it constructs a potential field matching the boundary normal component and evaluates $H_R = \\int_V (\\mathbf{A}+\\mathbf{A}_P)\\cdot(\\mathbf{B}-\\mathbf{P})\\,dV$. The fractional contribution of an incremental layer is estimated as $\\Delta H / H_{\\mathrm{total}}$, where $\\Delta H$ is the difference in helicity between two nested volumes; the paper explicitly notes this is a numerical estimate because relative helicity is not additive in subvolumes. Magnetic energy is the volume integral of $B^2/8\\pi$ over each subvolume. The 97% threshold is tied to the known approximate 3% spread among finite-volume helicity codes.","core_discovery":"The central claim is a height calibration for NLFFF extrapolation: for the SMFT active-region sample, cumulative relative magnetic helicity reaches 97% of its full-box value at about 81 Mm, and cumulative magnetic energy reaches the same fraction at about 49 Mm. Because finite-volume helicity methods agree with each other to within about 3%, the authors adopt the 97% retention level as the criterion for a sufficient extrapolation height. They therefore propose 81 Mm as a practical lower bound for helicity studies and 49 Mm for energy studies, with modest extensions permissible; going higher yields less than 3% extra accuracy while adding about 38% to computational cost. This result is presented as an empirical calibration that is consistent across three flux groups and only mildly sensitive to magnetogram resolution.","pith_inferences":["The slice-based calibration likely holds only if the upper layers of a full-box extrapolation are nearly independent of the lower boundary; a genuinely shorter extrapolation could alter the lower layers, so the 97% and 38% numbers should be rechecked by actually re-extrapolating at the proposed heights.","The 97% threshold is tied to the mutual 3% scatter among finite-volume helicity estimators, meaning the recommended height is an accuracy-matching convention rather than a physical transition height.","For eruption-focused studies, the missing 3% of helicity above 81 Mm may reside in exactly the high-reaching flux-rope structures that matter most, so the recommendation is safest for statistical budgets, not event analysis.","The observed resolution dependence suggests a simple operational rule: as magnetogram resolution degrades, raise the extrapolation height modestly; this could be tested directly by comparing 1, 2, and 4 arcsecond data for the same active regions."],"forward_implications":["Large-sample statistical studies of active-region helicity can cap NLFFF extrapolations at about 81 Mm and retain about 97% accuracy while saving roughly 38% of computational time.","Studies that need only magnetic energy can use an even lower cap of about 49 Mm with the same retention level.","The proposed heights are stable across active regions of different magnetic flux, so the calibration is not limited to one flux regime.","Lower-resolution magnetograms require slightly taller boxes (about 84 Mm at half resolution, 91 Mm at quarter resolution) to reach the same 97% helicity retention.","The 81 Mm cutoff is physically plausible because it encompasses typical filament heights and the heights at which quiescent filaments often begin to erupt."],"supporting_citations":[{"why":"Provides the Coulomb-gauge finite-volume method used to compute relative magnetic helicity in every subvolume.","marker":"Yang et al. (2013, 2018)"},{"why":"Supplies the NLFFF model and optimization formulation that underlies all the coronal extrapolations.","marker":"Wiegelmann and Sakurai (2012)"},{"why":"Defines the preprocessing procedure applied to the vector magnetograms before extrapolation.","marker":"Wiegelmann, Inhester, and Sakurai (2006)"},{"why":"Documents the approximate 3% agreement among finite-volume helicity methods, which is the basis for the 97% retention threshold.","marker":"Valori et al. (2016)"},{"why":"Introduces the relative magnetic helicity defined in terms of a reference potential field, the quantity computed throughout the paper.","marker":"Berger and Field (1984)"},{"why":"Provides the Ediv/E quality metric and its 0.05 threshold used to validate the reliability of the helicity calculations.","marker":"Thalmann et al. (2019)"},{"why":"Describes the calibration of SMFT polarization data across four decades, a prerequisite for the 150 magnetograms used here.","marker":"Su et al. (2024a)"},{"why":"Offers observational filament heights used to argue that an 81 Mm cutoff includes key eruptive structures.","marker":"Filippov (2013)"}],"fun_headline_variants":["97% of coronal helicity lies within 81 Mm, energy at 49 Mm","NLFFF height cut: 81 Mm for helicity, 49 Mm for energy, saves 38%","150 active regions: 81 Mm captures 97% of helicity and energy","Optimal extrapolation height: 81 Mm keeps 97% helicity, 38% faster"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The recommendation assumes that slicing a single 133 Mm extrapolation into nested subvolumes reproduces what a genuinely shorter extrapolation would give, even though NLFFF solutions depend on the computational box and its boundary treatment; the authors state that they do not re-extrapolate at shorter heights.","fun_headline_variants_meta":{"raw":{"variants":["97% of coronal helicity lies within 81 Mm, energy at 49 Mm","NLFFF height cut: 81 Mm for helicity, 49 Mm for energy, saves 38%","150 active regions: 81 Mm captures 97% of helicity and energy","Optimal extrapolation height: 81 Mm keeps 97% helicity, 38% faster"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00032,"raw_usage":{"total_tokens":1773,"prompt_tokens":883,"completion_tokens":890,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":499,"completion_tokens_details":{"reasoning_tokens":791}},"tokens_in":499,"tokens_out":890,"duration_ms":8310,"temperature":1.0,"reasoning_tokens":791,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T10:12:23.386577+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the extrapolations with actual box heights of 81 Mm and 49 Mm using the same preprocessing, code, and boundary conditions, and compare the resulting relative helicity and energy with the values from the 133 Mm box. If the shorter boxes systematically return less than 97% of the full-box quantities, or if the lower layers of the field change substantially when the ceiling is lowered, the proposed calibration does not apply to real extrapolations. A complementary check is to extrapolate a few cases to a taller box (e.g., 160 Mm) to confirm that the 133 Mm reference is itself converged.","supporting_citations":[],"review_version":1}