arxiv: 2604.27200 · v2 · submitted 2026-04-29 · ⚛️ physics.ao-ph

Recognition: no theorem link

Estimating the Kinetic Energy Spectrum from the Second-Order Velocity Structure Function using a Regularized Fitting Approach

Ayantika Bhattacharjee, Dhruv Balwada, Manuel Gutierrez-Villanueva, Shane Elipot, Spencer Jones

Pith reviewed 2026-05-12 01:46 UTC · model grok-4.3

classification ⚛️ physics.ao-ph

keywords kinetic energy spectrumvelocity structure functionocean turbulenceLagrangian observationsregularized inversionspectral estimationdrifter data

0 comments

The pith

A regularized optimization method recovers the kinetic energy spectrum from second-order structure functions derived from sparse Lagrangian data.

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

The paper develops a method to estimate the kinetic energy wavenumber spectrum by inverting its relationship to the second-order velocity structure function. Direct inversion is unstable because the structure function is a non-local weighted integral over all scales and because observations are limited by sampling and noise. The authors assume the spectrum consists of a finite number of segments each having its own slope and amplitude, then cast the inversion as an optimization problem that adjusts these parameters until the forward integral matches the observed structure function. Validation on idealized spectra shows exact recovery of the input parameters; tests on high-resolution model output confirm robustness to added noise; application to GLAD and LASER drifter data demonstrates the method works with real sparse Lagrangian observations.

Core claim

By representing the kinetic energy spectrum as a sum of a small number of power-law segments and solving an optimization problem whose objective is to reproduce the measured second-order structure function through the known forward integral relation, the procedure stably recovers the segment slopes and amplitudes even when direct numerical inversion of the relationship fails.

What carries the argument

The regularized optimization that parameterizes the spectrum as a finite collection of segments with adjustable slopes and amplitudes, then minimizes the mismatch between the forward-computed structure function and the observed one.

If this is right

The kinetic energy spectrum becomes obtainable from Lagrangian drifter trajectories without access to gridded Eulerian velocity fields.
Reconstruction accuracy persists when realistic levels of noise are present in the structure-function measurements.
Spectral diagnostics can be performed on any dataset that supplies pairwise velocity differences at varying separations.
Application to the GLAD and LASER drifter experiments produces usable kinetic energy spectra from surface ocean observations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same fitting procedure could be tested on atmospheric velocity data or on tracer variance structure functions to see whether the finite-segment assumption holds in other turbulence contexts.
If the ocean spectrum contains sharp transitions or continuous curvature that cannot be captured by a few power laws, the recovered parameters may represent an averaged approximation rather than the true distribution.
Systematic application across multiple drifter campaigns could map regional differences in the slopes of the recovered segments and link them to local forcing or stratification.
The optimization framework might be extended by allowing the number of segments itself to be chosen adaptively rather than fixed in advance.

Load-bearing premise

That the true kinetic energy spectrum consists of a small number of distinct power-law segments whose slopes and amplitudes can be stably recovered by optimization.

What would settle it

Apply the method to structure functions computed from a high-resolution ocean model whose true kinetic energy spectrum is independently known from a direct Fourier transform on the same gridded velocity field; large systematic differences between recovered and direct spectra would falsify the claim.

read the original abstract

Ocean turbulence plays a key role in shaping large-scale circulation, heat uptake, and biogeochemical processes. The kinetic energy (KE) wavenumber spectrum is a fundamental diagnostic, quantifying how KE is distributed across spatial scales. The second-order structure function -- computed from velocity differences between spatially separated observations -- provides a complementary measure, but unlike the KE spectrum, it reflects a non-local, weighted integral of KE over all scales. Analytic relationships link the two metrics, permitting forward and inverse transformations between them. However, recovering the KE spectrum from the structure function via the inverse relationship is highly sensitive to sampling limitations and numerical discretization errors. Here we propose a regularized approach in which the spectrum is assumed to consist of a finite number of segments with distinct slopes and amplitudes, and the inversion is formulated as an optimization problem. The approach is first validated in an idealized setting; for a number of idealized KE spectra with prescribed sets of spectral slopes and amplitudes, the corresponding structure functions are computed by numerically evaluating the forward relationship. These structure functions are then used to determine the underlying parameters using our proposed approach, which shows that we are able to perfectly recover the parameters and consequently the KE spectra. The method is further evaluated on high-resolution ocean model output, where it reconstructs the underlying spectra well even in the presence of noise. Finally, we apply the method to surface drifter observations (GLAD and LASER experiments). The results show that the framework enables estimation of the KE spectrum from sparse Lagrangian data, extending spectral diagnostics beyond gridded Eulerian measurements.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a workable way to turn structure functions from drifters into KE spectra by fitting a few power-law pieces, but the results hinge on that piecewise assumption holding.

