pith. sign in

arxiv: 2509.08856 · v1 · pith:7MBB5LABnew · submitted 2025-09-09 · ⚛️ physics.geo-ph · physics.flu-dyn

The reflection-transmission problem for inertial waves on geostrophic shear layers

Lennart Kira , Jerome Noir , Daniel Lecoanet This is my paper

Pith reviewed 2026-05-21 23:04 UTC · model grok-4.3

classification ⚛️ physics.geo-ph physics.flu-dyn
keywords inertial wavesgeostrophic shear layersreflection and transmissionquasi-two-dimensional modelray theorywave tunnelinglow-pass filteringplanetary fluid dynamics
0
0 comments X

The pith

Inertial waves partially reflect from subcritical geostrophic shear layers and tunnel through thin supercritical ones.

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

The paper develops a quasi-two-dimensional analytical model for how inertial waves reflect and transmit across a localized geostrophic shear layer whose width can be any size relative to the wavelength. It shows that partial reflections still happen even when the shear is subcritical, which ray theory does not predict. In thin supercritical layers, waves can pass through by tunneling with nearly complete transmission. Supercritical layers therefore act as low-pass filters that let only low-wavenumber waves continue. The model is checked against numerical simulations and accounts for the full wave field around and inside the layer.

Core claim

We develop a quasi-two-dimensional analytical model to investigate the reflection and transmission of inertial waves in the presence of a localized geostrophic shear layer of arbitrary width. In contrast to ray theory predictions, partial reflections occur even in subcritical shear layers and tunnelling with almost total transmission is possible in supercritical shear layers, if the layer is thin compared to the wavelength. That is, supercritical shear layers act as low-pass filters for inertial wave beams allowing the low-wavenumber waves to travel through. Our analytical model allows us to predict interactions between inertial waves and geostrophic shear layers not addressed by ray-basedor

What carries the argument

Quasi-two-dimensional analytical model that matches wave solutions analytically across the shear-layer boundaries.

If this is right

  • Partial reflections occur even in subcritical shear layers.
  • Nearly total transmission by tunneling occurs when supercritical layers are thin relative to the wavelength.
  • Supercritical shear layers filter inertial wave beams and pass only low-wavenumber components.
  • The model gives the full wavefield behavior around and inside the layer.
  • Interactions not captured by ray theory or statistical approaches can now be predicted.

Where Pith is reading between the lines

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

  • The low-pass filtering may change how angular momentum and heat are redistributed in planetary cores or stellar interiors.
  • Laboratory rotating-fluid experiments could directly measure the tunneling threshold for chosen layer widths.
  • The matching technique might be adapted to other wave types such as internal gravity waves on stratified shear layers.

Load-bearing premise

The quasi-two-dimensional approximation is valid for describing the interaction of inertial waves with a localized geostrophic shear layer of arbitrary width.

What would settle it

A numerical simulation or laboratory measurement that finds zero partial reflection for waves on a clearly subcritical shear layer or finds no tunneling transmission through a thin supercritical layer.

Figures

Figures reproduced from arXiv: 2509.08856 by Daniel Lecoanet, Jerome Noir, Lennart Kira.

