Pith. sign in

REVIEW 3 major objections 5 minor 10 references

Ekman Theory with Damping

T0 review · 3 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict A clean extension of Ekman theory whose observational validation is partly circular, but the framework is useful and deserves serious peer review. read the letter →

arxiv 2505.02068 v1 pith:TE2T6GGV submitted 2025-05-04 physics.ao-ph physics.geo-ph

Jiacheng Wu , Yonggang Liu , Rui Xin Huang , Jinhan Xie , Zhaoying Wang , Shaoqing Zhang This is my paper · ORCID
classification physics.ao-phphysics.geo-ph MSC 86A05
keywords Ekmanspiraltransportpumpinglineardampingnumberturbulentdissipationwindstressdivergenceoceanboundarylayer
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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.

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 (3)
  1. [Comparison with the observed Ekman spiral (main text) and SI Eqs. S17–S18] 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.
  2. [Diagnosing damping in the realistic global experiment; Fig. 5, SI Eqs. S19–S20] 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.
  3. [Discussion] 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.
minor comments (5)
  1. [Discussion] In the Discussion section, 'the influence of damping on Ekman transport and damping' should read 'Ekman transport and pumping'; the current phrasing is confusing.
  2. [Fig. 2 caption] 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.
  3. [Table S2] 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.
  4. [Eq. 4] 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.
  5. [Fig. 5] 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.

Circularity Check

2 steps flagged · score 6.0 of 10

Observational validation is circular: d_amp and d_rot are fit from the same spirals, then A_z, R, and wind stress are derived from them, so Fig. 2's 'better fit' is by construction; independent support comes mainly from MITgcm.

  1. fitted input called prediction [Main text 'Comparison with the observed Ekman spiral'; SI 'Inversing damp and drot from the observed Ekman spiral' (Eqs. S17-S18, S11, S4)]
    "The damping-enabled Ekman theory developed here allows us to inverse the depth scales d_amp and d_rot (fig. S1), as well as A_z and R, from the observations of u(z) and v(z) (Materials and Methods and SI Appendix). The wind stress in here is determined by substituting the observed velocity at a certain depth into Eq. 2 rather than from the observations directly... Given d_amp and d_rot, the value of A_z and d_Ek can be obtained from Eq. S11, with which the value of R can be obtained from Eq. S4."

    The quantities displayed as 'predicted by the new theory' in Fig. 2 are not predicted: d_amp and d_rot are linear-fit slopes of the same observed log-speed and phase profiles (Eqs. S17-S18), and A_z, d_Ek, R, and the wind-stress vector are algebraically derived from those fitted values (Eqs. S11, S4, S13). The orange curve is therefore the assumed functional form with parameters estimated from the very spiral it is compared against. The 'better fit' than the classic one-parameter curve is guaranteed by the addition of two depth scales and a normalization point; it cannot discriminate damping from any other mechanism producing d_amp < d_rot, and it does not independently measure R.

  2. self definitional [Eq. (3) and SI Eq. (S11); the inversion procedure]
    "Since μ<1 and θ<45°, d_amp is smaller than d_Ek, while d_rot is larger than d_Ek (Fig. 1), with d_amp d_rot = d_Ek^2. Therefore, the flattening of Ekman spiral as observed in the ocean is a natural result of the damping."

    In the observational test, d_Ek is not an independently measured quantity: Eq. S11 (d_amp d_rot = d_Ek^2) is used to define d_Ek, and hence A_z, from the fitted d_amp and d_rot. The ordering d_amp < d_Ek < d_rot then follows automatically from d_amp < d_rot by the geometric-mean construction, and the product relation d_amp d_rot = d_Ek^2 is an identity imposed by the fitting algebra rather than an independently verified prediction. The observed flattening ratio is an input, not an output.

full rationale

