pith. sign in

arxiv: 2606.21095 · v2 · pith:LK4AVKWKnew · submitted 2026-06-19 · ❄️ cond-mat.stat-mech · cond-mat.other· math.PR· physics.app-ph· physics.class-ph· physics.flu-dyn

Asymptotic hydrographs and anomalous dispersion in mass-conserving storage cascades

Henrique Santos Lima , M\'ark Honti , Bal\'azs S\'andor This is my paper

Pith reviewed 2026-07-01 07:28 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech cond-mat.othermath.PRphysics.app-phphysics.class-phphysics.flu-dyn
keywords q-exponential densityLévy stable lawsmass-conserving cascadesanomalous dispersionhydrographsstorage reservoirsnon-Gaussian spreading
0
0 comments X

The pith

Cascades of reservoirs with q-exponential waiting times produce asymptotic Lévy-stable hydrographs for 1<q<2.

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

The paper replaces the usual exponential waiting-time density with the q-exponential density from nonextensive statistics inside a chain of mass-conserving reservoirs. It derives an analytic large-cascade expression for the outflow rate and shows that this expression converges to a stable Lévy density, shifted by a Galilean transformation when q equals 5/3. A reader cares because the construction relies only on repeated convolution and strict mass conservation yet reproduces the non-Gaussian spreading normally obtained from fractional derivatives. The result supplies a direct probabilistic mechanism that turns nonextensive statistics into macroscopic anomalous dispersion.

Core claim

The main result is an analytical asymptotic expression for the outflow of a mass-conserving cascade driven by a q-exponential waiting-time kernel. In the critical case q=5/3 the large-cascade flow rate converges to a stable Lévy density whose time argument is shifted by a Galilean-type transformation. For the entire regime 1<q<2 the macroscopic dynamics are governed by α-stable Lévy laws. This establishes that anomalous non-Gaussian dispersion can emerge from pure mass-conserving convolutional chains without fractional derivatives.

What carries the argument

The q-exponential waiting-time kernel iterated through mass-conserving convolutions in a reservoir cascade.

If this is right

  • The outflow of a sufficiently long cascade is described by an α-stable Lévy density whose stability index depends on q.
  • For q=5/3 the Lévy density appears with an explicit Galilean shift in its argument.
  • Macroscopic dynamics throughout 1<q<2 follow stable Lévy laws rather than Gaussian diffusion.
  • Anomalous dispersion arises solely from the choice of waiting-time kernel and mass conservation.

Where Pith is reading between the lines

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

  • If the q-exponential kernel is realizable in physical storage systems, the same mechanism could generate Lévy hydrographs in hydrology or chemical engineering without invoking fractional operators.
  • Finite-cascade simulations could be compared directly to the analytic asymptote to measure the rate of convergence.
  • The Galilean shift at q=5/3 may produce observable time offsets in measured outflow peaks that distinguish this model from ordinary stable processes.

Load-bearing premise

The q-exponential density functions as a proper waiting-time kernel whose successive convolutions in the mass-conserving cascade admit the large-cascade asymptotic analysis presented.

What would settle it

A direct numerical convolution of many q-exponential densities followed by checking whether the resulting density matches the predicted shifted Lévy form within sampling error.

Figures

Figures reproduced from arXiv: 2606.21095 by Bal\'azs S\'andor, Henrique Santos Lima, M\'ark Honti.

