{"id":"58fb2984-68d7-4a75-a728-e105f0b375e3","arxiv_id":"2505.02068","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A damped Ekman model with a Rayleigh drag term reproduces observed flattened Ekman spirals and changes estimates of wind-driven transport and upwelling.","lead":"This paper adds one linear drag term to the classic ocean Ekman spiral equations and shows it can explain why real ocean currents decay faster with depth than they rotate. If the theory holds, standard estimates of wind-driven ocean transport and upwelling are too high by roughly 10 to 20 percent in key regions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ocean validation is circular: R is inverted from the same spiral shape the theory explains, so Fig. 2 cannot discriminate damping from other flattening mechanisms; an independent damping estimate is needed.","rationale":"The analytical model is internally consistent: solving Eq. 1 (A U_zz - (R+if)U = 0) gives d_amp = mu d_Ek/(sqrt(2) cos theta), d_rot = mu d_Ek/(sqrt(2) sin theta), and the product identity d_amp d_rot = d_Ek^2 follows from sin(2 theta) = mu^2; the transport and pumping formulae also follow under uniform R. I therefore do not object to the mathematical development. The load-bearing concern is attribution to a real oceanic damping. In the observational comparison, the same data provide d_amp, d_rot, R, Az, and the wind stress, so Fig. 2 is not a predictive test; the only genuinely independent checks are the MITgcm experiment (where the added term is explicitly a Laplacian damping, so support is limited) and the five-point xi-Da correlation, which is qualitative. The paper explicitly concedes the difficulty of determining R or Da. Since this concern is essentially the reader's weakest assumption, I agree with the conditional verdict: the theory is plausible and internally sound, but the observational case for turbulent dissipation as the cause of flattening is not yet established. A targeted momentum-balance diagnosis in a high-resolution simulation would settle whether the process is actually a linear damping.","tokens_in":16210,"tokens_out":12744,"duration_ms":169642,"concrete_test":"Run a targeted CESM (or LES) experiment at the five sites with high-frequency output and diagnose, from the full momentum balance, the residual force F = du/dt + (u.grad)u + f k x u - A_z d2u/dz2 - (horizontal mixing and forcing terms). Then form R_indep(z) = -(F.u)/|u|^2. If R_indep is positive and roughly depth-uniform in the Ekman layer, and its vertical or area average is within a factor of 2 of the R values in Table S3, then the linear-damping representation is supported; if not, the inferred R is an artifact of the inversion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is not the algebra of Eq. 2 (which is internally consistent) but the identification of R with a real oceanic process. In the observational validation, the damping number Da is not measured: d_amp and d_rot are first obtained from the slope of log speed and of phase angle versus depth (SI Eqs. S17-S18), then Az and R are computed from d_amp*d_rot = d_Ek^2 and R = A_z(1/d_amp^2 - 1/d_rot^2), and the wind stress magnitude and direction are fixed by substituting one observed velocity into Eq. 2. Thus R, Az, and the surface forcing are all inverted from the same spiral shape the theory is meant to explain. Any observed spiral with a decaying amplitude and a rotating direction can be accommodated by choosing R > 0; the improvement over the classic theory is therefore a measure of added parameters, not evidence that damping is the physical cause. The MITgcm experiment is the cleanest support, but there the added term in the model is literally a horizontal-Laplacian Rayleigh friction, so it confirms the math rather than the ocean mechanism. The CESM xi-Da correlation (Fig. 5) has only five points, and xi is not converted to an R value. The paper itself states that 'it is difficult to determine the value of R or Da' (Discussion), which is the soft spot.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an extension of the classic steady Ekman model by adding a linear Rayleigh damping term R to the momentum equation (Eq. 1), yielding an analytical solution (Eq. 2) with two depth scales, d_amp < d_Ek < d_rot and d_amp d_rot = d_Ek^2. From this solution the paper derives a reduced Ekman transport (Eq. 4) and a modified Ekman pumping that depends on both wind-stress curl and divergence (Eq. 5). The authors then invert d_amp and d_rot from five observed Ekman spirals, compare the new and classic solutions, test the model against an idealized MITgcm experiment and a CESM1.3 simulation, and estimate the effect of damping on zonally integrated transports and pumping, including a ~0.4 Sv (12%) reduction at 26.5°N for Da=0.2.","tokens_in":16490,"tokens_out":5098,"duration_ms":46082,"significance":"If correct, the analytic result is an elegant and useful extension of Ekman theory: it explains the universal flattening with a single additional parameter, predicts testable reductions in transport and a new divergence contribution to pumping, and it unifies several previously proposed mechanisms (e.g., complex viscosity) under one damping term. The derivation from Eq. 1 to Eqs. 2–5 is internally consistent, and the idealized MITgcm test with an independently estimated horizontal-friction damping is a genuine, falsifiable check of the mathematics. The main weaknesses are that the observational fits are not independent (R, Az, and the effective wind stress are all derived from the same spiral data via SI Eqs. S17–S18), the CESM support rests on only five points and a qualitative proxy ξ that is not converted to R, and the paper itself concedes that 'it is difficult to determine the value of R or Da.' As it stands, the manuscript establishes an internally consistent theory and one idealized numerical confirmation, but not the stronger claim that turbulent dissipation in the real ocean is the main cause of Ekman spiral flattening.","major_comments":[{"comment":"The observed-spiral validation is circular. In the SI, d_amp and d_rot are obtained by linear fits to ln|u| and arg(u) versus depth (Eqs. S17–S18); Az is then determined from d_amp d_rot = d_Ek^2 and R from R = Az(1/d_amp^2 − 1/d_rot^2), and the wind-stress magnitude and direction are fixed by substituting one observed velocity into Eq. 2. The theoretical curves in Fig. 2 are therefore constructed from the same data points they are compared with. Any observed spiral with monotonically decaying speed and rotating direction will be fitted by some R > 0, so the improvement over the classic theory measures the number of added parameters, not the physical correctness of damping. The statement that the new theory 'fits the observed data points well, much better than' the classic theory (main text, Comparison section) is therefore not evidence that damping is the cause of flattening. I ask the authors to either (i) obtain R from independent measurements (e.g., microstructure dissipation, momentum-budget residuals, or the model's own friction coefficient) and then compare predicted versus observed spirals without re-fitting, or (ii) explicitly re-frame the observational comparison as an inversion/parameterization exercise rather than a hypothesis test.","section":"Comparison with the observed Ekman spiral (main text) and SI Eqs. S17–S18"},{"comment":"The independent-support evidence from CESM1.3 is weak. Fig. 5B plots Da against ξ at five locations with no error bars and no significance test, and the SI itself states that converting ξ to R is not possible because the proportionality coefficient in Eq. S20 depends on the dissipation rate and the monthly output cannot yield R. With n = 5, a monotone relationship is suggestive at best; it does not establish that the inverted Da values (Table S3) are physically tied to the model's actual damping. Please provide a quantitative uncertainty analysis, additional locations or time periods, or a direct estimate of R from the model's momentum budget, and report a correlation coefficient or equivalent.","section":"Diagnosing damping in the realistic global experiment; Fig. 5, SI Eqs. S19–S20"},{"comment":"The paper's own caveat that 'it is difficult to determine the value of R or Da' (Discussion) is load-bearing. Since R is the only new physical parameter and the observational determination of it is circular, the central oceanographic claim — that damping is the main cause of the universal flattening — is not yet established. The idealized MITgcm experiment uses a horizontal-Laplacian friction term that is literally the R-term in the model, so it confirms the mathematics of Eq. 2 but not that ocean turbulence acts as a depth-uniform linear Rayleigh damping. The manuscript should be revised so that the abstract and conclusions either include a clear statement that the ocean evidence is consistent with, but does not uniquely identify, damping as the flattening mechanism, or the authors should supply the independent R estimate needed to close the loop.","section":"Discussion"}],"minor_comments":[{"comment":"In the Discussion section, 'the influence of damping on Ekman transport and damping' should read 'Ekman transport and pumping'; the current phrasing is confusing.","section":"Discussion"},{"comment":"The caption of Fig. 2 states that the coordinate axes differ among panels (east/north for A and D; downwind/crosswind for B, C, and E). This is not obvious from the figure itself; please add axis labels or an inset compass to each panel so the reader can interpret the spiral geometry without reading the caption closely.","section":"Fig. 2 caption"},{"comment":"Table S2 lists observed wind stress magnitudes τ_obs that differ considerably from the fitted τ_mag (e.g., 0.09 vs 0.127 N m⁻² at the coast of California; 0.08 vs 0.161 N m⁻² at Drake Passage). Since the text says the wind stress is obtained from Eq. 2 rather than from observations, please clarify whether these differences are expected given Stokes drift contamination or other uncertainties, and how the chosen data point for substitution was selected.","section":"Table S2"},{"comment":"Equation (4) in the main text appears to have typographical issues in the subscripts and superscripts; please ensure the formula is typeset correctly and all symbols are defined explicitly.","section":"Eq. 4"},{"comment":"Fig. 5 would benefit from error bars or colored markers distinguishing the five locations; with only five points the visual impression of proportionality is difficult to assess.","section":"Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for the journal and the analytic derivation is worth publishing if the causal claims are appropriately tempered. The main gatekeeping question is whether the authors can provide an independent estimate of R from real-ocean data or will instead revise the claims to acknowledge that the observational comparison is an inversion rather than a test. I lean toward major revision rather than rejection because the derivation and the idealized numerical test are sound, but the current observational evidence does not support the strong 'mainly due to damping' conclusion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read on this paper. The genuinely useful thing is the clean separation of the amplitude and rotation depth scales once you add a Rayleigh damping term to the Ekman equation. The product relation d_amp*d_rot = d_Ek^2 and the transport reduction by mu^2 are neat, and the pumping formula with the divergence term is a real consequence that people doing wind-driven diagnostics should know. The authors are honest that the damped equation itself isn't new; they cite McWilliams and Huckle and show equivalence to complex-viscosity treatments. The novelty is in the interpretation and in working out the consequences.\n\nWhat the paper does well is the analytic derivation and the MITgcm check. In that experiment R is estimated from Ah, not inverted from the spiral shape, and the damped solution fits the modeled spiral. That gives independent support for the mechanism within the model. The CESM correlation (Fig. 5) is a good idea but has only five points and xi is not converted to R, so it is weak.\n\nThe soft spot is the observational validation. The paper fits d_amp and d_rot from the same current profiles, then computes Az and R from those two scales, and even the wind stress is set by plugging one observed velocity into the solution. So Fig. 2 is closer to interpolation than to a predictive test. Any spiral with decaying amplitude and rotating direction can be accommodated with R > 0. That doesn't make the theory wrong, but it means the observations cannot discriminate damping from other flattening mechanisms. The paper admits it is difficult to determine R or Da, and some data points are excluded after the fact, which adds to the concern. The transport and pumping corrections for Da = 0.2 or 0.4 are scenario estimates, not predictions; the authors mostly present them that way.\n\nWho is this for? Anyone who uses Ekman theory to estimate meridional transports, especially at 26.5N with RAPID, or upwelling from wind stress. The formulas are simple enough to be used as an uncertainty check. It deserves serious peer review. A good referee should push for an independent estimate of R, maybe from microstructure dissipation or from a model where turbulent closure produces the damping explicitly, and for a cleaner statement of what is fit versus what is predicted. I would bring it to reading group, maybe for the methods section.","headline":"A clean extension of Ekman theory whose observational validation is partly circular, but the framework is useful and deserves serious peer review.","tokens_in":17071,"tokens_out":2828,"would_cite":true,"duration_ms":34936,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["86A05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Adding a linear damping term to the Ekman balance reproduces the observed flattening of ocean spiral currents and implies weaker wind-driven transport and pumping.","keywords":["Ekman spiral","Ekman transport","Ekman pumping","linear damping","damping number","turbulent dissipation","wind stress divergence","ocean boundary layer"],"falsifier":"At a wind front where the curl of the wind stress vanishes but the divergence is large, the theory predicts Ekman pumping driven only by the divergence term; repeated absence of the corresponding vertical velocity in observations would falsify the damping-enabled theory.","tokens_in":15983,"feed_emoji":"🌊","tokens_out":8797,"duration_ms":81650,"temperature":0.7,"pith_summary":"This paper sets out to explain why observed ocean Ekman spirals are always flatter than the classic theory predicts: the current speed decays with depth faster than its direction rotates. The proposed explanation is a single linear damping term added to the steady Ekman momentum balance, standing for turbulent dissipation, lateral friction, and related unresolved processes. With that one term the analytical solution changes in a testable way: the amplitude and rotation depth scales decouple, the surface deflection angle falls below 45 degrees, Ekman transport is reduced, and Ekman pumping gains a wind-stress divergence contribution. The paper compares the modified spiral with observations at five sites and with numerical simulations, and estimates that along 26.5 degrees north the zonally integrated Ekman transport shrinks by roughly 0.4 Sv (about 12 percent) under a modest damping. If the claim holds, classic Ekman-based estimates of wind-driven circulation need correcting.","feed_headline":"One damping term flattens ocean spirals, cuts wind-driven transport","feed_subtitle":"A single drag term makes modeled ocean currents match observed spirals and trims wind-driven transport by 12 percent.","key_machinery":"The central object is a linear Rayleigh damping term $R$ added to the steady Ekman momentum equation, expressed through the damping number $Da=R/|f|$. In the resulting solution the single classic Ekman depth $d_{Ek}=\\sqrt{2A_z/|f|}$ splits into two depth scales: $d_{amp}$ for the exponential decay of current speed and $d_{rot}$ for the rotation of the velocity vector, with $d_{amp} \\cdot d_{rot} = d_{Ek}^{2}$. This split is the mechanism that produces the flattening, and the same solution carries the transport factor $\\mu^{2}=(1+Da^{2})^{-1/2}$ and the extra wind-stress divergence term in Ekman pumping.","core_discovery":"The central claim is that adding a linear damping term $R$ to the steady Ekman equation yields a solution in which the amplitude decay scale $d_{amp}$ and the rotation scale $d_{rot}$ no longer coincide: $d_{amp} < d_{Ek} < d_{rot}$ and $d_{amp} \\cdot d_{rot} = d_{Ek}^{2}$. The paper calls this decoupling the natural mechanism behind the flattening of observed Ekman spirals. The same solution gives a surface deflection angle $\\theta = \\tfrac{1}{2}\\arctan(Da^{-1})$ below 45 degrees, a surface-speed reduction factor $\\mu = (1+Da^{2})^{-1/4}$, an Ekman transport reduced by $\\mu^{2}$ and no longer perpendicular to the wind, and an Ekman pumping that combines the wind-stress curl with a wind-stress divergence term. The paper supports this by fitting the new spiral to observations at five sites spanning different oceanic regimes, by reproducing flattened spirals in idealized and realistic simulations, and by showing that the inferred damping number tracks an independent dissipation-rate estimate.","pith_inferences":["Inference: because the damping number $Da=R/|f|$ grows toward the poles at fixed $R$, the same physical dissipation would produce stronger spiral flattening and larger transport corrections at high latitudes; this latitude dependence is testable with existing current-meter and drifter data.","Inference: the divergence term in the pumping rule means Ekman upwelling can occur under purely convergent winds even where the curl vanishes, so upwelling diagnostics based only on wind-stress curl may misattribute the forcing in frontal regions.","Inference: a direct turbulence measurement of $R$ at one of the five sites would turn the inversion from a fitting exercise into an independent test, since the paper's $R$ values are currently inferred from the same spirals the theory explains."],"forward_implications":["At any nonzero damping number, the spiral is flatter than classic: amplitude decays over $d_{amp} < d_{Ek}$ while direction rotates over $d_{rot} > d_{Ek}$, with $d_{amp} \\cdot d_{rot} = d_{Ek}^{2}$.","The surface current is reduced by the factor $\\mu=(1+Da^{2})^{-1/4}$ and turns through a deflection angle below 45 degrees.","Ekman transport is multiplied by $\\mu^{2}$ and rotates so it is no longer perpendicular to the wind stress.","Ekman pumping is generally weakened by $\\mu^{2}$ but gains a wind-stress divergence term, so it may locally strengthen even where the classic curl-only rule would predict no pumping.","Zonally integrated Ekman transport along 26.5 degrees north is about 0.4 Sv (12 percent) smaller when $Da=0.2$ than the classic theory predicts."],"supporting_citations":[{"why":"Defines the classic Ekman spiral and the baseline solution that the new theory modifies.","marker":"[1]"},{"why":"Supplies the observational evidence that observed current speed decays faster than rotation, the flattening phenomenon to be explained.","marker":"[3]"},{"why":"Documents unequal observed depth scales and a complex-viscosity interpretation that the paper recasts as linear damping.","marker":"[5]"},{"why":"Provides a stratified-Ekman alternative that the new theory is compared against at the observation sites.","marker":"[9]"},{"why":"Provides a time-dependent, Stokes-inclusive Ekman solution used as a comparison baseline.","marker":"[11]"},{"why":"Gives the Drake Passage comparison where the new theory is claimed to fit better than a geostrophic-shear-corrected classic model.","marker":"[24]"},{"why":"Supplies the second-order structure function method used to estimate dissipation independently from surface velocities.","marker":"[25]"}],"fun_headline_variants":["Damping flattens Ekman spiral, cuts wind-driven transport","New Ekman theory: damping explains spiral flatness","One drag term decouples Ekman scales, weakens flow","Damping trims Ekman transport by 12 percent","Ekman spiral flattening traced to turbulent damping"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that a single 'braking' force proportional to current speed and uniform with depth can stand in for all turbulence and friction; if the real braking is nonlinear, varies with depth, or differs from place to place, the predicted flattening, the $\\mu^{2}$ transport reduction, and the divergence term in pumping all change.","fun_headline_variants_meta":{"raw":{"variants":["Damping flattens Ekman spiral, cuts wind-driven transport","New Ekman theory: damping explains spiral flatness","One drag term decouples Ekman scales, weakens flow","Damping trims Ekman transport by 12 percent","Ekman spiral flattening traced to turbulent damping"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000213,"raw_usage":{"total_tokens":1399,"prompt_tokens":898,"completion_tokens":501,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":514,"completion_tokens_details":{"reasoning_tokens":419}},"tokens_in":514,"tokens_out":501,"duration_ms":5342,"temperature":1.0,"reasoning_tokens":419,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T01:03:49.778209+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At a wind front where the curl of the wind stress vanishes but the divergence is large, the theory predicts Ekman pumping driven only by the divergence term; repeated absence of the corresponding vertical velocity in observations would falsify the damping-enabled theory.","supporting_citations":[{"cited_title":"The wavy Ekman layer: Langmuir circulations, breaking waves, and Reynolds stress","cited_arxiv_id":null,"evidence_quote":"Defines the classic Ekman spiral and the baseline solution that the new theory modifies."},{"cited_title":"Direct observations of the Ekman balance at 10 N in the Pacific","cited_arxiv_id":null,"evidence_quote":"Supplies the observational evidence that observed current speed decays faster than rotation, the flattening phenomenon to be explained."},{"cited_title":"Observations of Ekman currents in the Southern Ocean","cited_arxiv_id":null,"evidence_quote":"Documents unequal observed depth scales and a complex-viscosity interpretation that the paper recasts as linear damping."},{"cited_title":"Kolmogorov constants for the second-order structure function and the energy spectrum","cited_arxiv_id":null,"evidence_quote":"Provides a stratified-Ekman alternative that the new theory is compared against at the observation sites."}],"review_version":1}