The core algebra (Eqs. 1-5 and SI) is internally consistent: adding a linear damping term to the steady Ekman equation produces two depth scales with d_amp d_rot = d_Ek^2, a reduced surface velocity, and transport/pumping formulas containing both wind-stress curl and divergence. This mathematical part is not circular; it is a closed-form solution under the stated damping assumption. The circularity enters in the observational validation. The depth scales d_amp and d_rot are not measured independently; they are obtained by linear fits to the same log-speed and phase-versus-depth profiles (Eqs. S17-S18), and A_z, d_Ek, R, and the wind-stress magnitude/direction are then computed from those fitted values (Eqs. S11, S4, S13). Consequently, Fig. 2 compares the assumed functional form with parameters fit from the same data, and the highlighted relation d_amp < d_Ek < d_rot is imposed by defining d_Ek as sqrt(d_amp d_rot). The paper itself concedes that 'it is difficult to determine the value of R or Da' in the Discussion. The CESM xi-Da correlation (Fig. 5) is an independent comparison, though based on only five points and with xi not converted into an R value, so it is suggestive rather than conclusive. The MITgcm experiment is the cleanest independent support because R is estimated directly from the model's explicit horizontal-mixing term rather than from the spiral fit, and it confirms the mathematics of the model though not the ocean mechanism. Overall, the theoretical derivation is self-contained, but the central oceanographic claim that damping explains the observed universal flattening rests on a fitted input called a prediction, giving partial circularity rather than full circularity.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The central derivation adds one free parameter R, plus site-fitted Az and assumed Da values, and relies on the assumption that all unresolved processes collapse to a depth-uniform linear drag. No new physical entities are introduced.

free parameters (3)
  • R (damping coefficient) = per site values in Table S2: (6.7, 4.2, 0.7, 8.8, 5.8) x 10^-5 s^-1
    Inverted from the observed or simulated Ekman spiral through the fitted depth scales damp and drot (SI Eqs. S4, S17-S18). This is the central parameter of the new theory and is not independently constrained by turbulence measurements in the observations.
  • Da (damping number, R/f) = 0.76, 0.52, 0.28, 0.70, 0.56 for the five sites (Table S2); assumed value 0.2 for transport and upwelling estimates
    Derived from the fitted R and local Coriolis parameter. The 0.2 value used for the RAPID and global estimates is an assumed scenario, not a measured value.
  • Az (vertical eddy viscosity) = per site values in Table S2: (2.9, 2.5, 8.0, 10.4, 2.9) x 10^-2 m^2 s^-1
    Obtained from the product of the fitted damp and drot (Eq. S11). In the observations it is not independently measured, so it is part of the same inversion that produces R.
assumptions (6)
  • domain assumption Steady, horizontally homogeneous Ekman balance with constant Az and depth-uniform damping R
    This is the starting model (main text Eq. 1, SI Eq. S5). Real ocean Ekman layers are unsteady, inhomogeneous, and have vertically varying mixing; the paper chooses uniform R after arguing that other processes can be lumped into it.
  • ad hoc to paper The net effect of turbulent dissipation, lateral friction, waves, and stratification can be represented as a linear Rayleigh damping term R U
    The paper postulates this lumped representation, for example in the Discussion: 'the major effect of various processes that have been proposed in the literature can be lumped into a linear damping.' No derivation from a turbulence closure is given, and R is not measured independently for the observations.
  • domain assumption Surface boundary condition of specified wind stress and vanishing velocity at depth
    Standard Ekman boundary conditions used in solving Eq. 2 (SI Eq. S6). The wind stress magnitude and direction are later inferred from the current profile instead of taken from independent wind observations.
  • domain assumption The second-order longitudinal structure function follows the two-thirds law in an inertial range
    Used to estimate the dissipation index xi in CESM1.3 (SI Eqs. S19-S20). The authors note that the Kolmogorov constant is uncertain and that the Gulf Stream site shows anisotropy, so xi is only a qualitative proxy.
  • domain assumption Horizontal variation of R is neglected when deriving Ekman pumping
    The SI states that 'ignoring the horizontal variation of R gives' Eq. S15, so the divergence term in Eq. 5 omits spatial gradients of the damping coefficient.
  • standard math Equivalence between complex vertical viscosity and linear damping
    The algebraic derivation in the last SI section shows that a complex Az maps to a real effective Az plus a damping R. This is a mathematical equivalence, not independent physical evidence.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Ekman Theory with Damping." pith.science (2026). https://pith.science/paper/TE2T6GGV

@misc{pith2026250502068,
  author       = {Pith},
  title        = {Pith review of: Ekman Theory with Damping},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TE2T6GGV}},
  note         = {Machine review of arXiv:2505.02068}
}
read the original abstract

The observed Ekman spirals in the ocean are always "flatter" than that predicted by the classic theory. We propose that the universal flattening of Ekman spiral is mainly due to the damping associated with turbulent dissipation. Analytical solutions and numerical simulations show convincingly a better fitting between the new theory and observations. Most importantly, the new theory indicates that the damping can lead to weakened Ekman transport and pumping, with the latter not only driven by the curl but also the divergence of wind stress. Under a modest damping, the Ekman transport along 26.5{\deg}N will be ~0.4 Sv (12%) smaller than that predicted by the classic theory. Hence, the damping due to turbulent dissipation can noticeably affect the wind-driven circulation in the upper ocean.