Figure 1
Figure 1. Figure 1: FIG. 1: Log-log comparison of the exact [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Exact numerical convolution of the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
read the original abstract

Sums of independent exponential random variables lead to the Erlang distribution, providing a direct probabilistic route from exponential waiting times to the integer-shape gamma law. This paper investigates how this classical construction changes when the exponential waiting-time density is replaced by the $q$-exponential density of nonextensive statistics. Our main result is an analytical asymptotic expression for the outflow of a mass-conserving cascade of reservoirs driven by a $q$-exponential waiting-time kernel. In the critical case $q=5/3$, the large-cascade flow rate converges to a stable L\'{e}vy density whose time argument is shifted by a Galilean-type transformation. This shifted L\'{e}vy law gives the asymptotic hydrograph of the cascade. We also found that for the entire regime $1<q<2$ the macroscopic dynamics are governed by $\alpha$-stable L\'{e}vy laws. This proves that anomalous non-Gaussian dispersion can emerge from pure mass-conserving convolutional chains without invoking fractional derivatives.

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

1 major / 1 minor

Summary. The paper replaces the exponential waiting-time density in a classical Erlang construction with the q-exponential density from nonextensive statistics. It derives an analytical asymptotic expression for the outflow hydrograph of a mass-conserving cascade of reservoirs under n-fold convolutions of this kernel. The central claim is that for 1 < q < 2 the large-n limit is governed by α-stable Lévy laws, with the specific case q = 5/3 yielding a Galilean-shifted Lévy density; this is presented as a demonstration that anomalous non-Gaussian dispersion arises from pure mass-conserving convolutional chains without fractional derivatives.

Significance. If the asymptotic analysis is rigorously established, the result supplies a concrete probabilistic route from q-exponentials to stable laws in a mass-conserving setting. It strengthens the link between nonextensive statistics and anomalous transport in storage cascades, offering an alternative to fractional-derivative models while preserving the normalization that guarantees mass conservation. The explicit connection to the domain of attraction of power-law tails (index 1/(q-1)) is a notable strength.

major comments (1)
  1. [Abstract] The abstract states that the large-cascade flow rate converges to a stable Lévy density with a Galilean-type shift at q = 5/3, but the explicit form of the shift and the steps of the asymptotic analysis (e.g., via characteristic functions or Laplace transforms of the n-fold convolution) are not supplied. Without these, it is impossible to verify that the claimed limit is indeed the stated shifted α-stable law rather than a different member of the domain of attraction.
minor comments (1)
  1. [Abstract] The normalization of the q-exponential kernel (ensuring integration to unity) should be stated explicitly, together with the support on which it is defined, to make the mass-conserving property immediate.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive evaluation of the work's significance and for the detailed comment. We address the single major comment below and will incorporate the suggested clarification.

read point-by-point responses

  1. Referee: [Abstract] The abstract states that the large-cascade flow rate converges to a stable Lévy density with a Galilean-type shift at q = 5/3, but the explicit form of the shift and the steps of the asymptotic analysis (e.g., via characteristic functions or Laplace transforms of the n-fold convolution) are not supplied. Without these, it is impossible to verify that the claimed limit is indeed the stated shifted α-stable law rather than a different member of the domain of attraction.

    Authors: We agree that the abstract would be strengthened by a concise indication of the analytical route. The manuscript derives the result via the characteristic function of the q-exponential kernel, shows that the n-fold convolution belongs to the domain of attraction of an α-stable law with α = 1/(q−1), and obtains the Galilean shift explicitly at q = 5/3 by centering the stable density. These steps appear in Sections 3 and 4. In the revised version we will expand the abstract by one or two sentences to state the characteristic-function approach and the explicit form of the shift, thereby making the verification path visible at the abstract level while preserving its brevity. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper derives the large-n asymptotic form of n-fold convolutions of the normalized q-exponential kernel directly from its power-law tail index 1/(q-1) and the domain-of-attraction properties of stable laws. This is a standard application of the generalized central limit theorem to a mass-conserving kernel whose integral equals unity by construction; the resulting shifted Lévy density for q=5/3 and α-stable laws for 1<q<2 follow as mathematical consequences without any parameter fitting to the target hydrograph, without self-citations invoked as load-bearing uniqueness theorems, and without renaming an external empirical pattern. The central claim that anomalous dispersion arises from convolutional chains is therefore self-contained in the convolution analysis itself.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the substitution of q-exponential densities into the classical Erlang construction and on the existence of large-cascade asymptotics for the resulting convolutional chain; no free parameters beyond the q value itself are identified in the abstract.

free parameters (1)
  • q
    The nonextensivity parameter q controls the waiting-time kernel and takes the specific critical value 5/3.
axioms (1)
  • domain assumption The outflow of the cascade is obtained by successive convolution of the q-exponential waiting-time densities under mass conservation.
    This modeling choice is the direct replacement for the classical exponential case described in the abstract.

pith-pipeline@v0.9.1-grok · 5734 in / 1176 out tokens · 40663 ms · 2026-07-01T07:28:14.176986+00:00 · methodology

Review history (2 revisions) →

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

39 extracted references · 39 canonical work pages · 3 internal anchors

  1. [1]

    Atx= 0: ˜O(0, s) = exp(0) = 1.(39)

    work page

  2. [2]

    Thus: ˜O(∞, s) = lim x→∞ exp − √s√D0 − s v x = 0.(40) Substituting Eqs

    Asx→ ∞: for any complex frequencysin the right half-plane where Re( √s)>0, the spatial decay of the anomalous diffusion term dominates the advective lin- ear growth. Thus: ˜O(∞, s) = lim x→∞ exp − √s√D0 − s v x = 0.(40) Substituting Eqs. (39) and (40) into the mass relation yields: ˜S(s) =− 1 s[0−1] = 1 s .(41) Applying the inverse Laplace transform, we f...

    work page 2024

  3. [3]

    J. E. Nash,The form of the instantaneous unit hydrograph, Proceedings of the General Assembly of Toronto, Interna- tional Association of Scientific Hydrology3, 114–121 (1957)

    work page 1957

  4. [4]

    J. E. Nash,Systematic determination of unit hydrograph parameters, Journal of Geophysical Research64, 111–115 (1959)

    work page 1959

  5. [5]

    J. C. I. Dooge,A general theory of the unit hydrograph, Jour- nal of Geophysical Research64, 241–256 (1959)

    work page 1959

  6. [6]

    Rodr´ ıguez-Iturbe and J

    I. Rodr´ ıguez-Iturbe and J. B. Vald´ es,The geomorphologic structure of hydrologic response, Water Resources Research 15, 1409–1420 (1979)

    work page 1979

  7. [7]

    Sz¨ oll˝ osi-Nagy,The discretization of the continuous linear cascade by means of state space analysis, Journal of Hydrol- ogy58, 223–236 (1982)

    A. Sz¨ oll˝ osi-Nagy,The discretization of the continuous linear cascade by means of state space analysis, Journal of Hydrol- ogy58, 223–236 (1982)

    work page 1982

  8. [8]

    Szilagyi,Discrete state-space approximation of the contin- uous Kalinin–Milyukov–Nash cascade of noninteger storage elements, Journal of Hydrology328, 132–140 (2006)

    J. Szilagyi,Discrete state-space approximation of the contin- uous Kalinin–Milyukov–Nash cascade of noninteger storage elements, Journal of Hydrology328, 132–140 (2006)

    work page 2006

  9. [9]

    V. P. Singh,Kinematic Wave Modeling in Water Resources: Environmental Hydrology(John Wiley & Sons, 1997)

    work page 1997

  10. [10]

    V. P. Singh and D. A. Woolhiser,A nonlinear kinematic wave model for watershed surface runoff, Journal of Hydrology31, 221–243 (1976)

    work page 1976

  11. [11]

    V. P. Singh and D. A. Woolhiser,Mathematical modeling of watershed hydrology, Journal of Hydrologic Engineering7, 270–292 (2002)

    work page 2002

  12. [12]

    Tsallis,Possible generalization of Boltzmann–Gibbs statis- tics, Journal of Statistical Physics52, 479–487 (1988)

    C. Tsallis,Possible generalization of Boltzmann–Gibbs statis- tics, Journal of Statistical Physics52, 479–487 (1988)

    work page 1988

  13. [13]

    Tsallis,Introduction to Nonextensive Statistical Mechan- ics: Approaching a Complex World(Springer, New York, 2023)

    C. Tsallis,Introduction to Nonextensive Statistical Mechan- ics: Approaching a Complex World(Springer, New York, 2023)

    work page 2023

  14. [14]

    E. P. Borges,A possible deformed algebra and calculus in- spired in nonextensive thermostatistics, Physica A: Statistical Mechanics and its Applications340, 95–101 (2004)

    work page 2004

  15. [15]

    Tsallis, S

    C. Tsallis, S. V. F. Levy, A. M. C. Souza, and R. Maynard, Statistical-mechanical foundation of the ubiquity of L´ evy dis- tributions in nature, Physical Review Letters75, 3589–3593 (1995)

    work page 1995

  16. [16]

    Umarov, C

    S. Umarov, C. Tsallis, and S. Steinberg,On a q-central limit theorem consistent with nonextensive statistical mechanics, 7 Milan Journal of Mathematics76, 307–328 (2008)

    work page 2008

  17. [17]

    Lutz,Anomalous diffusion and Tsallis statistics in an op- tical lattice, Physical Review A67, 051402 (2003)

    E. Lutz,Anomalous diffusion and Tsallis statistics in an op- tical lattice, Physical Review A67, 051402 (2003)

    work page 2003

  18. [18]

    H. S. Lima and E. M. F. Curado,Emergent heavy-tailed dis- tributions from a Markovian random walk, arXiv preprint arXiv:2605.22933 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  19. [19]

    S. M. Soares, L. Squillante, H. S. Lima, C. Tsallis, and M. de Souza,Universally nondiverging Gr¨ uneisen parameter at critical points, Physical Review B111, L060409 (2025)

    work page 2025

  20. [20]

    H. S. Lima and C. Tsallis,Ising chain: Thermal conduc- tivity and first-principle validation of Fourier’s law, Physica A: Statistical Mechanics and its Applications628, 129161 (2023)

    work page 2023

  21. [21]

    H. S. Lima, C. Tsallis, D. Eroglu, and U. Tirnakli,Fourier’s law breakdown for the planar-rotor chain with long-range in- teractions, Physica A: Statistical Mechanics and its Applica- tions685, 131314 (2026)

    work page 2026

  22. [22]

    H. S. Lima, R. M. Pereira, L. Moriconi, K. R. Sreenivasan, and C. Tsallis,Superstatistical Approach to Turbulent Circu- lation Fluctuations, arXiv preprint arXiv:2604.15277 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  23. [23]

    Nazari, D.-J

    B. Nazari, D.-J. Seo, and S. Kim,Analytical solution for non- linear hydrologic routing with general power-law storage func- tion, Journal of Hydrologic Engineering26, 04021016 (2021)

    work page 2021

  24. [24]

    A spectral model of power-law decay in natural and engineered systems

    B. S´ andor, M. Honti, and H. S. Lima,A spectral model of power-law decay in natural and engineered systems, arXiv preprint arXiv:2606.08342 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  25. [25]

    A. R. Plastino and A. Plastino,Non-extensive statistical me- chanics and generalized Fokker-Planck equation, Physica A 222, 347–354 (1995)

    work page 1995

  26. [26]

    Borland,Microscopic dynamics of the nonlinear Fokker- Planck equation: A phenomenological model, Physical Re- view E57, 6634–6642 (1998)

    L. Borland,Microscopic dynamics of the nonlinear Fokker- Planck equation: A phenomenological model, Physical Re- view E57, 6634–6642 (1998)

    work page 1998

  27. [27]

    Bologna, C

    M. Bologna, C. Tsallis, and P. Grigolini,Anomalous diffu- sion associated with nonlinear fractional derivative Fokker- Planck-like equation: Exact time-dependent solutions, Physi- cal Review E62, 2213–2218 (2000)

    work page 2000

  28. [28]

    F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark,NIST Handbook of Mathematical Functions(Cam- bridge University Press, 2010)

    work page 2010

  29. [29]

    Landau,On the energy loss of fast particles by ionization, J

    L. Landau,On the energy loss of fast particles by ionization, J. Phys. USSR8, 201–205 (1944)

    work page 1944

  30. [30]

    J. M. Chambers, C. L. Mallows, and B. W. Stuck,A method for simulating stable random variables, J. Amer. Statist. As- soc., vol.71, 340–344 (1976)

    work page 1976

  31. [31]

    J. P. Nolan,Numerical calculation of stable densities and distribution functions,Comm. Statist. Stochastic Models, vol.13, 759–774 (1997)

    work page 1997

  32. [32]

    Gorska and K

    K. Gorska and K. A. Penson,L´ evy stable two-sided distribu- tions: exact and explicit densities for asymmetric case, Phys. Rev. E, vol.83, 061125 (2011)

    work page 2011

  33. [33]

    Feller,An Introduction to Probability Theory and Its Ap- plications, Volume II(John Wiley & Sons, 1971)

    W. Feller,An Introduction to Probability Theory and Its Ap- plications, Volume II(John Wiley & Sons, 1971)

    work page 1971

  34. [34]

    V. M. Zolotarev,One-Dimensional Stable Distributions (American Mathematical Society, 1986)

    work page 1986

  35. [35]

    S. M. Duarte Queir´ os, L. G. Moyano, J. de Souza, and C. Tsallis,A nonextensive approach to the dynamics of financial observables, The European Physical Journal B55, 161–167 (2007)

    work page 2007

  36. [36]

    Beck,Dynamical foundations of nonextensive statistical mechanics, Physical Review Letters87, 180601 (2001)

    C. Beck,Dynamical foundations of nonextensive statistical mechanics, Physical Review Letters87, 180601 (2001)

    work page 2001

  37. [37]

    Beck and E

    C. Beck and E. G. D. Cohen,Superstatistics, Physica A: Sta- tistical Mechanics and its Applications322, 267–275 (2003)

    work page 2003

  38. [38]

    Sattin,Superstatistics from a different perspective, Physica A: Statistical Mechanics and its Applications338, 437–444 (2004)

    F. Sattin,Superstatistics from a different perspective, Physica A: Statistical Mechanics and its Applications338, 437–444 (2004)

    work page 2004

  39. [39]

    S. M. Duarte Queir´ os,On non-Gaussianity and dependence in financial time series: a nonextensive approach, Quantita- tive Finance5, 475–487 (2005)

    work page 2005