read the letter

The main takeaway is that this work offers a practical route to estimating kinetic energy spectra from the second-order structure functions that are easy to compute on sparse Lagrangian data like drifters. By casting the inversion as a regularized fit to a small number of piecewise power-law segments, they avoid the instability that comes with direct numerical inversion. They demonstrate exact recovery when the synthetic spectra are built exactly from the assumed form. On noisy output from an ocean model the reconstruction stays close to the true spectrum. The final step applies it to real surface drifter data from the GLAD and LASER experiments, which is where the value shows up for people who have lots of trajectory data but no gridded velocity fields. The approach is new in its specific formulation of the regularized segment fit. It does a solid job of showing the method works under the conditions it was designed for. The soft spot is the reliance on the spectrum being well approximated by a few straight segments in log-log space. When the actual spectrum has smooth transitions or more wiggles, the fit will still produce something, but it may not reflect the true distribution of energy. The paper's tests do not include cases that deliberately violate the piecewise assumption to measure the resulting error. Without that, it is hard to know how much the results on the drifter data are shaped by the model choice rather than the data. They also do not give much on how the number of segments is chosen or how results change with that choice. This paper is aimed at oceanographers and fluid dynamicists who analyze turbulence from drifter or float observations. A reader looking for a tool to get spectral information from such data will get a concrete method they can try. It is not revolutionary, but it is a useful incremental step that addresses a real practical problem. I would send it to peer review. The core idea is sound enough and the application is relevant, even if the validation could be more thorough on the limitations of the parameterization.

Referee Report

3 major / 0 minor

Summary. The paper claims to develop a regularized optimization approach for inverting the second-order velocity structure function to recover the kinetic energy wavenumber spectrum, under the assumption that the spectrum is composed of a finite number of segments with distinct power-law slopes and amplitudes. It demonstrates exact parameter recovery in idealized tests, good performance on noisy model data, and applies it to GLAD and LASER drifter observations to estimate spectra from Lagrangian data.

Significance. If valid, the method would allow spectral diagnostics of ocean turbulence from sparse drifter data, extending beyond Eulerian gridded measurements. Strengths include the exact recovery in forward tests matching the model and handling of noise in model output, providing a practical tool if the segment assumption holds for real spectra.

major comments (3)

[Idealized setting] The exact recovery of parameters is shown for idealized KE spectra with prescribed segments, but this test is circular as the data is generated from the same functional form assumed in the inversion; it does not demonstrate robustness to deviations from piecewise power laws.
[High-resolution ocean model output] The reconstruction is reported as good even with noise, but the manuscript lacks quantitative metrics such as error norms or direct comparisons to spectra computed from the model grid, and does not report sensitivity to the choice of number of segments, which is a key free parameter.
[Application to drifter observations] The extension to sparse Lagrangian data from GLAD and LASER is promising, but without tests on spectra with smooth transitions or additional features, it is unclear if the method recovers the true spectrum or fits the assumed segmented form.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and insightful comments. We address each major comment point by point below, providing the strongest honest defense of the manuscript while acknowledging where revisions are warranted to improve clarity and robustness.

read point-by-point responses

Referee: The exact recovery of parameters is shown for idealized KE spectra with prescribed segments, but this test is circular as the data is generated from the same functional form assumed in the inversion; it does not demonstrate robustness to deviations from piecewise power laws.

Authors: We agree that the idealized tests use spectra generated from the same segmented power-law form assumed by the inversion method. This is intentional as a first validation step to confirm that the regularized optimization can exactly recover parameters and spectra when the underlying assumption holds, thereby demonstrating that the forward-inverse relationship is numerically invertible without instability under ideal conditions. However, we acknowledge that this does not address robustness to model mismatch. In the revised manuscript we will add synthetic tests using spectra with smooth transitions (e.g., continuous -5/3 to -3 roll-off without sharp breaks) to quantify how the method approximates non-segmented forms. revision: yes
Referee: The reconstruction is reported as good even with noise, but the manuscript lacks quantitative metrics such as error norms or direct comparisons to spectra computed from the model grid, and does not report sensitivity to the choice of number of segments, which is a key free parameter.

