REVIEW 2 major objections 6 minor 18 references
Shear preserves plume reach while rewriting its spectral pathway
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
Imposed mean shear preserves the vertical reach and contact timing of finite-depth salt-finger plume forests while redistributing their spectral energy pathway.
T0 review reviewed 2026-07-09 challenge →
load-bearing objection A clean simulation study with one real gap: the sheared case never gets the threshold-sensitivity test that the paper's own logic demands. the 2 major comments →
Route survival and spectral modification of finite-depth salt-finger plume forests under imposed mean shear
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
The central discovery is a decoupling between finite-depth reach and spectral pathway in salt-finger plume forests. An imposed mean shear preserves the vertical reach and contact timing of a delayed mixed route while redistributing the planform spectral content that maintains that reach. Specifically, the sheared plume forest reaches the remote layers at the same times as the unsheared reference, but the intermediate spectral fraction drops to 0.530 times the unsheared value and the short-wave fraction rises to 1.278 times, showing that the plume forest can sustain the same vertical connection through a different flow organization.
What carries the argument
The machinery is a set of direct three-dimensional nonhydrostatic simulations of a two-layer salt-finger system at fixed density ratio, Prandtl number, and diffusivity ratio. The route is defined by coupled measures: active vertical width, distance-based contact timing with relaxation zones, planform spectral partition into broad, intermediate, and short-wave bands, and down-gradient salinity flux. A mixed-seed replicate defines a tolerance layer for natural variability. The shear perturbation is a single initial tanh velocity profile imposed at initialization and then allowed to evolve freely with the flow, testing whether the coupled system preserves or destroys the delayed route.
Load-bearing premise
The route-survival claim rests on a single shear amplitude, a single shear profile shape, a single density ratio, and a single mixed-route realization. The mixed route's proximity to the scalar-contact boundary means the contact-timing match could be fragile under small parameter changes.
What would settle it
If a different shear amplitude or a different initial roughness realization produced a case where the contact timing shifted away from the mixed-route values while the spectral partition stayed within tolerance, the decoupling claim would fail. The route-survival result specifically requires that reach and spectral pathway respond independently; finding them coupled under any perturbation would undermine the central claim.
If this is right
- If reach and spectral pathway are genuinely decoupled, transport parameterizations for double-diffusive staircases may need separate treatment of vertical connection efficiency and planform structure rather than a single flux law.
- The route-survival concept could be tested in laboratory salt-finger experiments with imposed shear layers, where contact timing and spectral content are independently measurable.
- If the decoupling holds across density ratios, it would suggest that oceanic salt-finger interfaces embedded in sheared currents can maintain their vertical exchange role while the horizontal structure of that exchange shifts.
- The spectral modification pattern, with reduced intermediate scales and elevated short-wave content, could serve as a diagnostic signature of shear interaction in field observations of thermohaline staircases.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript presents direct three-dimensional simulations of finite-depth salt-finger convection to test whether a 'route'—defined as a coupled state of contact timing, active vertical width, spectral branch, and scalar exchange—survives an imposed mean-shear perturbation. The author first establishes a no-shear 'route atlas' by varying only the interfacial roughness spectrum (low-mode, high-annulus, mixed, and a mixed-seed replicate), holding all other parameters fixed. The mixed route is then perturbed by an initial tanh mean shear profile (Eq. 1). The central finding is that the sheared run preserves finite-depth reach and contact timing (matching the unsheared mixed reference at t=57.75 for velocity and t=59.5 for salinity) while redistributing the spectral branch: broad fraction enhanced, intermediate fraction halved, and short-wave fraction elevated. The study is well-scoped, the metrics are clearly defined (Eqs. 2–7), and the mixed-seed tolerance layer provides an internal reproducibility check. The resolution support (Table 6) demonstrates that high-wavenumber tails are smaller on a finer grid, supporting the baseline.
Significance. The paper introduces a useful conceptual framework—'route survival with spectral modification'—that decouples finite-depth reach from spectral pathway in double-diffusive convection. This decomposition is physically motivated and could inform future comparisons of finite-depth mixing events where interfaces carry inherited roughness or are embedded in larger-scale motion. The simulation setup is clearly specified and reproducible (Oceananigans, fixed parameters, Zenodo archive [3]). The falsifiable prediction—that reach and spectral partition can respond independently to the same perturbation—is a concrete contribution. The threshold-sensitivity analysis of the no-shear atlas (Figure 3) is a commendable methodological detail that strengthens the internal validity of the route-family comparison.
major comments (2)
- §3.2, Figure 3: The paper's own logic creates a gap in the shear case. The manuscript demonstrates that the mixed route sits at the scalar-contact boundary—binary salinity contact occurs in only 2/4 tested thresholds for M and 0/4 for M' (§3.2)—and explicitly states that 'the route-family comparison is based on the quantities that are stable under a seed change, while the thresholded contact flags identify the narrow part of the response that is intrinsically sensitive.' However, Figure 3 includes L, M, M', and H but NOT M+Us. The sheared case's contact-timing match (t=57.75 for w, t=59.5 for salinity) is reported only for the 'standard threshold' (§2.2) without any threshold-sensitivity test. If M+Us salinity contact is similarly threshold-sensitive (e.g., contacting in only 1/4 or 2/4 thresholds), then the 'preserved contact timing' claim reduces to a coincidence of one threshold value
- §2.3, §3.1, Table 2: The route-connection index Rc and plume-coherence index Cp are used as primary ordering quantities for the route atlas (Table 2, Figure 1), but their definitions are not given in the manuscript. The reader cannot verify how Rc is computed or whether it is constructed from the other listed components (contact, vertical extent, area, spectral, probe). If Rc is a weighted composite of those components, the claim that routes 'separate along several axes' (§3.1) is partly circular. Please provide the formulas or explicit construction for Rc and Cp, or restructure the atlas comparison around the explicitly defined metrics (active width, contact time, spectral fractions, flux) from §2.2–2.3.
minor comments (6)
- §2.2: The 'standard threshold' value of α used for the main contact-time results is not stated. Please specify it.
- §2.1: The far-field 'relaxation zones' are mentioned but their vertical extent and sponge-layer formulation are not described. A brief statement of their geometry should be added.
- Table 1: The symbol M' (mixed-seed) is visually similar to a prime notation that could be confused with a derivative. Consider a more distinct label (e.g., M_seed).
- §3.6: The statement 'the intermediate fraction is nearly halved' is supported by the ratio 0.530, but the absolute values (0.0801 vs. 0.1510) should be stated alongside ratios throughout for clarity.
- Figure 8: The x-axis label 'w active width' and y-axis 'mean intermediate spectral fraction' are clear, but the time window over which the 'mean' is computed should be stated in the caption.
- References [1] and [2] are both by the present author and dated 2026. Their status (published, submitted, preprint) should be clarified for readers assessing novelty.
Circularity Check
No significant circularity; self-citations are contextual framing, not load-bearing inputs
full rationale
The paper's derivation chain is self-contained. The route metrics (active width in §2.2 Eqs. 2-4, spectral fractions in §2.3 Eq. 5, salt flux in §2.3 Eq. 6) are defined directly from simulation fields and computed from DNS output. The shear perturbation (Eq. 1) is an imposed initial condition with fixed parameters (U0=0.10, δs=2hi), not a fitted quantity. The route-survival test (§2.4) compares sheared vs. unsheared metrics that are independently computed — no 'prediction' reduces to a fit. The two self-citations ([1] Kalathoor 2026a, [2] Kalathoor 2026b) provide reference states and terminology ('route-memory language follows the spectral-memory framing in Kalathoor [2]'), but the route atlas itself is defined by the four simulations in this paper (Table 1), not imported from prior work. The central claim — that shear preserves reach while changing spectral pathway — emerges from comparing independently computed quantities (active widths within 1.5%/0.2% of mixed, spectral fractions shifted to 1.116x/0.530x/1.278x). No step in this chain is equivalent to its inputs by construction. The skeptic's threshold-sensitivity concern for the sheared case's contact timing is a robustness/correctness issue, not a circularity: the contact times are computed from simulation output via Eq. 4, not defined in terms of the survival claim. Score 1 reflects the minor self-citations that provide framing but do not load-bear the central result.
Axiom & Free-Parameter Ledger
free parameters (4)
- U0 =
0.10
- δs =
2*hi = 6
- α (threshold fraction) =
not explicitly stated
- Spectral band boundaries (broad, intermediate, short-wave) =
not explicitly stated
axioms (4)
- domain assumption Nonhydrostatic Boussinesq equations with salt-finger favorable stratification govern the finite-depth plume forest evolution.
- domain assumption The far-field relaxation zones adequately represent finite-depth boundaries without reflecting artifacts.
- ad hoc to paper A single mixed-seed replicate defines an adequate tolerance layer for route-family variability.
- ad hoc to paper The tanh initial shear profile is a representative perturbation for testing route survival.
invented entities (2)
-
Route (as a coupled state variable)
no independent evidence
-
Route-connection index Rc
no independent evidence
Cite this review
Pith. "Pith review of Route survival and spectral modification of finite-depth salt-finger plume forests under imposed mean shear." pith.science (2026). https://pith.science/paper/E7AVKLVC
@misc{pith2026260707012,
author = {Pith},
title = {Pith review of: Route survival and spectral modification of finite-depth salt-finger plume forests under imposed mean shear},
year = {2026},
howpublished = {\url{https://pith.science/paper/E7AVKLVC}},
note = {Machine review of arXiv:2607.07012}
}
abstract
Salt-finger plume forests in a finite layer can differ in strength and in the route by which interfacial activity becomes vertically connected. We use direct three-dimensional simulations to test whether such a route is a short-lived realization-specific transient or a persistent route family under an added mean-shear perturbation. The baseline route atlas holds density ratio, diffusivity ratio, Prandtl number, interface thickness, roughness amplitude, domain, and resolution fixed while varying the imposed interfacial roughness spectrum. Low-mode roughness forms a broad connecting endpoint, high-annulus roughness forms a localized route-memory endpoint, and mixed roughness forms a delayed scale-transfer route. A second mixed realization preserves continuous active-width, spectral, and transport measures after \(t=45\), with mean absolute differences of \(3.1\%\) in \(w\)-active width, \(1.6\%\) in salinity-active width, \(2.8\%\) in broad spectral fraction, and \(3.6\%\) in salt flux, while shifting the binary scalar-contact label. We then impose an initial tanh mean shear on the mixed route. The full-resolution shear case reaches \(t=60\) and preserves finite-depth reach: first velocity contact occurs at \(t=57.75\), first salinity contact occurs at \(t=59.5\), and both times match the unsheared mixed reference. The spectral branch is redistributed. At \(t=60\), the broad fraction is \(1.116\) times the mixed value, the intermediate fraction is \(0.530\) times the mixed value, and the short-wave fraction is \(1.278\) times the mixed value. In this finite-depth configuration, route survival means preserved reach and contact timing with a changed spectral pathway.
Reference graph
Works this paper leans on
-
[2]
Interfacial Spectral Memory as a State Variable for Finite-Depth Salt-Finger Exchange
Kalathoor SP (2026) Interfacial spectral memory as a state variable for finite-depth salt-finger exchange. https://doi.org/10.48550/arXiv.2606.26045, URL https:// doi.org/10.48550/arXiv.2606.26045, arXiv:2606.26045 17
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2606.26045 2026
-
[3]
Kalathoor SP (2026) Reduced data for route survival and spectral modification of finite-depth salt-finger plume forests under imposed mean shear. https://doi. org/10.5281/zenodo.21081376, URL https://doi.org/10.5281/zenodo.21081376, version 3
-
[4]
Journal of Marine Research 45(3):533–556
Kunze E (1987) Limits on growing, finite-length salt fingers: A richardson number constraint. Journal of Marine Research 45(3):533–556. https://doi.org/10.1357/ 002224087788326885
work page 1987
-
[5]
Radko T (2013) Double-Diffusive Convection. Cambridge University Press
work page 2013
-
[6]
Journal of Open Source Software 5(53):2018
Ramadhan A, Wagner GL, Hill C, et al (2020) Oceananigans.jl: Fast and friendly geophysical fluid dynamics on GPUs. Journal of Open Source Software 5(53):2018. https://doi.org/10.21105/joss.02018
-
[7]
Deep-Sea Research Part A 26(1):23–40
Schmitt RW (1979) The growth rate of super-critical salt fingers. Deep-Sea Research Part A 26(1):23–40. https://doi.org/10.1016/0198-0149(79)90083-9
-
[8]
Annual Review of Fluid Mechanics 26(1):255–285
Schmitt RW (1994) Double diffusion in oceanography. Annual Review of Fluid Mechanics 26(1):255–285. https://doi.org/10.1146/annurev.fl.26.010194.001351
-
[9]
SchmittRW,LedwellJR,MontgomeryET,etal(2005)Enhanceddiapycnalmixing by salt fingers in the thermocline of the tropical atlantic. Science 308(5722):685–
work page 2005
-
[10]
https://doi.org/10.1126/science.1108678
-
[11]
Silvestri S, Wagner GL, Hill C, et al (2023) Oceananigans.jl: A Julia library that achieves breakthrough resolution, memory and energy efficiency in global ocean simulations. https://doi.org/10.48550/arXiv.2309.06662, URL https://arxiv.org/ abs/2309.06662, arXiv:2309.06662
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2309.06662 2023
-
[12]
Direct Numerical Simulation of 3D Salt Fingers: From Secondary Instability to Chaotic Convection
Simeonov JA, Stern ME, Radko T (2009) Direct numerical simulation of 3d salt fingers: From secondary instability to chaotic convection. https://doi.org/10. 48550/arXiv.0910.2434, URL https://arxiv.org/abs/0910.2434, arXiv:0910.2434
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[13]
Stellmach S, Traxler AL, Garaud P, et al (2011) Dynamics of fingering convection. Part 2. The formation of thermohaline staircases. Journal of Fluid Mechanics 677:554–571. https://doi.org/10.1017/jfm.2011.99
-
[14]
Stern ME (1960) The salt-fountain and thermohaline convection. Tellus 12(2):172–
work page 1960
-
[15]
https://doi.org/10.3402/tellusa.v12i2.9378
-
[16]
Journal of Fluid Mechanics 35(2):209–218
Stern ME (1969) Collective instability of salt fingers. Journal of Fluid Mechanics 35(2):209–218. https://doi.org/10.1017/S0022112069001066
-
[17]
Traxler AL, Stellmach S, Garaud P, et al (2011) Dynamics of fingering convection. Part 1. Small-scale fluxes and large-scale instabilities. Journal of Fluid Mechanics 677:530–553. https://doi.org/10.1017/S002211201100098X 18
-
[18]
Annual Review of Fluid Mechanics 6(1):37–54
Turner JS (1974) Double-diffusive phenomena. Annual Review of Fluid Mechanics 6(1):37–54. https://doi.org/10.1146/annurev.fl.06.010174.000345
-
[19]
https://doi.org/10.48550/arXiv
Wagner GL, Silvestri S, Constantinou NC, et al (2025) High-level, high-resolution ocean modeling at all scales with Oceananigans. https://doi.org/10.48550/arXiv. 2502.14148, URL https://arxiv.org/abs/2502.14148, arXiv:2502.14148 19
work page internal anchor Pith review doi:10.48550/arxiv 2025
This paper was first reviewed by glm-5.2 on July 9, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.