Figures

Figures reproduced from arXiv: 2505.02068 by the authors.

Figure 1
Figure 1. Dependence of various variables on the damping number Da. (A) decay depth scales damp and drot normalized by the classic decay scale dEk. (B) flattening ratio damp/drot. (C) deflection angle of surface current. (D) the amplitude factor μ (solid line) in equation (2) [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 3
Figure 3. Influence of damping on Ekman transport. (A) zonally integrated Ekman transport over the ACC region for different damping number Da, averaged over the period 1991-2020; The wind stress used in the calculations is from GODAS (B) zonally integrated annual mean Ekman transport along 26.5°N, where RAPID array is installed, calculated using the classic Ekman theory (black) and damping-enabled theory (purple) with Da = 0.… view at source ↗
Figure 4
Figure 4. Influence of damping on Ekman pumping. (A) and (B) show the Ekman pumping (positive upward) calculated with Da being 0.0 (classic Ekman theory) and 0.2, respectively. (C) shows the difference between (B) and (A) (former minus latter). (D) is the same as (C) except that here the relative changes are shown. The wind stress used in the calculations is from GODAS [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figures from the paper (1 more)
Figure 5
Figure 5. Figure 5: The relationship between the damping (ξ) estimated utilizing turbulence theory and (A) the damping coefficient R or (B) the damping number Da. The meaning of ξ can be found from Materials and Methods and SI Appendix. 10 -5 10 -4 10 -2 10 -1 A Equatorial Pacific The Coa…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

10 extracted references · 10 canonical work pages

  1. [1]

    The wavy Ekman layer: Langmuir circulations, breaking waves, and Reynolds stress

    McWilliams JC, Huckle E, Liang J-H et al. The wavy Ekman layer: Langmuir circulations, breaking waves, and Reynolds stress. Journal of Physical Oceanography. 2012; 42(11): 1793-1816

    work page 2012

  2. [2]

    Wind-driven ocean currents and Ekman transport

    Price JF, Weller RA, Schudlich RR. Wind-driven ocean currents and Ekman transport. Science. 1987; 238(4833): 1534-1538

    work page 1987

  3. [3]

    Direct observations of the Ekman balance at 10 N in the Pacific

    Wijffels S, Firing E, Bryden H. Direct observations of the Ekman balance at 10 N in the Pacific. Journal of Physical Oceanography. 1994; 24(7): 1666-1679

    work page 1994

  4. [4]

    Direct evidence for an Ekman balance in the California Current

    Chereskin T. Direct evidence for an Ekman balance in the California Current. Journal of Geophysical Research: Oceans. 1995; 100(C9): 18261-18269

    work page 1995

  5. [5]

    Observations of Ekman currents in the Southern Ocean

    Lenn Y-D, Chereskin TK. Observations of Ekman currents in the Southern Ocean. Journal of Physical Oceanography. 2009; 39(3): 768-779

    work page 2009

  6. [6]

    Detecting and characterizing Ekman currents in the Southern Ocean

    Roach CJ, Phillips HE, Bindoff NL et al. Detecting and characterizing Ekman currents in the Southern Ocean. Journal of Physical Oceanography. 2015; 45(5): 1205-1223

    work page 2015

  7. [7]

    Observations of seasonal variation in the Ekman layer

    Schudlich RR, Price JF. Observations of seasonal variation in the Ekman layer. Journal of physical oceanography. 1998; 28(6): 1187-1204

    work page 1998

  8. [8]

    An unprecedented set of high‐resolution earth system simulations for understanding multiscale interactions in climate variability and change

    Chang P, Zhang S, Danabasoglu G et al. An unprecedented set of high‐resolution earth system simulations for understanding multiscale interactions in climate variability and change. Journal of Advances in Modeling Earth Systems. 2020; 12(12): e2020MS002298

    work page 2020

Show all 10 references
  1. [9]

    Kolmogorov constants for the second-order structure function and the energy spectrum

    Ni R, Xia K-Q. Kolmogorov constants for the second-order structure function and the energy spectrum. Physical Review E—Statistical, Nonlinear, and Soft Matter Physics. 2013; 87(2): 023002

    2013

  2. [10]

    The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers

    Kolmogorov AN. The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers. Proceedings of the Royal Society of London Series A: Mathematical and Physical Sciences. 1991; 434(1890): 9-13

    1991

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.