Authors: This observation is correct; the current text describes the model-output results qualitatively without providing error norms, grid-based Fourier comparisons, or sensitivity tests to segment number. We will revise the relevant section to include: (i) L2-norm and relative-error metrics between recovered and directly computed spectra, (ii) side-by-side plots of both spectra for the noisy and noise-free cases, and (iii) a sensitivity study varying the number of segments from 2 to 5, reporting changes in fit residual and recovered spectral slopes. revision: yes
Referee: The extension to sparse Lagrangian data from GLAD and LASER is promising, but without tests on spectra with smooth transitions or additional features, it is unclear if the method recovers the true spectrum or fits the assumed segmented form.

Authors: Because the true spectrum is unknown for the real drifter observations, direct validation against ground truth is not feasible. The segmented assumption is motivated by established turbulence theory (e.g., expected transitions between inertial and dissipation ranges). We will add a dedicated limitations paragraph discussing potential biases when the true spectrum deviates from piecewise power laws. In addition, we will include the same synthetic non-segmented tests mentioned above to illustrate the approximation behavior of the method on data that violate the assumption. revision: partial

Circularity Check

0 steps flagged

No significant circularity; central claim rests on independent validation

full rationale

The paper derives an analytic forward relationship between the second-order structure function and the KE spectrum, then proposes a regularized optimization that assumes the spectrum takes the form of a finite number of power-law segments with constant slopes and amplitudes. Idealized tests recover the input parameters exactly because the synthetic structure functions are generated from spectra matching that exact functional family, which is expected behavior for any consistent fitting procedure and does not constitute a self-definitional reduction. The method is then applied to high-resolution ocean model output (with added noise) and to real surface drifter observations from the GLAD and LASER experiments; these constitute external data sources whose spectra are not constructed from the paper's ansatz. No load-bearing self-citations, uniqueness theorems imported from prior author work, or renaming of known results appear in the derivation chain. The claim that the framework extends spectral estimation to sparse Lagrangian data therefore retains independent empirical content.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The method rests on the standard analytic integral relationship between the second-order structure function and the KE spectrum, plus the modeling choice that the spectrum can be represented by a small number of linear segments whose slopes and amplitudes are free parameters fitted by optimization.

free parameters (2)

number of segments
Chosen by the user or optimization procedure to represent the spectrum shape
slope and amplitude of each segment
Fitted via the regularized optimization to match the observed structure function

axioms (1)

domain assumption Analytic forward relationship linking structure function to integrated KE spectrum
Standard result in turbulence theory invoked for both forward computation and the inverse problem

pith-pipeline@v0.9.0 · 5598 in / 1261 out tokens · 52340 ms · 2026-05-12T01:46:20.298766+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

[1]

In the limit𝑟→ ∞,𝐽 0(𝑘𝑟) →0, so Eq

Derivation of KE Spectrum from the Second-Order Structure Function The forward relation between the second-order structure function,𝑆𝐹2(𝑟), and the kinetic energy spectrum,𝐸(𝑘), is 𝑆𝐹2(𝑟)=2 ∫ ∞ 0 𝐸(𝑘) [1−𝐽 0(𝑘𝑟) ] 𝑑𝑘,(1) where𝐽 0 is the zeroth-order Bessel function of the first kind. In the limit𝑟→ ∞,𝐽 0(𝑘𝑟) →0, so Eq. 1 reduces to, 𝑆𝐹2(∞)=2 ∫ ∞ 0 𝐸(𝑘)𝑑𝑘 ...

work page
[2]

The discrete wavenum- ber points used in the reconstruction fall either within a segment or on its boundaries, as determined by the algorithm

Piecewise Power-Law Representation of the Reconstructed Model Spectra Figure S2 shows the reconstructed kinetic energy spectra for the MITgcm llc4320 model for representative winter and summer months, with gray dashed vertical lines indicating the segment boundaries, which are uniformly spaced in logarithmic wavenumber space. The discrete wavenum- ber poi...

work page
[3]

Best-fit Error Distribution for GLAD and LASER𝑆𝐹 2 Figure S3 shows the distribution of best-fit errors and the optimal number of spectral segments for each of the 600 bootstrap realizations of the GLAD (left) and LASER (right) datasets. For both datasets,𝑆𝐹 2 is evaluated at𝑁=21 discrete separation distances, so the optimization is performed by varying th...

work page