{"id":"a159442a-6082-40f9-809c-f1ad78f5f207","arxiv_id":"2505.03723","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Simulated stellar flares and proton events on TRAPPIST-1e-like planets produce thermospheric cooling, mesospheric warming, ozone depletion, and up to 40 m/s wind enhancements in the middle atmosphere.","lead":"This paper uses a 3D climate model with interactive chemistry to simulate how stellar flares and energetic particles change the atmosphere of a tidally locked Earth-like exoplanet around TRAPPIST-1. It finds strong upper-atmosphere cooling, warming in the middle and lower atmosphere, and wind speed increases of up to 40 m/s at 30 to 50 km altitudes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The thermospheric-cooling claim is reported from a model whose top is ~10^-6 mbar, while the authors state that resolving flare photo-ionization/excitation needs a 300-400 km top; the cooling at 1.2e-5 mbar may be a boundary artifact.","rationale":"The reader's weakest_assumption—that the 150 km model top truncates the thermospheric chemistry and cooling relevant to the central claim—is also the most load-bearing concern I find. The claim is specifically about thermospheric cooling by NO and CO2, and the analysis level 1.223e-5 mbar is only a few scale heights below the 10^-6 mbar ceiling. The authors' own Section 4 statement that 300–400 km would be needed to resolve the relevant photo-ionization and photo-excitation processes provides direct in-text support for this concern. This is a correctness risk, not a disagreement with community consensus: WACCM's finite vertical domain is a model configuration choice, and the paper's own limitation statement concedes the issue. The middle-atmosphere warming and the 40 m/s wind enhancement are less directly affected, so the overall 'transient stellar emissions materially force climate' conclusion may survive in weakened form; hence CONDITIONAL remains the appropriate verdict rather than REJECT. I considered the outcome-based region selection for the wind-speed profile as an alternative concern, but that affects a secondary regional result and is easier to remedy by spatial multiple-testing corrections, whereas the model-top issue strikes at the thermospheric component that is showcased in the abstract. I also credit the paper for archiving data and scripts on Zenodo, which makes the proposed higher-top or vertical-profile checks feasible, but I did not treat the archive itself as verification of the thermospheric mechanism in the absence of a rerun.","tokens_in":20897,"tokens_out":6771,"duration_ms":76679,"concrete_test":"Rerun the Extreme flare case with WACCM-X, or with the model top extended to 300–400 km, using the same flare spectra, proton fluence, and planetary parameters; then compare the 8-year median temperature and NO VMR at 1.223e-5 mbar and 0.008 mbar with the present run. If the ~100–150 K thermospheric cooling and the NO enhancement do not persist when the upper boundary is moved well above the cooling layer, the reported thermospheric cooling is a model-top artifact and the headline claim should be restricted to middle-atmosphere effects.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central thermospheric result—abrupt NO/CO2 radiative cooling, with anomalies of 100–150 K at 1.223e-5 mbar (Fig. 2a; Table 2)—comes from a model whose ceiling is ~10^-6 mbar (~150 km), leaving only a few pressure scale heights of resolved atmosphere above the analysis level. In Section 4 the authors state explicitly that resolving the flare-associated photo-ionization and photo-excitation processes would require a model top of at least 300–400 km. This is not a peripheral caveat: the claimed mechanism is thermospheric, and the physical processes that set NO abundance, CO2 cooling, and thermal balance in the thermosphere are truncated at the upper boundary. The cooling maximum at the uppermost resolved levels could therefore be a numerical boundary response or a missing-source artifact rather than a robust radiative signal. No radiative-budget decomposition is presented to show that NO or CO2 emissions dominate the cooling, so the attribution is inferred from correlated NO enhancements. The middle-atmosphere warming and wind anomalies are less exposed to this objection, but the headline thermospheric claim is not supported by the configured model without a higher-top verification.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript uses the 3D Whole Atmosphere Community Climate Model (WACCM) with interactive photochemistry to simulate the response of an Earth-like, tidally locked TRAPPIST-1e-like planet to transient stellar UV flares and stellar energetic particle events. Four scenarios are compared: Quiescent, Moderate, Active, and Extreme, using flare statistics from the MUSCLES survey and proton spectra scaled from the 2003 Halloween event. The authors report 100–150 K cooling in the upper atmosphere attributed to NO and CO2 radiative emissions, 30–50 K warming in the middle and lower atmosphere associated with enhanced N2O and H2O, wind speed increases up to 40 m s−1 at 30–50 km altitude, and ozone depletion. They argue that transient stellar emissions materially alter the climate, chemistry, and circulation of synchronously rotating rocky exoplanets and affect their observability.","tokens_in":21135,"tokens_out":5778,"duration_ms":51617,"significance":"If the results hold, this would be one of the first systematic 3D chemistry-climate investigations of flare-driven atmospheric effects on temperate tidally locked exoplanets, with implications for interpreting phase curves and transmission spectra. The paper's strengths include the use of an observationally anchored stochastic flare generator, external input spectra from MUSCLES rather than fitted simulation outputs, a publicly archived dataset, and a sensible effort to apply non-parametric significance tests. However, the central thermospheric cooling claim is not yet supported by the configured model, and the statistical design as implemented does not cleanly isolate the forced response.","major_comments":[{"comment":"The caption of Figure 1 states that the Extreme and Moderate temporal medians are taken only for days when flares occur, whereas the Quiescent medians are not restricted in this way. This conditions the two groups on different subsets of the time series, so the anomaly and the Wilcoxon rank-sum test may reflect the flare-day selection rather than the atmospheric response to stellar activity. Please recompute all medians, anomalies, and significance tests on identically conditioned samples (for example, all simulated days for every scenario) or explicitly justify why the conditioning is equivalent. The same issue may propagate to the boxplot anomalies in Figure 2 and the Appendix figures if they use the same conditioning.","section":"Figure 1 caption; §3.1"},{"comment":"The headline thermospheric cooling of 100–150 K is reported at 1.223×10^-5 mbar, yet the model top is about 10^-6 mbar (~150 km), leaving only about one pressure scale height of resolved atmosphere above the analysis level. The manuscript itself states in Section 4 that resolving flare-associated photo-ionization and photo-excitation processes would require a model top of 300–400 km. With the cooling maximum so close to the upper boundary, the signal could be a boundary artifact or a missing-source effect rather than a robust radiative response. In addition, the abstract attributes the cooling to radiative emissions from NO and CO2, but the results show NO enhancements only and no CO2 diagnostic or radiative-budget decomposition. Please provide a higher-top verification or a radiative flux decomposition, or reframe the claim to the middle and lower thermosphere and soften the 'abrupt thermospheric cooling' attribution.","section":"§3.1, §4, Fig. 2a, Table 2"},{"comment":"Each of the four scenarios is a single 8-year realization, but the Wilcoxon rank-sum and Kruskal-Wallis tests are applied to daily output across many grid points as if these were independent samples. Atmospheric state variables are strongly autocorrelated in space and time, so the effective sample size is far smaller than the number of grid-cell-days, and the reported p<0.01 overstates the confidence that the scenarios differ. Please report the effective sample size or ensemble spread, and phrase the significance claims accordingly.","section":"§2.2, §3.2"},{"comment":"The altitude nomenclature is internally inconsistent: Table 2 labels 1.223×10^-5 mbar as the lower thermosphere, while Figure 1a and Figure 2 call the same level the middle atmosphere; Table 2 labels 0.008 mbar as the mesosphere, while Figure 3 calls it the middle atmosphere. Since the paper's central claim is about thermospheric cooling, these labels must be harmonized so readers can tell which results are thermospheric, mesospheric, or stratospheric. This is not merely cosmetic; it determines whether the abstract's 'thermospheric cooling' claim is actually tested.","section":"Table 2 vs. Figures 1–3"}],"minor_comments":[{"comment":"The caption says the panels show NO and OH, but the axes in panels (e) and (f) are labeled N2O; please correct the caption to match the displayed quantities.","section":"Figure 2 caption"},{"comment":"The Flare Rate row in Table 1 is listed as 3–6 day^-1, while the text specifies 6 day^-1 for Moderate and Extreme and 4 day^-1 for Active; please unify these numbers.","section":"Table 1 vs. §2.1"},{"comment":"The phrase 'after the ≈63 yr mark' should refer to months, not years, given the 8-year simulation and the month axis in Figure 3.","section":"§3.3"},{"comment":"There is a typo '10^-6 mba' that should be '10^-6 mbar'.","section":"§4"},{"comment":"The Data Availability statement says 'scrips'; this should be 'scripts'.","section":"§5"},{"comment":"There is a typo 'distributuon' in the first paragraph; it should be 'distribution'.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the data availability is a genuine strength. The central finding would be more robust if the middle-atmosphere results were presented as the headline, with the thermospheric cooling explicitly contingent on a higher model top and a radiative-budget verification. The statistical conditioning issue in Figure 1 and the single-realization design should be fixed before publication. I would not reject the paper, but the load-bearing thermospheric claim needs substantial additional support or reframing."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is the first 3D chemistry-climate model study of how transient stellar emissions (UV flares and SEPs) actually drive climate and circulation on a tidally locked Earth-like planet. The middle-atmosphere findings—30–50 K warming from N2O and H2O, wind anomalies near 40 m/s at 30–50 km, ozone-mediated variability in moderate-activity regimes—are new, plausible, and a genuine step past the 1D photochemistry canon. The paper also does the field a service by naming its own limits: model-top height, synthetic flare spectra, a single proton spectrum, and daily flare cadence are all flagged in the discussion. Data and scripts are on Zenodo, which is reproducible and welcome.\n\nNow the soft spots, in order of how much they matter. The headline thermospheric cooling (100–150 K at ~1e-5 mbar via NO/CO2 emission) is the weakest result, and the stress-test is right to point at the model top. WACCM here reaches ~10^-6 mbar (~150 km); the authors themselves say resolving flare photoionization and photoexcitation would need 300–400 km. With only a few scale heights above the analysis level, that cooling could be a boundary artifact, and no radiative budget decomposition is shown to prove NO/CO2 actually dominate the balance. I would not stake a claim on it without a higher-top run, even a 1D one.\n\nTwo smaller issues. The Figure 1 caption says the Extreme and Moderate medians are taken only on days when flares occur, while the Quiescent median is not restricted that way—this biases the anomalies upward. And each scenario is a single stochastic realization; the Wilcoxon tests show within-run separation, not run-to-run robustness. The 40 m/s wind profile region was also chosen post hoc for showing the largest anomaly, so I'd read that number as an upper bound. None of this is fatal; the paper is an honest pilot study and reads like one.\n\nThe citation pattern and external inputs look clean: flare statistics from MUSCLES, proton fluence from the literature, nothing fitted to output.\n\nWho gets value: M-dwarf atmosphere modelers, and observers planning transmission or direct-imaging campaigns who need to know whether flare-driven weather can corrupt retrievals. I'd send it to review, with the model-top issue and the flare-day-only medians as required major revisions. The middle-atmosphere story should survive; the thermospheric claim needs either a higher-top test or a downgrade to tentative.","headline":"First 3D chemistry-climate study of how flares and stellar protons affect climates of tidally locked exo-Earths; the middle-atmosphere results are the real news, but the headline thermospheric cooling needs a higher model top before I'd trust it.","tokens_in":734,"tokens_out":1102,"would_cite":true,"duration_ms":52785,"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":"Transient stellar flares and particle storms materially alter the climate and circulation of tidally locked Earth-like planets around M dwarfs, cooling the thermosphere through NO and CO2 emission while warming the middle and lower…","keywords":["stellar flares","stellar energetic particles","tidally locked exoplanets","M-dwarf habitability","3D chemistry-climate modeling","thermospheric cooling","exoplanet atmospheric dynamics","ozone variability"],"falsifier":"Run the same flare scenarios with the model top raised to 300 to 400 km and explicit ion chemistry: if the NO/CO2 radiative cooling no longer produces the roughly 100 K thermospheric anomaly, that claim is falsified. A complementary test is observational: monitor a known flaring M-dwarf planet through repeated transits and look for flare-correlated N2O and H2O enhancements and O3 depletions; their absence would falsify the chemical-climate coupling.","tokens_in":20691,"feed_emoji":"🌍","tokens_out":10612,"duration_ms":98661,"temperature":0.7,"pith_summary":"Using a three-dimensional global chemistry-climate model, this paper asks whether transient stellar emissions—ultraviolet flares and stellar energetic particles—do more than alter photochemistry on tidally locked Earth-like exoplanets. The authors simulate a TRAPPIST-1e-like planet and find that recurring flares cool the upper atmosphere by roughly 100 K through radiative emission by NO and CO2, while warming the mesosphere and lower atmosphere as N2O and water vapor accumulate. They also report that intense flares accelerate nightside winds by up to 40 m/s at 30 to 50 km altitude and can push water vapor toward moist-greenhouse abundances. If right, flaring activity is not a side note but a primary driver of climate, chemistry, and observable variability for planets around active low-mass stars.","feed_headline":"Flares cool an exo-Earth's upper air and heat its lower air","feed_subtitle":"A 3-D chemistry-climate model shows NO and CO2 cooling the thermosphere while N2O and H2O warm the layers below.","key_machinery":"The central object is a three-dimensional general circulation model with interactive photochemistry applied to a tidally locked, Earth-composition exo-Earth forced by stochastic flare spectra and proton precipitation. The mechanism that carries the result is the coupling of proton-driven odd-nitrogen and odd-hydrogen production with radiative transfer: NO and CO2 become thermospheric coolants, N2O and H2O become lower-atmosphere warmers, and ozone depletion or recovery modulates the temperature response on monthly to annual timescales.","core_discovery":"On its own terms, the paper's central claim is that time-dependent stellar emissions produce a vertically split thermal response in a synchronously rotating, Earth-composition planet: stellar protons and UV drive odd-nitrogen and odd-hydrogen chemistry that radiatively cools the thermosphere (up to roughly 100 to 150 K in the most extreme runs), while the same events raise middle- and lower-atmosphere temperatures by tens of kelvin through increased infrared absorbers N2O and H2O. Ozone acts as a modulating knob: under moderate flaring it can increase, but under extreme flaring it is eroded, and the temperature response tracks this variability. The simulations also show a dynamical consequence—horizontal winds on the nightside at 30 to 50 km altitude increasing by up to about 40 m/s—and water vapor occasionally reaching a mixing ratio above $10^{-3}$, connecting transient stellar events to both long-term habitability and short-term observability.","pith_inferences":["The model's 150 km top likely truncates the very thermospheric chemistry responsible for the cooling signal; extending the top to 300 to 400 km with explicit ion chemistry could make the reported NO/CO2 cooling stronger or shift its altitude, while the middle-atmosphere warming should be more robust.","Because the authors use a single 9000 K blackbody flare spectrum and note that real M-dwarf flares are harder in the FUV and NUV, the balance between ozone production and destruction—and hence the temperature response—could shift with spectral shape; rerunning with observationally anchored flare spectra is a direct test.","On Earth, geomagnetic shielding moderates energetic particle effects; if tidally locked planets have weak dynamos, protons could penetrate deeper and amplify the lower-atmosphere greenhouse response this paper reports.","An observational discriminator: a flaring M-dwarf planet should show correlated enhancements of N2O and H2O and depletion of O3 in spectra taken weeks after large flares; the absence of such a correlation would challenge the modeled coupling."],"forward_implications":["Climate simulations of M-dwarf habitable-zone planets that omit transient stellar emissions will misrepresent both the thermal structure and the circulation, since the paper finds persistent upper-atmosphere cooling and middle-atmosphere warming that depend on flare frequency and energy.","Observed spectra will be time-variable: repeated flares change N2O, H2O, and O3 abundances, so transmission and emission observations averaged over many transits must account for flare-modulated composition rather than a single quiescent state.","Habitability assessments should include flare statistics: moderate flaring may raise greenhouse gases while extreme flaring erodes ozone, so surface UV dose and surface temperature depend on the flare frequency distribution rather than mean stellar XUV alone.","The reported wind-speed enhancements of up to 40 m/s at 30 to 50 km altitude could alter day-night chemical transport and the dynamical regime of the middle atmosphere, with consequences for how species are mixed and observed.","Extreme flare scenarios can push water vapor into the moist-greenhouse regime (mixing ratio above $10^{-3}$), implying that actively flaring hosts may drive their planets toward faster water loss than quiescent hosts."],"supporting_citations":[{"why":"Establishes the exoplanet application of the chemistry-climate model used here, including the Earth-like surface and land-mask assumptions.","marker":"Chen et al. 2019"},{"why":"Prior modeling showing ozone depletion and greenhouse-gas enhancement under stellar proton events; provides the photochemical baseline this study extends to three dimensions.","marker":"Chen et al. 2021"},{"why":"Reference photochemistry result for a single large flare's ozone impact on an Earth-like planet, which this work generalizes to recurrent three-dimensional forcing.","marker":"Segura et al. 2010"},{"why":"One-dimensional simulations of repeated flares eroding ozone motivate the need for a global climate-chemistry treatment.","marker":"Tilley et al. 2019"},{"why":"Proposal that stellar proton events enhance greenhouse gases such as N2O, the mechanism this paper tests in a three-dimensional climate model.","marker":"Airapetian et al. 2016"},{"why":"Stochastic flare generator used to produce the ultraviolet flare spectra and frequencies driving the simulations.","marker":"Loyd et al. 2018"},{"why":"Ultraviolet survey of low-mass stars whose spectra and variability underpin the flare input model.","marker":"France et al. 2016"},{"why":"Solar proton event parameterization connecting proton flux to atmospheric ion-pair production, used to inject energetic particle effects.","marker":"Jackman et al. 2008"},{"why":"Prescription linking a flare's SiIV energy to peak proton fluence, used to set the stellar energetic particle inputs.","marker":"Youngblood et al. 2017"}],"fun_headline_variants":["Flares chill exo-Earth upper air, warm lower air","Stellar flares cool exo-Earth thermosphere, warm lower atmosphere","Flares split exo-Earth climate: cooling aloft, heating below","Exo-Earth thermal split: flares chill top, warm bottom","Flares drive opposing temperature shifts in exo-Earth atmospheres"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reported thermospheric cooling depends on the model resolving NO and CO2 radiative emission with an atmospheric top near 150 km, but the authors note that flare photo-ionization and photo-excitation require a top at 300 to 400 km; if the truncation removes the relevant thermospheric chemistry, that part of the claim collapses.","fun_headline_variants_meta":{"raw":{"variants":["Flares chill exo-Earth upper air, warm lower air","Stellar flares cool exo-Earth thermosphere, warm lower atmosphere","Flares split exo-Earth climate: cooling aloft, heating below","Exo-Earth thermal split: flares chill top, warm bottom","Flares drive opposing temperature shifts in exo-Earth atmospheres"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001002,"raw_usage":{"total_tokens":4277,"prompt_tokens":1022,"completion_tokens":3255,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":3161}},"tokens_in":638,"tokens_out":3255,"duration_ms":23285,"temperature":1.0,"reasoning_tokens":3161,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:45:12.472824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same flare scenarios with the model top raised to 300 to 400 km and explicit ion chemistry: if the NO/CO2 radiative cooling no longer produces the roughly 100 K thermospheric anomaly, that claim is falsified. A complementary test is observational: monitor a known flaring M-dwarf planet through repeated transits and look for flare-correlated N2O and H2O enhancements and O3 depletions; their absence would falsify the chemical-climate coupling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Ultraviolet survey of low-mass stars whose spectra and variability underpin the flare input model."},{"cited_title":"2008, Atmospheric Chemistry and Physics, 8, 765","cited_arxiv_id":null,"evidence_quote":"Solar proton event parameterization connecting proton flux to atmospheric ion-pair production, used to inject energetic particle effects."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prescription linking a flare's SiIV energy to peak proton fluence, used to set the stellar energetic particle inputs."}],"review_version":1}