{"id":"073dbb9f-3051-4ed0-a0d6-4d4cb07a7470","arxiv_id":"2608.12457","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A 130,000-simulation climate survey of the nearby exoplanet Gl 514 b finds it is most often either globally frozen or ice-free, with polar ice caps or an ice belt in only about 1.3% of cases, and requires a CO2 pressure of 7.25-9.5 bar to become temperate.","lead":"Researchers used a fast climate model to run over 130,000 simulations of the nearby super-Earth Gl 514 b, varying its orbit, spin, and carbon dioxide levels. They found that most simulated versions are either fully frozen or ice-free, with only about 1.3% showing Earth-like polar ice caps, and that a thick CO2 atmosphere of 7.25 to 9.5 bar would be needed for a livable surface.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 7.25–9.5 bar CO2 habitability window and the 1.27% partial-ice fraction rest on the WK97 OLR scheme used at 7–11 bar without validation in that regime.","rationale":"The paper is a careful, honest parameter exploration: VPLanet is open source, the model has been compared through the FILLET intercomparison, the authors explicitly warn that the 1.27% fraction is not a probability, and they include a heat-capacity sensitivity study (Figure 11) showing the global mean temperature remains temperate across a wide range of C_W and C_L. These are real strengths. The most load-bearing concern I find is not the prior-weighting issue, which the paper itself disclaims in the Discussion, but the unvalidated use of the WK97 OLR scheme at 7–11 bar CO2. The abstract's pCO2 range and the climate-state fractions are direct outputs of that scheme, and the paper only shows OLR comparisons up to 5 bar before selecting WK97. Because the \"inferred snowball\" and \"inferred runaway greenhouse\" categories are defined by the WK97 temperature bounds, the very large aggregate fractions in Table 5 are sensitive to the OLR parameterization's boundary behavior. A line-by-line radiative transfer check at the relevant pressures would either confirm the WK97 extrapolation and leave the headline intact, or reveal a quantitative shift in the habitability window and ice-state fractions. The reader's CONDITIONAL verdict already captures the need for such validation, so I do not propose changing it; my concern is more specific than the reader's broad 'model fidelity' worry and is therefore recorded as partial agreement.","tokens_in":17154,"tokens_out":11120,"duration_ms":104288,"concrete_test":"Run a line-by-line radiative transfer calculation (e.g., LBLRTM, petitRADTRANS, or a validated narrowband code) for a CO2–N2–H2O atmosphere under Gl 514's stellar spectrum at pCO2 = 7, 8, 9, and 10 bar, with surface temperatures spanning 200–360 K, and compare the resulting OLR and top-of-atmosphere albedo against the WK97 prescription. If WK97 deviates by more than about 10 W/m2 in the 7–10 bar range, rerun the Set A, C, D, and E sweeps with a corrected OLR scheme and check whether the 7.25–9.5 bar habitability window and the 1.27% partial-ice fraction persist. If the corrected OLR shifts the window or changes the fractions by more than a few percent, the abstract's central claims are not robust to the OLR model choice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims—that pCO2 in the range 7.25–9.5 bar permits habitability and that only 1.27% of simulations yield polar ice caps or an ice belt—are generated by the Williams & Kasting (1997) OLR scheme in POISE. Section 2.1 selects WK97 because it \"behaves realistically at higher pCO2 values\" and is computationally cheap, but the manuscript provides no quantitative validation of WK97 against line-by-line radiative transfer or a 3D GCM in the 7–11 bar CO2 regime that Sets A, C, D, and E actually use. The OLR curves shown in Figure 2 are only at 280 ppm and 5 bar, and the model's temperature limits derive from a second-order polynomial fit (Appendix B of Williams & Kasting 1997). Moreover, the model's upper/lower temperature bounds define the large \"inferred snowball\" and \"inferred runaway greenhouse\" categories that dominate the aggregate statistics, so those fractions are partly artifacts of the OLR model's boundary behavior. The paper also acknowledges (Section 2.2 and Discussion) that the recalibrated model places the northern sea ice line at 47.7°N rather than the observed 72°N, indicating an ice bias. If WK97's OLR is inaccurate at 7–10 bar CO2—for example, by tens of W/m2 near the greenhouse limit—the required CO2 range for habitability and the snowball/ice-free/partial-ice fractions would shift. This is the load-bearing physical assumption behind the abstract's headline numbers.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models the climate of the super-Earth Gl 514 b, a high-eccentricity (e = 0.45) habitable-zone planet around an M dwarf, using the one-dimensional energy balance model POISE within VPLanet. The authors recalibrate six free parameters (ice/water/land albedos, water/land heat capacities, and a diffusion coefficient) against modern Earth's global mean temperature, OLR, and four latitudinal ice lines, then run 133,316 simulations across five sets (A-E) that vary pCO2, eccentricity, obliquity, precession angle, land fraction, and land distribution. They report that the majority of simulations yield ice-free or snowball states, that about 1.27% of simulations produce polar ice caps or an ice belt, and that a pCO2 range of 7.25-9.5 bar permits surface habitability. The paper concludes that Gl 514 b could be a temperate, potentially habitable planet if it possesses a multi-bar CO2 atmosphere.","tokens_in":17462,"tokens_out":5541,"duration_ms":48693,"significance":"If the quantitative results hold, the paper provides useful guidance for direct-imaging follow-up of Gl 514 b and a systematic, low-cost climate classification of a high-eccentricity habitable-zone planet. The study's strengths include the use of the open-source VPLanet package, a large and clearly described parameter sweep, a transparent calibration procedure against Earth observables, and an explicit heat-capacity sensitivity test (Figure 11). The authors also honestly acknowledge several limitations, including the model's ice bias and the non-probabilistic nature of the reported fractions. However, the headline numbers (1.27% partial-ice fraction, 7.25-9.5 bar CO2 window) rest on the Williams & Kasting (1997) OLR scheme applied at CO2 pressures far above the validation shown in the paper, so the significance is contingent on that scheme's behavior in the 7-11 bar regime.","major_comments":[{"comment":"The WK97 OLR scheme is used for simulations with pCO2 between 7 and 11 bar, but the only validation shown (Figure 2) is at 280 ppm and 5 bar. The model's temperature boundaries (-83.15 C to 86.85 C), which come from the WK97 polynomial fit, define the 'inferred snowball' and 'inferred runaway greenhouse' categories that dominate Table 5. Because the abstract's headline fractions and the CO2 habitability window are extracted from simulations that frequently hit these boundaries, the authors should demonstrate that WK97 is accurate at 7-11 bar (e.g., against line-by-line radiative transfer or a GCM) or, at minimum, show robustness of the main results to the choice among the WK97, Haqq-Misra, and Kadoya-Tajika OLR schemes.","section":"Section 2.1 and Figure 2"},{"comment":"The calibration uses 'reasonable but arbitrary' uncertainties for the ice-line observables and reproduces the northern sea ice line at 47.7 N instead of the adopted 72 N. This poor fit indicates a systematic ice bias, which the authors acknowledge in Section 4. Since one of the paper's main products is the 1.27% partial-ice fraction, the authors should assess how this calibration bias affects that fraction, for example by rerunning the classification under alternative ice-line weights or by reporting how the fraction changes if the sea-ice constraint is relaxed.","section":"Section 2.2 and Table 3"},{"comment":"The abstract's phrase 'most likely' overstates an unweighted parameter sweep. The paper explicitly states that the 1.27% fraction 'is not intended to represent the actual probability' (Section 4), and no priors are placed over pCO2, obliquity, eccentricity, or land fraction. The abstract should be reworded to reflect the sampled parameter space (e.g., 'in the sampled parameter space, the majority of simulations...') or the authors should provide a formal weighting based on observational constraints such as the eccentricity posterior from Damasso et al. (2022).","section":"Abstract and Section 4"},{"comment":"The paper never defines the quantitative criterion for 'habitable' or 'temperate' used to classify simulations and to derive the pCO2 range 7.25-9.5 bar. Please state the temperature threshold (or set of thresholds) used to designate a simulation as habitable, and clarify whether the reported range corresponds to global mean temperature, surface liquid-water conditions, or something else.","section":"Sections 2.3 and 3"}],"minor_comments":[{"comment":"The rows 'OLR Model NC79 WK97' and 'pCO2 [ppm] - 280' are confusing because NC79 does not use pCO2; clarify that the pCO2 column applies only to the WK97 calibration.","section":"Table 3"},{"comment":"The text says that 230 out of 15,555 Set A simulations finished with polar ice caps (1.48%), while Table 5 gives a total polar-cap fraction of 1.05% across all sets; please state explicitly whether 'polar ice caps' in the Set A text includes single and bipolar caps and reconcile the numbers in the text or table.","section":"Section 3, Set A results"},{"comment":"The caption mentions inferred snowball and inferred runaway greenhouse states from out-of-bounds temperatures; adding the -83.15 C and 86.85 C boundaries to the figure itself would help readers interpret the classification.","section":"Figure 4"},{"comment":"The statement that none of Forster et al. (2021), Walsh et al. (2017), de la Mare (2009), Edinburgh & Day (2016), or SCAR (2024) include uncertainties should be supported by citing the specific tables or sections checked, since at least some of these sources report ranges or spread measures.","section":"Section 2.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the authors have done a thorough parameter exploration with a standard, open-source EBM. My main concern, as reflected in the major comments, is that the abstract's quantitative claims are more confident than the underlying OLR validation and the unweighted sweep support. I would recommend requesting either validation of the WK97 scheme at 7-11 bar or a substantial softening of the headline statements, plus an explicit operational definition of habitability. The paper is not fatally flawed, but it needs these load-bearing points addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe quick take: this is an honest, thorough parameter sweep, but the two headline results—1.27% partial ice coverage and a 7.25–9.5 bar CO2 window for habitability—are softer than the abstract suggests. They come from a 1D energy balance model using the Williams & Kasting (1997) OLR scheme in a pressure regime where that scheme is not validated, and the \"most likely\" language in the abstract is not backed by any weighting of the sampled parameters.\n\nWhat is actually new: this is the first full EBM exploration of Gl 514 b covering obliquity, precession, eccentricity, land fraction, and land distribution, with 133,000 simulations. Adding the WK97 OLR option to POISE is a useful extension, and the recalibration against modern Earth geography (mean temperature, OLR, four ice lines) is a reasonable step beyond the old uniform-land calibration. The paper is candid about its own caveats, and the code is open source. The figures are clear, and the comparison with Biasiotti et al. (2024) helps the reader see what is added.\n\nWhere the soft spots are, in proportion: the stress-test concern is real. The paper shows WK97 OLR curves only at 280 ppm and 5 bar, yet the runs use 7–11 bar. There is no check against line-by-line radiative transfer or a GCM in that regime. Since the WK97 temperature bounds (below -83.15°C, above 86.85°C) define the \"inferred snowball\" and \"inferred runaway greenhouse\" categories, which together make up 73% of all runs, the aggregate statistics are partly artifacts of the OLR model's boundary behavior. The calibration also reproduces the northern sea ice line poorly (47.7°N vs the target 72°N) and yields a water-to-land heat capacity ratio of roughly 300. These are real issues; the heat capacity sensitivity test is reassuring, and the authors openly acknowledge the ice bias. The abstract's \"most likely\" overreaches, though the Discussion correctly says the 1.27% fraction is not a true probability.\n\nWho benefits: anyone planning direct-imaging observations of Gl 514 b or designing GCM follow-up. The paper is a competent, reproducible study that deserves peer review. I would send it to a good referee with a request that the high-pCO2 regime of WK97 be addressed, and that the abstract be aligned with the paper's own caveats. I would accept with major revision.","headline":"A serious, honest EBM parameter sweep whose headline numbers are softer than the abstract implies—worth refereeing, but the high-pCO2 regime and the probability framing need work.","tokens_in":18078,"tokens_out":3407,"would_cite":true,"duration_ms":27139,"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":"Climate simulations of the nearby super-Earth Gl 514 b find the surface is most likely either globally frozen or ice-free, with only about 1.27% of modeled cases showing polar ice caps or an ice belt.","keywords":["exoplanet climate","energy balance model","habitable zone","eccentric orbit","carbon dioxide atmosphere","ice coverage","Gl 514 b","M dwarf planets"],"falsifier":"A future direct-imaging observation of Gl 514 b that measures a photometric phase curve or color signature consistent with partial ice coverage (polar caps or an ice belt) in a configuration the model labels snowball or ice-free would falsify the model's climate-state distribution. A simpler near-term check is a 3D general circulation model run at the same 8.375-bar, low-obliquity parameters: reproducing polar caps would support the EBM, while yielding an ice-free or globally glaciated surface would falsify it.","tokens_in":16911,"feed_emoji":"🪐","tokens_out":5778,"duration_ms":45957,"temperature":0.7,"pith_summary":"This paper tries to establish what the climate of Gl 514 b, a nearby super-Earth on a highly eccentric orbit, is likely to be. Using a one-dimensional energy balance model, the authors explore the wide range of orbital, rotational, atmospheric, and surface properties permitted by current observations. They find the planet is most likely in one of two extreme states—a snowball or an ice-free world—and that a partial pressure of CO2 between about 7.25 and 9.5 bar would be needed to keep its surface temperate. The result matters because Gl 514 b is a prime target for future direct-imaging telescopes, and knowing whether its surface is partly ice-covered affects the photometric signal observers will see.","feed_headline":"Gl 514 b likely frozen or ice-free, not partly icy","feed_subtitle":"Climate simulations put only 1.27% of cases in polar ice caps or an ice belt.","key_machinery":"The load-bearing tool is a one-dimensional seasonal energy balance model (the POISE module of VPLanet), which solves for latitudinal temperature and albedo along the seasonal cycle and includes dynamic ice sheets that depress the lithosphere. For this study the model was recalibrated against modern Earth's geography, global mean temperature, outgoing longwave radiation, and all four latitudinal ice lines, and it adopts the Williams & Kasting (1997) outgoing longwave radiation scheme, which makes the greenhouse effect of CO2 explicit. The model is fast enough to allow a sweep of over 130,000 parameter combinations, and the recalibration provides the bridge from Earth-calibrated physics to an exoplanet with many bars of CO2.","core_discovery":"The central claim is that Gl 514 b's surface is probably not partially ice-covered under most plausible configurations. In more than 130,000 simulations spanning eccentricities out to 0.9, obliquities from 0 to 90 degrees, CO2 partial pressures from 7 to 11 bar, land fractions, and land distributions, the planet ends in either a snowball state, an ice-free state, or (in about 1.27% of cases) a state with polar ice caps or an ice belt. A CO2 partial pressure in the range 7.25–9.5 bar permits a temperate, potentially habitable surface, with a specific configuration at 8.375 bar and low obliquity producing an Earth-like polar-cap climate. The authors stress that the 1.27% figure is not a true probability but the fraction in the parameter space they explored.","pith_inferences":["The 1.27% partial-ice fraction is likely an upper-bound estimate because the recalibrated model predicts more ice than Earth actually has (northern sea ice line at 47.7N versus the observed 72N); a model biased toward ice would undercount ice-free states, so the true partial-ice probability could be even lower.","The paper's 'most likely' wording is a statement about the explored parameter grid, not a Bayesian posterior; folding in plausible observational priors—such as a preference for low obliquity from tidal damping—could shift the balance between snowball and ice-free outcomes.","A 3D general circulation model run at the 8.375-bar, low-obliquity configuration that does produce polar caps would support the EBM result, while one that yields a fully ice-free or globally glaciated surface would indicate the EBM's heat-capacity treatment overestimates ice stability.","An ice belt, though rare at 0.22% of simulations, would produce a distinctive double-peaked photometric phase curve, giving direct-imaging light curves a clean way to distinguish the rare partial-ice states from the common ice-free and snowball states."],"forward_implications":["Direct-imaging surveys of Gl 514 b should most often see either a high-albedo, fully ice-covered disk or a low-albedo ocean world, rather than a mottled partial-ice pattern.","A habitable surface on Gl 514 b would require a massive CO2 atmosphere (7.25–9.5 bar), implying strong greenhouse forcing; such thick atmospheres are physically plausible around M dwarfs.","Improved measurements of the orbital eccentricity will sharpen climate predictions, because the model puts the planet into a snowball state for eccentricities below about 0.45 when other parameters are fixed.","Land distribution is not a critical factor for habitability, but land concentrated at the poles or equator can change ice coverage patterns and thus alter the photometric signature of the planet.","The dominance of snowball and ice-free outcomes suggests Gl 514 b may occupy a bistable climate regime, where future tidal evolution of the rotation rate or obliquity could flip the planet between the two states."],"supporting_citations":[{"why":"Supplies the measured stellar and planetary properties of Gl 514 b, including eccentricity, semi-major axis, period, and minimum mass.","marker":"Damasso et al. 2022"},{"why":"Provides the outgoing longwave radiation and top-of-atmosphere albedo scheme used in POISE, which sets the CO2 greenhouse response and the temperature bounds for snowball and runaway-greenhouse classification.","marker":"Williams & Kasting 1997"},{"why":"Describes the POISE energy balance model and its original calibration, which this paper recalibrates against modern Earth geography and ice lines.","marker":"Deitrick et al. 2018b"},{"why":"Provides the dynamic ice-sheet and lithosphere depression formulation used to simulate ice growth and retreat.","marker":"Huybers & Tziperman 2008"},{"why":"Supplies the 1850–1900 northern ice-line reconstructions used as calibration targets for the model.","marker":"Walsh et al. 2017"},{"why":"Provides the ROCKE-3D version 2.0 land and ocean distribution data used to set modern Earth geography in the recalibration.","marker":"Tsigaridis et al. 2025"},{"why":"Supplies the observed outgoing longwave radiation value (239±3 W/m2) used as a calibration target.","marker":"Forster et al. 2021"},{"why":"Prior climate modeling of Gl 514 b that this paper extends by considering a fuller range of rotational, surface, and orbital parameters.","marker":"Biasiotti et al. 2024"},{"why":"Establishes the expected behavior of ice belts and polar caps at high obliquity in POISE, guiding the interpretation of Gl 514 b's ice states.","marker":"Wilhelm et al. 2022"}],"fun_headline_variants":["Gl 514 b: snowball, bare, or rarely partly icy","Only 1.27% of Gl 514 b climates yield polar ice","Sims: Gl 514 b rarely has partial ice coverage","Gl 514 b mostly frozen or ice-free, not partly icy","Gl 514 b climate: extreme states dominate, rare ice caps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that a one-dimensional energy balance model calibrated to modern Earth's ice lines, with the Williams-Kasting radiation scheme and its recalibrated heat capacities, correctly captures Gl 514 b's climate under high eccentricity and CO2 partial pressures of many bars.","fun_headline_variants_meta":{"raw":{"variants":["Gl 514 b: snowball, bare, or rarely partly icy","Only 1.27% of Gl 514 b climates yield polar ice","Sims: Gl 514 b rarely has partial ice coverage","Gl 514 b mostly frozen or ice-free, not partly icy","Gl 514 b climate: extreme states dominate, rare ice caps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000554,"raw_usage":{"total_tokens":2676,"prompt_tokens":1021,"completion_tokens":1655,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":637,"completion_tokens_details":{"reasoning_tokens":1562}},"tokens_in":637,"tokens_out":1655,"duration_ms":12602,"temperature":1.0,"reasoning_tokens":1562,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:07:52.373746+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future direct-imaging observation of Gl 514 b that measures a photometric phase curve or color signature consistent with partial ice coverage (polar caps or an ice belt) in a configuration the model labels snowball or ice-free would falsify the model's climate-state distribution. A simpler near-term check is a 3D general circulation model run at the same 8.375-bar, low-obliquity parameters: reproducing polar caps would support the EBM, while yielding an ice-free or globally glaciated surface would falsify it.","supporting_citations":[{"cited_title":"M., & Kasting, J","cited_arxiv_id":null,"evidence_quote":"Provides the outgoing longwave radiation and top-of-atmosphere albedo scheme used in POISE, which sets the CO2 greenhouse response and the temperature bounds for snowball and runaway-greenhouse classification."},{"cited_title":"2008, Paleoceanography, 23, doi: https://doi.org/10.1029/2007PA001463","cited_arxiv_id":null,"evidence_quote":"Provides the dynamic ice-sheet and lithosphere depression formulation used to simulate ice growth and retreat."},{"cited_title":"E., Fetterer, F., Stewart, J","cited_arxiv_id":null,"evidence_quote":"Supplies the 1850–1900 northern ice-line reconstructions used as calibration targets for the model."},{"cited_title":"S., Aleinov, I., et al","cited_arxiv_id":null,"evidence_quote":"Provides the ROCKE-3D version 2.0 land and ocean distribution data used to set modern Earth geography in the recalibration."},{"cited_title":"2021, in Climate Change 2021: The Physical Science Basis, ed","cited_arxiv_id":null,"evidence_quote":"Supplies the observed outgoing longwave radiation value (239±3 W/m2) used as a calibration target."},{"cited_title":"2022, The Planetary Science Journal, 3, 13, doi: 10.3847/PSJ/ac3b61","cited_arxiv_id":null,"evidence_quote":"Establishes the expected behavior of ice belts and polar caps at high obliquity in POISE, guiding the interpretation of Gl 514 b's ice states."}],"review_version":1}