{"id":"69187350-37ce-4086-b380-8e8869d00ae8","arxiv_id":"2501.11667","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Including CO2 ice albedo in a 1D energy balance model raises the starlight required for F-dwarf planets to exit a snowball state by about 29 percent relative to water-ice-only models.","lead":"This paper simulates how carbon dioxide frost on the surface of an eccentric exoplanet changes the amount of starlight needed to melt a global ice cover. The model suggests that planets around hotter F-type stars need about 29 percent more starlight to thaw when CO2 frost is included, which affects which exoplanets are judged habitable.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The cold-start CO2-to-water albedo switch is described two incompatible ways (Section 2: switch at 271/273 K; Section 3 Eq. 4: switch at 131.06 K); the 29% F-dwarf thaw-threshold claim depends on which rule was actually implemented.","rationale":"The reader's weakest assumption already flags the Section 2 vs. Section 3 contradiction about when CO2 ice albedo is replaced by water ice albedo in cold-start runs. My review agrees that this is the most load-bearing concern because the headline 29% increment is a cold-start, F-dwarf, e = 0.5 result, and the magnitude of that increment is controlled by how long the very high CO2 albedo is applied while the planet warms. The rest of the paper's qualitative conclusions, such as eccentric planets deglaciating more easily and larger CO2 grains reducing thaw thresholds, are more robust and are supported by the model's internal logic and by broad agreement with previous EBM studies. I do not see an error that would invalidate the central qualitative claim, but the quantitative headline cannot be fully trusted until the albedo-switch ambiguity is resolved. Since the reader already reached CONDITIONAL on essentially this basis, my stress-test does not move the verdict.","tokens_in":26060,"tokens_out":5026,"duration_ms":50889,"concrete_test":"Obtain or reconstruct the EBM code and rerun the Figure 8 cold-start simulation for the F2V star at e = 0.5 with 200 micron CO2 ice under both albedo-switch rules: (A) Section 2 rule, CO2 ice retained until T reaches 271 K (ocean) / 273 K (land); (B) Section 3 Eq. (4) rule, switch to water-ice albedo at T > 131.06 K. Compare the orbit-averaged flux at which global water ice cover is lost. If the reported 29% value changes by more than a few percentage points between the two runs, the headline claim is implementation-dependent and should be reported together with the exact switch rule and the code.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim, that adding the CO2 ice-albedo parameterization makes F-dwarf planets require 29% more orbit-averaged flux to thaw out of global water ice cover (Abstract; Fig. 8), rests on a cold-start model rule that the paper describes inconsistently. Section 2 states that cold-start planets are initialized everywhere with CO2 ice and that 'as the instellation increases, the albedo values of CO2 ice are decreased to that of water ice when the temperature reaches 271 K in the ocean and 273 K on land.' Section 3, Eq. (4), instead specifies a switch to water-ice albedo as soon as the zonal temperature exceeds 131.06 K, and the prose following Eq. (4) repeats the 131.06 K switch. These two rules give very different behavior: under the Section 2 rule, the very reflective CO2 albedo (0.96 for F2V, 200 micron grain, Table 2) remains in effect across the entire water-ice melting branch, strongly raising the deglaciation threshold; under the Eq. (4) rule, CO2 albedo acts only below 131 K and the planet spends most of the thaw on the water-ice branch, so the incremental 29% could be much smaller. Since no code is provided and the manuscript does not report which rule was used for Figure 8, the headline percentage is not uniquely defined by the published text. This is not a stylistic inconsistency; it changes the radiative forcing applied during the part of the simulation that produces the headline result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a modified one-dimensional energy balance model (North and Coakley 1979) that adds a wavelength-dependent CO2 ice surface albedo parameterization, using CO2 ice spectra from Hansen (1997) weighted by host-star spectral energy distributions. The model is applied to warm-start and cold-start climate hysteresis calculations for Earth-like planets around F2V, G2V, K2V, and M3V stars at eccentricities 0, 0.5, and 0.9, with a nominal 200 micron CO2 ice grain size and 400 ppmv CO2. The headline result is that cold-start F-dwarf planets require 29% more orbit-averaged flux to exit global water ice cover when CO2 ice albedo is included than with a water-ice-only albedo parameterization, with a 30% analogue for circular orbits; eccentric planets require smaller additional flux because of periastron heating. Additional results concern CO2 condensation thresholds, hysteresis widths, runaway greenhouse onset, and a strong sensitivity to CO2 grain size. The model is validated against Earth and Mars, and the discussion acknowledges several idealizations, including dust-free pure CO2 ice, a single CO2 abundance, and the absence of a dynamic OLR response.","tokens_in":26432,"tokens_out":18224,"duration_ms":188791,"significance":"The qualitative finding that CO2 ice albedo can materially affect snowball deglaciation on planets orbiting F and G stars is timely and would be a useful addition to EBM-based habitability studies. The work is not circular: the CO2 ice spectra, water ice albedos, and OLR parameterizations are externally sourced, no target threshold is fitted, and the Earth/Mars comparisons are independent. The grain-size sensitivity shown in Figure 14 is a genuinely informative result. However, the central quantitative claim is not uniquely determined by the text because the cold-start CO2-to-water albedo switch is described in two incompatible ways, and the 29% value is a point estimate without a sensitivity range around the acknowledged dust, coverage, and CO2-abundance idealizations. These issues affect the headline number rather than the basic premise, so the manuscript needs revision before the quantitative claims can be accepted.","major_comments":[{"comment":"The cold-start CO2-to-water albedo switch is described in two incompatible ways. §2 states that 'the albedo values of CO2 ice are decreased to that of water ice when the temperature reaches 271 K in the ocean and 273 K on land.' §3, Eq. (4), instead specifies CO2 ice albedo for T ≤ 131.06 K and water ice albedo for T > 131.06 K, and the prose following Eq. (4) repeats the 131.06 K switch. These rules are not equivalent: the F2V 200 micron CO2 albedo is 0.960 while the 50% snow/blue-marine water ice albedo is 0.536 (Table 2), so under the §2 rule the high CO2 albedo is applied across the entire water-ice melting branch, whereas under Eq. (4) it is active only below 131 K. The cold-start deglaciation threshold in Figure 8, which is the basis for the abstract's '29% more orbit-averaged flux' claim, depends directly on which rule was coded. Since no code is provided and the manuscript does not state the implemented rule, the headline percentage is not reproducible from the published text. Please specify the implemented algorithm, rerun the affected simulations if necessary, and restate the affected percentages.","section":"§2 and §3, Eq. (4)"},{"comment":"The 29% and 30% headline values are single point estimates. The paper's own Discussion acknowledges that Mars-like CO2 ice is dust-contaminated with albedo around 0.6, that pure dust-free CO2 ice is assumed, and that only 400 ppmv CO2 is tested; it also treats the CO2 ice surface albedo as an effective TOA albedo. Given that Figure 14 shows the deglaciation flux varying by roughly a factor of 7 between the 2 micron and 2000 micron grain sizes, the F-dwarf result is likely to be similarly sensitive to the assumed dust content, surface coverage, and CO2 abundance. Please add quantitative sensitivity tests, such as a dust-mixed CO2 albedo or partial surface coverage, or explicitly bracket the range within which the 29% and 30% values hold.","section":"Abstract, Fig. 8, §5"},{"comment":"Different sections report inconsistent numbers for the same G-dwarf, e = 0.5 cold-start grain-size calculations: §4.4 gives 340.6% (2 micron) and 46.19% (2000 micron) of the orbit-averaged solar constant for entry into the ice-free state, while §5 states 295% and 40%. Similarly, §4.3 says that e = 0.9 planets require 'considerably less' instellation for the runaway greenhouse transition than e = 0.5 planets, but the listed values are 126.18% versus 121.24%, which is more, not less. Please audit all reported percentages and make the abstract, results, and discussion mutually consistent.","section":"§4.4 and §5"}],"minor_comments":[{"comment":"Please specify the temperature units in Eq. (5); as printed, A + BT with A = 203.3 and B = 2.08 gives roughly 866 W m^-2 at 319 K if T is in kelvin, so the relationship to the 300 W m^-2 cap needs clarification.","section":"Eq. (5)"},{"comment":"The K2V and M3V rows of Table 2 contain repeated '2 [micron]' entries and misaligned water-ice albedo labels; please reformat the table so that the CO2 ice and water ice columns are unambiguous.","section":"Table 2"},{"comment":"The Mars validation reports only that the model yielded a global mean temperature of 211 K; please provide the comparison value and its uncertainty, and state specifically which aspects of the CO2 ice parameterization are constrained by this comparison.","section":"§3.1"},{"comment":"The statement that when CO2 condenses 'the entire atmosphere, including water vapor, has condensed' is imprecise: at 400 ppmv CO2, CO2 condenses at 131.06 K, while water vapor would already have condensed at much higher temperatures. Please rephrase to describe the assumed layer structure accurately.","section":"§2 and §5"},{"comment":"Please correct typographical errors such as '9GCMs' in Section 1, 'with can eccentricity' in Section 5, and the corrupted inequality symbol in Eq. (5).","section":"General"},{"comment":"The data availability statement lists spectra and data products but no model code; providing the EBM code, or at least a detailed decision-tree description of the cold-start albedo logic, would resolve the ambiguity raised in the first major comment.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's qualitative conclusions are likely salvageable, but the internal contradiction about the cold-start albedo-switch rule must be resolved before the 29% claim can be evaluated. I also recommend asking the authors to audit all numerical values in the results and discussion and to make the code or a precise algorithmic description available. The topic is within the journal's scope, and the paper is a reasonable contribution to EBM-based exoplanet habitability studies once these issues are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper adds a CO2 ice-albedo parameterization to the Palubski/Shields 1D EBM and applies it to eccentric planets around F, G, K, and M stars, varying CO2 grain size. That is genuinely new: previous CO2 condensation studies were mostly 3D GCMs for circular orbits or Mars, and the earlier EBM work did not include CO2 ice albedo. The model reproduces Earth's ice line and Mars' mean temperature, and the paper is honest about its assumptions (pure dust-free CO2, 400 ppmv, no carbon cycle, no dynamic OLR). The qualitative result—that CO2 ice albedo raises cold-start deglaciation thresholds, especially for F/G stars—is plausible and likely robust.\n\nThe problem is the cold-start albedo switch. Section 2 says CO2 ice albedo stays until water ice melts at 271/273 K. Section 3, Eq. (4), says the model switches to water ice albedo as soon as zonal temperature exceeds 131.06 K. These are different physics. Under the first rule the planet thaws with an albedo of ~0.96 for most of the branch; under the second, CO2 albedo disappears almost immediately and the planet spends most of the thaw on water ice. The headline \"29% more flux for F-dwarf planets\" is produced by one of these rules, but the paper does not say which, and no code is provided. That is not a stylistic slip; it changes the forcing that generates the central number.\n\nThe other soft spot is uncertainty quantification. Percentages like 29% and 30% are presented without error bars or sensitivity tests around the switch rule, the fixed 400 ppmv CO2, and the assumption of pure CO2 ice. The authors do flag many of these, so the issue is not neglect; it is that the model's headline numbers are not yet uniquely defined by the published text.\n\nI think this deserves a serious referee. The community working on eccentric-planet habitability will want this parameterization, and the validation gives it some credibility. But the referee should require (1) a consistent statement of the cold-start switch, ideally with code or at least a clear statement of which rule was implemented, and (2) a sensitivity test showing the 29% under both switch temperature rules and, ideally, with a few CO2 pressure values. If the 29% collapses under the 131 K rule, the abstract needs to change.\n\nWho this is for: EBM modelers and anyone estimating snowball deglaciation thresholds for eccentric exoplanets. Worth engaging seriously, but the authors need to fix the ambiguity before I trust the quantitative headline.\n\nRecommendation: send to peer review, and hold the authors to the clarification.","headline":"A useful incremental EBM extension whose headline 29% F-dwarf deglaciation number is ambiguous because the cold-start CO2-to-water albedo switch is described two incompatible ways.","tokens_in":27151,"tokens_out":4152,"would_cite":true,"duration_ms":43522,"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":"Adding CO2 ice albedo to a climate model raises the stellar flux F-dwarf planets need to deglaciate by 29%.","keywords":["exoplanet habitability","CO2 ice albedo feedback","energy balance model","eccentric planets","snowball state","stellar spectral energy distribution","ice-grain size"],"falsifier":"Run a 3D general circulation model of an eccentric F-dwarf planet at e=0.5 with 400 ppm CO2 that predicts surface CO2 frost coverage, dust contamination, and sublimation: if deglaciation requires much less than 29% extra orbit-averaged flux, the pure-ice parameterization overstates the feedback, and if it requires more, the feedback is stronger.","tokens_in":25823,"feed_emoji":"❄️","tokens_out":11992,"duration_ms":117597,"temperature":0.7,"pith_summary":"The paper is trying to establish that carbon dioxide ice condensing on the surfaces of cold or eccentric exoplanets is radiatively important for habitability, not a minor detail. It adds a CO2 ice-albedo parameterization to a one-dimensional energy balance model and finds that, compared with a traditional water-ice-only treatment, an F-dwarf planet needs 29% more orbit-averaged stellar flux to deglaciate from a water snowball, and that F-dwarf planets on circular orbits need 30% more than M-dwarf planets. The physical reason is spectral: F-dwarf light is concentrated at visible and ultraviolet wavelengths where CO2 ice is highly reflective, while M-dwarf light is infrared where CO2 ice absorbs. The paper also finds that eccentric planets escape global ice cover more easily than circular-orbit planets because of intense periastron heating, and that the size of CO2 ice grains strongly controls how much extra flux is needed. If these results hold, models that ignore CO2 ice overestimate the ease with which bright-star planets can become habitable again after global glaciation.","feed_headline":"F-dwarf planets need 29% more starlight to thaw with CO2 ice","feed_subtitle":"Adding CO2 ice albedo to climate models changes when eccentric planets thaw, reshaping habitability estimates.","key_machinery":"The load-bearing device is the CO2 ice-albedo parameterization: a temperature-threshold rule inside the energy balance model that replaces a latitude band's albedo with the SED-weighted reflectivity of pure CO2 ice, computed from laboratory spectra for grain sizes of 2, 20, 200, and 2000 μm, whenever the surface temperature falls below 131.06 K. Cold-start runs initialize the planet with CO2 ice everywhere and switch back to water ice as the planet warms. The parameterization carries the argument by raising the planet's albedo exactly in the cold state where ice-albedo feedback is strongest, which widens the climate hysteresis loop and pushes the deglaciation threshold to higher instellation, most strongly for F-, G-, and K-dwarf hosts whose stellar spectra overlap CO2 ice's high reflectivity.","core_discovery":"The paper finds that once a planet's zonal temperature falls below 131.06 K, the condensation point of 400 ppm CO2, its surface albedo should be replaced by the spectrally weighted reflectivity of pure CO2 ice. In cold-start simulations, F-dwarf planets require 29% more orbit-averaged flux to thaw out of global water ice cover than simulations using only the water ice-albedo parameterization; at zero eccentricity, F-dwarf planets require 30% more orbit-averaged flux than M-dwarf planets to exit a water snowball. The same physics makes the deglaciation threshold strongly grain-size dependent: a G-dwarf planet at eccentricity 0.5 needs 46.19% of the orbit-averaged solar constant with 2000 μm CO2 ice grains but 340.6% with 2 μm grains. Eccentric planets behave differently from circular ones: the two orders of magnitude higher flux received at periastron melts ice, so eccentric planets can leave a snowball with a smaller relative increase in instellation and sustain liquid-water or water-ice conditions over a wider range of orbit-averaged instellation.","pith_inferences":["Inference: The clean-ice numbers are likely the strongest version of the effect; real CO2 frost would carry dust and could be patchy, and lowering the albedo would shrink the 29% and 30% threshold shifts.","Inference: The same albedo parameterization implies an observable phase-curve signature: a cold eccentric planet would brighten at visible and near-ultraviolet wavelengths near apoastron as CO2 frost forms, which future reflected-light observations could test.","Inference: Grain metamorphosis over an orbit could make the deglaciation threshold time-dependent rather than fixed; fresh 2 μm frost at apoastron coarsening to larger grains would weaken the CO2 ice-albedo feedback before periastron."],"forward_implications":["Habitability studies that use only a water ice-albedo parameterization underestimate the instellation needed for F-, G-, and K-dwarf planets to deglaciate; for F-dwarf planets the shortfall is 29%.","Eccentric planets can remain in warm or partially ice-covered states across a broader range of orbit-averaged instellation because periastron heating melts ice, so the traditional circular-orbit snowball threshold may misclassify them.","CO2 ice-grain size is a first-order climate parameter: going from 2 μm to 2000 μm grains changes the cold-start deglaciation flux for a G-dwarf planet at eccentricity 0.5 from 340.6% to 46.19% of the orbit-averaged solar constant.","CO2 ice can condense on Earth-like planets at moderate eccentricity (0.5) and on low-instellation circular orbits, so the parameterization is relevant beyond highly eccentric cases.","Warm-start, ice-free planets are largely insensitive to CO2 ice-grain size because water ice controls their albedo, meaning the CO2 ice effect matters most after global glaciation has occurred."],"supporting_citations":[{"why":"It supplies the 1D energy balance model framework, Eq. 1, that is modified here.","marker":"North and Coakley (1979)"},{"why":"It supplies the CO2 ice reflectance spectra for 2 to 2000 μm grains used in the new albedo parameterization.","marker":"Hansen (1997)"},{"why":"It provides the ultraviolet-to-near-infrared CO2 ice absorption data behind the grain-size dependence of deglaciation thresholds.","marker":"Hansen (2005)"},{"why":"It provides the SED-dependent bond albedos for water ice, land, and ocean that define the baseline simulations.","marker":"Shields et al. (2013)"},{"why":"It supplies the water ice albedo spectra for snow, blue marine ice, and their mixture used for the water-ice state.","marker":"Joshi and Haberle (2012)"},{"why":"It provides the runaway-greenhouse OLR parameterization and the eccentric-planet habitability framework this study extends.","marker":"Palubski et al. (2020)"},{"why":"It provides the linearized OLR coefficients A and B used in the EBM's energy balance.","marker":"Spiegel et al. (2010)"},{"why":"It establishes CO2 condensation as a limit to deglaciation of Earth-like planets, the effect this paper parameterizes for eccentric planets.","marker":"Turbet et al. (2017)"}],"fun_headline_variants":["CO2 ice forces 29% more starlight for F-dwarf thaw","Eccentric planets exit snowballs with less instellation boost","Grain size shifts thaw threshold: 2000 μm vs 2 μm CO2 ice","F-dwarf snowball exit needs 30% more flux than M-dwarf","CO2 ice albedo reshapes eccentric planet habitability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative thresholds assume that when the surface temperature falls below 131.06 K the entire surface becomes covered with pure, dust-free CO2 ice whose reflectivity is used as the planet's overall reflectivity, and that cold-start planets keep that CO2 ice until they reach water-ice melting temperatures; patchy, dusty, or earlier-subliming frost would shrink the 29% and 30% effects.","fun_headline_variants_meta":{"raw":{"variants":["CO2 ice forces 29% more starlight for F-dwarf thaw","Eccentric planets exit snowballs with less instellation boost","Grain size shifts thaw threshold: 2000 μm vs 2 μm CO2 ice","F-dwarf snowball exit needs 30% more flux than M-dwarf","CO2 ice albedo reshapes eccentric planet habitability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00023,"raw_usage":{"total_tokens":1568,"prompt_tokens":1119,"completion_tokens":449,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":735,"completion_tokens_details":{"reasoning_tokens":350}},"tokens_in":735,"tokens_out":449,"duration_ms":4940,"temperature":1.0,"reasoning_tokens":350,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:59:28.863341+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a 3D general circulation model of an eccentric F-dwarf planet at e=0.5 with 400 ppm CO2 that predicts surface CO2 frost coverage, dust contamination, and sublimation: if deglaciation requires much less than 29% extra orbit-averaged flux, the pure-ice parameterization overstates the feedback, and if it requires more, the feedback is stronger.","supporting_citations":[],"review_version":1}