{"id":"8212515e-fbb0-46c2-950f-60983c49be22","arxiv_id":"2412.19357","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For Earth-like planets with Earth's annual sunlight, rotation period controls habitable land area; habitability falls from about 70% to 20% when day length exceeds about 20 days.","lead":"This paper simulates 93 Earth-like planets with different spins, tilts, and orbits, and trains a statistical model to predict which ones keep land warm and wet enough for life. The main finding is that day length dominates: planets rotating once per 20 days or slower lose most of their habitable land.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Monthly-average habitability metric may artifactually create or shift the ~20-day drop; diurnal cycle unresolved in the critical rotation regime.","rationale":"The reader identified monthly averaging as the weakest assumption, and I agree that it is the most load-bearing concern. The paper's headline claim is a quantitative statement about the location and magnitude of a habitability drop, and that drop is computed from a metric whose temporal averaging is inconsistent across the rotation-period range being compared: diurnal cycles are resolved for Prot ≥ 64 days but averaged over for Prot < ~30 days. The acknowledged limitation in §6.1 directly applies to the transition region, so the '~20 days' boundary and the '~70% to ~20%' numbers could be artifacts. This is more fundamental than the sparse eccentricity sampling, which is a comparative claim that the test set partially validates, and more direct than the lack of internal-variability propagation, which affects error bars rather than the central tendency. I considered whether the biased test residuals (11 of 14 outside 1σ underestimates at Prot < 64 days) indicate emulator miscalibration, but the magnitudes are modest and the overall pattern is roughly consistent with the stated uncertainties; that would not change the qualitative conclusion. The proposed re-run with higher-cadence output directly tests whether the monthly metric drives the drop; if the shift is small, the paper's claim stands, and if large, the claim needs to be reframed around the chosen growing-season definition. Because the authors already disclose the limitation and the verdict is CONDITIONAL, my assessment does not move the verdict; it sharpens the condition that should be satisfied before the quantitative boundary is used in follow-up work.","tokens_in":33395,"tokens_out":7976,"duration_ms":74462,"concrete_test":"Re-run a subsample of ~10 ROCKE-3D simulations spanning the critical region (Prot ≈ 8, 16, 24, 32, 48, 64 days; obliquity 0–60°; e = 0–0.2) with 6-hourly or daily output for the final 30 years. Recompute fH using the same 0–100°C and 300 mm/yr thresholds but with temporal averaging at 6-hourly, daily, and monthly scales. Determine the rotation period at which the emulated fH falls below ~0.4 and the plateau values for each averaging choice. If the boundary shifts by more than a factor of ~2 (e.g., from ~20 days to ≥40 days) or the plateau changes by more than 15%, the '~20 days' drop is an artifact of the monthly averaging choice, and the central claim requires re-qualification.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—the sharp drop in fH between ~20 and ~32 days—rests on a habitability metric computed from monthly mean temperatures (Equation 2). For rotation periods shorter than the monthly output interval (~30 days), the diurnal temperature cycle is averaged away, so the metric cannot distinguish between a cell whose monthly mean is below 0°C and one that is above freezing for a substantial fraction of the month. The authors acknowledge in §6.1 that the diurnal cycle is only resolved for Prot ≥ 64 days and that faster rotators would give different fH values if shorter averaging were used. The claimed transition at ~20 days lies exactly in this unresolved regime. For a planet with Prot = 20–30 days, a land cell can experience daytime temperatures above 0°C for ~10–15 consecutive days while its monthly mean remains below 0°C; the metric counts such a cell as uninhabitable, even though the warm interval may exceed the 'growing season' the authors invoke. The defense in §6.1 (that a few days is too short for a biosphere) is not quantified and does not address warm intervals of a week or more. If fH were recomputed from daily or 6-hourly outputs, the slow-rotation values could rise substantially, shifting the boundary to longer rotation periods and undermining the headline statement that rotation period is the overarching primary driver.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 93 ROCKE-3D general circulation model simulations of Earth-like planets, varying rotation period, obliquity, orbital eccentricity, and longitude of periastron via Latin Hypercube Sampling. A Gaussian-process emulator is trained on 46 runs and tested on an independent 46-run sample to map a combined temperature-precipitation climate habitability metric fH defined over land from monthly model outputs. The central claim is that, for eccentricities up to 0.225, rotation period is the primary driver of fH, with a sharp drop from ~70% to ~20% for rotation periods beyond ~20 days, and that obliquity is a secondary factor for fast rotators. The paper also compares its metric against previous work by Jansen et al. (2019) and He et al. (2022) and discusses the physical mechanisms behind the rotation-period dependence.","tokens_in":33653,"tokens_out":12605,"duration_ms":106295,"significance":"If the central claim holds, the paper provides a practically important result for exoplanet habitability: a single spin parameter, rather than orbital shape, may set the climate habitability of Earth-like planets in the habitable zone, with direct implications for target selection and follow-up observations. The methodological framework—LHS sampling of a four-dimensional parameter space, GP emulation with a held-out test set, and public code and data—is a strength and makes the study reproducible. The paper also gives a careful comparison to prior habitability metrics and honestly acknowledges several limitations, including the unresolved diurnal cycle. However, the headline quantitative finding (the ~20-day drop and its magnitude) is sensitive to the monthly averaging used to define fH, a sensitivity the authors themselves flag in Section 6.1; this needs to be addressed before the central claim can be taken at face value.","major_comments":[{"comment":"The central claim of a sharp drop in habitability between ~20 and ~32 days is computed from monthly-mean temperatures, but the authors concede in §6.1 that the diurnal cycle is only resolved for rotation periods of 64 days and above, and that shorter averaging would change fH. For Prot = 20–30 days, the sunlit interval lasts ~10–15 consecutive days, which is not the 'few days' that the §6.1 growing-season argument dismisses; the monthly-mean indicator in Eq. (2) counts a cell as uninhabitable for an entire month whose mean temperature is below 0°C even if the cell is above freezing for a substantial fraction of that month. The claimed location and magnitude of the drop are therefore not robust to the averaging interval, and the abstract's '~20 days' and '~70% to ~20%' are overstated as physical findings. I recommend either recomputing fH from daily or 6-hourly outputs for a subset of runs spanning Prot = 16–64 days to test robustness, or explicitly qualifying the abstract and conclusions as applying to the monthly-mean metric.","section":"§6.1, Eq. (2)"},{"comment":"Table 1 does not list rotation periods for the test set, and the text in §5.4 gives test Case 16 Prot = 4.59 days and test Case 27 Prot = 21.1 days, which are identical to the training-set Prot values for the same case numbers. As printed, this makes the independence of the test sample in the primary driver dimension unverifiable, and if the test set reuses the training Prot values, the held-out validation in §5.4 would not test interpolation in the rotation-period direction. Please add the missing test Prot column or otherwise clarify how the two LHS samples are independent in this dimension.","section":"Table 1, §5.4"},{"comment":"Of the 14 test points that deviate from the emulator prediction by more than 1σ, 11 are underestimates and all occur at rotation periods shorter than 64 days. The authors state that the 14/46 count is consistent with Gaussian expectations, but they do not test the sign asymmetry; a two-tailed binomial test on 11 of 14 deviations being one-signed has p ≈ 0.06, which is marginal. If the emulator systematically underpredicts fH in the fast-rotation regime, the contrast between fast and slow rotators—and hence the quantitative '70% to 20%' drop—could be exaggerated. Please discuss this asymmetry and, if possible, test for bias.","section":"§5.4, Fig. 13"}],"minor_comments":[{"comment":"The precipitation indicator Iprec assigns a 1 to a cell for the entire year if the annual total exceeds 300 mm, regardless of whether the precipitation falls in a single month; the paper notes this in §4.2, but the text could more explicitly state that fprec is insensitive to seasonal concentration of rainfall.","section":"§3.2, Eq. (2)"},{"comment":"The caption states that the global fractional averages agree closely with Jansen et al. (2019) up to 32 days; the plotted symbols for the different studies are not always easy to distinguish in gray-scale, so different marker shapes would improve readability.","section":"Fig. 1"},{"comment":"The sentence noting that 'the faster rotation periods would also have different fractional habitabilities if the calculation was made from shorter averaging periods' is an understatement that could be quantified; a simple sensitivity estimate for Prot = 20–30 days would help the reader judge the magnitude of the effect.","section":"§6.1"},{"comment":"The RBF length scales for e cos φp are reported as ≳10^3 for fT and H, which effectively means the emulator is flat in that dimension; the table could note that such large values indicate the dimension is uninformative, reinforcing the paper's claim of a weak eccentricity dependence.","section":"Table 2"},{"comment":"The phrase 'for rotation periods greater than ~20 days, habitability drops significantly' should be qualified with 'for the monthly-mean metric used here' or similar, given the acknowledged dependence of fH on the averaging interval.","section":"Abstract"},{"comment":"The emulator's predictive variance is zero at training points because the GCM is treated as deterministic; the paper states this, but it is worth reminding the reader that the reported uncertainties exclude GCM internal variability, which is partially validated by the test-set residuals but not propagated into the headline fH values.","section":"§3.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and presents a substantial, reproducible simulation campaign. The main concern is the temporal-averaging issue: the central ~20-day drop may be an artifact of using monthly means in a regime where the diurnal cycle is unresolved. The authors' own §6.1 admission makes this load-bearing, so I recommend revision rather than rejection. Please also check the original Table 1: the test-set rotation periods appear missing in the manuscript version I reviewed, and if they are in fact identical to the training-set values, the validation design would need serious rethinking."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Rotation-period dominance is real in these models, and the eccentricity null result is the genuinely new piece. The held-out LHS test is the best part: 46 training runs, 46 independent test runs, residuals consistent with the emulator's error bars. That is how a parameter survey should be done. The comparison with Jansen et al. and He et al. is careful, and the emulator code is archived.\n\nThe soft spot is temporal resolution. Everything is computed from monthly means, and §6.1 admits the diurnal cycle is only resolved for Prot ≥ 64 days. The claimed ~20-day drop sits exactly in the unresolved regime. The stress-test concern lands: for Prot = 20–30 days, a land cell can be above freezing for a week or more while its monthly mean stays below 0°C, and the metric counts it uninhabitable. The §6.1 defense — a few days is not a growing season — does not address warm intervals of a week or more, which is the relevant length for those rotation periods. So the quantitative boundary is metric-dependent. The qualitative message that rotation is primary would likely survive a daily-resolution check, but the paper should test it or quantify the growing-season timescale before asserting the sharp transition.\n\nTwo smaller issues: the GCM outputs themselves are not archived, only the emulator code, and internal variability is not propagated into fH, so the emulator uncertainties are interpolation-only. Neither is fatal.\n\nBottom line: a solid, useful parameter study for HWO target prioritization, honestly written and with a real held-out test. It deserves peer review. A good referee will ask for a daily or 6-hourly sensitivity run across the 20–30 day rotation range and for archived GCM outputs. I would not quote the 20-day threshold without the metric caveat.","headline":"A solid but metric-sensitive map of 4-D habitability: rotation dominates, eccentricity doesn't matter much, but the ~20-day drop is an artifact of monthly averaging that needs a sensitivity check.","tokens_in":34212,"tokens_out":3791,"would_cite":true,"duration_ms":35172,"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":"For Earth-like planets, one number — rotation period — appears to control climate habitability, with the fraction of warm, wet land collapsing from roughly 70% to 20% once days grow longer than about 20 Earth days.","keywords":["habitability metric","rotation period","obliquity","orbital eccentricity","longitude of periastron","exoplanet climate","general circulation models","Gaussian process emulation"],"falsifier":"Take the existing 64-day and 128-day rotation runs, the models where the diurnal cycle is actually resolved, and recompute the temperature metric against daily output instead of monthly means. If daily-thresholded f_T rises above the monthly value by more than the run-to-run scatter, then the sharp habitability collapse near 20 days is inflated by temporal averaging, and the cliff's true location and steepness would need re-measuring.","tokens_in":33201,"feed_emoji":"🪐","tokens_out":9803,"duration_ms":82643,"temperature":0.7,"pith_summary":"Earth-like planets in the habitable zone can have wildly different climates depending on how fast they spin and how the spin is oriented, and this paper tries to map out which configurations keep land warm and wet enough for life as we know it. The authors ran 93 three-dimensional climate models of Earth analogs spanning rotation period, obliquity, orbital eccentricity, and longitude of periastron, then used statistical emulation to turn those scattered runs into a smooth prediction of habitability across all four parameters. Their central claim is that rotation period dominates: for days longer than about 20 Earth days, the fraction of land that stays above freezing with adequate rain drops from roughly 70% to 20%, driven by cooler land temperatures. Obliquity matters only as a secondary lever on fast rotators, and eccentricity up to 0.225 barely registers. If correct, this singles out one measurable quantity — planet spin — as the key to predicting which Earth analogs are worth observing.","feed_headline":"Day length, not orbit shape, sets Earth-analog habitability","feed_subtitle":"93 climate models show warm, wet land shrinking from ~70% to ~20% once days stretch past 20.","key_machinery":"The load-bearing object is the fractional climate habitability metric f_HZ (denoted H in the body), defined as the land-area-weighted fraction of grid cells whose monthly mean surface temperature lies between 0 and 100 °C and whose yearly precipitation reaches 300 mm, with the temperature indicator applied before time averaging so that slow-rotator diurnal cycles count only where they persist for roughly a month. Around that metric the paper builds a four-dimensional parameter study: Latin Hypercube Sampling distributes 92 non-Earth GCM runs evenly across rotation period, obliquity, eccentricity, and longitude of periastron, and a Gaussian process regressor with a radial-basis-function kernel interpolates habitability between the runs, with the eccentricity pair recast as (e cos φ_p, e sin φ_p) so the kernel sees a Cartesian space. The emulator is validated against a held-out test set of 46 runs, and the rotation-period gradient it recovers, steepest between 16 and 32 days, is the mechanism that carries the paper's main claim.","core_discovery":"The paper's claim, stated on its own terms, is this: for planets with Earth's surface conditions and annual insolation, the fraction of land that is both above freezing and wet enough for life is set overwhelmingly by rotation period, over the full explored ranges of obliquity (0–90°), eccentricity (up to 0.225), and longitude of periastron. The emulated habitability map shows a cliff near 20 days of rotation: at shorter periods roughly 60–80% of land meets the temperature-and-precipitation criterion, while no model with a rotation period longer than 32 days shows average habitability above 0.25, no matter how the other three parameters are set. The same map assigns obliquity a factor-of-two influence, largest at intermediate obliquity, but only for rotation periods under about 20 days, and assigns eccentricity an effect too weak for the emulator to attribute with confidence.","pith_inferences":["If habitability were judged on a growing-season or diurnal timescale instead of monthly means, the ~20-day cliff would likely soften and shift toward slower rotation: the monthly metric writes off short warm spells by construction, and the paper only resolves the diurnal cycle at rotation periods of 64 days and beyond.","The weak eccentricity signal is at least partly a consequence of averaging: the metric integrates over full orbits and over all land, so seasonally extreme but locally habitable intervals are washed out, and the emulator's one hint of an eccentricity dependence rests on a single training point.","Holding the same sampling-and-emulation pipeline to parameters this paper fixes — ocean fraction, atmospheric CO2, continental layout — would show whether the rotation-dominated result transfers beyond the single Earth-analog configuration modeled here."],"forward_implications":["Mission target selection should treat rotation period as a primary filter: a planet with a ~10-day spin and low-to-intermediate obliquity is far more likely to expose temperate, rainy land than a ~30-day rotator on an otherwise identical orbit.","Time-series observations that constrain rotation would directly test the predicted habitability break, since current transit and radial-velocity surveys constrain orbital elements but not spin.","The habitability cliff is a land-temperature effect, not a rainfall effect: oceans retain near-100% habitable fractions out to 128-day rotations, so ocean-dominated worlds would stay habitable by this metric even where continents freeze.","The training-and-test emulator design makes four-dimensional parameter sweeps of expensive climate models affordable, offering a template for extending GCM-based habitability studies beyond the rotation-and-orbit parameters explored here."],"supporting_citations":[{"why":"Supplies the temperature-only fractional habitability metric and the rotation-period peak near 16 days that this paper's metric and emulator are benchmarked against.","marker":"Jansen et al. (2019)"},{"why":"Adds the obliquity dimension and the temperature-plus-precipitation habitability condition that H adopts in modified form.","marker":"He et al. (2022)"},{"why":"Introduced the fractional habitability concept — the fraction of planetary area warm enough for liquid surface water, integrated over an orbit.","marker":"Spiegel et al. (2008)"},{"why":"Is the ROCKE-3D general circulation model paper, the code used for all 93 simulations.","marker":"Way et al. (2017)"},{"why":"Provides the Gaussian process regression formalism the emulator is built on.","marker":"Rasmussen & Williams (2006)"},{"why":"Frames emulation as a way to model expensive simulators with built-in uncertainty, the approach the paper follows.","marker":"O'Hagan (2006)"},{"why":"Supplies the terrestrial precipitation dataset used to justify the 300 mm/yr cutoff in the precipitation metric.","marker":"Willmott & Matsuura (2018)"}],"fun_headline_variants":["Rotation, not orbit, controls Earth-like habitability","Day length rules: 70% to 20% land habitability drop","Climate cliff: 20-day rotation kills Earth-analog life zones","Spin dominates: obliquity matters only for short days","93 GCMs: rotation period sets habitable land fraction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole map rests on monthly-averaged model output being a fair unit of habitability, even though the daily temperature cycle is resolved only for rotation periods of 64 days and longer, so the reported fractions — and the location of the ~20-day drop — would change if shorter warm spells were counted as habitable.","fun_headline_variants_meta":{"raw":{"variants":["Rotation, not orbit, controls Earth-like habitability","Day length rules: 70% to 20% land habitability drop","Climate cliff: 20-day rotation kills Earth-analog life zones","Spin dominates: obliquity matters only for short days","93 GCMs: rotation period sets habitable land fraction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000244,"raw_usage":{"total_tokens":1592,"prompt_tokens":1065,"completion_tokens":527,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":681,"completion_tokens_details":{"reasoning_tokens":441}},"tokens_in":681,"tokens_out":527,"duration_ms":4755,"temperature":1.0,"reasoning_tokens":441,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:40:09.244407+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the existing 64-day and 128-day rotation runs, the models where the diurnal cycle is actually resolved, and recompute the temperature metric against daily output instead of monthly means. If daily-thresholded f_T rises above the monthly value by more than the run-to-run scatter, then the sharp habitability collapse near 20 days is inflated by temporal averaging, and the cliff's true location and steepness would need re-measuring.","supporting_citations":[{"cited_title":"E., & Williams, C","cited_arxiv_id":null,"evidence_quote":"Provides the Gaussian process regression formalism the emulator is built on."},{"cited_title":"J., & Matsuura, K","cited_arxiv_id":null,"evidence_quote":"Supplies the terrestrial precipitation dataset used to justify the 300 mm/yr cutoff in the precipitation metric."}],"review_version":1}