Imagine you're averaging GPS tracks of a hundred runners circling a track. Each runner's path is slightly different — some cut corners, some drift wide. You average their positions at each time step and get a "mean runner." But here's the problem: that mean runner's path quietly shrinks inward. The track gets smaller. Nobody actually ran a shorter loop, but the averaging process itself created a phantom contraction. That's the exact bug in generalised Lagrangian mean (GLM) theory that this paper fixes. The committed claim: by replacing the standard GLM mean flow map with its closest volume-preserving approximation — identified via polar factorization from optimal-transport theory — the authors eliminate the spurious divergence that has plagued Lagrangian averaging of incompressible flows since Andrews & McIntyre introduced GLM in 1978. This is a methodological fix to a known, well-documented defect, not a discovery of a new physical phenomenon. The architecture is clean and mathematically principled. Polar factorization decomposes any map into a volume-preserving part and a gradient-of-convex part — the same decomposition that Brenier proved in optimal transport. The authors keep only the volume-preserving factor. Their numerical implementation builds on recently developed on-the-fly Lagrangian-mean algorithms using exponential and Butterworth filters, which handle the time-averaging step. The computational overhead sits on top of existing simulation pipelines for 2D incompressible and shallow-water models. On the ladder: the baseline here is standard GLM itself. The paper doesn't claim to beat GLM at predicting anything new — it claims to remove an artifact. The comparison is direct: Lagrangian-mean fields computed with and without the volume-preservation constraint, shown side by side. The key visual result is that vortices no longer shrink under averaging. There's no external benchmark because the problem is internal consistency, not competitive accuracy against another method. Integrity is reasonable for a methods paper in mathematical fluid dynamics. The validation is same-team numerical simulation on two model problems (2D incompressible Euler and shallow-water equations). These are standard test beds, not cherry-picked. The comparison against unconstrained GLM is the right one — that's the thing they're fixing. But there's no independent replication, no community benchmark suite, and the demonstrations are 2D only. The milestone question is where it gets concrete: the obvious next number is 3D. Everything demonstrated here is two-dimensional. Real ocean and atmosphere models run in three dimensions, where the optimal-transport computation (finding the closest volume-preserving map) becomes substantially harder — both algorithmically and in compute cost. If the method scales to 3D turbulent flow simulations at resolutions used in operational climate models, it becomes a practical tool. If the 3D optimal-transport step is prohibitively expensive, this stays a theoretical contribution. The successor experiment the authors did not run is 3D. My honest read: this is option (c) — they're saving it for the next paper. The mathematical framework generalizes to any dimension, the 2D proof of concept is complete, and extending the optimal-transport numerics to 3D is a substantial enough project to warrant its own publication. The shallow-water demonstration hints at geophysical relevance without committing to the full 3D atmosphere/ocean problem.