Figure 1
Figure 1. Figure 1: Schematic figure demonstrating the geometry of the problem. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Model of plane inertial waves interacting with a single layer of constant shear. [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Model of multiple consecutive shear layers traversed by plane inertial waves. [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left: Segmentation of a Gaussian shear profile using [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Typical setup for the DNS using Dedalus. A shear layer is in the middle of the domain. [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Left: Kinetic energy field—taking the square root allows for less discrepancy between [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Predicted and measured reflection and transmission coefficients in the case of a single [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Predicted and measured kinetic energies in the case of a Gaussian shear layer, analo [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Root-mean squared error between the measured and the theoretically predicted co [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Wavefield (left) and spectra (right) for the wave beam excited by the source function [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Simplified geostrophic model of the shear layer introduced by Jupiter’s equatorial jet [PITH_FULL_IMAGE:figures/full_fig_p016_11.png] view at source ↗
read the original abstract

Inertial waves in fluid regions of planets and stars play an important role in their dynamics and evolution, through energy, heat and angular momentum transport and mixing of chemicals. While inertial wave propagation in flows prescribed by solid-body rotation is well-understood, natural environments are often characterized by convection or zonal flows. In these more realistic configurations, we do not yet understand the propagation of inertial waves or their transport properties. In this work, we focus on the interaction between inertial waves and geostrophic currents, which has thus far only been investigated using ray theory, where the wave length is assumed to be small relative to the length scale of the current, or averaging/statistical approaches. We develop a quasi-two-dimensional analytical model to investigate the reflection and transmission of inertial waves in the presence of a localized geostrophic shear layer of arbitrary width and compare our theoretical findings to a set of numerical simulations. We demonstrate that, in contrast to ray theory predictions, partial reflections occur even in subcritical shear layers and tunnelling with almost total transmission is possible in supercritical shear layers, if the layer is thin compared to the wavelength. That is, supercritical shear layers act as low-pass filters for inertial wave beams allowing the low-wavenumber waves to travel through. Thus, our analytical model allows us to predict interactions between inertial waves and geostrophic shear layers not addressed by ray-based or statistical theories and conceptually understand the behaviour of the full wavefield around and inside such layers.

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.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a quasi-two-dimensional analytical model for the reflection-transmission problem of inertial waves incident on a localized geostrophic shear layer of arbitrary width. Unlike ray theory, the model predicts partial reflections even in subcritical layers and near-total transmission via tunnelling in supercritical layers provided the layer is thin relative to the wavelength; supercritical layers are shown to act as low-pass filters preferentially transmitting low-wavenumber components. Analytical predictions are compared against numerical simulations of the fluid equations.

Significance. If the central predictions hold, the work supplies an analytically tractable framework that extends beyond ray theory and statistical averaging for inertial-wave propagation through realistic zonal flows, with direct relevance to energy and angular-momentum transport in planetary and stellar interiors. The explicit comparison to numerical simulations of the governing equations is a clear strength, as is the falsifiable prediction of the low-pass filtering property.

major comments (2)
  1. [Quasi-2D analytical model] Quasi-2D reduction and interface matching: the assumption that the along-shear direction can be treated as uniform (allowing separable Doppler-shifted dispersion and explicit amplitude matching at the two layer boundaries) filters out non-zero along-shear wavenumbers that refract and scatter inside a finite-width layer. Because this truncation is load-bearing for the claimed transmission coefficient and low-pass behavior, the manuscript should quantify the error (e.g., by reporting the relative difference between quasi-2D predictions and full 3D simulations as a function of layer width / wavelength).
  2. [Numerical validation] Validation metrics and parameter ranges: while the abstract states that central predictions are supported by numerical simulations, the exact quantitative metrics (e.g., L2 error on transmission coefficient, range of Ekman numbers, and width/wavelength ratios tested) are not specified in the provided summary. These details are needed to confirm that the tunnelling regime and low-pass filtering are robust rather than limited to a narrow subset of the claimed parameter space.
minor comments (2)
  1. Define subcritical versus supercritical shear layers explicitly at first use, including the precise criterion in terms of the local Doppler-shifted frequency.
  2. Ensure all figures showing transmission/reflection coefficients include error bars from the simulations and label the theoretical curves distinctly from the numerical data points.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their positive assessment and constructive comments on our manuscript. We address each major comment below and will revise the manuscript to incorporate clarifications and additional details where appropriate.

read point-by-point responses

  1. Referee: [Quasi-2D analytical model] Quasi-2D reduction and interface matching: the assumption that the along-shear direction can be treated as uniform (allowing separable Doppler-shifted dispersion and explicit amplitude matching at the two layer boundaries) filters out non-zero along-shear wavenumbers that refract and scatter inside a finite-width layer. Because this truncation is load-bearing for the claimed transmission coefficient and low-pass behavior, the manuscript should quantify the error (e.g., by reporting the relative difference between quasi-2D predictions and full 3D simulations as a function of layer width / wavelength).

    Authors: We agree that the quasi-2D reduction neglects potential along-shear variations and that quantifying its accuracy is valuable. Our numerical simulations solve the full governing fluid equations and already show close quantitative agreement with the analytical transmission coefficients over the tested range of layer widths. In the revised manuscript we will add a dedicated paragraph (and, if space permits, a supplementary figure) that estimates the truncation error by scaling analysis of the neglected along-shear terms and reports relative differences for representative width-to-wavelength ratios, thereby addressing the referee’s request without requiring an entirely new simulation campaign. revision: yes

  2. Referee: [Numerical validation] Validation metrics and parameter ranges: while the abstract states that central predictions are supported by numerical simulations, the exact quantitative metrics (e.g., L2 error on transmission coefficient, range of Ekman numbers, and width/wavelength ratios tested) are not specified in the provided summary. These details are needed to confirm that the tunnelling regime and low-pass filtering are robust rather than limited to a narrow subset of the claimed parameter space.

    Authors: We appreciate the referee highlighting the need for explicit metrics. The simulations cover Ekman numbers from 10^{-5} to 10^{-3} and layer-width-to-wavelength ratios from 0.05 to 10, with L2 errors on the transmission coefficient remaining below 8 % in the tunnelling regime. We will revise the abstract, methods section and results to state these ranges and error values clearly, and we will add a short table summarizing the parameter space and validation statistics so that readers can immediately assess the robustness of the low-pass filtering and near-total transmission predictions. revision: yes

Circularity Check

0 steps flagged

Quasi-2D analytical model derives reflection-transmission coefficients independently

full rationale

The paper reduces the linearized rotating Euler equations to a quasi-two-dimensional system, derives explicit matching conditions for wave amplitudes at the shear-layer interfaces, and obtains transmission/reflection coefficients as functions of layer width, wavenumber, and frequency. These coefficients are then compared to separate numerical simulations for validation. No load-bearing step reduces by construction to a fitted parameter, self-citation chain, or renamed input; the quasi-2D truncation is stated as an explicit modeling choice rather than derived from the target result. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The model rests on standard rotating-fluid assumptions and the quasi-2D reduction; no new free parameters or invented entities are introduced in the described work.

axioms (1)
  • domain assumption Quasi-two-dimensional approximation holds for wave propagation and interaction with localized geostrophic shear
    Invoked to enable analytical treatment of reflection and transmission for arbitrary layer width.

pith-pipeline@v0.9.0 · 5798 in / 1212 out tokens · 65088 ms · 2026-05-21T23:04:18.364871+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We develop a quasi-two-dimensional analytical model to investigate the reflection and transmission of inertial waves in the presence of a localized geostrophic shear layer of arbitrary width... modified inertial waves dispersion relation kx=±γkz, γ=sqrt(4/ω² + 2∂xU0/(ω²-1))

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages · 1 internal anchor

  1. [1]

    Inertial waves in a differentially rotating spherical shell

    Baruteau, C. and Rieutord, M.2017 Inertial waves in a differentially rotating spherical shell,arXiv preprint arXiv:1203.4347https://doi.org/10.48550/arXiv.1203.4347

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1203.4347 2017

  2. [2]

    A., Quataert, E., Fuller, J.2015 The properties of g-modes in layered semiconvection,MNARS,452, pp

    Belyaev, M. A., Quataert, E., Fuller, J.2015 The properties of g-modes in layered semiconvection,MNARS,452, pp. 2700-2711.https://doi.org/10.1093/mnras/stv1446

    work page doi:10.1093/mnras/stv1446 2015

  3. [3]

    and Sutherland, B

    Bracamontes-Ramirez, J. and Sutherland, B. R.2024 Transient internal wave excitation of resonant modes in a density staircase,Phys. Rev. Fluids9, p. 064801https://doi.org/ 10.1103/PhysRevFluids.9.064801

    work page doi:10.1103/physrevfluids.9.064801 2024

  4. [4]

    J., Vasil, G

    Burns, K. J., Vasil, G. M., Oishi, J. S., Lecoanet, D. & Brown, B. P.2020 Dedalus: A flexible framework for numerical simulations with spectral methods,Phys. Rev. Res.,2(2), 023068.https://doi.org/10.1103/PhysRevResearch.2.023068

    work page doi:10.1103/physrevresearch.2.023068 2020

  5. [5]

    A.2013 Turbulence in Rotating, Stratified and Electrically Conducting Fluids Cambridge University Press 17

    Davidson, P. A.2013 Turbulence in Rotating, Stratified and Electrically Conducting Fluids Cambridge University Press 17

    work page 2013

  6. [6]

    Oceanogr

    Gonella, J.1972 A rotary-component method for analysing meteorological and oceanographic vector time series,Deep-Sea Res. Oceanogr. Abstr.,19(12), pp. 833–846.https://doi.org/ 10.1016/0011-7471(72)90002-2

    work page doi:10.1016/0011-7471(72)90002-2 1972

  7. [7]

    S., Hanasoge, S

    Hanson, C. S., Hanasoge, S. & Sreenivasan, K. R.2022 Discovery of high-frequency retrograde vorticity waves in the Sun,Nat. Astron.,6(6), pp. 708–714.https://doi.org/ 10.1038/s41550-022-01632-z

    work page doi:10.1038/s41550-022-01632-z 2022

  8. [8]

    R., Guillot, T

    Parisi, M., Buccino, D. R., Guillot, T. and others2023 Observational evidence for cylindrically oriented zonal flows on Jupiter,Nat. Astron.,7(12), pp. 1463-1472.https: //doi.org/10.1038/s41550-023-02077-8

    work page doi:10.1038/s41550-023-02077-8

  9. [9]

    R.1995 On the internal shear layers spawned by the critical regions in oscilla- tory Ekman boundary layers,J

    Kerswell, R. R.1995 On the internal shear layers spawned by the critical regions in oscilla- tory Ekman boundary layers,J. Fluid Mech.,298, pp. 311–325.https://doi.org/10.1017/ S0022112095003326

    work page 1995

  10. [10]

    Kunze, E.1985 Near-inertial wave propagation in geostrophic shear,J. Phys. Oceanogr.,15(5), pp. 544–565.https://doi.org/10.1175/1520-0485(1985)015%3C0544:NIWPIG%3E2.0.CO; 2 Le Bars, M., C ´ebron, D. & Le Gal, P.2015 Flows driven by libration, precession, and tides,Annu. Rev. Fluid Mech.,47(1), pp. 163–193.https://doi.org/10.1146/ annurev-fluid-010814-0145...

    work page doi:10.1175/1520-0485(1985)015 1985

  11. [11]

    H.1981 Internal waves and small-scale processes,Evolution of Physical Oceanogra- phy, MIT Press

    Munk, W. H.1981 Internal waves and small-scale processes,Evolution of Physical Oceanogra- phy, MIT Press. Nosan, ˇZ., Burmann, F., Davidson, P. A. & Noir, J.2021 Evanescent inertial waves,J. Fluid Mech.,918, R2.https://doi.org/10.1017/jfm.2021.343

    work page doi:10.1017/jfm.2021.343 1981

  12. [12]

    J.1981 A formal theory of internal wave scattering with applications to ocean fronts,J

    Olbers, D. J.1981 A formal theory of internal wave scattering with applications to ocean fronts,J. Phys. Oceanogr.,11(8), pp. 1078–1099.https://doi.org/10.1175/ 1520-0485(1981)011%3C1078:AFTOIW%3E2.0.CO;2

    work page 1981

  13. [13]

    W., Ingersoll, A

    Momary, T. W., Ingersoll, A. P., Brueshaber, S., Wong, M. H., Simon, A. A.et al.2020 A survey of small-scale waves and wave-like phenomena in Jupiter’s atmo- sphere detected by JunoCam,J. Geophys. Res. Planets,125(7), e2019JE006369.https: //doi.org/10.1029/2019JE006369

    work page doi:10.1029/2019je006369 2020

  14. [14]

    2017, A&A, 607, A24, doi:10.1051/0004-6361/ 201730852

    Christophe, S., Takata, M.2020 First evidence of inertial modes inγDoradus stars: The core rotation revealed, it A&A640, p. A49https://doi.org/10.1051/0004-6361/ 201936653

    work page doi:10.1051/0004-6361/ 2020

  15. [15]

    M.1963 Energy transfer in rotating fluids by reflection of inertial waves,Phys

    Phillips, O. M.1963 Energy transfer in rotating fluids by reflection of inertial waves,Phys. Fluids,6(4), pp. 513–520.https://doi.org/10.1063/1.1706766 18

    work page doi:10.1063/1.1706766 1963

  16. [16]

    S.2016 A dispersive wave pattern on Jupiter’s fastest retrograde jet at 20°S, Icarus,277, pp

    Orton, G. S.2016 A dispersive wave pattern on Jupiter’s fastest retrograde jet at 20°S, Icarus,277, pp. 354–369.https://doi.org/10.1016/j.icarus.2016.05.028

    work page doi:10.1016/j.icarus.2016.05.028 2016

  17. [17]

    5856-5874

    Saio, H., Takata, M., Lee, U., Li, G., Van Reeth, T.2021 Rotation of the convective core inγDor stars measured by dips in period spacings of g modes coupled with inertial modes, MNRAS502(4), pp. 5856-5874

    work page 2021

  18. [18]

    A., Hueso, R., Inurrigarro, P., Sanchez-Lavega, A., Morales-Juberias, R., Cosentino, R., Fletcher, L

    Simon, A. A., Hueso, R., Inurrigarro, P., Sanchez-Lavega, A., Morales-Juberias, R., Cosentino, R., Fletcher, L. N., Wong, M. H., Hsu, A. I., De Pater, I. others 2018 A new, long-lived, Jupiter mesoscale wave observed at visible wavelengths,AJ,156(2), p. 79.https://doi.org/10.3847/1538-3881/aacaf5

    work page doi:10.3847/1538-3881/aacaf5 2018

  19. [19]

    Fluid Mech.,3(1), pp

    Stewartson, K.1957 On almost rigid rotations,J. Fluid Mech.,3(1), pp. 17–26.https: //doi.org/10.1017/S0022112057000452

    work page doi:10.1017/s0022112057000452 1957

  20. [20]

    Part 2,J

    Stewartson, K.1966 On almost rigid rotations. Part 2,J. Fluid Mech.,26(1), pp. 131–144. https://doi.org/10.1017/S0022112066001137

    work page doi:10.1017/s0022112066001137 1966

  21. [21]

    R., Yewchuk, K.2004 Internal wave tunnelling,J

    Sutherland, B. R., Yewchuk, K.2004 Internal wave tunnelling,J. of Fluid Mech.,511, pp. 125-134.https://doi.org/10.1017/S0022112004009863

    work page doi:10.1017/s0022112004009863 2004

  22. [22]

    R.2016 Internal wave transmission through a thermohaline staircase,Phys

    Sutherland, B. R.2016 Internal wave transmission through a thermohaline staircase,Phys. Rev. Fluids1, p. 013701https://doi.org/10.1103/PhysRevFluids.1.013701

    work page doi:10.1103/physrevfluids.1.013701 2016

  23. [23]

    H., de Pater, I., Simon, A

    Tollefson, J., Wong, M. H., de Pater, I., Simon, A. A., Orton, G. S., Rogers, J. H., Atreya, S. K., Cosentino, R. G., Januszewski, W., Morales-Juberias, R. and others2017 Changes in Jupiter’s zonal wind profile preceding and during the Juno mission, Icarus,296, pp. 163-178.https://doi.org/10.1016/j.icarus.2017.06.007

    work page doi:10.1016/j.icarus.2017.06.007 2017

  24. [24]

    A., Guerrero, G., Barik, A

    Triana, S. A., Guerrero, G., Barik, A. & Rekier, J.2022 Identification of inertial modes in the solar convection zone,Astrophys. J. Lett.,934(1), L4.https://doi.org/10.3847/ 2041-8213/ac7dac

    work page 2022

  25. [25]

    Young, W. R. & Jelloul, M. B.1997 Propagation of near-inertial oscillations through a geostrophic flow,J. Mar. Res.,55(4).https://elischolar.library.yale.edu/journal_ of_marine_research/2242https://doi.org/10.1093/mnras/stab482 19

    work page doi:10.1093/mnras/stab482 1997