{"id":"692ee9f7-2bd7-4836-abb6-6f66c042ad57","arxiv_id":"1908.02661","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Fast-rotating Earth-like planets should show larger equator-to-pole temperature contrasts when they spin faster or have thinner atmospheres, with bulk lapse rate peaking at baroclinic criticality.","lead":"This paper derives simple scaling relations for the equator-to-pole temperature contrast and bulk lapse rate of fast-rotating terrestrial exoplanet atmospheres using baroclinic criticality theory. The scalings are tested with general circulation model simulations and could help interpret future space telescope observations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (5)'s Held-Larichev closure sets every exponent in Eqs. (6), (13), and (14), but the GCM suite never tests it directly, and all tested runs sit near ξ≈1, leaving the claimed baroclinically unstable regime outside the validated closure.","rationale":"I agree with the reader's weakest_assumption: the Held-Larichev closure in Equation (5) is the sole source of the power-law exponents in Equation (6), and the subsequent equations are algebraic consequences of that scaling. The reader's CONDITIONAL verdict is appropriate because the concern is concrete and testable but not immediately fatal: within the tested near-critical regime, Figure 2's agreement with Equation (6) provides indirect support, and the paper explicitly restricts its claims to baroclinically unstable rotation rates. I would not change the verdict because the central argument is internally coherent and the concern is principally about extrapolation rather than an internal inconsistency. The proposed direct diagnostic test would settle whether the closure holds across the parameter sweep. Other caveats, such as model-derived tropopause heights and the fixed Δ_hθ_eq taken from a single albedo, are real but secondary: they weaken the strength of the validation without undermining the scaling logic itself.","tokens_in":11585,"tokens_out":16250,"duration_ms":190487,"concrete_test":"Use the existing GCM output to directly diagnose the eddy diffusivity D = <v'θ'>/|∂θ/∂y|, the deformation radius L_d, and the criticality ξ (as in Equation 15) for each rotation and pressure run. Then plot D against β(ξL_d)^3 across all runs. If the ratio D/[β(ξL_d)^3] is not constant to within a factor of about 2 over Ω∈[0.5,8]Ω_Earth and p∈[0.25,4] bar, the closure in Equation (5) is not supported and the exponents should be re-derived. A complementary run at Ω=16Ω_Earth or p=0.125 bar would test whether Equation (6) continues to hold in a more strongly supercritical regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The scaling exponents in Equations (6), (13), and (14) all inherit the closure D_eddy ~ β(ξ L_d)^3 from Equation (5), introduced in Section 2.1. This is a geostrophic-turbulence result (Held & Larichev 1996) that assumes an inverse energy cascade arrested at the Rhines scale and a Rhines scale proportional to ξL_d. Equation (5) is never tested directly in the paper. The GCM comparison in Figure 2 tests Equation (6) as a whole, so agreement there cannot distinguish the closure from compensating errors in τ_rad, H, or L_d. Moreover, the simulated parameter sweep clusters near ξ≈1, the regime where the closure is most plausible; the paper's broader claim of applicability throughout the baroclinically unstable regime requires extrapolating to larger Ω and different pressures with no direct validation. If the true eddy diffusivity has a different dependence on ξL_d—for example, a weaker scaling in strongly supercritical turbulence—the predicted Ω, p, and H exponents in Equation (6) shift, and the central predictions of Equations (13) and (14) change accordingly.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript derives scaling relations for the baroclinic criticality parameter, equator-to-pole potential temperature contrast, and bulk lapse rate of terrestrial exoplanet atmospheres. Building on Jansen & Ferrari (2013) and the Held-Larichev eddy-diffusivity closure, the authors obtain Equation (6) for the criticality parameter as a function of rotation rate, surface pressure, tropopause height, radius, and gravity, and then Equations (13) and (14) for the temperature contrast and lapse rate. The predictions are compared to ExoCAM GCM simulations with rotation rates from 0.0625 to 8 Earth values and surface pressures from 0.25 to 4 bars. The paper reports broad agreement for rotation rates greater than about one third of Earth's rotation rate and discusses observational tests of the scalings. Limitations of the theory for slowly rotating and tidally locked planets are acknowledged.","tokens_in":11850,"tokens_out":3731,"duration_ms":41657,"significance":"If the scalings hold, they provide compact, observationally relevant predictions for how the circulation of Earth-like exoplanets depends on basic planetary parameters. The paper connects a well-established body of geophysical fluid dynamics to exoplanet observations, and the GCM comparison is a reasonable and transparent test. The authors are honest about key inputs, such as the Earth-normalized criticality and the model-derived tropopause height and albedo, and they explicitly delimit the regime of applicability. The main weakness is that the validation is concentrated near the marginally critical regime, so the claim of applicability throughout the baroclinically unstable regime is not yet strongly supported.","major_comments":[{"comment":"The exponents in Equations (6), (13), and (14) all inherit the Held-Larichev closure D_eddy ~ beta (xi L_d)^3 from Equation (5). The GCM comparisons test the integrated scaling rather than this closure directly, and the simulated runs lie mostly at or near xi ~ 1 (as noted in the discussion of Figure 4). Consequently, the paper's claim of applicability 'throughout the baroclinically unstable regime' rests largely on extrapolation to strongly supercritical states (xi >> 1) that the suite does not sample. I recommend adding a simulation or reanalysis at higher rotation rate or lower pressure that pushes xi well above 1, or substantially softening the scope claim to the near-critical regime.","section":"Section 2.1, Equations (5)-(6) and Figures 2-4"},{"comment":"The pressure scaling is not a fully independent prediction because the predicted criticality uses the tropopause height H taken from the GCM experiments (the text states this explicitly). Since H enters Equation (6) with an exponent -3/5 and itself varies with surface pressure, the agreement between theory and model in Figure 2 partly reflects the model's H rather than a pure a priori prediction. The disclosure is appreciated, but the discussion should distinguish which parts of the pressure dependence are predicted and which are diagnosed, and should quantify the sensitivity to the choice of H (e.g., using a fixed H would change the predicted exponent to -1/5).","section":"Section 4.1, Figure 2 (right panel)"},{"comment":"The fixed equilibrium equator-to-pole contrast Delta_h_theta_eq = 242 K is derived using an albedo of 0.55 taken from the Earth-like simulation, and the manuscript notes that albedo changes due to sea ice are ignored. For the 2 and 4 bar cases, which are ice-covered, the change in Delta_h_theta_eq is not accounted for and is invoked as the reason for the discrepancy. This is an acknowledged limitation, but it directly affects the pressure scaling that is a central target of the paper. Please estimate how Delta_h_theta_eq changes with the simulated albedo across the suite, so readers can judge how much of the apparent pressure trend in Figure 3 is attributable to the fixed Delta_h_theta_eq assumption.","section":"Section 4.2, Equation (13) and Figure 3"}],"minor_comments":[{"comment":"The word 'baroclincally' is a typo and should be 'baroclinically'.","section":"Abstract"},{"comment":"The phrase 'there is no tunable parameter in our equations' is misleading because the theory is anchored by setting xi = xi_Earth for Earth-like parameters and by fixing Delta_h_theta_eq with an albedo taken from the model. These are inputs rather than tunable fits, but a more precise statement would acknowledge them as calibration choices.","section":"Section 4.2, paragraph 1"},{"comment":"The caption refers to 'the arrow' but the figure does not visibly include an arrow in the rendered version; please either include the arrow in the figure or remove the reference.","section":"Figure 2 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its inputs and limitations, and the derivation is internally consistent. The main concern is that the validation does not cover the strongly supercritical regime that the title and abstract emphasize. I would encourage the authors to either add a high-rotation or low-pressure simulation that clearly places them in the xi >> 1 regime, or to explicitly reframe the conclusions as applying to the near-critical regime. The pressure-comparison caveat about using GCM tropopause heights should also be elevated from a sentence to a more prominent discussion of what the test actually constrains."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper gives the exoplanet community something useful: explicit power-law scalings for equator-to-pole temperature contrast and bulk lapse rate as functions of rotation rate, surface pressure, tropopause height, radius, and gravity, built on Jansen & Ferrari's baroclinic criticality framework. The derivation is clear and honest, and the comparison to ExoCAM GCM simulations across rotation rates and surface pressures shows the theory captures the qualitative trends and much of the quantitative behavior. That is a genuine new contribution, not just a repackaging of existing results.\n\nThe soft spots are real but not fatal. The stress-test note is right that the Held-Larichev closure in Eq. (5) is the load-bearing assumption behind every exponent, and the paper never tests D_eddy directly. What the GCM comparison does test is the end-to-end scaling, so agreement there does not isolate the closure from compensating errors in tau_rad or L_d. Also, the simulated runs sit close to xi ~ 1 (roughly 0.7 to 3), so the claim that the theory works throughout the \"baroclinically unstable regime\" is a modest extrapolation to strongly supercritical cases. I would like to see a few faster-rotation runs or a direct look at the eddy diffusivity, but this is a limitation to flag, not a reason to reject.\n\nTwo other caveats are worth naming. The surface-pressure comparison uses tropopause heights taken from the GCM itself, which makes that test partly circular as an independent prediction. And the 2 bar and 4 bar runs are in an ice-covered state, so the fixed Delta_h_theta_eq = 242 K with albedo 0.55 misses the albedo change; the theory cannot capture the step-like contrast shift in those runs. Both are addressable and the authors are upfront about the ice-coverage issue.\n\nOverall, the central argument holds up. The scalings are derived transparently, the GCM evidence is reasonably supportive, and the observational implications are concrete. Who gets value: exoplanet atmospheric modelers, observers planning phase-curve or retrieval campaigns, and anyone needing quick estimates of climate contrasts for Earth-like planets. This deserves a serious referee, and I would send it to review rather than desk reject. My recommendation: accept after a moderate revision that addresses the closure caveat and clarifies the limits of the pressure comparison.","headline":"A clean extension of baroclinic criticality theory to exoplanets, with GCM support for the main trends; the main caveat is that the closure is only validated indirectly, with runs clustering near xi ~ 1.","tokens_in":12390,"tokens_out":1791,"would_cite":true,"duration_ms":22755,"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":"This paper derives two scaling equations that predict the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmospheres from rotation rate, surface pressure, tropopause height, radius, and gravity, and…","keywords":["baroclinic criticality","terrestrial exoplanet atmospheres","equator-to-pole temperature contrast","bulk lapse rate","atmospheric heat transport","general circulation models","ExoCAM","rotation rate scaling"],"falsifier":"Run a GCM at $16\\Omega_\\oplus$ and 1 bar: the theory predicts $\\xi$ about three times Earth's value and an equator-to-pole potential temperature contrast near 218 K (for $\\Delta_h \\theta_{eq} = 242$ K), versus about 121 K at Earth's rotation; if the simulated contrast does not approach that value, or its growth with rotation deviates from the predicted $\\Omega^{2/5}$ trend, the central scaling fails.","tokens_in":11399,"feed_emoji":"🌍","tokens_out":7882,"duration_ms":76954,"temperature":0.7,"pith_summary":"The paper asks how fast-rotating terrestrial exoplanet atmospheres transport heat poleward, and derives closed-form scaling laws for two observable quantities: the equator-to-pole potential temperature contrast and the bulk lapse rate. The laws follow from balancing baroclinic eddy heat fluxes against radiative relaxation, building on a prior theory of baroclinic criticality. The paper compares the scalings against ExoCAM general circulation model simulations spanning rotation rates 0.0625 to 8 times Earth's and surface pressures 0.25 to 4 bars, and finds broad agreement wherever baroclinic instability is active. If the theory is right, future observations could use thermal phase curves and spectra to estimate a planet's rotation rate and surface pressure from its temperature structure.","feed_headline":"Spin rate sets exoplanet equator-pole temperature contrast","feed_subtitle":"Criticality theory matches GCM runs; future telescopes could read spin and pressure from thermal maps.","key_machinery":"The central object is the baroclinic criticality parameter $\\xi = s a / H$, where $s$ is the slope of mid-latitude isentropes, $a$ is planetary radius, and $H$ is the smaller of the tropopause height and the scale height. The argument runs through three steps: the isentropic slope is set by the distance eddies diffuse heat in one radiative relaxation time, $s \\sim H/\\sqrt{\\tau_{\\mathrm{rad}} D_{\\mathrm{eddy}}}$; the eddy diffusivity obeys the geostrophic-turbulence closure $D_{\\mathrm{eddy}} \\sim \\beta (\\xi L_d)^3$ from Held & Larichev (1996), with $L_d$ the Rossby deformation length; and the horizontal and vertical eddy heat fluxes are related along isentropes, giving a constraint that separates the equator-to-pole and vertical contrasts. With the radiative relaxation timescale scaling $\\tau_{\\mathrm{rad}} \\propto p/(gT^3)$, this machinery converts an abstract measure of instability into observable temperature contrasts.","core_discovery":"The central claim is that the baroclinic criticality parameter $\\xi$—a measure of how slanted mid-latitude isentropes are, roughly the ratio of equator-to-pole to surface-to-tropopause potential temperature contrast—controls both heat transport and lapse rate, and that it scales as $\\xi \\propto (\\Omega/\\Omega_\\oplus)^{2/5}(p/p_\\oplus)^{-1/5}(H/H_\\oplus)^{-3/5}(a/a_\\oplus)^{3/5}(g/g_\\oplus)^{-1/10}$. Combining this with the assumption that eddy heat flux is directed along isentropes yields equations (13) and (14), which predict the equator-to-pole potential temperature contrast and the bulk lapse rate (defined as the smaller of the surface-to-tropopause contrast and the contrast over one scale height) with no free parameters once the radiative-equilibrium contrast is fixed. The GCM comparisons show the equator-to-pole contrast increases with rotation rate and decreases with surface pressure, while the bulk lapse rate varies weakly and peaks near $\\xi \\approx 1$. The paper argues this confirms that baroclinic instabilities, not just radiative-convective balance, set the temperature structure of fast-rotating terrestrial exoplanets.","pith_inferences":["A retrieved equator-to-pole contrast alone would not pin down rotation rate or surface pressure, because the scalings depend on the combination $\\xi \\sim \\Omega^{2/5} p^{-1/5}$ together with tropopause height, radius, and gravity; independent estimates of pressure or tropopause height would be needed to break the degeneracy.","Because the theory fixes the radiative-equilibrium contrast $\\Delta_h \\theta_{eq}$ at one value, ice-albedo feedbacks are an unmodeled pathway; coupling the lapse-rate prediction to a sea-ice model might explain the anomalous ice-covered 2-4 bar cases noted in the paper.","The predicted $\\Omega^{2/5}$ exponent could be checked before exoplanet spectroscopy is ready, using existing dry dynamical-core GCMs or rotating-tank experiments that vary rotation rate over a wider range than the present suite.","If the eddy-diffusivity closure survives, the same criticality balance may give testable predictions for the phase-curve amplitudes of fast-rotating gas giants as functions of their rotation rate and radiative timescale."],"forward_implications":["At fixed incident stellar flux, faster rotation widens the equator-to-pole potential temperature contrast, so fast rotators should have colder poles and larger sea-ice cover, an effect that could appear in orbital-phase albedo variations.","The bulk lapse rate is only weakly dependent on planetary parameters near Earth-like values and peaks near $\\xi \\approx 1$; for rotation rates above about $8\\Omega_\\oplus$, further spin-up should reduce it rather than increase it.","Spectroscopic retrievals that constrain the vertical temperature profile and tropopause height can directly test the predicted lapse rate, while full-phase thermal light curves can constrain the equator-to-pole contrast.","The theory applies only to planets with active baroclinic instability, roughly those with rotation periods of about three days or shorter; slowly rotating and tidally locked planets require separate treatments.","Extending the same reasoning to warm Jupiters and warm Neptunes is plausible, and the predicted increase of equator-to-pole contrast with rotation rate matches previous simulations of warm Jupiters."],"supporting_citations":[{"why":"Supplies the core criticality framework relating isentropic slope to the eddy heat flux balance, and the separate scalings for horizontal and vertical potential temperature contrasts that become equations (13) and (14).","marker":"Jansen & Ferrari (2013)"},{"why":"Provides the geostrophic-turbulence closure $D_{\\mathrm{eddy}} \\sim \\beta (\\xi L_d)^3$ that determines the exponents in the criticality scaling of equation (6).","marker":"Held & Larichev (1996)"},{"why":"Gives the radiative relaxation timescale scaling $\\tau_{\\mathrm{rad}} \\propto p/(gT^3)$ used to convert the criticality relation into explicit planetary-parameter scalings.","marker":"Showman & Guillot (2002)"},{"why":"The GCM simulation suite (same ExoCAM setup, extended here to faster rotation) is the dataset against which the scaling theory is compared.","marker":"Komacek & Abbot (2019)"},{"why":"Defines the ExoCAM aquaplanet configuration, including the N2-H2O atmosphere and slab ocean, used for all simulations.","marker":"Kopparapu et al. (2017)"},{"why":"Motivates taking $H$ as the minimum of the tropopause height and the scale height in the criticality and temperature-contrast scalings.","marker":"Chai & Vallis (2014)"}],"fun_headline_variants":["Baroclinic criticality sets exoplanet heat transport","Exoplanet temperature contrast scales with spin and pressure","New scaling predicts exoplanet equator-pole temperature gap","Theory matches GCMs for terrestrial exoplanet atmospheres","Spin and pressure control exoplanet thermal structure"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions inherit the assumption that eddy mixing follows the quasi-geostrophic turbulence law $D_{\\mathrm{eddy}} \\sim \\beta (\\xi L_d)^3$, with the Rhines scale tied to the deformation radius; if atmospheres with very different rotation or pressure mix according to a different law, the derived rotation and pressure exponents would not hold.","fun_headline_variants_meta":{"raw":{"variants":["Baroclinic criticality sets exoplanet heat transport","Exoplanet temperature contrast scales with spin and pressure","New scaling predicts exoplanet equator-pole temperature gap","Theory matches GCMs for terrestrial exoplanet atmospheres","Spin and pressure control exoplanet thermal structure"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000187,"raw_usage":{"total_tokens":1381,"prompt_tokens":1049,"completion_tokens":332,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":665,"completion_tokens_details":{"reasoning_tokens":252}},"tokens_in":665,"tokens_out":332,"duration_ms":3760,"temperature":1.0,"reasoning_tokens":252,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:38:20.840553+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a GCM at $16\\Omega_\\oplus$ and 1 bar: the theory predicts $\\xi$ about three times Earth's value and an equator-to-pole potential temperature contrast near 218 K (for $\\Delta_h \\theta_{eq} = 242$ K), versus about 121 K at Earth's rotation; if the simulated contrast does not approach that value, or its growth with rotation deviates from the predicted $\\Omega^{2/5}$ trend, the central scaling fails.","supporting_citations":[{"cited_title":"& Ferrari, R","cited_arxiv_id":null,"evidence_quote":"Supplies the core criticality framework relating isentropic slope to the eddy heat flux balance, and the separate scalings for horizontal and vertical potential temperature contrasts that become equations (13) and (14)."},{"cited_title":"& Larichev, V","cited_arxiv_id":null,"evidence_quote":"Provides the geostrophic-turbulence closure $D_{\\mathrm{eddy}} \\sim \\beta (\\xi L_d)^3$ that determines the exponents in the criticality scaling of equation (6)."},{"cited_title":"& Guillot, T","cited_arxiv_id":null,"evidence_quote":"Gives the radiative relaxation timescale scaling $\\tau_{\\mathrm{rad}} \\propto p/(gT^3)$ used to convert the criticality relation into explicit planetary-parameter scalings."},{"cited_title":"& Abbot, D","cited_arxiv_id":null,"evidence_quote":"The GCM simulation suite (same ExoCAM setup, extended here to faster rotation) is the dataset against which the scaling theory is compared."},{"cited_title":"& Vallis, G","cited_arxiv_id":null,"evidence_quote":"Motivates taking $H$ as the minimum of the tropopause height and the scale height in the criticality and temperature-contrast scalings."}],"review_version":1}