{"id":"e23c705e-adac-4559-82e7-f53e34129ccf","arxiv_id":"2607.06245","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":4,"one_line_summary":"Stable layers in Jupiter's interior modify the background density, creating a degeneracy with the vertical wind profile that prevents unique determination of jet depth from gravity data alone.","lead":"This paper shows that stable stratified layers in Jupiter's interior change the planet's gravity field in ways that mimic deeper jet streams, meaning the vertical structure of Jupiter's jets cannot be uniquely determined from gravity data alone. A smart generalist should read it because it quantifies a fundamental degeneracy in interpreting Juno spacecraft data, showing that what we thought we knew about how deep Jupiter's winds penetrate depends on an assumption about the深层热","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The degeneracy claim is mathematically sound, but its quantitative magnitude hinges on a near-maximal stratification amplitude (k=0.9∇ad) that is never varied; the qualitative claim survives but the quantitative results are illustrative only.","rationale":"The reader correctly identified the most load-bearing concern: the artificial prescription of the stable layer with k=0.9∇ad, fixed for all configurations, is the weakest link in the quantitative argument. The reader's assessment that the qualitative degeneracy claim is well-supported while the quantitative details are 'illustrative rather than definitive' is accurate. The mathematical argument in Eq. 9 is sound — both terms on the RHS do contribute to the dynamical density, and the optimization experiments (Figs. 7–9) demonstrate that different (ρ_stat, Q) combinations can reproduce the observed harmonics. The concern is not about correctness of the framework but about the sensitivity of the quantitative results to an unconstrained parameter. The reader's secondary concerns (mass conservation, missing error bars for SL profiles) are valid but less load-bearing: mass non-conservation primarily affects even harmonics that are explicitly omitted, and the visual fits in Fig. 8 suggest the SL solutions are viable even without formal error bars. The GCM/ERA5 comparison sections (§4–5) are indeed tangential — they contextualize the degeneracy but do not strengthen the core result. The CONDITIONAL verdict with MODERATE confidence is appropriate. The paper makes a legitimate contribution by explicitly demonstrating the ρ_stat–Q degeneracy in the gravity inversion problem, but the quantitative magnitude of the effect remains undetermined without varying k. No verdict adjustment is needed.","tokens_in":23131,"tokens_out":3795,"duration_ms":245819,"concrete_test":"Re-run the constrained optimization (Eq. 10, §3.2) for the most extreme SL case (1–10^5 bar) with k = 0.1, 0.3, 0.5, and 0.9 × ∇ad. For each k, record the optimal wind decay scale H and the cost function value L (Eq. B1). If the optimal H shifts by less than ~20% across the full k range, the quantitative claim is robust to stratification strength. If H changes substantially (e.g., the wind adjustment halves when going from k=0.9 to k=0.3), the quantitative results should be explicitly framed as scaling with the unconstrained stratification amplitude, and the headline claim about 'substantial' modification should be qualified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central qualitative claim — that a degeneracy exists between ∂ρ_stat/∂z and ∂Q/∂z in the thermal wind balance (Eq. 9) — is mathematically straightforward and holds for any k > 0. However, the paper's quantitative claims ('shallow, extensive stable layers substantially modify the background density, requiring more rapid decay of zonal winds') depend critically on the choice k = 0.9∇ad (Eq. 1, §2.1). This value is near the maximum physically possible subadiabatic departure, meaning the density enhancement within the SL is maximized. The paper explicitly states that k is 'fixed for all stable layers considered' to isolate location/extent effects, but never tests sensitivity to k itself. If the real stratification is weaker (e.g., k = 0.1–0.3∇ad, consistent with some estimates cited in §2.1), the density modification shrinks proportionally, the required wind adjustment diminishes, and the 'substantial' modification becomes modest. The paper is honest that 'the actual Brunt–Väisälä frequency within Jupiter's interior remains poorly constrained,' but does not frame its results as an upper bound, which they effectively are at k = 0.9∇ad. Additionally, error bars are shown only for the adiabatic profile (Fig. 6, §3.1), so the reader cannot assess whether SL profiles achieve statistically equivalent fits or merely visually acceptable ones. The mass non-conservation (§2.1) is a lesser concern because the paper omits the even harmonics that would be most affected, and the odd harmonics arise entirely from the dynamical density.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript investigates how subadiabatic stable layers (SLs) in Jupiter's interior affect the inference of zonal wind vertical structure from Juno gravity harmonics. The authors construct temperature-pressure profiles with prescribed stable layers (Eq. 1), derive density profiles using the Chabrier et al. (2019) EOS, and compute gravitational harmonics via thermal wind balance (Eq. 9). The central result is a degeneracy: increasing ∂ρ_stat/∂z (from a stable layer) produces a similar gravitational signature to increasing ∂Q/∂z (from deeper winds), meaning that the jet stream structure cannot be uniquely determined without independent constraints on internal stability. The paper demonstrates this by showing that shallow, extensive SLs require more rapid wind decay to match Juno observations, and by comparing the resulting envelope of wind profiles with GCM outputs and Earth's eddy-driven jets.","tokens_in":23351,"tokens_out":1413,"duration_ms":126094,"significance":"The degeneracy identified in Eq. 9 — that ∂ρ_stat/∂z and ∂Q/∂z enter the thermal wind balance symmetrically — is a clear and important conceptual point for the interpretation of Juno gravity data. The systematic exploration of 15 SL configurations (Fig. 6) and the demonstration that free-decay wind inversions remain feasible under non-adiabatic profiles (Fig. 8) are valuable. The comparison with GCM jet structures and ERA5 data (Fig. 11) provides useful physical context. The falsifiable prediction that winds must retain ~10% amplitude to ~3000 km depth (§3.4) is a concrete, testable constraint. However, the quantitative results depend on an unvaried parameter (k = 0.9∇_ad), which limits the strength of the quantitative claims.","major_comments":[{"comment":"§2.1, Eq. (1): The amplitude parameter k = 0.9∇_ad is fixed for all stable layers and never varied. This value is near the maximum physically possible subadiabatic departure, so the density enhancement within the SL is maximized. The paper's quantitative claims ('shallow, extensive stable layers substantially modify the background density, requiring more rapid decay of zonal winds') scale with k. If the real stratification is weaker (e.g., k = 0.1–0.3∇_ad, consistent with some estimates the authors cite in §2.1), the density modification and required wind adjustment diminish proportionally. The qualitative degeneracy claim holds for any k > 0, but the magnitude of the effect is load-bearing for the paper's central quantitative conclusion. At minimum, the authors should (i) run one or two cases with smaller k to demonstrate sensitivity, or (ii) explicitly reframe the results as an upper-~","section":null},{"comment":"§3.1, Fig. 6: Model uncertainty bars are shown only for the adiabatic profile, not for the SL profiles. The reader therefore cannot assess whether SL profiles achieve statistically equivalent fits to Juno data or merely visually acceptable ones. Since the paper's argument involves comparing how different SL configurations shift the harmonics relative to Juno's error bars, the absence of uncertainty estimates for the SL cases leaves the comparison incomplete. Adding at least wind-measurement-derived uncertainties for the extreme SL case (1–10^5 bar) would address this.","section":null},{"comment":"§2.1: The profiles omit a core and do not conserve mass ('inserting SLs does not strictly conserve mass'). The authors argue this is acceptable because they omit the low-order even harmonics (J_2, J_4, J_6, J_8) most affected by deep structure. However, the mass non-conservation could also affect the high-order harmonics that ARE used in the fits (J_10–J_20), particularly for the most extensive SLs. A quantitative estimate of the mass deficit/excess introduced by the extreme SL (1–10^5 bar) and its direct contribution to J_10–J_20 would strengthen the argument that this omission is benign.","section":null}],"minor_comments":[{"comment":"Abstract: 'inadicating' should be 'indicating'.","section":null},{"comment":"Fig. 4 y-axis label: ';-P profiles with SL' appears to be a rendering artifact; should read 'ρ-P profiles with SL'. Same issue in Fig. 5.","section":null},{"comment":"§2.1: The statement 'we neglect the effects of mean molecular weight gradients and latent heat release due to condensation' could briefly note that mean molecular weight gradients from helium rain are relevant at megabar pressures, which overlaps with some of the deeper SL configurations considered.","section":null},{"comment":"Fig. 11: The multiple y-axes (left, right-orange, rightmost-pink) make this figure difficult to parse. A clearer indication of which axis corresponds to which curves, or splitting into panels, would improve readability.","section":null},{"comment":"§3.3: The statement about super-adiabatic layers producing 'localized density deficits' and requiring 'slower decay of the zonal winds' is an interesting prediction but is stated qualitatively without any supporting calculation. A brief note that this is a conjecture by analogy, or a simple confirmation, would suffice.","section":null},{"comment":"Appendix B, Eq. (B1): The weight matrix w_ij is described as the inverse of the covariance matrix, but it would help to state explicitly whether off-diagonal covariance terms between harmonics are included, as this affects the optimization.","section":null},{"comment":"§2.1: The reference temperature T' = T_1bar (P/1bar)^0.25 uses an adiabatic scaling to evaluate ∇_ad at the SL midpoint, but the actual temperature at the midpoint will differ from this once the SL is imposed. A brief note on whether this inconsistency matters for the results would be helpful.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The core mathematical argument (Eq. 9 degeneracy) is sound and the qualitative claim survives any k > 0. The main concern is that the quantitative framing overstates the effect given the fixed near-maximal k. This is fixable by either a sensitivity test or reframing as an upper bound, which is why I recommend minor revision rather than major. The paper is a good fit for MNRAS given its methodology and relevance to Juno science. The citation pattern appears normal for this subfield."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. The central conceptual point — the degeneracy between ∂ρ_stat/∂z and ∂Q/∂z in the thermal wind balance — is acknowledged as valid, and we are grateful that the referee recognizes its importance. Below we address each major comment in turn. We agree that all three points warrant revision and will incorporate them in the revised manuscript.","responses":[{"response":"The referee is correct that the quantitative magnitude of the density modification, and hence the required wind adjustment, scales with k. We chose k = 0.9∇_ad as a deliberate end-member: it maximizes the subadiabatic departure and thus brackets the largest possible effect of a stable layer at a given location and extent. This choice was made to explore the full parameter space, but we agree that without sensitivity tests the reader cannot assess how the results scale for weaker stratification. In the revised manuscript we will: (i) run two additional cases with k = 0.1∇_ad and k = 0.3∇_ad for the most extreme SL configuration (1–10^5 bar), demonstrating that the density perturbation and required wind adjustment scale approximately linearly with k while the qualitative degeneracy persists for any k > 0; (ii) add a paragraph in §2.1 and §3.2 explicitly stating that k = 0.9∇_ad represents an upper bound on the effect, and that for weaker stratification the wind adjustment diminishes proportionally; and (iii) reframe the abstract and conclusions to clarify that the quantitative results correspond to the maximum-stratification case. We note that even for k = 0.1∇_ad, the degeneracy identified in Eq. 9 remains conceptually important — it is the symmetry between ∂ρ_stat/∂z and ∂Q/∂z that matters, not the amplitude — but we agree the paper should make clear that the magnitude of the wind adjustment is an upper limit.","revision_made":"yes","referee_comment":"§2.1, Eq. (1): k = 0.9∇_ad is fixed and never varied. Quantitative claims scale with k. Should run smaller-k cases or reframe as upper bound."},{"response":"This is a fair point. The uncertainty bars in Fig. 6 are derived from cloud-level wind measurement errors and are straightforward to propagate to the SL cases, since the same wind field is used. We omitted them from the SL curves for visual clarity, but we agree this leaves the comparison incomplete. In the revised manuscript we will add wind-measurement-derived uncertainty bands for at least the extreme SL case (1–10^5 bar), and ideally for all 15 SL profiles (as a shaded envelope if individual error bars are too cluttered). This will allow the reader to directly assess whether the SL profiles that deviate most from the adiabatic case remain within or exceed the observational uncertainties relative to Juno's error bars.","revision_made":"yes","referee_comment":"§3.1, Fig. 6: Model uncertainty bars shown only for the adiabatic profile, not for SL profiles. Cannot assess whether SL profiles achieve statistically equivalent fits."},{"response":"We agree that a quantitative estimate of the mass anomaly introduced by the SLs, and its direct contribution to the harmonics used in our fits, would strengthen the argument. We will add this calculation in the revised manuscript. Specifically, we will: (i) compute the integrated mass deficit/excess introduced by the extreme SL (1–10^5 bar) relative to the adiabatic profile, expressed both as an absolute mass and as a fraction of Jupiter's total mass; (ii) estimate the direct contribution of this mass anomaly to J10–J20 by computing the static (non-dynamical) harmonic contribution from the density difference between the SL and adiabatic profiles; and (iii) compare this contribution to Juno's measurement uncertainties on J10–J20 to demonstrate that it is subdominant. We expect this contribution to be small because the SL density anomalies are confined to low-pressure regions (≲1 Mbar) where the mass is a small fraction of the total, and the high-order harmonics weight the outer layers more heavily but the density anomalies are modest in absolute terms. However, we acknowledge that for the most extensive SLs this should be verified quantitatively rather than asserted, and we will do so. If the contribution turns out to be non-negligible for the extreme cases, we will note this as a limitation and discuss how it could be compensated by adjustments to the deeper density structure (e.g., core mass or heavy-element distribution) without affecting the dynamical harmonics that are the focus of this work.","revision_made":"yes","referee_comment":"§2.1: Profiles omit a core and do not conserve mass. Mass non-conservation could affect J10–J20 used in the fits. Need quantitative estimate of mass deficit/excess for extreme SL and its contribution to J10–J20."}],"tokens_in":22838,"tokens_out":1370,"duration_ms":56026,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The core result here is straightforward and correct: in the thermal wind balance (Eq. 9), the term ∂ρ_stat/∂z from a stable layer and the term ∂Q/∂z from wind decay contribute similarly to the dynamical density, so multiple wind profiles can reproduce the same Juno gravity harmonics when non-adiabatic layers are present. This is a genuine and useful point that prior gravity-inversion work has not explicitly addressed. The mathematical argument is clean, the numerical experiments confirm it, and the paper is honest about its limitations throughout. The systematic sweep over 15 stable-layer configurations is a real strength — it shows clearly that shallow, extensive layers matter while deep ones do not, which follows directly from where density is temperature-sensitive versus pressure-controlled. The experiment forcing jet decay within the stable layer (§3.4) is also well done: it shows that winds must retain ~10% amplitude down to ~3000 km regardless, which is a concrete, falsifiable constraint. The GCM and ERA5 comparisons (§4–5) are tangential to the main argument but do illustrate that no existing model uniquely matches the gravity solutions, which supports the non-uniqueness claim. The main soft spot is real but bounded. The stratification amplitude k = 0.9∇_ad is near the maximum physically possible subadiabatic departure, and it is never varied. The paper says this explicitly — k is fixed to isolate location and extent effects — but does not frame the quantitative results as an upper bound, which they effectively are. If the real stratification is weaker (say k = 0.1–0.3∇_ad, consistent with some cited estimates), the density modification and required wind adjustment shrink proportionally. The qualitative degeneracy survives for any k > 0, but the claim that shallow extensive layers 'substantially modify' the background density is specific to this near-maximal choice. A sensitivity test varying k would have addressed this directly. Two minor issues: error bars are shown only for the adiabatic profile (Fig. 6), so the reader cannot assess whether SL profiles achieve statistically equivalent fits or merely visually acceptable ones; and the profiles omit a core and do not conserve mass, though this is defensible since the paper deliberately excludes the even harmonics that would be most affected. The stress-test concern about k is valid and is the most important point a referee should raise. The reader's 'conditional' verdict is slightly too cautious — the degeneracy result is mathematically sound and the paper earns its conclusion. But the reader is right that the quantitative details should be treated as illustrative. This paper is for researchers working on Juno gravity inversion and giant-planet interior dynamics. It deserves a serious referee who should ask for (1) a k-sensitivity test or at minimum explicit framing of results as upper bounds, and (2) error bars on the SL profile fits. With those addressed, this is a solid contribution to the subfield.","headline":"The degeneracy between stable-layer density structure and wind decay in Jupiter's gravity inversion is real and clearly demonstrated; the quantitative magnitude depends on a near-maximal stratification amplitude that is never varied.","tokens_in":24034,"tokens_out":697,"would_cite":true,"duration_ms":78111,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Stable layers on Jupiter deepen the mystery of how far its jet streams reach","keywords":[],"falsifier":"If future measurements of Jupiter's actual Brunt–Väisälä frequency (e.g., from polar cyclone depths or seismology) show stratification strengths far weaker than k = 0.9∇_ad, the density perturbations from stable layers would be much smaller, the degeneracy with wind depth would shrink, and the adiabatic wind-decay solutions would remain approximately valid.","tokens_in":23294,"feed_emoji":"🌀","tokens_out":894,"duration_ms":117746,"temperature":0.7,"pith_summary":"This paper argues that if Jupiter has stable, non-convective layers in its outer envelope — as several recent observations suggest — then the vertical profile of its jet streams cannot be uniquely determined from gravity data alone. The core mechanism is a mathematical degeneracy in the thermal wind balance (Eq. 9): the gravity signal produced by increasing the radial gradient of the background density (via a stable layer) is nearly identical to the signal produced by increasing the radial gradient of the wind decay function (deeper winds). Because Juno's gravity harmonics measure the combined product of these two quantities, a stable layer that steepens the density gradient forces the winds to decay faster to produce the same observed signal. The paper demonstrates this by constructing 15 density profiles with prescribed stable layers at various depths and pressures, computing the resulting gravitational harmonics, and then optimizing the wind decay profile to match Juno measurements. Shallow, extensive stable layers (e.g., 1–10⁵ bar) substantially increase the local density, requiring winds to drop to less than half their cloud-level strength within the top 100 km — far faster than in adiabatic models. Deep stable layers (above ~1 Mbar) have negligible effect on density and thus on the gravity-wind degeneracy. The paper also tests a physically motivated scenario where jets decay within the stable layer itself, finding that this generally fails to reproduce Juno data because winds must retain at least ~10% of their cloud-level amplitude down to ~3000 km depth. Comparisons with Jupiter GCMs and Earth's eddy-driven jets show that no existing dynamical model produces a decay profile that uniquely matches the gravity-inversion envelope, especially once stable layers are admitted.","feed_headline":"Jupiter's jet depth can't be pinned down if stable layers exist","feed_subtitle":"Stable atmospheric layers mimic the gravity signal of deep winds, creating a degeneracy that Juno data alone cannot break.","key_machinery":"Thermal wind balance (Eq. 7–9) on a giant planet; prescribed subadiabatic stable layers (Eq. 1) modifying the static density profile; gravitational harmonic computation (Eq. 4–6) comparing modeled vs. Juno-observed values; constrained (Eq. 10) and free (Eq. 11–12) wind-decay optimization (Eq. B1) to fit odd and high-order harmonics.","core_discovery":"The paper identifies a direct mathematical degeneracy between two terms in the thermal wind equation (Eq. 9): the vertical gradient of static density (∂ρ_stat/∂z, modified by stable layers) and the vertical gradient of the wind decay function (∂Q/∂z, controlled by how deep the winds penetrate). Both terms enter the dynamical density anomaly identically, so increasing one can compensate for decreasing the other. This means that a stable layer that steepens the density profile produces the same gravitational harmonic signature as deeper-penetrating winds. Consequently, the vertical structure of Jupiter's jet streams — a key target of the Juno mission — cannot be uniquely inferred from gravityh","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Stable layers obscure the true depth of Jupiter's jet streams","Jupiter's jet depth remains hidden due to stable layer degeneracy","Juno gravity data cannot uniquely resolve Jupiter's jet stream depth","Stable layers weaken constraints on Jupiter's deep jet stream structure","Jupiter's jet stream depth is unconstrained without internal stability data"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The stable layers are prescribed with an artificial temperature gradient (Eq. 1) using a fixed amplitude parameter k = 0.9∇_ad rather than derived from opacity or convection calculations, and the authors note that the actual stratification strength within Jupiter remains poorly constrained. The entire quantitative result — how much the winds must adjust — depends on this arbitrary choice. If the real stratification is much weaker, the degeneracy shrinks toward the adiabatic (","fun_headline_variants_meta":{"raw":{"variants":["Stable layers obscure the true depth of Jupiter's jet streams","Jupiter's jet depth remains hidden due to stable layer degeneracy","Juno gravity data cannot uniquely resolve Jupiter's jet stream depth","Stable layers weaken constraints on Jupiter's deep jet stream structure","Jupiter's jet stream depth is unconstrained without internal stability data","A stable layer creates a mathematical blind spot for Jupiter's jet depth"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":1113,"prompt_tokens":599,"completion_tokens":514,"prompt_tokens_details":null},"tokens_in":599,"tokens_out":514,"duration_ms":20357,"temperature":1.0,"reasoning_tokens":462,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T12:04:34.276600+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If future measurements of Jupiter's actual Brunt–Väisälä frequency (e.g., from polar cyclone depths or seismology) show stratification strengths far weaker than k = 0.9∇_ad, the density perturbations from stable layers would be much smaller, the degeneracy with wind depth would shrink, and the adiabatic wind-decay solutions would remain approximately valid.","supporting_citations":[],"review_version":1}