{"id":"008260db-aacf-430a-95b2-bb891c77cccc","arxiv_id":"2509.02759","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Jupiter's upper-troposphere vertical motions inferred from the transformed Eulerian mean circulation have upwelling in zones and downwelling in belts, matching the observed cloud structure and resolving the apparent contradiction with Eulerian mean dynamics.","lead":"This paper argues that Jupiter's belt-zone cloud pattern, with upwelling in bright zones, is compatible with the momentum-flux constraints once the flow that transports tracers is recognized as the transformed Eulerian mean circulation, which includes an eddy-driven Stokes drift. Using observed winds and an eddy diffusion closure, the authors infer upwelling in zones and downwelling in belts, and a global circulation model supports the picture.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"TEM sign depends on the unmeasured constant-K downgradient PV closure; a GCM-based local-K diagnostic can test whether this assumption is load-bearing.","rationale":"The paper does something strong: it uses TEM theory and a PV-mixing closure to resolve a real observational contradiction, and the derivation from Eq. (10) to Eq. (12) is internally consistent under the stated QG assumptions. The acknowledged limitations—constant K, neglect of stretching, and transport barriers—are exactly where the sign of the vertical velocity could fail. The reader identified the downgradient closure as the weakest assumption; my concrete test sharpens this by using the GCM's own eddy fluxes to compute a local K rather than assuming one. The GCM support is real evidence but does not close the issue, because the GCM's TEM/Eulerian magnitude ratio differs from observations, so the GCM is not a direct check of the observational K estimate. Therefore the CONDITIONAL verdict stands unchanged.","tokens_in":18709,"tokens_out":5087,"duration_ms":57186,"concrete_test":"Using the published Young et al. (2019a,b) GCM output, at ~200 mbar and |φ| > 10°, diagnose the local eddy diffusivity K(φ) = −v′q′ / ∂y q̄ from the model's own eddy PV flux and mean PV gradient. Then reconstruct v̄* two ways: (a) with the constant observational K = 6×10^6 m²/s and (b) with the diagnosed K(φ), and compare the implied w̄* = −(H/2)(1/(a cosφ)) ∂φ(v̄* cosφ) patterns. If the zone/belt sign pattern of w̄* is identical in both reconstructions, the constant-K assumption is not load-bearing; if it changes, the central reconciliation is conditional on the closure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The crux is the chain in §2.2: Eq. (10) gives f v̄* = −v′q′, and the downgradient closure Eq. (11), v′q′ ∼ −K ∂y q̄ with constant K, turns this into v̄* ≈ (K/f) ∂y q̄ (Eq. 12). The headline claim that w̄* is positive in zones and negative in belts follows from the convergence/divergence of this v̄*. Neither K nor v′q′ is directly measured on Jupiter; the observed eddy momentum flux divergence S constrains only the barotropic part of v′q′, not the heat-flux part entering the TEM correction. The authors explicitly assume K varies on scales large compared with the jets and neglect transport barriers at jet cores, which can suppress K precisely where ∂y q̄ is enhanced. If K has local minima at prograde jet cores, the 'more strongly poleward in prograde jets' pattern weakens or reverses, and the inferred sign of w̄* in zones/belts need not hold. The GCM is suggestive but not a direct validation of the observational closure: it produces the same sign pattern while its TEM and Eulerian meridional flows are comparable in magnitude, not an order apart, so its effective K differs from the assumed value.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the apparent contradiction between two inferences about Jupiter's upper troposphere: eddy momentum flux convergence implies an Eulerian mean meridional circulation with downwelling in zones and upwelling in belts, while cloud and haze structure suggests upwelling in zones and downwelling in belts. The authors argue that the relevant tracer-transport circulation is the transformed Eulerian mean (TEM), not the Eulerian mean. Using the observed zonal wind profile, a downgradient potential vorticity flux closure v'q' ~ -K ∂y q̄ with a constant eddy diffusivity K, and neglecting the stretching contribution to the mean PV gradient, they infer a TEM meridional flow that is generally poleward and an order of magnitude stronger than the Eulerian mean. Its convergence/divergence pattern then implies TEM upwelling in zones and downwelling in belts, resolving the contradiction. The arguments are supplemented by GCM simulations from Young et al. (2019a,b), which qualitatively reproduce the sign pattern of the TEM vertical velocity.","tokens_in":18980,"tokens_out":7357,"duration_ms":90906,"significance":"If the central claim is correct, it is an important step: it would replace the apparent contradiction with a circulation whose vertical branches match the observed cloud structure and which is the relevant mean flow for ammonia and other tracers. The paper has real strengths: the scaling argument is transparent, it makes explicit use of observed Cassini eddy momentum fluxes and zonal winds, and the authors candidly acknowledge the main caveats, including the uncertain eddy diffusivity and the limitations of the GCM. The GCM data are made available. The result is also falsifiable in principle through future tracer or heat-flux observations. However, the chain from Eqs. (10)-(12) is the load-bearing part, and it currently rests on an unmeasured closure; the GCM provides only qualitative support. The manuscript is therefore promising but not yet conclusive.","major_comments":[{"comment":"The headline claim—TEM vertical motion is positive in zones and negative in belts and an order of magnitude stronger than the Eulerian mean—rests entirely on the downgradient closure v′q′ ∼ −K∂y q̄ with constant K. Neither v′q′ nor K is measured on Jupiter. If K has local minima at prograde jet cores, as the paper acknowledges from transport-barrier studies, the product K∂y q̄ need not peak there; the convergence/divergence of v̄* and hence the sign of w̄* can change in individual bands. The GCM does not provide the missing test: its TEM and Eulerian meridional flows are comparable in magnitude, and its eddy velocities are a factor ~3 smaller, so its effective K differs from the 6×10^6 m^2 s^-1 assumed in Fig. 1. Please report local-K diagnostics from the GCM, or an observational constraint on v′q′, before the central sign claim is considered established.","section":"§2.2, Eqs. (10)–(12)"},{"comment":"The GCM is the only independent check of the scaling, yet Fig. 4 shows that the simulated TEM and Eulerian mean meridional flows are of similar magnitude, not an order apart. The text attributes this to weaker eddy amplitudes and smaller S, but that means the model does not reproduce the assumed K or the claimed dominance of the Stokes-drift contribution. To validate the observational inference, the model should be diagnosed to show whether the sign pattern of w̄* follows from the same balance (Eq. 10) and to identify why the magnitude ratio differs. Without this, the GCM supports qualitative sign consistency but not the 'order of magnitude' part of the central claim.","section":"§3.1.1, Fig. 4"},{"comment":"The potential vorticity gradient in Eq. (8) contains a stretching term that is discarded on the authority of Read et al. (2006)'s ~10% estimate. This is defensible as a first-order scaling, but it becomes load-bearing where the barotropic part ∂y(f+ζ) is small, e.g., near retrograde jet cores and band edges, where the sign of ∂y q̄ may be controlled by the stretching term. The paper should either quantify the stretching contribution from the GCM temperature field at the levels shown in Fig. 1, or restrict the sign claim to latitudes where the barotropic term clearly dominates.","section":"§2.2, Eq. (8)"}],"minor_comments":[{"comment":"The eddy diffusivity is written as K ∼ V L ∼ 6 × 10^6 m s^-2; the units should be m^2 s^-1. The same typo appears in the Fig. 1 caption.","section":"§2.2, after Eq. (12)"},{"comment":"In the list of tracers, 'hydrogen sulfilde (H2S)' should be 'hydrogen sulfide (H2S)'.","section":"§3"},{"comment":"The S/I = 0.05 threshold is described as chosen for continuity with the standard TEM at 200 mbar. Since Fig. 5 extends the residual circulation to depth using this choice, a brief sensitivity test or a statement of how the deep pattern changes with the threshold would aid robustness.","section":"§3.1.1, Eqs. (13)–(15)"},{"comment":"The four panels use two vertical axes with different scales and line styles; a small legend or explicit axis labels for the orange and dashed-blue curves in each panel would improve readability.","section":"Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The validation GCM and several of the theoretical foundations come from the same research group (Schneider and Liu 2009; Liu and Schneider 2010; Young et al. 2019a,b). This is not a disqualifier, but it means the independent confirmation of the central closure is weaker than it appears. An editor may wish to seek a referee who can assess the TEM/potential-vorticity closure in detail."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper genuinely does something new. It shows the old contradiction between cloud-implied upwelling in Jupiter's zones and the Eulerian-mean downwelling inferred from momentum flux convergence can be dissolved once you look at the transformed Eulerian mean (TEM) circulation rather than the Eulerian mean. Using the observed zonal wind and a constant-K downgradient PV diffusion closure, the TEM meridional flow is more strongly poleward in prograde jets, which yields TEM upwelling in zones and downwelling in belts—about an order of magnitude stronger than the Eulerian mean. That specific inference is new relative to the earlier TEM work by Gierasch et al. (1986), Lee & Kaspi (2021), and Duer et al. (2021), and it's a useful, clean idea.\n\nThe paper is also honest in method. The scaling argument in §2.2 is transparent, the limitations are listed rather than buried, and the GCM testbed is used the way it should be: the authors note that the simulated eddy amplitudes are too weak and that the model's TEM and Eulerian meridional flows are comparable in magnitude, not an order apart. That is the right way to present a supporting model.\n\nThe soft spot is the load-bearing closure in Eq. (11). The sign of the TEM vertical velocity follows directly from multiplying the absolute vorticity gradient by a positive constant K. K is not measured on Jupiter, and the transport-barrier literature says K can drop sharply at prograde jet cores—exactly where ∂y q̄ is largest. If K has local minima there, the 'strongly poleward in prograde jets' pattern weakens or reverses, and with it the inferred sign of w*. The authors acknowledge this in prose but do not test it. They also drop the stretching term on the grounds that it's about 10%, which is plausible but not checked. The GCM gives same-sign support, but because its effective K evidently differs from the assumed constant, it does not verify the quantitative 'order of magnitude stronger' claim. These are real limitations, though they should not be overstated: the paper is a testable scaling hypothesis, and at present it is the most coherent picture we have.\n\nRead this if you work on giant-planet belt/zone dynamics, ammonia transport, or TEM theory applied to planets. It deserves a serious referee, and I'd send it out—with a request that the authors stress-test the K closure, ideally using a latitude-dependent diffusivity or a local-K diagnostic from the GCM.","headline":"A clean TEM-scaling argument that likely resolves Jupiter's upwelling-downwelling contradiction, but the sign and magnitude rest on a constant-K PV diffusion closure that the paper doesn't yet stress-test.","tokens_in":19494,"tokens_out":4105,"would_cite":true,"duration_ms":46860,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Jupiter's bright zones are upwelling after all — once the eddies' Stokes drift is counted in, the paper argues, the cloud bands and the dynamics agree.","keywords":["Jupiter","transformed Eulerian mean","atmospheric circulation","potential vorticity","belts and zones","eddy fluxes","ammonia clouds","general circulation model"],"falsifier":"Measure, at the cloud level and outside the equatorial jet, the meridional eddy flux of potential temperature (equivalently, the baroclinic part of the Eliassen-Palm flux). If that flux is too small for the Stokes-drift term to dominate the Eulerian mean meridional flow, or if its sign pattern is not what downgradient diffusion of potential vorticity implies, the claimed TEM upwelling in zones would not hold. A more direct falsification would be a cloud-tracking or tracer-inversion estimate of vertical velocity showing downwelling in the bright zones.","tokens_in":18603,"feed_emoji":"🪐","tokens_out":7566,"duration_ms":76558,"temperature":0.7,"pith_summary":"Jupiter's cloud bands have long posed a puzzle: the bright zones appear to be upwelling and the dark belts downwelling, yet the momentum balance of the observed eddies implies the opposite sense of vertical motion. This paper argues that the contradiction dissolves when tracers are tracked by the transformed Eulerian mean (TEM) flow — the Eulerian mean plus an eddy-driven Stokes drift — rather than by the Eulerian mean alone. Using potential vorticity gradients implied by the observed winds and a mixing-length diffusivity, the authors estimate TEM upwelling in zones and downwelling in belts, about an order of magnitude stronger than the Eulerian mean vertical velocity. A global circulation model of Jupiter's upper troposphere reproduces the pattern. If right, the work reconciles dynamics with cloud appearance and identifies turbulent eddies, not the mean flow, as the main agents of vertical tracer transport.","feed_headline":"Jupiter's bright zones upwell once eddy drift is counted","feed_subtitle":"Counting the eddies' Stokes drift turns belt-zone downwelling into upwelling, resolving a long-standing contradiction.","key_machinery":"The load-bearing object is the transformed Eulerian mean (TEM) circulation, an approximation of the Lagrangian mean flow that adds a Stokes-drift term built from meridional eddy heat fluxes to the Eulerian mean. Its key identity, f v̄* = −v′q′, sets the TEM meridional flow by the eddy potential vorticity flux rather than the eddy momentum flux divergence; the TEM vertical flow follows from continuity. Combined with the downgradient diffusive closure v′q′ ∼ −K ∂yq̄, with diffusivity K estimated from observed eddy scales, and with Jupiter's mean potential vorticity gradient being barotropically dominated, this makes the TEM meridional flow converge in belts and diverge in zones — hence upwelli","core_discovery":"The paper's central claim is that the apparent contradiction between Jupiter's cloud bands and its inferred vertical motions is an artifact of averaging. Nearly conserved tracers such as ammonia are advected not by the Eulerian mean flow but by the transformed Eulerian mean (TEM) circulation, which adds the eddies' Stokes drift. From the potential vorticity structure implied by the observed jets and a downgradient diffusive closure for the eddy potential vorticity flux, the TEM meridional flow is poleward almost everywhere, strongest in prograde jets; continuity then gives TEM upwelling in zones and downwelling in belts, with magnitudes about an order of magnitude stronger than the Eulerian","pith_inferences":["If the TEM picture is right, the vertical motion that shapes Jupiter's ammonia cloud bands is largely a tracer-transport (Lagrangian) motion rather than a mass motion, and the Eulerian circulation could be nearly its reverse; a direct test would be measuring the meridional eddy heat flux at cloud level, since a small Eliassen-Palm contribution would make the TEM and Eulerian circulations coincide ","The same argument, with adjusted parameters, should apply to Saturn and other giant planets whose banded cloud structure is also anticorrelated with the Eulerian-mean circulation inferred from eddy momentum fluxes; the mechanism predicts TEM upwelling in their bright zones as well.","The constant-diffusivity assumption yields a specific, testable relation between jet curvature and vertical velocity: upwelling should be strongest where the mean potential vorticity gradient is most positive (prograde jet cores) and downwelling where it is most negative; cloud-top ammonia contrast maps could be compared quantitatively with this pattern.","Because the TEM vertical velocity is an order of magnitude larger than the Eulerian one, global tracer budgets for Jupiter — for example, ammonia supply from depth or lightning occurrence in belts — may need to be re-evaluated in terms of Lagrangian-mean transport rather than Eulerian mean motion."],"forward_implications":["The apparent conflict between eddy-momentum-based circulation inferences and cloud appearance in Jupiter's belts and zones is resolved: the eddy-driven (Stokes drift) component of the flow reverses the sign of vertical motion relative to the Eulerian mean.","Transient tracers such as ammonia in the upper troposphere are transported by a circulation whose upwelling regions coincide with the bright zones where ammonia ice clouds form.","The TEM vertical velocity below the cloud layer is about an order of magnitude larger than the Eulerian mean value, so turbulent eddies, not the mean meridional circulation, dominate vertical tracer transport in the upper troposphere.","The circulation is driven from the top down: shallow upper-tropospheric eddy fluxes force overturning cells that close in a deep frictional layer, consistent with the downward-control principle.","Simulations that reproduce the observed jet structure and energetics produce the same belt-zone pattern of TEM vertical motion, making the inference more than a scaling argument."],"supporting_citations":[{"why":"Supplies the observed zonal wind profile, eddy momentum flux divergence, and eddy velocity scales used to construct both the Eulerian and TEM estimates.","marker":"Salyk et al., 2006"},{"why":"Provides the transformed Eulerian mean / Lagrangian-mean theory on which the central argument rests.","marker":"Andrews and McIntyre (1976, 1978)"},{"why":"Gives the leading-order zonal momentum balance and the energetic constraint that eddy momentum fluxes are shallow, grounding the downward-control closure.","marker":"Schneider and Liu (2009)"},{"why":"The downward-control principle that lets the shallow eddy forcing be converted into deep vertical motion of the TEM circulation.","marker":"Haynes et al., 1991"},{"why":"Observed potential vorticity structure of Jupiter's upper troposphere, the basis for treating the mean PV gradient as barotropically dominated.","marker":"Read et al., 2006"},{"why":"Formulates the Eliassen-Palm flux divergence that connects the TEM momentum balance to the eddy potential vorticity flux.","marker":"Edmon et al., 1980"},{"why":"Provides the modified residual circulation used to extend TEM estimates below the stably stratified upper troposphere.","marker":"Held and Schneider (1999)"},{"why":"The global circulation model simulations used to test the scaling arguments and reproduce the ammonia ice band structure.","marker":"Young et al. (2019a,b)"},{"why":"The observation that eddy momentum fluxes converge in prograde jets, which sets up the original contradiction.","marker":"Ingersoll et al., 1981"},{"why":"The baroclinic eddy generation mechanism and energetics justifying the shallow eddy driving of Jupiter's jets.","marker":"Liu and Schneider (2010)"}],"fun_headline_variants":["Stokes drift flips Jupiter’s downwelling into upwelling","Jupiter’s zones upwell when eddy drift is included","Transformed Eulerian motion solves Jupiter band paradox","Eddy drift reveals upwelling in Jupiter’s bright zones","Counting Stokes drift reconciles Jupiter cloud bands"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"Everything hinges on the assumption that eddy potential vorticity fluxes are downgradient and diffusive, with a single constant eddy diffusivity; if real Jupiter fluxes are counter-gradient, or if the jet cores act as strong transport barriers that shrink the diffusivity there, the computed sign of the vertical motion in zones and belts could change.","fun_headline_variants_meta":{"raw":{"variants":["Stokes drift flips Jupiter’s downwelling into upwelling","Jupiter’s zones upwell when eddy drift is included","Transformed Eulerian motion solves Jupiter band paradox","Eddy drift reveals upwelling in Jupiter’s bright zones","Counting Stokes drift reconciles Jupiter cloud bands"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000467,"raw_usage":{"total_tokens":2164,"prompt_tokens":740,"completion_tokens":1424,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":1353}},"tokens_in":484,"tokens_out":1424,"duration_ms":11733,"temperature":1.0,"reasoning_tokens":1353,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:24:44.048307+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, at the cloud level and outside the equatorial jet, the meridional eddy flux of potential temperature (equivalently, the baroclinic part of the Eliassen-Palm flux). If that flux is too small for the Stokes-drift term to dominate the Eulerian mean meridional flow, or if its sign pattern is not what downgradient diffusion of potential vorticity implies, the claimed TEM upwelling in zones would not hold. A more direct falsification would be a cloud-tracking or tracer-inversion estimate of vertical velocity showing downwelling in the bright zones.","supporting_citations":[{"cited_title":", author Ingersoll, A.P","cited_arxiv_id":null,"evidence_quote":"Supplies the observed zonal wind profile, eddy momentum flux divergence, and eddy velocity scales used to construct both the Eulerian and TEM estimates."},{"cited_title":", author McIntyre, M.E","cited_arxiv_id":null,"evidence_quote":"Provides the transformed Eulerian mean / Lagrangian-mean theory on which the central argument rests."},{"cited_title":", author Liu, J.J","cited_arxiv_id":null,"evidence_quote":"Gives the leading-order zonal momentum balance and the energetic constraint that eddy momentum fluxes are shallow, grounding the downward-control closure."},{"cited_title":", author Marks, C.J","cited_arxiv_id":null,"evidence_quote":"The downward-control principle that lets the shallow eddy forcing be converted into deep vertical motion of the TEM circulation."},{"cited_title":", author Gierasch, P.J","cited_arxiv_id":null,"evidence_quote":"Observed potential vorticity structure of Jupiter's upper troposphere, the basis for treating the mean PV gradient as barotropically dominated."},{"cited_title":", author Hoskins, B.J","cited_arxiv_id":null,"evidence_quote":"Formulates the Eliassen-Palm flux divergence that connects the TEM momentum balance to the eddy potential vorticity flux."},{"cited_title":", author Beebe, R.F","cited_arxiv_id":null,"evidence_quote":"The observation that eddy momentum fluxes converge in prograde jets, which sets up the original contradiction."},{"cited_title":", author Schneider, T","cited_arxiv_id":null,"evidence_quote":"The baroclinic eddy generation mechanism and energetics justifying the shallow eddy driving of Jupiter's jets."}],"review_version":1}