{"id":"b1bf0fea-b4b4-48c2-b954-644928076154","arxiv_id":"2507.23493","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Vorticity jet interfaces in beta-plane turbulence show sub-Gaussian height fluctuations, fat-tailed speeds, rough non-differentiable time series, and multifractal structure that smoothens as the zonostrophy parameter grows.","lead":"By tracking the wavy boundaries between fast and slow bands in computer-simulated planetary jets, the authors measure how these boundaries wiggle over time. The wiggles are mostly small and smooth, but occasionally jump; they have a rough, fractal structure that depends on how strongly rotation organizes the jets, which matters for how these flows spread particles and pollutants.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Interface-tracking rule is only visually validated; since every reported statistic is computed on those tracked contours, the quantitative laws could be artifacts of the tracking algorithm.","rationale":"The reader's weakest-assumption analysis and mine converge on the same point: the identification of vorticity interfaces with local maxima of u is the load-bearing step, and it is validated only by visual inspection (animations, Fig. 2c). This is not a matter of disagreeing with the field's consensus; it is an internal methodological gap. The paper's own language ('It appears', 'seem') is honest, but honesty does not supply the missing quantitative validation. Everything that follows—height field, speeds, spectra, structure functions, and multifractal spectra—is a statistic of the tracked contours, so any tracking bias propagates to all reported exponents. The consistency between N=1024 and N=2048 is useful evidence of numerical convergence, but it does not validate the physical correspondence between velocity maxima and vorticity jumps, since both resolutions use the same tracking criterion. The proposed test is deliberately concrete: replace the tracking rule with an independent, physically motivated definition and check whether the conclusions survive. Because the concern is real but does not by itself force rejection—it could be resolved by quantitative validation—the appropriate verdict remains conditional, as the reader concluded. No change to the reader's verdict is needed.","tokens_in":11823,"tokens_out":4388,"duration_ms":50771,"concrete_test":"On the same DNS snapshots, track C_alt(x,t) by an independent interface definition—e.g., the y-position of the maximum of |∂ω/∂y| at fixed x, or the zero of ∂²u/∂y² closest to each velocity maximum. Quantify the separation between C and C_alt; then recompute P(h), P(vC), the four spectra, S2/ζ2, and D_q using C_alt. If any reported exponent or tail shape shifts by more than the statistical uncertainty (bootstrap over sectors/snapshots), the central laws are tracking artifacts; if they coincide, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on identifying the vorticity interface C(x,t) with the y-location of the local maximum of u(x,y,t) at each x,t (Sec. II). The paper's own justification is observational: 'It appears that the vorticity interfaces correspond to the local maxima...' and 'quite remarkable to notice from an animation.' No quantitative comparison against an independent interface definition is made. The entire statistical characterization—height PDFs, speed distributions, f^-2/f^-5/f^-5/3/f^-7/2 spectra, structure-function exponents ζ2, and multifractal D_q—is a set of functionals of C(x,t). If the algorithm occasionally locks onto a transient eddy maximum, jumps between neighbouring maxima as a jet wobbles, or fails to maintain contour identity when jets merge, then the supposed roughness (ζ2<2), the shock-like events, and the fat-tailed speeds are properties of the tracking rule rather than of the physical vorticity boundaries. The paper does not state how multiple candidate maxima are resolved, how the time-averaged profile Cbar in Eq. (3) is defined when a contour vanishes, or how merging/splitting is handled, so these failure modes are not excluded.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies fluctuating vorticity interfaces in forced-dissipative barotropic beta-plane turbulence by direct numerical simulation. The authors identify each interface with the local maximum of the zonal velocity u(x,y,t) along y at fixed x, define a height field h(x,t) as the deviation from the time-averaged contour, and then characterize its statistics: height PDFs with sub-Gaussian tails, heavy-tailed PDFs of fluctuation speeds, power-law frequency spectra of sectorial mean and variance, a second-order temporal structure function S2(τ) with a short-time exponent ζ2<2 that approaches 2 with increasing zonostrophy parameter Rβ, and a multifractal spectrum Dq that narrows as Rβ increases. The authors interpret these results as quantitative statistical laws for jet-interface fluctuations and suggest relevance for transport of microplastics and other particles.","tokens_in":12060,"tokens_out":3140,"duration_ms":36073,"significance":"If the reported exponents are robust, this would be a useful systematic characterization of interface fluctuations in a classical geophysical turbulence model, and the paper adds to the relatively sparse literature on fluctuating jet boundaries. The DNS setup is standard and clearly described: N=1024 and 2048 pseudo-spectral runs, hyperviscosity, narrow-band forcing, five values of β, and 1000 statistically independent snapshots for statistics. The paper is also honest that the tracking rule is motivated by visual inspection. The main novelty lies in treating the jet edge as a tracked random curve and extracting scaling laws; however, because every reported statistic is a functional of the tracked contour, the lack of a quantitative validation of the tracking rule is a load-bearing weakness. The manuscript would be strengthened substantially by error bars and fit ranges for all exponents, a quantitative comparison with independent interface definitions, and a clarification of how contour identity is maintained under merging or disappearance.","major_comments":[{"comment":"The entire analysis rests on identifying the vorticity interface with the local maximum of u(x,y,t) at each x and t. The paper's justification is visual and qualitative: 'It appears that the vorticity interfaces correspond to the local maxima...' and 'quite remarkable to notice from an animation.' No quantitative comparison is made against an independent definition, such as maxima of |∂ω/∂y| or level sets of the zonally filtered vorticity. If the algorithm occasionally locks onto a transient eddy maximum, jumps between neighboring maxima during jet wobble, or fails when a contour shrinks and disappears, then the reported sub-Gaussian PDFs, spectral exponents, ζ2<2, and multifractal widths may be properties of the tracking rule rather than of the physical vorticity boundaries. Please provide a quantitative validation, at least for selected snapshots and Rβ values, and state how multiple candidate maxima are resolved and how contour identity is maintained over time.","section":"Section II, contour identification"},{"comment":"The height field h(x,t)=C(x,t)−Cbar(x) assumes that C(x,t) is a single-valued, continuously trackable function of x and t. The manuscript does not state what happens when a jet interface weakens, merges with another jet, or disappears locally, nor how Cbar(x) is computed in regions where the contour is absent. Without such a procedure, the time average in Eq. (3) is not well defined for fragmented contours, and the PDFs and structure functions may be biased by arbitrary splicing. Please specify the tracking and identity-maintenance algorithm explicitly.","section":"Equation (3) and Section II"},{"comment":"The abstract states that the authors calculate 'the moments of the time-increments of the interfacial height fluctuations,' but the paper reports only the second-order structure function S2(τ) in Eq. (5) and Fig. 5(a). No higher-order moments are shown, so the claim of roughness and the discussion of intermittency are based on a single moment. Please either add higher-order structure functions or revise the abstract and the corresponding statements to refer specifically to the second-order moment.","section":"Abstract and Section V (structure functions)"},{"comment":"The reported spectral exponents f^-2, f^-5, f^-5/3, and f^-7/2, as well as the structure-function exponents ζ2 and ξ2≈1.1, are presented without fit ranges, confidence intervals, or error bars. In particular, the f^-5 and f^-7/2 regimes are asserted from log-log plots with a short frequency range, and the text gives no criterion for what counts as a scaling range. Please provide a table of all exponents with the fitted frequency/time ranges, the number of independent samples, and an uncertainty estimate (for example, bootstrap over sectors or over independent time blocks).","section":"Fig. 4(d,e) and Fig. 5(a)"},{"comment":"The multifractal analysis defines H(x,y,t)=|h(x,y,t)|+0.001 and then computes partition functions over segments of the interfacial contour. There are two problems. First, h is a function of x and t, so the notation H(x,y,t) is inconsistent and the coarse-graining sum over x is ambiguous if H is meant to depend on y. Second, the additive constant 0.001 is a free parameter that strongly affects Zq(l) for negative q, where small values of H dominate; the reported broadening of Dq for negative q may therefore be an artifact of this offset. Please clarify the definition, test sensitivity to the offset, and consider using a measure that is not contaminated by an arbitrary additive constant.","section":"Section VI, multifractal analysis"},{"comment":"Several central claims—ζ2→2 with increasing Rβ, the narrowing of the multifractal spectrum, and the approach of the speed PDF to Gaussian—are described qualitatively rather than numerically. The paper would be significantly more convincing if, for each Rβ, the authors reported the measured ζ2, ξ2, the width of the singularity spectrum, and a measure of the tail of P(vC), together with their statistical uncertainties. This would allow the reader to assess the Rβ dependence claimed in the text.","section":"General reporting of exponents"}],"minor_comments":[{"comment":"There are several typographical errors: 'multifractactal' in the last line of the introduction, 'fluctuationspeeds' two sentences later, and 'Ekmann' for Ekman in the model paragraph. Please proofread the text.","section":"Introduction"},{"comment":"The time average in Eq. (3) is written as an integral over T, but the practical computation uses 1000 discrete snapshots. Please clarify whether the average is over the same snapshots used for the PDFs and spectra, and whether the integration is approximated by a sum.","section":"Section II, Eq. (3)"},{"comment":"The caption states Rβ = 7.5 (β = 50) for all three panels, but the text and Fig. 1 present β=10 as Rβ=6.3. Please check that the caption in Fig. 2 is consistent with the actual parameters of the snapshot.","section":"Fig. 2 caption"},{"comment":"The partition function is written with a sum over i from 0 to Nl−1, but the text uses Hl,i = Σ_{x=li}^{l(i+1)} H(x,y); the lower limit should be li or l(i) with a clearly defined integer index. Please rewrite this equation with unambiguous notation.","section":"Section VI, Eq. for partition function"}],"recommendation":"major_revision","confidential_remarks":"The core idea is interesting and the simulation campaign is adequately described, but the central measurement object—the tracked interface—is validated only visually and several reported exponents lack error bars and fit ranges. These issues are fixable within the scope of the manuscript. I also recommend that the authors explicitly discuss the relation to their previous work on multifractal interfaces in active suspensions (ref. [37]), since the multifractal methodology appears to be identical and the novelty here is the application to beta-plane jets."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a competent, clearly written DNS study that gives the first systematic statistical characterization of vorticity interfaces in forced-dissipative barotropic beta-plane turbulence, including R_beta dependence. The object is new: no previous work in the cited literature tracks these vorticity shear zones and measures their height PDFs, spectra, structure functions, and multifractal spectra. The DNS setup is standard and well reported, with two resolutions and clear parameters. The concrete exponents — f^-2, f^-5, f^-5/3, f^-7/2, zeta2 < 2 approaching 2 with R_beta, and the narrowing multifractal spectrum — are useful targets for stochastic models of jet meandering.\n\nThe soft spot is load-bearing, and it is one the authors flag as observational. The interface C(x,t) is identified with the local maximum of the zonal velocity u(x,y,t) in y at each x, justified by \"it appears\" and by animations. Every statistic in the paper is computed on those tracked contours. The paper does not explain how multiple candidate maxima are resolved, what happens when a jet weakens or vanishes, or how merging and splitting are handled. If the tracker locks onto transient eddy maxima or jumps between jets, the roughness, shock-like events, and fat-tailed speeds are properties of the rule, not of the physical vorticity boundaries. That is a serious concern, not cosmetic.\n\nOther issues are minor. Exponents come without error bars or fit ranges; the abstract says \"moments\" while only the second-order structure function is computed; no code or data are released, which matters because the observable is defined by an unpublished algorithm. All fixable.\n\nThe circularity concern from the stress-test does not really land: nothing is fitted to reproduce exponents, and the method borrowing from earlier work is legitimate. The paper is honest about its evidence base.\n\nWho benefits: researchers modeling zonal jet interfaces, particle transport in planetary atmospheres, and stochastic contour dynamics. A serious referee should ask for quantitative validation of the tracking rule against an independent interface definition, plus error bars and convergence checks. I would send it to review, not desk reject; the contribution is clear and useful, and the validation gap is fixable.\n\nRecommendation: peer review, major revision.","headline":"Systematic, well-written statistics of jet-interface fluctuations in beta-plane turbulence, but the load-bearing tracking rule is only visually validated and the exponents lack error bars; worth a serious referee.","tokens_in":12542,"tokens_out":3256,"would_cite":false,"duration_ms":33171,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["47.27.-i"],"model":"deepseek-v4-flash","headline":"This paper establishes that the vorticity interfaces bounding zonal jets in barotropic beta-plane turbulence are rough, multifractal fluctuating contours with sub-Gaussian heights, fat-tailed speeds, and specific power-law spectra.","keywords":["zonal jets","beta-plane turbulence","vorticity interfaces","interface tracking","multifractal","height field statistics","structure functions","zonostrophy parameter"],"falsifier":"Extract interfaces from the same DNS with an independent definition, for example level sets of the vorticity field or ridges of $|\\nabla \\omega|$, and compare the contours and the resulting height statistics with the velocity-maximum contours. If the two definitions frequently diverge, especially during jet merging or strong eddy activity, or if the reported exponents $\\zeta_2$, $\\xi_2$, and the multifractal widths change materially under the alternative definition, the central claim is not robust.","tokens_in":11625,"feed_emoji":"🌀","tokens_out":6192,"duration_ms":62697,"temperature":0.7,"pith_summary":"The paper develops a way to track the sharp interfaces between zonal jets in forced-dissipative barotropic $\\beta$-plane turbulence by following the local maxima of the eastward zonal velocity, and then works out the statistical laws of those interfaces. It claims that the interfacial height fluctuations are sub-Gaussian, while the fluctuation speeds are fat-tailed, becoming more Gaussian as the zonostrophy parameter increases. It further claims that the time series of the mean and variance of the interfacial height have power-law spectra with exponents $f^{-2}$, and $f^{-5}$, $f^{-5/3}$, $f^{-7/2}$ respectively, and that the temporal height field is non-differentiable, with a short-time second-order structure exponent $\\zeta_2 < 2$ that approaches 2 with increasing $R_\\beta$. A multifractal analysis shows the interfaces have a broad singularity spectrum that narrows as $R_\\beta$ grows. If these results hold, they give quantitative fluctuation laws for jet interfaces in planetary and oceanic flows.","feed_headline":"Turbulent jet edges obey rough, multifractal fluctuation laws","feed_subtitle":"Simulations track jet interfaces to reveal sub-Gaussian heights, fat-tailed speeds, and power-law spectra","key_machinery":"The central object is the tracked contour $C(x,t)$, defined operationally as the set of local maxima of the eastward zonal velocity when $u$ is plotted against $y$ at each fixed $x$, and then the height field $h(x,t)=C(x,t)-\\bar{C}(x)$ measures its deviation from the time-averaged profile. This scalar height field carries the entire analysis: its probability density function, the speed $v_C = \\lim_{\\delta t \\to 0} \\delta h/\\delta t$, sectorial mean $M_{h,s}$ and variance $V_{h,s}$ spectra, the second-order temporal structure function $S_2(\\tau)$, and the partition-function multifractal spectrum built from $H=|h|+0.001$. The Rossby time $\\tau_w = 1/(L_\\epsilon \\beta)$ normalizes time increments, and the zonostrophy parameter $R_\\beta = L_{Rh}/L_\\epsilon$ controls the strength and smoothness of the jets.","core_discovery":"The interface between successive vorticity shear zones, tracked as the locus of local maxima of the zonal velocity $u(x,y,t)$ at fixed $x$, is not a smooth material line but a rough, irregular contour whose fluctuations carry turbulence-like statistics. The height field $h(x,t)=C(x,t)-\\bar{C}(x)$ has a Gaussian core with sub-Gaussian tails; the speed field $v_C = dh/dt$ has heavy tails; the height variance spectrum shows the $\\beta$-plane turbulence power laws $f^{-5}$, $f^{-5/3}$, and $f^{-7/2}$, while the sectorial mean height spectrum decays as $f^{-2}$; the second-order temporal structure function violates differentiability at short times ($\\zeta_2 \\neq 2$) with $\\zeta_2 \\to 2$ for larger $R_\\beta$; and the interfaces are multifractal, with the singularity range narrowing as $R_\\beta$ increases.","pith_inferences":["Beyond the paper: the same tracking rule should be applied to westward jets or to both velocity extrema to test whether the reported exponents are specific to eastward jet flanks or universal to all jet interfaces.","Beyond the paper: a quantitative validation of the velocity-maximum tracking against an independent vorticity-based interface definition, such as level sets of $\\omega$ or ridges of $|\\nabla \\omega|$, would settle whether the statistical laws belong to the physical interfaces or to the tracking rule.","Beyond the paper: because interfaces induce lateral shear, the reported statistics suggest a testable extension in Lagrangian transport: microplastic-like tracers in the same DNS should show enhanced stretching and preferential sampling near interfaces with larger $R_\\beta$."],"forward_implications":["If the claim is correct, the fluctuating edges of zonal jets have a complete statistical description: sub-Gaussian heights, fat-tailed speeds, and specific spectral exponents, which can serve as quantitative benchmarks for reduced models of planetary jet dynamics.","The variance spectrum with $f^{-5}$, $f^{-5/3}$, and $f^{-7/2}$ regimes indicates that interfacial height variance tracks the underlying beta-plane turbulence energy spectrum, so interfacial variance measurements could act as a proxy for the turbulence at these scales.","The short-time non-differentiability ($\\zeta_2 < 2$) rules out smooth stochastic contour models and would need to be reproduced by any stochastic differential equation intended to model jet-interface evolution.","The systematic trend of $\\zeta_2 \\to 2$ and the narrowing singularity spectrum as $R_\\beta$ increases show that stronger, more persistent jets host smoother, less intermittent interfaces, a distinction relevant to comparing terrestrial and gas-giant zonal flows."],"supporting_citations":[{"why":"Supplies the zonostrophy parameter $R_\\beta$ and the beta-plane jet regime whose range $6.3$ to $7.5$ the simulations target.","marker":"[10]"},{"why":"Provides the narrow-band random forcing scheme used to drive the statistically steady turbulent state.","marker":"[53]"},{"why":"Provides the exponential time-differencing integration method used for the direct numerical simulations.","marker":"[56]"},{"why":"Establishes the structure of zonal jets in geostrophic turbulence that motivates interpreting velocity maxima as jet interfaces.","marker":"[60]"},{"why":"Gives the reference for fluctuating interface variance spectra and multifractal analysis that the present study contrasts with binary active-fluid interfaces.","marker":"[37]"},{"why":"Supplies the partition-function multifractal procedure applied to the interfacial height field.","marker":"[2]"}],"fun_headline_variants":["Jet interfaces: sub-Gaussian heights, fat-tail speeds, power laws","Rough jet edges: heavy-tailed speeds and multifractal structure","Beta-plane jets: interfaces are rough, multifractal, non-differentiable","Jet interface statistics: sub-Gaussian heights and heavy-tailed speeds","Rough, multifractal jet interfaces: sub-Gaussian heights, fat-tail speeds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the vorticity interfaces are exactly the local maxima of the eastward zonal velocity at each $x$, a correspondence checked by eye through animations rather than by quantitative comparison with jumps in vorticity; if this identification fails when jets wobble, merge, or eddies pass, every reported exponent belongs to the tracking rule, not to the physical interface.","fun_headline_variants_meta":{"raw":{"variants":["Jet interfaces: sub-Gaussian heights, fat-tail speeds, power laws","Rough jet edges: heavy-tailed speeds and multifractal structure","Beta-plane jets: interfaces are rough, multifractal, non-differentiable","Jet interface statistics: sub-Gaussian heights and heavy-tailed speeds","Rough, multifractal jet interfaces: sub-Gaussian heights, fat-tail speeds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000792,"raw_usage":{"total_tokens":3462,"prompt_tokens":893,"completion_tokens":2569,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":509,"completion_tokens_details":{"reasoning_tokens":2471}},"tokens_in":509,"tokens_out":2569,"duration_ms":21787,"temperature":1.0,"reasoning_tokens":2471,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T10:41:36.746700+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Extract interfaces from the same DNS with an independent definition, for example level sets of the vorticity field or ridges of $|\\nabla \\omega|$, and compare the contours and the resulting height statistics with the velocity-maximum contours. If the two definitions frequently diverge, especially during jet merging or strong eddy activity, or if the reported exponents $\\zeta_2$, $\\xi_2$, and the multifractal widths change materially under the alternative definition, the central claim is not robust.","supporting_citations":[{"cited_title":"Galperin, S","cited_arxiv_id":null,"evidence_quote":"Supplies the zonostrophy parameter $R_\\beta$ and the beta-plane jet regime whose range $6.3$ to $7.5$ the simulations target."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the narrow-band random forcing scheme used to drive the statistically steady turbulent state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the exponential time-differencing integration method used for the direct numerical simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the structure of zonal jets in geostrophic turbulence that motivates interpreting velocity maxima as jet interfaces."},{"cited_title":"Turbulence-Induced Fluctuating Interfaces in Heterogeneously-Active Suspensions","cited_arxiv_id":"2502.16443","evidence_quote":"Gives the reference for fluctuating interface variance spectra and multifractal analysis that the present study contrasts with binary active-fluid interfaces."},{"cited_title":"Mukherjee, S","cited_arxiv_id":null,"evidence_quote":"Supplies the partition-function multifractal procedure applied to the interfacial height field."}],"review_version":1}