{"id":"29c1df21-ff64-4a88-b99d-3f660577fb7d","arxiv_id":"2607.29261","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A stress-free boundary condition at the base of the Ekman layer forces the modeled current to grow below that depth, so the condition is internally inconsistent unless the equations change there.","lead":"This paper shows that a common way of closing wind-driven ocean current models—setting the stress to zero at the bottom of the Ekman layer—forces the computed current to strengthen again below that depth, contradicting the idea that the wind's influence dies out there. It then argues that a no-slip condition on the actual seafloor gives the standard right-angle transport rule to very high accuracy in deep water.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 1 is correct, but the 'necessarily unphysical' conclusion assumes Eq. (1) holds below -D; if the deep ocean is geostrophic or the model domain ends at -D, the reverse spiral below -D is an artifact of extending the frictional ODE.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing concern: the physical interpretation relies on the same linear Ekman ODE (1) remaining valid below -D down to the seafloor. Theorem 1 is a correct consequence of integration by parts and is not in doubt. The central weakness is that the paper extends a mathematical result about a single f-plane frictional ODE into a universal physical statement about stress-free boundary conditions. In real ocean settings, and in many of the cited truncated-layer models, z=-D is the base of the turbulent boundary layer and the interior is geostrophic/inviscid; under those dynamics the reverse spiral need not appear. This is a scope limitation rather than a mathematical error, so the reader's CONDITIONAL verdict is appropriate. The additional WKB slip about |U_D|<|U(0)| is real but peripheral; it does not change the conclusion that the stress-free theorem is sound within its stated ODE model. No verdict adjustment is needed.","tokens_in":7788,"tokens_out":13912,"duration_ms":162353,"concrete_test":"Build a two-layer model: for -D<z<0 keep Eq. (1) with U'(-D)=0; for -H<z<-D replace (1) with the steady geostrophic balances f v = ρ^{-1}p_x and f u = -ρ^{-1}p_y, with a smooth prescribed pressure p and continuity of U and zero shear at -D. Check whether a bounded interior solution exists with the same far-field/no-slip condition. If yes for generic p, the growing reverse spiral is an artifact of extending the frictional ODE; if no, the paper's conclusion is supported. This can be done analytically in the constant-coefficient limit and numerically for a KPP-type m(z) profile.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equations (9)-(10) prove the minimum/reverse-spiral statement within the single-ODE model (7); that mathematical part is sound. The load-bearing step is the transfer from the ODE to the ocean: after imposing U'(-D)=0, the paper reads r'<0 on (-H,-D) as a physical 'increase beneath the Ekman layer' and calls the stress-free condition unphysical. This is only forced if Eq. (1), the steady f-plane turbulent-viscosity balance with no horizontal pressure gradient, remains valid through the abyssal water column. The paper asserts this by taking m to the molecular value below -D (WKB section, line before Eq. (11)), but the deep ocean is typically in geostrophic balance; the momentum balance there includes pressure-gradient terms, not just (mU')' = ifρU. In a two-layer formulation where z=-D is the base of the turbulent layer and the interior is geostrophic, a stress-free condition at -D does not force a growing reverse spiral: the interior flow is set by the large-scale pressure gradient and can match U(-D) with zero stress. Thus 'necessarily unphysical' is an overgeneralization. The paper should either restrict the conclusion to models that solve (1) below -D or supply a physical justification that the deep ocean obeys (1). A separate, smaller issue is the WKB section's appeal to Theorem 1 to assert |U_D|<|U(0)| in the no-slip case; Theorem 1 assumes stress-free at -D, so that step is formally invalid, though the exponential bound in Eq. (12) is so small that this slip does not affect the no-slip orthogonality conclusion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the f-plane Ekman equations with depth-dependent density and eddy viscosity, in complex notation (Eq. (7)). Its main theorem proves that if a stress-free boundary condition is imposed at an intermediate depth z = -D, then the speed |U| decreases with depth down to -D, attains a strict minimum there, and increases again below -D; the direction angle behaves analogously. The authors interpret this 'reverse Ekman spiral' as making the mid-depth stress-free condition unphysical. They then use a WKB approximation to show that, for a deep ocean with no-slip bottom boundary and with eddy viscosity falling to the molecular value below the Ekman depth, the bottom velocity derivative is exponentially small, so the Ekman transport is nearly orthogonal to the wind stress. The paper concludes that the no-slip bottom condition is the physically correct choice and that choosing an a priori Ekman depth with a stress-free condition is incorrect.","tokens_in":8129,"tokens_out":8242,"duration_ms":87228,"significance":"The mathematical core of the paper is sound and elegant. Theorem 1 is a parameter-free identity: multiplying Eq. (7) by the complex conjugate of U and integrating gives a signed relation that forces the qualitative behaviour of r(z) and phi(z). This is a useful and nontrivial caution for modelers who solve the full-column f-plane Ekman equation with a stress-free condition at an intermediate depth. The WKB estimate in the second half is also constructive and gives an explicit, quantitative bound on the bottom stress in the deep-ocean no-slip case. However, the paper's headline conclusion—that the stress-free condition is 'necessarily unphysical'—goes beyond the theorem: it requires the additional, unproved premise that the same linear frictional ODE (1) remains valid below the Ekman layer all the way to the seafloor. If that premise is not supplied or the claim is not restricted, the central message overreaches.","major_comments":[{"comment":"The inference from Theorem 1 to 'necessarily unphysical' is load-bearing and not justified. The theorem is a conditional statement about solutions of the ODE (7) on the full interval (-H,0), and the proof integrates from -D to z in (-H,0). The conclusion r'(z)<0 for z<-D is a mathematical consequence only when Eq. (1) is assumed to hold below -D. In the ocean, the abyssal flow is typically geostrophic, with horizontal pressure gradients that are absent from (1). Moreover, in many practical implementations of the stress-free condition at the base of the Ekman layer, the model domain ends at z=-D and no solution is defined below that depth; Theorem 1 does not apply to that truncated problem. The paper should either restrict its conclusion to full-column f-plane frictional models that solve (1) down to -H, or provide a physical justification that (1) describes the deep ocean. Without this,","section":"Abstract and Section 'Main Result' (Eq. (7))"}],"minor_comments":[{"comment":"The proof says 'multiply the first equation in (7) by U(z)', but Eq. (8) contains |U'(s)|^2 and |U(s)|^2, which is only consistent with multiplication by the complex conjugate \\bar U(z). The subsequent extraction of the real and imaginary parts in (9)-(10) is correct with \\bar U, but the notation should be fixed.","section":"Proof of Theorem 1, Eq. (8)"},{"comment":"The sentence 'a stress-free condition at -D analogous to (4)' should refer to Eq. (5), not Eq. (4), since (4) is the no-slip condition.","section":"Introduction, paragraph after Eq. (4)"},{"comment":"The statement 'due to Theorem 1, we have |U_D|<|U(0)|' is formally invalid: Theorem 1 assumes the stress-free condition U'(-D)=0, whereas the WKB scenario imposes no-slip at -H. The same monotonicity can be obtained directly from the integral identity applied to the no-slip problem, or the remark can be omitted, as the WKB estimate does not depend on it.","section":"WKB section, before Eq. (11)"}],"recommendation":"major_revision","confidential_remarks":"The mathematical identity in Theorem 1 is correct and worth publishing, and the WKB estimate is a useful quantitative add-on. The problem is the interpretive leap: the paper's title and abstract claim that the mid-depth stress-free condition is incorrect for Ekman flows, whereas the theorem only addresses solutions of the single-ODE frictional model extended below the Ekman layer. I would support acceptance after a revision that either narrows the claim to that setting or supplies a physical justification for the deep-ocean validity of Eq. (1). If the authors are unwilling to qualify the conclusion, I would not recommend publication in its current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's real contribution is a simple but useful theorem: in the f-plane Ekman model with positive depth-dependent eddy viscosity and density, imposing U'(-D)=0 forces |U| to decrease down to -D, hit a minimum, and then increase again below it. The proof is a one-line integration by parts after multiplying by the conjugate, and it's correct and self-contained. As far as the citations go, the specific statement is new, though it's an immediate corollary of the technique in the authors' companion paper [6]. That's fine.\n\nThe paper also contains a WKB estimate showing that if you instead impose no-slip at the bottom and the ocean is deep, U'(bottom) is exponentially small, so the Ekman transport is nearly orthogonal to the wind. The estimate is standard but competently carried out.\n\nWhere I part company with the authors is the headline conclusion that the stress-free condition is \"necessarily unphysical.\" The theorem says: within the single ODE (1), if U'(-D)=0, then |U| grows below -D. That is a statement about that ODE on an interval that extends below -D. If you are modeling the whole water column with the same turbulent-viscosity balance, then yes, the flow does not vanish below -D. But if the deep ocean is in geostrophic balance, or if the model domain ends at -D and the interior is handled differently, then the stress-free condition does not force any reverse spiral. The paper never acknowledges that this transfer requires an assumption about the deep-ocean momentum balance; it simply asserts that below -D the eddy viscosity drops to its molecular value and (11) applies. That is not enough. So the physically strong claim is not supported; the mathematically precise claim is sound.\n\nThere's also a small slip in the WKB section: to bound |U_D|, the authors cite Theorem 1, but that theorem assumes stress-free at -D, which does not hold in the no-slip case. The exponential estimate is so strong that the slip doesn't change the conclusion, but it should be fixed.\n\nOverall, the paper is a good candidate for publication after revision. The core theorem is a genuine contribution. The physical narrative needs to be more careful: either restrict the \"unphysical\" claim to models that solve (1) below -D, or justify that balance. I'd send it to a fluid-dynamics refereed journal, expecting the authors to tone down the interpretation and fix the WKB citation.","headline":"A clean, correct theorem about the Ekman ODE, but the 'unphysical' label depends on an unstated assumption that the same ODE holds below the chosen depth.","tokens_in":8655,"tokens_out":5314,"would_cite":true,"duration_ms":47678,"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":"This paper proves that imposing a stress-free boundary condition at the base of the Ekman layer forces the modeled current to intensify beneath that depth, a 'reverse Ekman spiral' that contradicts the assumption that wind-driven flow is ne","keywords":["Ekman layer","stress-free boundary condition","Ekman transport","wind-driven currents","eddy viscosity","WKB approximation","reverse Ekman spiral","ocean modelling"],"falsifier":"An observed current profile in the open ocean that exhibits a shear-free point at mid-depth while the speed continues to decay beneath it, and where the local dynamics are independently verified to be governed by the linear Ekman balance (1), would contradict the paper's central claim.","tokens_in":7598,"feed_emoji":"🌊","tokens_out":8511,"duration_ms":75122,"temperature":0.7,"pith_summary":"The paper targets a boundary condition commonly used in ocean models: setting the vertical shear of the horizontal velocity to zero at an intermediate depth -D, taken as the base of the wind-driven Ekman layer. The authors prove that any solution of the f-plane Ekman equations with this condition must have its speed minimum exactly at -D, with the current strengthening again below -D (a 'reverse Ekman spiral') and the wind-current angle decreasing again. This contradicts the very premise of the condition, namely that wind-driven flow is negligible below -D. They then show, via a WKB estimate for slowly varying deep density, that the physically motivated no-slip condition at the true seafloor gives Ekman transport nearly perpendicular to the wind when the ocean is deep, with an exponentially small error.","feed_headline":"Stress-free Ekman boundary forces reverse current spiral","feed_subtitle":"A common boundary condition forces modeled Ekman currents to strengthen at depth, not vanish.","key_machinery":"The central mechanism is an integral identity obtained by multiplying the complex Ekman equation by U and integrating from -D to z: this yields m(z)|U|'(z)|U|(z) = ∫ m|U'|² and m(z)φ'(z)|U|² = f∫ρ|U|², both evaluated from -D. Because the integrands are positive (m,ρ>0), these identities force |U| to have a minimum and φ (relative to the wind) to have a maximum at z=-D. The second piece of machinery is a WKB ansatz, U(z)≈κ(z)^{-1/2} exp(±∫κ), with κ = sqrt(ifρ/m), valid when |ρ'/ρ| ≪ |κ|; this is used to estimate the bottom shear U'(-H) under a no-slip condition and show it is exponentially small for deep oceans.","core_discovery":"On its own terms, the paper establishes Theorem 1: for the Ekman system (mU')' = ifρU with a stress-free condition U'(-D)=0 at an intermediate depth -D (with m>0, ρ>0), the modulus |U| is strictly decreasing from the surface to -D and strictly increasing from -D toward the seafloor; in the Northern Hemisphere the angle between the wind stress and the current varies in the opposite sense, with a maximum at -D. Consequently the flow below -D is not negligible but grows, which is the 'reverse Ekman spiral' the authors deem unphysical. The argument is a two-line integration by parts: multiplying the ODE by U and integrating yields sign-definite integral expressions for r'(z) and φ'(z). The paper","pith_inferences":["A natural extension is to replace the stress-free condition with a matching condition to a geostrophic interior; the theorem suggests any derivative-zero condition at finite depth will misrepresent the deep flow unless the interior dynamics differ from the Ekman balance.","The result implies that coarse-resolution models that parameterize the Ekman layer with a fixed depth and shear-free base will systematically misplace the level of no motion, with consequences for the vertical structure of the simulated currents.","The WKB estimate suggests a testable quantitative prediction: departures of the Ekman transport angle from 90° should be negligible in the deep open ocean but become significant where H-D is small (shelves and coastal regions), consistent with the coastal circulation literature.","One could test the theorem directly in a laboratory or numerical setting: solve (1)–(2)–(6) and measure |U| below -D; the predicted increase is independent of the viscosity profile, so any observed decay below -D would indicate that the Ekman balance itself breaks down there."],"forward_implications":["The stress-free condition at an intermediate depth is internally inconsistent with the assumption that the flow is negligible below that depth; any model that imposes it will generate spurious sub-Ekman currents.","Ocean models should not impose an artificial Ekman depth as a boundary; instead the Ekman depth should emerge from the vertical structure of the eddy viscosity (e.g., a KPP-style profile).","With no-slip at the true seafloor, the Ekman transport is nearly orthogonal to the wind in deep water, matching field observations, with an error that decays exponentially in (H-D)/δ.","The theorem is robust to the details of m(z) and ρ(z): it requires only positivity, so the inconsistency is structural rather than an artifact of a particular profile.","The orthogonality of Ekman transport and wind is recovered under no-slip even with depth-dependent density, contrary to what one might expect from the variable-density transport formula."],"fun_headline_variants":["Mid-depth stress-free Ekman boundary leads to reverse spiral","Stress-free Ekman boundary at depth yields reverse current growth","Mid-depth stress-free condition makes Ekman current deepen","Ekman flow: mid-depth stress-free boundary gives retrograde spiral","Why stress-free Ekman boundary at depth is wrong: spiral reverses"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The argument that the deepening current below -D is unphysical assumes the same Ekman equation—linear, with positive eddy viscosity and density—holds all the way from -D down to the seafloor; in the real ocean the deep flow is often geostrophic or otherwise outside that balance, in which case the reverse spiral may not be realized and the boundary condition might merely be an imperfect but not hopeless truncation.","fun_headline_variants_meta":{"raw":{"variants":["Mid-depth stress-free Ekman boundary leads to reverse spiral","Stress-free Ekman boundary at depth yields reverse current growth","Mid-depth stress-free condition makes Ekman current deepen","Ekman flow: mid-depth stress-free boundary gives retrograde spiral","Why stress-free Ekman boundary at depth is wrong: spiral reverses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000643,"raw_usage":{"total_tokens":2758,"prompt_tokens":669,"completion_tokens":2089,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":413,"completion_tokens_details":{"reasoning_tokens":2005}},"tokens_in":413,"tokens_out":2089,"duration_ms":12935,"temperature":1.0,"reasoning_tokens":2005,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T10:32:23.249649+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An observed current profile in the open ocean that exhibits a shear-free point at mid-depth while the speed continues to decay beneath it, and where the local dynamics are independently verified to be governed by the linear Ekman balance (1), would contradict the paper's central claim.","supporting_citations":[],"review_version":1}