{"id":"2b7b4655-330a-4775-8453-7d42cd07d0ed","arxiv_id":"2607.29057","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A 3D temperature parameterization derived from energy-balance timescales, fit to JWST phase-curve spectra of WASP-121b, recovers a day–night asymmetric atmosphere with a nightside methane enhancement and hints of magnetic drag.","lead":"The paper presents a new analytical 3D temperature model for hot Jupiters, built from an approximate radiative–advective–diffusive balance, and a data-driven method to extract phase-resolved spectra from JWST time series. Applied to WASP-121b, the model retrieves a strongly day–night asymmetric thermal structure, with pressure-dependent hotspot offsets that favor magnetic drag.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The kinematic 1D operator (Eq. 4) is an unrestricted shape model: without a test that fitted (ε,ψ) match the GCM's actual transport timescales, the physical interpretation and magnetic-drag comparison are unsupported.","rationale":"The reader's weakest_assumption correctly identifies the 1D kinematic operator as the load-bearing assumption. My independent reading of the paper confirms this. The central claim—that the retrieval recovers physical transport timescales and that the comparison favors a magnetically damped GCM—depends on the interpretation of ε and ψ as ratios of radiative to advective/diffusive timescales. The paper itself labels the model 'kinematic' and acknowledges the simplification, but then uses the retrieved ε(p),ψ(p) as physical diagnostics (e.g., 'dynamical transport strengthens with pressure') and as a discriminator between magnetic and non-magnetic GCMs. The GCM validation only demonstrates that the analytic function can reproduce large-scale temperature morphology, not that the fitted control parameters correspond to the GCM's true timescale ratios. Without that correspondence, the retrieval could be fitting shape parameters that mimic a 3G GCM's temperature field for reasons unrelated to magnetic drag (e.g., missing opacity, latitudinal degeneracy, or the assumed α(p) logistic). The proposed test—comparing fitted (ε,ψ) against actual GCM diagnostics—would directly establish or refute the physical interpretability. An injection-recovery test from GCM phase curves would additionally test the full pipeline. Both are computationally feasible with the existing GCM output, so the concern is addressable. No fatal flaw is identified; the correct action is to require this validation (or an explicit caveat) before accepting the dynamical interpretation, consistent with the reader's CONDITIONAL verdict.","tokens_in":33551,"tokens_out":4815,"duration_ms":51749,"concrete_test":"Using the Roth et al. (2024) GCM output, compute the zonal-mean zonal wind \\bar u(p), the eddy momentum flux, and the local radiative timescale τ_rad(p) via the Newtonian relaxation approximation. Form τ_adv = 2πR_p/\\bar u (or L/\\bar u) and τ_diff from the diffusion coefficient. Then fit the Eq. 6 parameterization to the GCM temperature field (as in Fig. 5), using the same retrieval-style optimizer, and compare the fitted ε(p) and ψ(p) and the implied hotspot offset to these directly computed timescale ratios. If the fitted values lie outside the GCM's physical values by more than the posterior uncertainty, the parameterization is an ad hoc shape model and the magnetic-drag interpretation is unsupported. A complementary injection-recovery test: generate a synthetic phase curve from the 0G GCM and run the DNChem retrieval; if it still prefers the 3G GCM on the basis of hotspot offset, the","verdict_should_be":"UNCHANGED","load_bearing_attack":"The parameterization is derived by reducing the full energy equation (Eq. 2) to a 1D zonal operator with constant advective and diffusive timescales (Eq. 4). This drops vertical entropy advection, meridional transport, and the thermodynamic feedback between winds and temperature. The dimensionless ratios ε=τ_rad/τ_adv and ψ=τ_rad/τ_diff are therefore not independently measured; they are fitted shape parameters that absorb any missing physics. The paper's central dynamical conclusion—that the small, pressure-dependent hotspot offsets favor a 3G magnetic GCM over 0G—rests on interpreting these effective parameters as physical transport timescales. The GCM validation (Figs. 5–6) only shows that the analytic function can mimic a GCM temperature field; it does not test whether the fitted ε,ψ correspond to τ_rad/τ_adv and τ_rad/τ_diff computed from the GCM's actual winds and radiative relaxation. If they do not, the retrieved ε(p),ψ(p) profiles (Fig. 11) and the GCM comparison (Figs. 10,12) are not evidence about atmospheric drag; they only indicate which analytic shape fits the spectra. This is the essential soft spot because the headline claims about dynamics and magnetic drag depend on it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a framework for extracting and interpreting multidimensional thermal structure from JWST phase curves, applied to the NIRSpec G395H observations of WASP-121b. It combines an 'eclipse normalisation' method that derives phase-resolved spectra by dividing by a stellar template from secondary eclipse, with an analytic 3D temperature parameterisation obtained by reducing the thermodynamic energy equation to a 1D zonal advection–diffusion balance with Newtonian cooling (Eq. 4). The retrieval, run with petitRADTRANS columns and NAUTILUS nested sampling, favours a model with independent dayside and nightside elemental abundances (DNChem). The preferred retrieval shows a strong day–night contrast, dayside and partial nightside thermal inversions, limb temperature asymmetry, and pressure-dependent hotspot offsets increasing from ~4° to ~9°. The authors argue that the confined dayside hot region and small offsets favour a magnetically damped 3 G GCM over a non-magnetic GCM, and that the enhanced nightside CH4 indicates disequilibrium chemistry.","tokens_in":33951,"tokens_out":6955,"duration_ms":74185,"significance":"If the dynamical interpretation survives scrutiny, the framework would be a valuable, computationally efficient tool for multidimensional retrievals. The paper has several concrete strengths: the source code is publicly available, the eclipse-normalised spectra are compared explicitly with standard light-curve fits, the radiative-transfer approximation is benchmarked against PICASO, and the retrieval uses a standard Bayesian framework with quantitative model comparison. The DNChem preference and the nightside CH4 interpretation are interesting and plausibly robust within the tested model family. However, the central dynamical claims rest on interpreting the fitted parameters ε and ψ as physical transport timescales, and that interpretation is not validated. The GCM comparisons are fits rather than predictions, and no test is offered that the retrieved (ε, ψ) correspond to the actual τ_rad/τ_adv and τ_rad/τ_diff in the GCMs. The paper is therefore useful and publishable in principle, but the dynamical conclusions currently outrun the evidence provided.","major_comments":[{"comment":"The reduction from Eq. (2) to Eq. (4) is a heuristic replacement of the full dynamical transport by a constant-coefficient 1D operator, as the paper itself acknowledges by calling it a 'kinematic model'. But the paper later interprets the fitted ε and ψ as physical timescale ratios and uses them to argue for magnetic drag. Because vertical entropy advection, meridional transport, and thermodynamic feedbacks are absorbed into these effective parameters, the retrieved ε(p), ψ(p) are shape parameters. The paper does not test whether they match the GCM's actual τ_rad/τ_adv and τ_rad/τ_diff. A concrete diagnostic would be to compute these timescale ratios from the GCM wind and radiative fields and compare them with the retrieved profiles, or to run synthetic retrievals on GCM spectra to see whether the true transport timescales are recovered. Without such a test, the magnetic-drag conclusion","section":"§3.1, Eq. (4); §4.3, Figs. 11–12; §5.1"},{"comment":"The abstract and §3.2 state that the parameterisation 'reproduces' GCM temperature fields. However, the validation is a least-squares fit of the parameterisation to each GCM field; agreement therefore measures expressive flexibility, not physical fidelity. The paper should either provide an out-of-sample test (e.g., fit one GCM and predict another, or withhold a pressure/longitude region) or compare the fitted (ε, ψ) with independently diagnosed GCM transport timescales. Without this, Figs. 5–6 do not establish that Eq. (4) captures the dynamical content claimed in the title and conclusions.","section":"§3.2, Figs. 5–6"},{"comment":"The monotonic increase of ε and ψ with pressure is interpreted as 'dynamical heat redistribution becomes progressively more efficient at depth' (Fig. 11 caption) and as 'radiative relaxation is faster than advective and diffusive transport at low pressures' (§4.3). This attribution is not justified because ε=τ_rad/τ_adv depends on τ_rad, which in real hot-Jupiter atmospheres increases strongly with pressure. The same retrieved ε(p) profile could result from a radiative-timescale gradient with constant advection. Since τ_rad is not independently constrained—it is absorbed into the dimensionless ratios—the pressure dependence in Fig. 11 alone does not constrain the pressure dependence of transport. The text should be revised to treat ε(p), ψ(p) as phenomenological, or the analysis should include an independent radiative-timescale prior or diagnostic.","section":"§4.3, Fig. 11"}],"minor_comments":[{"comment":"The phrase 'derived from radiative, advective and diffusive energy balance' overstates the status of Eq. (4). The paper itself calls it a kinematic model; the abstract and title should reflect this, e.g., 'inspired by' or 'parameterised following' a simplified energy balance.","section":"Abstract and §3.1"},{"comment":"The text says the hotspot offset is 'well approximated' by tan⁻¹(ε/(1+ψ)) in the strong-radiative regime, but retrieved ε and ψ reach values well above unity at depth. Please state explicitly whether the reported offsets come from the numerical solution of Eq. (14) or from the approximation, and verify the approximation against the full solution at the retrieved parameters.","section":"§3.1, Eqs. (13)–(14)"},{"comment":"The description of the linear detrending step is vague: it should be stated more precisely how the linear function fitted to in-eclipse data is extrapolated to out-of-eclipse phases and propagated into the uncertainties. This is important for reproducing the eclipse-normalisation pipeline.","section":"§2.2"},{"comment":"The bottom panel labelled 'Diff. (σ)' is described as a histogram in the caption but the axes and units are not explained. Please clarify what is plotted.","section":"Fig. 3"},{"comment":"The Cloud model is reported with ΔlnZ = −17.05. This is weak-to-moderate evidence and should be described with appropriate caution; the text currently says it is 'disfavoured' without qualifying the strength.","section":"Table 2"},{"comment":"Several instances read 'W ASP-121b' with an unwanted space, and 'JWSTandAriel' is missing a space. Also, the data DOI is given incompletely as 'doi:10.17909' and should be completed.","section":"References / typesetting"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is likely publishable after major revision. The eclipse-normalisation method and the analytic parameterisation are useful contributions, and the retrieval analysis is careful in many respects. The main concern is that the headline dynamical conclusions—especially the magnetic-drag interpretation—overreach the kinematic model. I would ask the authors to add the proposed validation (or substantially soften the dynamical claims). I see no grounds for rejection; the central retrieval methodology is sound and the chemical findings are interesting."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth your time if you work on phase-curve retrievals. The analytic temperature field (Eqs. 6–7) is a real addition: a closed-form 3D structure with a handful of parameters, able to reproduce the large-scale morphology of GCMs. The eclipse-normalisation method is simple and validated against standard light-curve fitting, and the WASP-121b retrieval is careful, with a sensible model comparison table. The thermal-structure results—strong day–night contrast, dayside and partial nightside inversions, limb asymmetry—are plausible and consistent with prior work.\n\nThe main soft spot is the physical interpretation. Eq. 4 is a kinematic reduction, not a derivation from the primitive equations. The paper says as much, but the abstract and some phrasing (\"derived from radiative, advective and diffusive energy balance\") oversell it. The fitted ε and ψ are effective shape parameters, and the GCM validation only shows the functional form can mimic a GCM field, not that the fitted values correspond to actual τ_rad/τ_adv and τ_rad/τ_diff from the GCM's winds. That matters for the interpretation of Figure 11. However, the magnetic-drag comparison mostly uses the hotspot offset, which is a direct property of the retrieved temperature field, so that part does not collapse. The comparison is qualitative and the paper correctly hedges with Rayleigh drag. Still, a stronger check—computing transport timescales from the GCM and comparing with fitted ε,ψ—would substantially upgrade the dynamical claims.\n\nMinor issues: some DNChem abundances sit at the prior boundary ([Si/H]_n), the disequilibrium-chemistry conclusion is indirect (a phenomenological day/night abundance split), and the 1D column radiative transfer could bias the retrieved structure; the paper acknowledges the last point.\n\nNet: the framework is a useful addition to the retrieval toolbox, and the WASP-121b analysis is solid enough. It deserves peer review. I'd push for revisions that temper the \"derived\" language, explicitly define ε,ψ as effective parameters, and add the GCM timescale test if they can. The magnetic-drag claim should remain a qualitative comparison.","headline":"A genuinely useful analytic 3D temperature parameterization and a clean data-extraction method, with a dynamical interpretation that overshoots what the kinematic model can support.","tokens_in":34427,"tokens_out":4690,"would_cite":true,"duration_ms":46892,"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":"The paper claims that a five-parameter analytic temperature law, derived from radiative-advective-diffusive balance, lets a JWST phase curve of WASP-121b be read as a three-dimensional thermal structure.","keywords":["exoplanet atmospheres","hot Jupiters","phase curves","atmospheric retrieval","three-dimensional temperature structure","WASP-121b","radiative-advective-diffusive balance","JWST"],"falsifier":"Generate a synthetic JWST-like phase curve from a general circulation model that includes strong vertical entropy advection, run the same retrieval, and compare the recovered ε(p), ψ(p), and hotspot offsets with the GCM's true values; any systematic mismatch demonstrates that the kinematic operator is absorbing unrepresented transport rather than measuring it.","tokens_in":33476,"feed_emoji":"🪐","tokens_out":5746,"duration_ms":61569,"temperature":0.7,"pith_summary":"This paper tries to establish that the large-scale three-dimensional temperature field of a strongly irradiated, tidally locked hot Jupiter can be captured by a closed-form analytic solution to a one-dimensional longitudinal energy-balance equation, with just a few physically meaningful timescale ratios as controls. If true, it gives retrievals a cheap, physically interpretable bridge between disk-integrated phase-curve spectra and 3D atmospheric structure, without running full general circulation models inside the retrieval loop. Applied to the JWST/NIRSpec G395H phase curve of WASP-121b, the framework yields a consistent 3D thermal map: a strong day–night contrast, dayside and partial nightside thermal inversions, limb temperature asymmetry, and pressure-dependent hotspot offsets increasing from about 4° to 9°. The paper further argues that the small offsets favor a magnetically damped circulation over a non-magnetic one, and that the spectra require distinct dayside and nightside chemical states, pointing to disequilibrium chemistry.","feed_headline":"Hotspot drift points to magnetic drag on WASP-121b","feed_subtitle":"The hotspot shifts from 4 to 9 degrees with depth, and day and night sides carry different chemistry.","key_machinery":"The engine is the kinematic energy-balance equation (Eq. 6): ε ∂T/∂θ − ψ ∂²T/∂θ² = T_eq − T, where ε = τ_rad/τ_adv and ψ = τ_rad/τ_diff are the ratios of the radiative relaxation timescale to the advective and effective diffusive timescales. This linear ODE in longitude has a closed-form piecewise solution (Eq. 7) that gives the full T(θ, ϕ, p) field once the substellar and antistellar temperature profiles, T_d(p) and T_n(p), and a longitudinal sharpness parameter α(p) are supplied. The hotspot longitude is defined by ∂T/∂θ = 0 and reduces to θ_h ≈ tan⁻¹(ε/(1+ψ)) in the strong-radiative limit. This single balance is what lets the observable phase-curve morphology be translated into physical","core_discovery":"The central claim is that the three-dimensional temperature field of a strongly irradiated, tidally locked giant planet can be written analytically as the solution of a steady one-dimensional longitudinal balance between Newtonian radiative relaxation, zonal advection, and effective diffusion, with the ratios ε = τ_rad/τ_adv and ψ = τ_rad/τ_diff as control parameters. This closed-form temperature law, combined with a data-driven eclipse-normalisation that reads phase-resolved spectra directly from the time series, lets a Bayesian retrieval recover large-scale thermal structure from a JWST/NIRSpec G395H phase curve. For WASP-121b the retrieval finds a pronounced day–night contrast, dayside an","pith_inferences":["If the kinematic balance holds, the retrieved vertical profiles of ε(p) and ψ(p) could be mapped onto pressure-dependent zonal wind speeds, giving a direct, testable link between phase-curve retrievals and high-resolution Doppler wind measurements.","The eclipse-normalisation extraction is model-free at the reduction stage; a natural extension is to feed the same phase-resolved spectra into high-resolution cross-correlation analyses, potentially breaking degeneracies that disk-integrated photometry alone leaves unresolved.","The recovered day/night chemical contrast is better read as a diagnostic of vertical mixing and chemical quench levels than as a literal elemental abundance discontinuity; coupled kinetic models could test whether nightside vertical transport timescales are short enough to maintain the CH4 excess.","Because the paper itself notes that 1D-column radiative transfer may bias retrieved gradients, a testable extension is to recompute the favoured WASP-121b model with full 3D radiative transfer and check whether the limb asymmetry and hotspot offsets persist."],"forward_implications":["If the parameterisation is correct, phase-curve retrievals can recover the large-scale 3D thermal structure of hot Jupiters at a fraction of the cost of full GCM-based retrieval, making multidimensional characterisation routine for JWST and Ariel phase curves.","For WASP-121b, the retrieved structure implies a pressure-dependent dynamical regime: radiation dominates the upper atmosphere while advection and diffusion become progressively more important at depth.","The small, pressure-dependent hotspot offsets are more consistent with a magnetically damped circulation (3 G case) than a non-magnetic one, although the paper stresses that Rayleigh drag or other damping cannot be excluded.","The strong Bayesian preference for independent dayside and nightside elemental abundances (ΔlnZ ≈ 65 over shared abundances) indicates that a single globally uniform chemical composition is inadequate within the tested model family, implying real chemical inhomogeneity.","The retrieved thermal structure implies inhomogeneous cloud condensation conditions, with the cooler nightside and morning limb favouring condensates while the hottest dayside remains largely cloud-free."],"fun_headline_variants":["Hotspot drift on WASP-121b hints at magnetic drag","JWST resolves 3D thermal map of WASP-121b","Phase curve yields analytic 3D atmosphere for WASP-121b","WASP-121b's hotspot shift reveals magnetic influence","New model reads JWST phase curves in 3D"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise, stated in Section 3.1, is that all dynamical transport—including vertical entropy advection—can be folded into a single one-dimensional longitudinal advection–diffusion operator with constant timescales; if vertical or meridional transport contributes significantly, the retrieved ε, ψ, and hotspot offsets are effective quantities and the drag interpretation could be biased.","fun_headline_variants_meta":{"raw":{"variants":["Hotspot drift on WASP-121b hints at magnetic drag","JWST resolves 3D thermal map of WASP-121b","Phase curve yields analytic 3D atmosphere for WASP-121b","WASP-121b's hotspot shift reveals magnetic influence","New model reads JWST phase curves in 3D"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000132,"raw_usage":{"total_tokens":1019,"prompt_tokens":846,"completion_tokens":173,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":84}},"tokens_in":590,"tokens_out":173,"duration_ms":2805,"temperature":1.0,"reasoning_tokens":84,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T14:25:34.476842+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Generate a synthetic JWST-like phase curve from a general circulation model that includes strong vertical entropy advection, run the same retrieval, and compare the recovered ε(p), ψ(p), and hotspot offsets with the GCM's true values; any systematic mismatch demonstrates that the kinematic operator is absorbing unrepresented transport rather than measuring it.","supporting_citations":[],"review_version":1}