pith. sign in

arxiv: 2605.29554 · v1 · pith:BSEWWUKRnew · submitted 2026-05-28 · 🧮 math.AP · cs.NA· math.NA

Water-at-Rest Equilibrium Stability Analysis of a first-moment Shallow Water Exner Moment Model with Sediment Entrainment and Deposition: Extended Technical Report

Pith reviewed 2026-06-29 06:47 UTC · model grok-4.3

classification 🧮 math.AP cs.NAmath.NA
keywords shallow water equationsExner equationsediment entrainmenthyperbolicityequilibrium manifoldmoment modeltransport closure
0
0 comments X

The pith

The SWEMED1 sediment model reaches a fully-settled water-at-rest equilibrium yet remains only weakly hyperbolic there because of its transport closure.

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

This paper derives the first-moment Shallow Water Exner Moment model with sediment entrainment and deposition, called SWEMED1. It establishes that the complete source term admits a fully-settled water-at-rest equilibrium manifold. The model is shown to be only weakly hyperbolic at this manifold, which rules out Yong's structural stability framework, yet linear spectral analysis and numerical tests reveal no instability. A fast-slow scaling of the source term produces a new suspended water-at-rest equilibrium manifold that is likewise weakly hyperbolic, and the paper traces the remaining obstruction to the transport closure while outlining how new closures might be built.

Core claim

The full source term of the SWEMED1 model possesses a fully-settled water-at-rest equilibrium manifold at which the system is only weakly hyperbolic. Linear spectral analysis and numerical results indicate no instability at this point. Introduction of a fast-slow scaling yields a suspended water-at-rest equilibrium manifold of different structure that is still only weakly hyperbolic, with the persistent obstruction linked directly to the transport closure of the model.

What carries the argument

The transport closure inside the SWEMED1 model, which generates the weak hyperbolicity at both the fully-settled and suspended water-at-rest equilibrium manifolds.

If this is right

  • Linear spectral analysis and numerical experiments confirm the absence of instability at the equilibrium despite weak hyperbolicity.
  • New transport closures can be constructed that preserve the physical entrainment and deposition terms while producing stronger hyperbolicity.
  • The fast-slow scaling supplies a systematic route to alternative equilibrium manifolds with altered structure.
  • The obstruction identified in the transport closure points toward a concrete direction for deriving models with improved analytical properties.

Where Pith is reading between the lines

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

  • Models built with the suggested new closures might permit direct application of structural stability frameworks to prove long-time behavior.
  • Numerical schemes that exactly preserve the identified equilibrium manifold could be designed once the hyperbolicity issue is resolved.
  • The same weak-hyperbolicity pattern may appear in other moment-based sediment models and could be diagnosed by the same fast-slow scaling technique.

Load-bearing premise

The transport closure chosen for the SWEMED1 model is the root cause of the persistent weak hyperbolicity, and alternative closures exist that would remove the obstruction while keeping the physical content of the entrainment and deposition terms unchanged.

What would settle it

Derivation of an alternative transport closure for the SWEMED1 model followed by explicit computation of the eigenvalues of the Jacobian matrix at the water-at-rest equilibrium to check whether they become real and distinct.

read the original abstract

We derive the first-moment Shallow Water Exner Moment model with sediment entrainment and deposition (SWEMED1) and show that the full source term has a fully-settled water-at-rest equilibrium manifold. We prove that the model is only weakly hyperbolic at this equilibrium, which prevents the use of Yong's structural stability framework. However, a linear spectral analysis and numerical results do not indicate instability. Based on numerical results, we introduce a fast-slow scaling of the source term, and for the fast limit, we derive a new suspended water-at-rest equilibrium manifold, which has a different structure but is still only weakly hyperbolic. Our results show that the remaining obstruction is linked to the transport closure of the SWEMED1, and we give a constructive direction for the derivation of new closures leading to models with more desirable analytical properties.

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

0 major / 1 minor

