{"id":"248a453b-1d66-4aad-b1b7-1cef7e2d7f4e","arxiv_id":"1908.09315","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Many dark H-alpha fibrils around magnetic network are statistically associated with preceding rapid blue and red excursions, suggesting they are cooling aftermaths of type II spicule heating events.","lead":"Dark fibrils in the Sun's hydrogen-alpha images may be cooling trails left behind by small, rapid heating explosions called type II spicules. The paper uses statistical correlations in high-resolution images to argue these trails are common and explain the Sun's web-like chromosphere.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No null model for co-location of recurrent events leaves the causal RBE-to-fibril link unproven","rationale":"The reader's weakest_assumption identifies precisely the load-bearing concern: the delayed dark-core darkening at RBE locations could arise from recurring but independent features that share the same network patches, rather than from a causal chain initiated by the RBE heating event. My review of the full text confirms that this is the least secure link in the argument. The paper's own discussion (Sect. 4) emphasizes rapid RBE recurrence and contamination by dynamic fibrils and RREs, and the delay histograms in Fig. 14 are interpreted qualitatively without any permutation, bootstrap, or surrogate-data control. The spatial masks (Fig. 10) and resolution-degradation test (Fig. 7) are thoughtful and partially address alternative explanations, but they do not test the temporal coincidence against chance. The darkest-pixel construction in Figs. 12–14 is particularly vulnerable because it selects extreme values per pixel over long sequences; under independent recurrent activity at a fixed location, the time separation between the deepest wing and deepest core events will be influenced by the event rate and the exclusion gap, potentially mimicking a causal 2–5 min peak. The claimed fill fractions (26%/19% and the summed 43%/27%) are presented as evidence of direct association, but without a null expectation these numbers are not interpretable as excess fractions. This concern is load-bearing because the paper's headline conclusion—that type II spicules are a major agent in producing dark network fibrils—would lose its foundation if the dark-dark associations are consistent with chance co-location. I therefore agree with the reader's conditional verdict. The paper is valuable as a suggestive observational study, and the proposed null test is feasible with the existing data; until such a test is done, the causal interpretation should remain conditional. My verdict is UNCHANGED because the reader already assigned CONDITIONAL, and my analysis reinforces rather than redirects that assessment.","tokens_in":34447,"tokens_out":3013,"duration_ms":33678,"concrete_test":"Perform a temporal randomization null test on data A and B: for each network pixel, destroy the temporal link between the blue-wing and core sequences by cyclically shifting the H-alpha core time series (or the core sequence) by a random lag exceeding the RBE lifetime (e.g., 10–30 min), independently per trial, while preserving spatial structure and the marginal distribution of events. Recompute the Fig. 14 upper delay histograms and the 2–5 min fill fractions for ~1000 shifts to build the null distribution. If the observed peak height and the 26%/19% and 43%/27% fill fractions fall above the 95th percentile of the null, the causal association is supported; if they fall within the null, the co-location confounder is not excluded. A complementary control is to repeat the Fig.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central causal claim—that the 2–5 min delayed H-alpha core darkening at RBE locations is produced by the preceding RBE heating event—rests on dark-dark association statistics (Figs. 8–12, 14) that are never compared against a null model. The paper explicitly notes that RBEs recur every ~84 s at fixed network locations (Sekse et al. 2013a) and that dynamic fibrils and RREs occupy the same pixels (Sect. 4). Under a null model in which RBE occurrence and core-darkening events are independent but spatially anchored to the same network patches, the darkest-wing/darkest-core selection in Fig. 12 and the [−2,+2] min exclusion gap in Fig. 14 could produce a spurious 2–5 min peak: a pixel with several recurrent events will have its deepest wing and deepest core events separated by an interval drawn from the event-rate distribution, and the gap biases the nearest allowed separation toward +2–5 min. The 26%/19% fill fractions and the '43%/27% directly associated' claim in Sect. 5 therefore lack a demonstrated excess over chance. Without a permutation or time-shift control, the conclusion that type II spicules are a major agent in fibril production is not yet established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes two high-resolution SST/CRISP H-alpha image sequences of quiet-Sun network bordering a coronal hole, building on the earlier RBE/RRE catalogs of Rouppe van der Voort et al. (2009) and the single-event case of Rutten & Rouppe van der Voort (2017). It uses time-delay scatter diagrams, RBE-aligned stacked cutouts, darkest-pixel maps, delay histograms, and time-averaged comparisons to search for statistical associations between blue-wing RBE/RRE darkenings and dark H-alpha core or red-core features appearing a few minutes later. The authors interpret the resulting 2–5 min dark-dark associations as cooling, recombining gas outlining the tracks of intermittent type II spicule heating events, and conclude that type II spicules are a major agent in producing dark fibrils around network. The analysis is purely observational and correlation-based; no model is fitted and no machine-checked proofs or code are claimed, though the IDL tools and SHOWEX browser are referenced.","tokens_in":34612,"tokens_out":5859,"duration_ms":64825,"significance":"If the causal association is established, the paper would significantly change the standard picture of H-alpha network fibrils by identifying them as classifiable aftermaths of small-scale heating events, and it would strengthen the case that nonequilibrium hydrogen opacity, rather than static temperature structure, controls H-alpha fibril visibility. The paper has notable strengths: it reuses an independently defined RBE skeleton list, presents two complementary datasets, includes a resolution-degradation control (Fig. 7), separates network from internetwork areas (Fig. 10), and makes a specific, falsifiable prediction about the delay between wing and core darkening. However, the central claim currently rests on visual inspection of scatter contours and histograms without a statistical null model, so the significance of the result is conditional on additional controls.","major_comments":[{"comment":"The central claim that the 2–5 min delayed H-alpha core darkening at RBE locations is caused by the preceding RBE heating event is not yet established because the delay statistics are never compared with a null model. Section 4 itself notes that RBEs recur every ~84 s at fixed network locations and that dynamic fibrils and RREs occupy the same pixels; under an independence-with-co-location null, a pixel with several recurrent events will have its deepest wing and deepest core events separated by an interval drawn from the event-rate distribution, and the [−2,+2] min exclusion gap biases the nearest allowed separation toward exactly the +2–5 min range emphasized in Fig. 14. The 26%/19% fill fractions and the 43%/27% \"directly associated\" percentages in Sect. 5 therefore do not demonstrate an excess over chance. I request a permutation or time-shift control, for example randomizing event times at fixed pixels or correlating darkest core maps with darkest wing maps shifted by a large fixed offset, with the resulting null distributions reported.","section":"Section 5; Figs. 12–14"},{"comment":"The \"darkest value per pixel\" constructs are affected by a common-seeing selection that is not controlled. The authors state that selecting the darkest value sequence-wide implies selecting moments of good seeing for both the wing and the core samples; because seeing quality varies on timescales comparable to the delays, two independent darkenings at the same pixel could both be selected preferentially on good-seeing frames and produce a spurious 2–5 min association even in the absence of a physical link. The resolution-degradation test in Fig. 7 addresses instantaneous scatter spurs, but not the sequence-wide extremal selection used in Figs. 9–12. Please add a control using seeing-matched random time pairs or restrict the extremal analysis to a fixed good-seeing subset.","section":"Figs. 9, 10, 12; Sect. 3"},{"comment":"The term \"highly significant\" is used repeatedly for the dark-dark spurs, but the paper provides no statistical significance measure. The Pearson and quadrant correlation coefficients are descriptive, and the Strous-format contours are binned densities without uncertainties. Given the strong spatial and temporal autocorrelation of the data, the apparent significance of the spurs cannot be judged visually. At minimum, the authors should report bootstrap confidence intervals on the quadrant coefficients or compare the spur strengths against the same quantities computed from time-shifted image pairs.","section":"Sect. 5; Figs. 8–11"}],"minor_comments":[{"comment":"The word \"fibrilar\" appears throughout; the standard spelling is \"fibrillar.\" This should be corrected in the abstract, text, and figure captions.","section":"Abstract and throughout"},{"comment":"The y-axis label \"px/bin\" in the lower histograms is ambiguous; the caption defines the axis as a fraction of the network area, so the label should read \"pixels per bin as fraction of network area\" or similar.","section":"Fig. 14"},{"comment":"The sentence containing \"mean RBE recurrence time of of 84 s\" contains a duplicated \"of\" and should be corrected.","section":"Sect. 5"}],"recommendation":"major_revision","confidential_remarks":"The missing null model is the key issue; the authors can likely fix it with permutation or time-shift controls, and the paper's physical mechanism is plausible. I recommend requiring a quantitative null-model control before the causal language is accepted, and asking the authors to soften conclusions such as \"type II spicules are a major agent\" if the excess over chance is not demonstrated. No concerns about duplicate publication or scope: the paper is self-contained within an established series."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this if you care about what makes the H-alpha chromosphere fibrilar. The paper gives the first systematic statistical evidence that network fibrils often appear 2–5 minutes after RBE/RRE launch events, and it does so with two high-quality SST data sets, well-documented scatter tools, and a healthy sense of how confused these scenes are. The resolution-degradation test and the internetwork/network split are good ideas, and the discussion of nonequilibrium hydrogen overopacity is physically cogent. I learned from the darkest-pixel maps even if I am not fully convinced by the delay histograms.\n\nThe load-bearing claim — that the delayed dark cores are caused by the preceding RBE, not merely co-located with it — is not established. The paper never compares its delay peaks to a null model. Given that RBEs recur every ~84 s in the same network patches, the [−2,+2] min exclusion gap in Fig. 14 can produce an apparent peak just outside the gap by construction: a pixel with two independent recurrent events will have its furthest-spaced pair separated by just over two minutes. The fact that the upper histograms show a negative-delay peak as well, which the authors interpret as recurrence, is exactly the signature of co-location without causality. The 26%/19% fill fractions are not scatters around a known baseline; without a permutation or time-shift control, calling them significant is premature. The simultaneous-vs-delayed comparison and the resolution test are helpful, but they do not substitute for a control that randomizes event times.\n\nThe soft spots are real but not fatal to the paper's value. The interpretations are appropriately hedged, the authors list competing agents, and the physics section is worth reading even if the statistics are not decisive. Who should read it: anyone working on type II spicules, H-alpha diagnostics, or chromosphere modeling; the observational phenomenology is the contribution, not the final causal proof.\n\nI would send this to a serious referee. It deserves scrutiny and could be strengthened with a simple time-shift control; the analysis is reproducible enough that a referee could request it. My own verdict would be conditional acceptance pending that control.","headline":"A serious, data-rich case for the contrail-fibril idea that stops short of proof because the statistics lack any null control for recurrent co-location.","tokens_in":35216,"tokens_out":2406,"would_cite":true,"duration_ms":27237,"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":"Many dark H-alpha fibrils spreading from network regions are cooling contrails left by rapid type II spicule heating events, according to statistical time-delay analysis of the H-alpha chromosphere.","keywords":["solar chromosphere","H-alpha fibrils","type II spicules","rapid blue-shifted excursions","rapid red-shifted excursions","non-equilibrium hydrogen ionization","network fibrils","time-delay statistics"],"falsifier":"Take the darkest-pixel delay maps and shuffle the RBE occurrence times per pixel while keeping the same spatial mask; if the 2 to 5 minute positive-delay peak survives the shuffle at the same strength, the dark-dark association is not causal. A direct spatial test is also decisive: the delayed dark fibrils should lie along the continuation of the RBE track and retract toward the network root, so darkening at RBE pixels with no track-aligned extension, or appearing simultaneously with the RBE, would argue against the contrail interpretation.","tokens_in":34181,"feed_emoji":"🌞","tokens_out":10790,"duration_ms":98954,"temperature":0.7,"pith_summary":"This paper tries to establish that many of the slender dark H-alpha fibrils spreading outward from solar network regions are cooling aftermaths of small, rapid heating events, not pre-existing cool structures. Using two high-resolution H-alpha sequences, the authors show that the moments when a pixel is darkest in the blue wing, marking a rapid blue-shifted excursion (RBE) and the on-disk signature of a type II spicule, are followed 2 to 5 minutes later by the moment when that same pixel is darkest in the H-alpha line core. They interpret this statistical dark-dark association as a causal chain: the heating event ionizes hydrogen along its track, and the subsequent recombination and slow non-equilibrium settling of the $n=2$ level of hydrogen overproduce H-alpha opacity, so the cooling track appears as a dark fibril. If true, the fibrilar appearance of the chromosphere around network is largely a map of recent small-scale heating, and the paper concludes that type II spicules are a major agent in producing network dark fibrils.","feed_headline":"Dark H-alpha fibrils are cooling contrails of hot spicule launches","feed_subtitle":"Time-delay maps show dark fibrils appear minutes after rapid spicule heating, as cooling hydrogen recombines.","key_machinery":"The load-bearing mechanism is intermittent hydrogen ionization followed by delayed non-equilibrium H-alpha overopacity: a brief heating event, such as an RBE or RRE type II spicule or a smaller unresolved agent, ionizes hydrogen along a thin track, and as the gas cools the $n=2$ population controlling H-alpha extinction decays slowly because collisional detailed balance in Ly-alpha is slow at low temperature, so the track stays dark for minutes. The statistical instrument is the time-delay scatter analysis in Strous format, which compares blue-wing RBE darkening with line-core and red-core darkening at delays from 0 to 12 minutes; it is supplemented by darkest-pixel composite images and per-pixel delay histograms that isolate the 2 to 5 minute positive-delay peak and the recurrence of darkest RBE instances in 3 to 5 minute groups.","core_discovery":"The central discovery claim is that the dark H-alpha fibrils around network are, to a large extent, the cooling aftermath of type II spicules, rapid, short-lived chromospheric jets seen on the disk as rapid blue-shifted and red-shifted excursions (RBEs and RREs), and of similar but smaller unresolved heating agents. Along a spicule track, hydrogen is heated and partially or fully ionized; when the gas cools, recombination together with slow collisional settling in the Ly-alpha line leaves the $n=2$ level of hydrogen overpopulated for minutes, producing a strong non-equilibrium increase in H-alpha extinction. That opacity excess makes the cooling track appear as a dark fibril, the chromospheric analogue of an airplane contrail. The evidence is statistical: in time-delay scatter diagrams, darkest-pixel maps, and delay histograms from two high-resolution H-alpha sequences, the darkest blue-wing RBE moments at a pixel are preferentially followed within 2 to 5 minutes by the darkest line-core or red-core darkening at that pixel, and the same locations build up dark-dark association in time-averaged images. On this basis the paper concludes that type II spicules represent a major agent in the production of dark fibrils around network.","pith_inferences":["If the causal chain holds, the dark fibril itself could be used as a high-cadence detector of prior heating history, with its darkness and decay time constraining the peak temperature and ionization fraction of the spicule that produced it.","A forward-modeling test would close the loop: simulate a single RBE track with time-dependent hydrogen kinetics, synthesize the H-alpha spectrum, and check whether a dark fibril appears 2 to 5 minutes later with the observed orientation and retraction.","The same time-delay statistics applied to diagnostics without Ly-alpha-controlled opacity, such as Ca ii 8542 Å or millimeter continua, could test whether the delayed darkening is specific to H-alpha's non-equilibrium mechanism or shared by other lines.","If the contrail interpretation survives, tracking the apparent motion and decay of dark fibrils could map the full trajectories and recurrence statistics of type II spicules more completely than the short-lived blue-wing detections, potentially informing coronal mass and energy budgets."],"forward_implications":["If the claim is right, a large share of the slender dark fibrils around network are chronologically ordered traces of recent heating: a spicule launch first, a dark fibril minutes later, often retracting toward the network root.","Because RBEs live only about half a minute while their dark aftermaths persist for minutes, many more dark fibrils than spicules will be visible at any instant, resolving the apparent mismatch between sparse RBE detections and dense fibril scenes.","H-alpha fibril darkness becomes a history diagnostic rather than a temperature diagnostic: it records gas that was recently hot enough to ionize hydrogen, so modeling network fibrils requires time-dependent non-equilibrium hydrogen opacity and the local Ly-alpha radiation environment.","The chromosphere around network would have to be understood as a three-dimensional, time-dependent, intermittently heated domain rather than a slowly varying plane-parallel layer, with much of the heating occurring below current resolution.","The same delayed-opacity argument could extend to other H-alpha features such as Ellerman bombs, surges, filaments, and prominences, as the paper suggests."],"supporting_citations":[{"why":"Establishes that H-alpha rapid blue-shifted excursions are on-disk type II spicule manifestations and supplies the RBE skeleton list used as the spatial basis for the time-delay statistics.","marker":"Pub A"},{"why":"Provides the single propagating-heating-event, contrail-fibril, return-aftermath case that motivates the hypothesis, and supplied the data B H-alpha sequence.","marker":"Pub B"},{"why":"Supplies RBE lifetimes (about half a minute) and recurrence rate (every 84 seconds) used to interpret the confusion, apparent longevity, and recurrence in the delay statistics.","marker":"Sekse et al. (2013a)"},{"why":"Identifies rapid red-shifted excursions and their mutual exclusion with RBEs, used to explain asymmetries and contamination in the scatter diagrams.","marker":"Sekse et al. (2013b)"},{"why":"Provides the fundamental 1D hydrodynamic simulation showing that slow Ly-alpha-controlled decay of the hydrogen $n=2$ population produces non-equilibrium H-alpha opacity after heating.","marker":"Carlsson & Stein (2002)"},{"why":"Extends the simulation to 2D MHD and demonstrates large $n=2$ overpopulations and H-alpha overextinction in cooling wakes, the physical basis for dark contrail fibrils.","marker":"Leenaarts et al. (2007)"},{"why":"Shows that type II spicule tips reach high temperatures and often full hydrogen ionization, justifying the premise that these events ionize hydrogen along their tracks.","marker":"De Pontieu et al. (2011)"}],"fun_headline_variants":["Dark H-alpha fibrils are cooling spicule contrails","Fibrils trace cooling gas from spicule launches","New view: spicule heat creates dark fibrils","Network fibrils are recombination afterglows","Spicules leave dark H-alpha cooling tracks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The delayed 2 to 5 minute darkening at RBE locations is caused by the RBE heating event itself, rather than by independent recurrent events that happen to occupy the same network pixels.","fun_headline_variants_meta":{"raw":{"variants":["Dark H-alpha fibrils are cooling spicule contrails","Fibrils trace cooling gas from spicule launches","New view: spicule heat creates dark fibrils","Network fibrils are recombination afterglows","Spicules leave dark H-alpha cooling tracks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000443,"raw_usage":{"total_tokens":2239,"prompt_tokens":940,"completion_tokens":1299,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":1225}},"tokens_in":556,"tokens_out":1299,"duration_ms":9964,"temperature":1.0,"reasoning_tokens":1225,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:15:19.654524+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the darkest-pixel delay maps and shuffle the RBE occurrence times per pixel while keeping the same spatial mask; if the 2 to 5 minute positive-delay peak survives the shuffle at the same strength, the dark-dark association is not causal. A direct spatial test is also decisive: the delayed dark fibrils should lie along the continuation of the RBE track and retract toward the network root, so darkening at RBE pixels with no track-aligned extension, or appearing simultaneously with the RBE, would argue against the contrail interpretation.","supporting_citations":[{"cited_title":"& Stein, R","cited_arxiv_id":null,"evidence_quote":"Provides the fundamental 1D hydrodynamic simulation showing that slow Ly-alpha-controlled decay of the hydrogen $n=2$ population produces non-equilibrium H-alpha opacity after heating."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extends the simulation to 2D MHD and demonstrates large $n=2$ overpopulations and H-alpha overextinction in cooling wakes, the physical basis for dark contrail fibrils."}],"review_version":1}