{"id":"7adfcf86-4a19-4b00-a4e4-14573e65216b","arxiv_id":"2602.06473","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Starting from Navier-Stokes on a rotating sphere, the authors derive an asymptotic Ekman-layer model with latitude-dependent Coriolis force and nonlinear wind stress, prove existence and uniqueness of the leading-order solution, and compute surface deflection angles for five eddy viscosity profiles.","lead":"This paper derives a simplified mathematical model for wind-driven ocean currents using the full spherical geometry of the Earth, rather than the flat-plane approximation used in the classic 1905 Ekman theory. The model recovers the known spiral shape of these currents and predicts how much the surface current is deflected from the wind direction, for several realistic profiles of turbulent viscosity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Flat-surface assumption is internally inconsistent: the non-dimensional kinematic boundary condition (2.14) appears mis-scaled by a factor ε, so the leading-order solution need not satisfy w=0 at a flat surface, undermining the derivation of the model.","rationale":"The reader correctly identified the flat-free-surface assumption as the weakest point. My stress test sharpens this: the issue is not only that waves and Stokes drift are neglected; the flat-surface condition is internally inconsistent with the leading-order continuity equation as derived. The printed equation (2.14) appears to have an extra ε factor, and the paper then uses that to deduce h_0 is an invariant. With the correct scaling, h_0=0 would impose w_0(0)=0, which is not enforced by the constructed solution. This affects the paper's central claim that (3.20) is the asymptotic limit of the Navier–Stokes free-surface problem, although the ODE itself and Theorem 4.1 could survive if the surface condition is reformulated as a rigid-lid problem with an unknown surface pressure. I therefore do not move the verdict away from the reader's CONDITIONAL: the mathematical core is plausibly salvageable, but the derivation and physical interpretation require correction. The spurious factor in (5.20) and the abstract/body mismatches are real but secondary compared to this boundary-condition scaling question.","tokens_in":43871,"tokens_out":31142,"duration_ms":276258,"concrete_test":"Substitute (2.10) into (2.6) and verify whether (2.14) should read w = ∂h/∂t + ... rather than εw = ... . Then, using the constant-viscosity solution (5.5), take a non-uniform wind such as W_w(φ)=A e^{iφ} with A>0 and compute w_0(0) = -(1/cosθ)∫_{z_0}^0 [∂_φ u_e + ∂_θ(v_e cosθ)] ds at a fixed latitude (e.g., θ=45°). If w_0(0)≠0, the flat-surface boundary condition is violated and the derived model is not a closed free-surface limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Under the stated scaling (2.10) with h'=D'h and w'=ε U'w, the dimensional kinematic condition (2.6) should become w = ∂h/∂t + u/((1+εz)cosθ)∂h/∂φ + v/(1+εz)∂h/∂θ, not εw = ... as printed in (2.14). The extra ε causes the expansion in §3 to yield (3.5) as an invariant condition for h_0; the correct leading-order condition would be w_0 = Dh_0/Dt. If one then sets h_0=0 (a flat surface), the correct condition forces w_0(0)=0. However, the paper does not impose this on the leading-order solution. From (3.14) and the definition of W after (3.20), w_0(0) = -(1/cosθ)∫_{z_0}^0 [∂_φ u_e + ∂_θ(v_e cosθ)] ds, which is generically nonzero for a spatially varying wind. Thus the flat-surface solution violates the free-surface boundary condition. Alternatively, if the flat surface is read as a rigid lid, the pressure at the lid should be an unknown P_0, not the atmospheric pressure P_s fixed by (3.12). Either way, the asymptotic reduction is not a consistent free-surface limit. Since the model (3.20) and its deflection-angle predictions rest on this reduction, this is a load-bearing gap, independent of the waves/Stokes-drift limitation acknowledged in the Discussion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives an asymptotic model for wind-driven Ekman currents from the Navier–Stokes equations in rotating spherical coordinates, using the small thin-shell parameter ε and a small Rossby number R. The leading-order problem is reduced to the complex ODE (m(z)W')' = 2i sinθ W with a nonlinear wind-stress boundary condition at z=0 and a vanishing condition at the Ekman depth z0. The authors prove existence and uniqueness of the solution, show it forms a classical Ekman spiral, derive a formula for the surface deflection angle, and provide explicit Bessel-function solutions for five eddy-viscosity profiles. They also formulate the first-order correction in R and give a priori bounds. The central claim is that this reduced model is a consistent spherical-coordinate generalization of the classical Ekman theory and that its deflection-angle predictions agree well with observations.","tokens_in":44282,"tokens_out":13531,"duration_ms":188041,"significance":"If the derivation were fully consistent, the paper would be a valuable contribution: it treats a classical problem without tangent-plane approximations, gives rigorous existence/uniqueness results for a nonlinear boundary-value problem, and supplies explicit analytic solutions for several physically motivated eddy-viscosity profiles. The ODE analysis in Section 4 and the explicit computations in Section 5 appear careful and are likely correct as standalone mathematics. However, the derivation from the free-surface problem contains a scaling error in the kinematic boundary condition, which undermines the claim that (3.20) is the correct leading-order model of a flat free surface. The comparison with observations is also qualitative, since the layer depth z0 and wind magnitude |Ww| are chosen ad hoc rather than fitted systematically.","major_comments":[{"comment":"The non-dimensional kinematic boundary condition is mis-scaled. With w' = ε U' w and h' = D' h, substituting into (2.6) gives w = ∂h/∂t + u/((1+εz)cosθ)∂h/∂φ + v/(1+εz)∂h/∂θ, not εw = .... The extra ε propagates into the expansion in §3, so (3.5) should read w0 = ∂h0/∂t + u0/cosθ ∂h0/∂φ + v0 ∂h0/∂θ on z=h0, not an invariant condition. Since the paper then sets h0=0, the correct condition would impose w0=0 at z=0. The solution obtained from (3.20) and the formula for w0 after (3.20) does not in general satisfy this, because w0(0) is determined by the horizontal divergence of the ageostrophic field. Thus the flat-surface reduction is not a consistent asymptotic limit of the free-surface problem. This is a load-bearing issue for the derivation of the central model (3.20), independent of the wave/Stokes-drift limitation acknowledged in §6. Please correct the scaling or explicitly reformulate","section":"§2.3, Eq. (2.14); §3.1, Eq. (3.5)"},{"comment":"The mathematical analysis of the reduced problem (3.20) is sound, but the physical interpretation of that problem is affected by the free-surface inconsistency. If the flat surface is instead interpreted as a rigid lid, the leading-order pressure should not be fixed to the atmospheric pressure P_s; an unknown lid pressure would enter the geostrophic balance and hence alter Ww. The paper does not provide this alternative interpretation. The central claim that (3.20) is the leading-order free-surface model therefore needs either a corrected derivation or a clearly stated change of modeling assumptions.","section":"§3.1 and §4, in relation to the model (3.20)"}],"minor_comments":[{"comment":"The abstract states that the fluid has 'depth-varying density' and mentions 'three explicit eddy viscosity profiles,' whereas the body assumes constant density above the thermocline and actually treats five profiles. Please align the abstract with the body.","section":"Abstract"},{"comment":"The parameter R = U'/(Ω'R') is called the 'inverse Rossby number,' but by the standard definition it is the Rossby number itself (up to a factor of 2 depending on convention). Please correct the terminology.","section":"§2.3"},{"comment":"There is a typo: 'One more setting the coefficients to zero' should likely read 'Again setting the coefficients to zero.'","section":"§3.2"},{"comment":"The notation for the Bessel-function coefficients is inconsistent: both c2 and ec2 are used without always clarifying the relationship c2 = ec2 c1. Standardizing this notation would improve readability.","section":"§5.2, §5.3"},{"comment":"The statement that the deflection-angle results are 'remarkably consistent with observations' is somewhat strong, because z0 and |Ww| are chosen by hand in the considered ranges and no systematic comparison or uncertainty quantification is given. Please temper the claim or state explicitly that the comparison is qualitative.","section":"§6"}],"recommendation":"major_revision","confidential_remarks":"The main concern is the mis-scaled kinematic boundary condition, which affects the derivation of the central model. If the authors cannot repair the free-surface scaling or convincingly switch to a rigid-lid formulation, the paper's central derivation would be invalid. However, the ODE analysis and explicit solutions have independent mathematical value, so I do not think outright rejection is necessary yet; a revision should be asked to resolve this load-bearing issue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the Puntini–Roberti–Stefanescu paper. The mathematical core is real: they set up a double asymptotic expansion in the thin-shell parameter and Rossby number, derive the leading-order Ekman problem with the nonlinear wind-stress boundary condition, prove existence and uniqueness for arbitrary depth-dependent eddy viscosity, obtain deflection-angle bounds and the vanishing deflection at the equator, and give explicit Bessel-function solutions for five eddy-viscosity profiles. The proofs of Theorems 4.1 and 4.2 look correct; the polynomial argument for uniqueness is clean. This is a legitimate extension of the Constantin–Johnson spherical-coordinate program and deserves to be in the literature.\n\nThere is a problem in the derivation that the authors have missed. The non-dimensional kinematic boundary condition (2.14) is wrong by a factor of ε. From their own scaling, w' = ε U' w and h' = D' h, so both sides of the dimensional condition carry one power of ε. Dividing gives w = ∂h/∂t + u/((1+εz)cosθ)∂h/∂φ + v/(1+εz)∂h/∂θ, not εw = ... as printed. At leading order they then get an invariant condition for h0, and when they set h0 = 0, they do not impose w0 = 0 at the surface. The vertical velocity recovered from continuity is generically nonzero at z = 0 for a spatially varying wind, so the flat-surface solution violates the free-surface condition. If the flat surface is meant as a rigid lid, the pressure at the lid should be an unknown, not the imposed Ps. Either way, the asymptotic reduction is not a consistent free-surface limit. This is load-bearing for the claim that the model is derived from first principles. The ODE (3.20) might still be a good heuristic model, but the derivation needs reworking.\n\nOther issues are minor by comparison. The statement that the results are \"remarkably consistent with observations\" is qualitative — the deflection-angle comparisons are visual, not a quantitative fit, and the spread across profiles is wide, so the claim is stronger than the evidence. Equation (5.20) has a spurious a_k in the coefficient; the explicit solutions don't appear to use it, but it could mislead. The abstract complaint in the reader's version (depth-varying density) does not apply to the manuscript text, which consistently says constant density.\n\nI would send this to peer review, but I would want the authors to fix the boundary-condition scaling and either add a quantitative comparison or tone down the observational claims. The paper is worth the referees' time.","headline":"Nice spherical-coordinate Ekman framework with solid proofs, but the free-surface boundary condition is mis-scaled by ε, leaving the flat-surface model without a consistent asymptotic justification.","tokens_in":44755,"tokens_out":6460,"would_cite":true,"duration_ms":55190,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q30","35Q35","76D05","76U60","86A05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Starting from the Navier–Stokes equations in rotating spherical coordinates, this paper derives an asymptotic model for wind-driven ocean drift currents that retains full spherical geometry, proves the leading-order solution is a unique Ekm","keywords":["wind-driven ocean currents","Ekman spiral","asymptotic expansion","spherical coordinates","eddy viscosity","surface deflection angle","Rossby number","thin-shell parameter"],"falsifier":"At a mid-latitude site with a known, roughly depth-uniform eddy viscosity and a steady wind near 10 m/s, measure the angle between the surface current and the wind. If the observed deflection differs systematically from the value predicted by formula (4.14) with λ computed from the viscosity profile and ϱ the unique positive root of (4.8), beyond measurement uncertainty, the flat-surface leading-order model would be ruled out.","tokens_in":43734,"feed_emoji":"🌊","tokens_out":5510,"duration_ms":54138,"temperature":0.7,"pith_summary":"Starting from the Navier–Stokes equations in rotating spherical coordinates, the paper derives an asymptotic model for wind-driven ocean drift currents that avoids the classical tangent-plane (f-plane) approximation. The derivation uses a double expansion in two small parameters: the ratio of the Ekman depth to Earth's radius, and the Rossby number. The leading-order result is a linear complex-valued ODE for the ageostrophic current with a nonlinear wind-stress boundary condition; the paper proves this problem has a unique solution forming a classical Ekman spiral for any depth-dependent eddy viscosity. From the solution it derives a closed-form surface deflection angle, and with five explicit viscosity profiles it produces angles consistent with observed values, including both greater and less than 45°.","feed_headline":"Wind-drift deflection angles reproduced by spherical model","feed_subtitle":"No flat-plane approximation needed; observed 45° variations follow from depth-dependent eddy viscosity.","key_machinery":"The key object is the complex-valued boundary value problem (3.20), obtained by splitting the horizontal velocity into geostrophic and ageostrophic components after the double asymptotic expansion in the thin-shell parameter and the Rossby number. The nonlinear wind-stress boundary condition is the only nonlinearity in the leading-order problem; it is what ties the solution's amplitude to the wind. Existence and uniqueness are proved by writing W=cY with Y the solution of the linear homogeneous initial value problem, then solving a quartic for the surface-current magnitude. The surface deflection angle then follows from the complex ratio λ=Y'(0)/Y(0), giving formula (4.14) and the bounds (4.","core_discovery":"The central claim is that the leading-order boundary value problem (m(z)W')' = 2i W sinθ, W'(0)=C|Ww−W(0)|(Ww−W(0)), W(z0)=0 has a unique solution for any positive viscosity profile m(z), and that solution is a classical Ekman spiral: its magnitude grows monotonically toward the surface and its phase rotates monotonically with depth, in the sense sign(θ)ϑ'>0. The proof reduces the nonlinear boundary condition to a fourth-order polynomial whose unique positive root fixes the surface current magnitude, and the deflection angle is then determined by the single complex number λ=Y'(0)/Y(0) associated with the linear homogeneous solution. This yields explicit deflection angles for five viscosity p","pith_inferences":["A testable extension: the model's deflection-angle predictions could be compared against a dataset of simultaneous wind and drifter measurements stratified by latitude, wind speed and estimated eddy viscosity profile; if observed angles fall outside the range spanned by the five profiles, the flat-surface and no-Stokes-drift assumptions would be implicated.","The authors' claim that the scaling forces the vertical velocity to be order ε smaller than horizontal, thereby excluding surface waves, suggests that including a non-flat free surface would require a different scaling regime; this could be a way to incorporate Coriolis–Stokes forcing without abandoning the spherical geometry.","Since the leading-order system is essentially steady with advective time scaling, an unsteady extension with a smaller time scale would include near-inertial oscillations; the same complex-variable framework and quartic-based deflection formula may extend to that regime."],"forward_implications":["The model can describe large-scale wind-drift currents without the f-plane approximation, extending Ekman theory to flows whose horizontal extent is comparable to the Earth's radius.","The surface deflection angle formula (4.14) yields predictions for any depth-dependent eddy viscosity without solving the full problem, and the bounds (4.15) approximate the angle directly from λ without solving the quartic.","For viscosity profiles decaying with depth, the model predicts deflection angles exceeding 45°, while profiles increasing with depth yield angles below 45°; the piecewise linear profile (increasing then decreasing) yields angles below 45°, suggesting the increasing part dominates.","The Ekman transport angle deviates from 45° for finite-depth layers, approaching 45° as the layer depth increases, recovering classical results in the infinite-depth limit.","The first-order correction in the Rossby number is governed by a linear problem whose solution is estimated via logarithmic matrix norms, providing a priori bounds in terms of the leading-order solution."],"fun_headline_variants":["Spherical model reproduces wind-drift angles","New asymptotic model for large-scale ocean drift","Uniqueness proved for Ekman spiral solutions","Wind-driven currents curve correctly on a sphere","Spherical coordinates fix wind-drift deflection"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The model assumes a flat free surface and no surface waves, so wave-induced Stokes drift and Coriolis–Stokes forcing are absent; if these wave effects materially change the near-surface momentum balance for typical winds, the predicted current directions will not correspond to the observed currents.","fun_headline_variants_meta":{"raw":{"variants":["Spherical model reproduces wind-drift angles","New asymptotic model for large-scale ocean drift","Uniqueness proved for Ekman spiral solutions","Wind-driven currents curve correctly on a sphere","Spherical coordinates fix wind-drift deflection"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00022,"raw_usage":{"total_tokens":1309,"prompt_tokens":799,"completion_tokens":510,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":543,"completion_tokens_details":{"reasoning_tokens":442}},"tokens_in":543,"tokens_out":510,"duration_ms":5830,"temperature":1.0,"reasoning_tokens":442,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T03:54:38.264761+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At a mid-latitude site with a known, roughly depth-uniform eddy viscosity and a steady wind near 10 m/s, measure the angle between the surface current and the wind. If the observed deflection differs systematically from the value predicted by formula (4.14) with λ computed from the viscosity profile and ϱ the unique positive root of (4.8), beyond measurement uncertainty, the flat-surface leading-order model would be ruled out.","supporting_citations":[],"review_version":1}