Summary. The manuscript derives the first-moment Shallow Water Exner Moment model with sediment entrainment and deposition (SWEMED1). It establishes that the full source term admits a fully-settled water-at-rest equilibrium manifold, proves that the system is only weakly hyperbolic at this equilibrium (preventing application of Yong's structural stability framework), and shows via linear spectral analysis together with numerical results that no instability is present. A fast-slow scaling of the source term is introduced; in the fast limit a new suspended water-at-rest equilibrium manifold is derived that remains weakly hyperbolic. The obstruction is linked to the transport closure, and a constructive direction for deriving improved closures is indicated.

Significance. If the derivations and numerical evidence hold, the work supplies a concrete stability analysis for a moment closure in a sediment-transport system, isolating the source of weak hyperbolicity while demonstrating that linear stability is retained. The explicit construction of two distinct equilibrium manifolds, the combination of spectral and numerical checks, and the forward-looking suggestion for alternative closures constitute genuine technical contributions to the study of hyperbolic systems with relaxation source terms in mathematical geophysics.

minor comments (1)
  1. The abstract is lengthy; a shorter version that still states the main claims and the constructive direction would improve readability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their detailed summary of the manuscript, for recognizing its technical contributions to the stability analysis of moment closures in sediment-transport systems, and for recommending acceptance.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper derives the SWEMED1 model from first principles, identifies the water-at-rest equilibrium manifold via the source term structure, proves weak hyperbolicity at equilibrium using standard hyperbolic PDE techniques, performs linear spectral analysis, and applies a fast-slow scaling to obtain a new manifold. Each step operates on the explicitly stated model equations and employs external mathematical tools (spectral analysis, scaling limits) without reducing any claimed result to a fitted parameter, self-definition, or unverified self-citation. The final interpretive link to the transport closure is presented as a suggested direction for future closures rather than a completed claim that loops back to the paper's own inputs. The analysis is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis depends on the validity of the first-moment closure used to obtain SWEMED1 and on standard background results from hyperbolic PDE theory; no free parameters or new physical entities are introduced in the abstract.

axioms (2)
  • domain assumption The first-moment closure for the shallow-water Exner system with entrainment and deposition yields a well-defined system of PDEs
    The entire analysis begins from the derived SWEMED1 model obtained via this closure.
  • standard math Standard results on hyperbolicity and spectral stability for systems of conservation laws apply to the linearized equations at equilibrium
    The proof of weak hyperbolicity and the linear spectral analysis rely on these background theorems.

pith-pipeline@v0.9.1-grok · 5692 in / 1572 out tokens · 40858 ms · 2026-06-29T06:47:46.124443+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

29 extracted references · 24 canonical work pages

  1. [1]

    Hydrological sciences journal 54(1), 147–159 (2009)

    Afzalimehr, H., Rennie, C.: Determination of bed shear stress in gravel-bed rivers using boundary- layer parameters. Hydrological sciences journal 54(1), 147–159 (2009). DOI 10.1623/hysj.54.1. 147

  2. [2]

    SIAM Journal on Scientific Com- puting 25(6), 2050–2065 (2004)

    Audusse, E., Bouchut, F., Bristeau, M.O., Klein, R., Perthame, B.: A fast and stable well-balanced scheme with hydrostatic reconstruction for shallow water flows. SIAM Journal on Scientific Com- puting 25(6), 2050–2065 (2004). DOI 10.1137/S1064827503431090

  3. [3]

    derivation and numerical validation

    Audusse, E., Bristeau, M.O., Perthame, B., Sainte-Marie, J.: A multilayer saint-venant system with mass exchanges for shallow water flows. derivation and numerical validation. ESAIM: Mathemat- ical Modelling and Numerical Analysis 45(1), 169–200 (2011). DOI 10.1051/m2an/2010036

  4. [4]

    Journal of Computational Physics 364, 209–234 (2018)

    Bonaventura, L., Fernández-Nieto, E.D., Garres-Díaz, J., Narbona-Reina, G.: Multilayer shallow water models with locally variable number of layers and semi-implicit time discretization. Journal of Computational Physics 364, 209–234 (2018). DOI 10.1016/j.jcp.2018.03.017

  5. [5]

    i: Formulation and numeri- cal analysis

    Bradford, S.F., Katopodes, N.D.: Hydrodynamics of turbid underflows. i: Formulation and numeri- cal analysis. Journal of Hydraulic Engineering 125(10), 1006–1015 (1999). DOI 10.1061/(ASCE) 0733-9429(1999)125:10(1006)

  6. [6]

    Journal of Scientific Computing 85, 1–40 (2020)

    Bürger, R., Fernández-Nieto, E.D., Osores, V .: A multilayer shallow water approach for polydis- perse sedimentation with sediment compressibility and mixture viscosity. Journal of Scientific Computing 85, 1–40 (2020). DOI 10.1007/s10915-020-01334-6 24 A. Parvin et al

  7. [8]

    Advances in Water Resources 132, 103387 (2019)

    Cantero-Chinchilla, F.N., Castro-Orgaz, O., Khan, A.A.: Vertically-averaged and moment equa- tions for flow and sediment transport. Advances in Water Resources 132, 103387 (2019). DOI 10.1016/j.advwatres.2019.103387

  8. [9]

    Numerical Mathematics: Theory, Methods and Applications 16(4), 1087–1126 (2023)

    Del Grosso, A., Castro Díaz, M.J., Chalons, C., Morales de Luna, T.: On lagrange-projection schemes for shallow water flows over movable bottom with suspended and bedload transport. Numerical Mathematics: Theory, Methods and Applications 16(4), 1087–1126 (2023). DOI 10.4208/nmtma.OA-2023-0082

  9. [11]

    Water 13(9), 1254 (2021)

    Elgamal, M.: A moment-based Chezy formula for bed shear stress in varied flow. Water 13(9), 1254 (2021). DOI 10.3390/w13091254

  10. [12]

    Sitzungsberichte der Akademie der Wissenschaften in Wien, Mathematisch- Naturwissenschaftliche Klasse 134(2a), 165–203 (1925)

    Exner Ewarten, F.M.: Über die Wechselwirkung zwischen Wasser und Geschiebe in Flüssen. Sitzungsberichte der Akademie der Wissenschaften in Wien, Mathematisch- Naturwissenschaftliche Klasse 134(2a), 165–203 (1925)

  11. [13]

    Journal of Statistical Physics 162(2), 457–486 (2016)

    Fan, Y ., Koellermeier, J., Li, J., Li, R., Torrilhon, M.: Model reduction of kinetic equations by operator projection. Journal of Statistical Physics 162(2), 457–486 (2016). DOI 10.1007/ s10955-015-1384-9

  12. [14]

    Journal of Geophysical Research: Oceans 98(C3), 4793–4807 (1993)

    Garcia, M.H., Parker, G.: Experiments on the entrainment of sediment into suspension by a dense bottom current. Journal of Geophysical Research: Oceans 98(C3), 4793–4807 (1993). DOI 10.1029/92JC02404

  13. [15]

    Communications in Computational Physics 30(3), 903– 941 (2021)

    Garres-Díaz, J., Castro Díaz, M.J., Koellermeier, J., Morales de Luna, T.: Shallow water moment models for bedload transport problems. Communications in Computational Physics 30(3), 903– 941 (2021). DOI 10.4208/cicp.OA-2020-0152

  14. [16]

    Communications on Applied Mathematics and Computation 5(2), 885–922 (2023)

    González-Aguirre, J.C., González-Vázquez, J.A., Alavez-Ramírez, J., Silva, R., Vázquez-Cendón, M.E.: Numerical simulation of bed load and suspended load sediment transport using well- balanced numerical schemes. Communications on Applied Mathematics and Computation 5(2), 885–922 (2023). DOI 10.1007/s42967-021-00162-1

  15. [18]

    Mathematical Methods in the Applied Sciences 45(10), 6459–6480 (2022)

    Huang, Q., Koellermeier, J., Y ong, W.A.: Equilibrium stability analysis of hyperbolic shallow wa- ter moment equations. Mathematical Methods in the Applied Sciences 45(10), 6459–6480 (2022). DOI 10.1002/mma.8180

  16. [19]

    Communications in Computational Physics 28(3), 1038–1084 (2020)

    Koellermeier, J., Rominger, M.: Analysis and numerical simulation of hyperbolic shallow water moment equations. Communications in Computational Physics 28(3), 1038–1084 (2020). DOI 10.4208/cicp.OA-2019-0065

  17. [20]

    Com- munications in Computational Physics 25(3), 669–702 (2019)

    Kowalski, J., Torrilhon, M.: Moment approximations and model cascades for shallow flow. Com- munications in Computational Physics 25(3), 669–702 (2019). DOI 10.4208/cicp.OA-2017-0263

  18. [21]

    Communications in Computational Physics 29(1), 80–110 (2020)

    Meng, X., Hoang, T.T.P ., Wang, Z., Ju, L.: Localized exponential time differencing method for shallow water equations: Algorithms and numerical study. Communications in Computational Physics 29(1), 80–110 (2020). DOI 10.4208/cicp.OA-2019-0214

  19. [22]

    Jour- nal of Differential Equations 155(1), 89–132 (1999)

    Y ong, W.A.: Singular perturbations of first-order hyperbolic systems with stiff source terms. Jour- nal of Differential Equations 155(1), 89–132 (1999). DOI 10.1006/jdeq.1998.3584

  20. [23]

    China Water and Power Press, Beijing (1993)

    Zhang, R., Xie, J.: Sedimentation research in China: Systematic selections. China Water and Power Press, Beijing (1993)

  21. [24]

    Journal of Hy- drology 570, 647–665 (2019)

    Zhao, J., Özgen-Xian, I., Liang, D., Wang, T., Hinkelmann, R.: A depth-averaged non-cohesive sediment transport model with improved discretization of flux and source terms. Journal of Hy- drology 570, 647–665 (2019). DOI 10.1016/j.jhydrol.2018.12.059

  22. [25]

    Communications in Mathematical Sciences 15(3), 609–633 (2017)

    Zhao, W., Y ong, W.A., Luo, L.S.: Stability analysis of a class of globally hyperbolic moment system. Communications in Mathematical Sciences 15(3), 609–633 (2017). DOI 10.4310/CMS. 2017.v15.n3.a2 Water-at-Rest Equilibrium Stability of SWEMED1 25

  23. [26]

    Communications in Mathematical Sciences 11(2), 547–571 (2013)

    Cai, Z., Fan, Y ., Li, R.: Globally hyperbolic regularization of Grad’s moment system in one dimensional space. Communications in Mathematical Sciences 11(2), 547–571 (2013). DOI 10.4310/CMS.2013.v11.n2.a12

  24. [27]

    Numerical Mathematics: Theory, Methods and Applications 10(2), 255–277 (2017)

    Di, Y ., Fan, Y ., Li, R., Zheng, L.: Linear stability of hyperbolic moment models for Boltzmann equation. Numerical Mathematics: Theory, Methods and Applications 10(2), 255–277 (2017). DOI 10.4208/nmtma.2017.s04

  25. [28]

    Commu- nications on Applied Mathematics and Computation 6(4), 2155–2195 (2024)

    Scholz, U., Kowalski, J., Torrilhon, M.: Dispersion in shallow moment equations. Commu- nications on Applied Mathematics and Computation 6(4), 2155–2195 (2024). DOI 10.1007/ s42967-023-00325-2

  26. [29]

    Advances in Water Resources 34(8), 980–989 (2011)

    Cordier, S., Le, M.H., Morales de Luna, T.: Bedload transport in shallow water models: Why splitting (may) fail, how hyperbolicity (can) help. Advances in Water Resources 34(8), 980–989 (2011). DOI 10.1016/j.advwatres.2011.05.002

  27. [30]

    ESAIM: Mathematical Modelling and Numerical Analysis 51(1), 115–145 (2017)

    Fernandez-Nieto, E.D., Morales de Luna, T., Narbona-Reina, G., Zabsonre, J.D.: Formal deduc- tion of the Saint-Venant–Exner model including arbitrarily sloping sediment beds and associated energy. ESAIM: Mathematical Modelling and Numerical Analysis 51(1), 115–145 (2017). DOI 10.1051/m2an/2016018

  28. [31]

    Advances in Water Resources140, 103575 (2020)

    Gonzalez-Aguirre, J.C., Castro, M.J., Morales de Luna, T.: A robust model for rapidly varying flows over movable bottom with suspended and bedload transport: Modelling and numerical ap- proach. Advances in Water Resources140, 103575 (2020). DOI 10.1016/j.advwatres.2020.103575

  29. [32]

    Journal of Nuclear Medicine 20(2), 162–164 (1979)

    Prince, J.R.: Comments on equilibrium, transient equilibrium, and secular equilibrium in serial radioactive decay. Journal of Nuclear Medicine 20(2), 162–164 (1979). Appendices A Derivation of the complete reference system A.1 Mapping of the mass balance From the divergence-free conditions ∇ · u = 0, we have 𝜕𝑥𝑢 + 𝜕𝑧 𝑤 = 0. (82) After multiplying (82) by ...