pith. sign in

arxiv: 2605.18785 · v1 · pith:6UG5VMCMnew · submitted 2026-05-07 · ⚛️ physics.soc-ph · cs.SY· eess.SY

A Novel Urban Flood Dynamical System Model and a Corresponding Nonstandard Finite Difference Method

Yongfu Tian , Shan Ding , Guofeng Su , Jianguo Chen This is my paper

Pith reviewed 2026-05-20 23:44 UTC · model grok-4.3

classification ⚛️ physics.soc-ph cs.SYeess.SY
keywords urban flooddynamical system modelnonstandard finite differencecellular automatawater conservationpositivity preservationHEC-RAS comparison
0
0 comments X

The pith

An urban flood dynamical system model paired with a conservation nonstandard finite difference method keeps solutions positive and conserves water while matching physical conditions and HEC-RAS results within millimeters.

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

The paper constructs an Urban Flood Dynamical System Model that treats the city as a cellular grid where cell types and simple water motion rules can be chosen to fit real urban data. A sufficient condition on those rules guarantees that the model's solutions respect the broad physical laws of water movement. To evolve the system numerically without breaking positivity or conservation, the authors introduce a first-order conservation nonstandard finite difference scheme that shares the same fixed points as the continuous model. Validation against an exact solution confirms the scheme's properties, and a full urban flood experiment produces distances of roughly 2 mm and 0.02 mm from HEC-RAS outputs with relative errors of 7.5 percent and 0.06 percent.

Core claim

The UFDSM lets users define cell types and water motion or distribution rules tailored to actual city layouts, subject to a sufficient condition that aligns discrete solutions with macroscopic water physics. The accompanying first-order conservation nonstandard finite difference algorithm preserves positivity of water depths, exact conservation of total water volume, and the fixed-point set of the underlying dynamical system, properties verified both analytically and against HEC-RAS urban simulations that differ by only a few millimeters.

What carries the argument

The first-order conservation nonstandard finite difference algorithm applied to the UFDSM, which enforces positivity and water conservation while inheriting the dynamical system's fixed points.

If this is right

  • The model can be coupled with additional hydrological processes such as infiltration or drainage networks without altering its core structure.
  • Data assimilation becomes straightforward because the scheme maintains the same equilibrium states as the continuous model.
  • Urban planners can trade off rule simplicity against computational cost while still satisfying the physical consistency condition.
  • The same nonstandard discretization approach extends directly to other cellular-automaton-based environmental models that require positivity and conservation.

Where Pith is reading between the lines

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

  • Coupling the UFDSM with real-time rainfall or sensor networks could enable operational flood forecasting at lower computational cost than full hydrodynamic codes.
  • The positivity-preserving property may reduce the need for artificial clipping or limiters when the model is embedded inside larger climate or infrastructure simulations.
  • Testing the sufficient condition on rules derived from high-resolution lidar or satellite imagery would show how much urban detail can be retained before accuracy degrades.

Load-bearing premise

The chosen simple water motion and distribution rules satisfy the paper's sufficient condition so the dynamical system solutions remain consistent with the macroscopic physical laws of water movement.

What would settle it

Run the nonstandard finite difference scheme on a closed test domain with an exact analytical solution; if total water volume changes by more than machine precision or any cell depth becomes negative, the conservation and positivity claims fail.

read the original abstract

Urban flood disaster is one of the most serious natural disasters. Numerous flood simulation models have been proposed and relatively matured. However, two major challenges persist: excessive simplification of the city system and high computational complexity. To break these limitations, this paper develops an Urban Flood Dynamical System Model (UFDSM) based on the concept of the Cellular Automata Urban Flood Model. This model allows flexible customization of cell types and selection of water motion or distribution rules based on actual urban environments to incorporate as much the urban system data as possible. The water motion and distribution rules can be simple, which could reduce the computational complexity, but not arbitrary. So, a sufficient condition is provided so that solutions of dynamical system align with macroscopic physical conditions governing water movement. Then, to preserve the evolutionary properties of the UFDSM, we propose a first-order conservation nonstandard finite difference algorithm. This numerical method ensures positive solutions and conservation of water while maintaining the same fixed-point characteristics as the dynamical system. And, this numerical method is validated by comparing it with an analytical solution.Furthermore, to verify the applicability of our model, we performed an urban flood simulation experiment and compared it to HEC-RAS. There is approximately a 2mm discrepancy in distance dp' and 0.02mm discrepancy in distance d2' , with the relative distance Rp about 7.5% and the relative distance R2 approximately 0.06%. Additionally, the proposed model is easily coupled with other hydrological processes and facilitates data assimilation, thereby offering promising practical applications.

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. This manuscript presents the Urban Flood Dynamical System Model (UFDSM), a cellular automata-based approach for simulating urban floods that permits customization of cell types and selection of water motion and distribution rules. A sufficient condition is derived to ensure that the dynamical system's solutions are consistent with macroscopic physical conditions for water movement. To numerically solve the system while preserving key properties, a first-order conservation nonstandard finite difference (NSFD) algorithm is proposed, which maintains positivity, water conservation, and the same fixed points as the continuous system. The NSFD scheme is validated against an analytical solution, and the overall model is tested in an urban flood scenario with quantitative comparison to the HEC-RAS model, reporting discrepancies of approximately 2 mm in dp' and 0.02 mm in d2', with relative distances of about 7.5% and 0.06%, respectively.

Significance. If the central claims hold, the work offers a flexible and computationally lighter alternative to established flood models by allowing incorporation of urban data through customizable rules while providing NSFD-based guarantees on positivity and conservation. The small reported discrepancies versus HEC-RAS and the potential for coupling with other hydrological processes and data assimilation represent practical strengths. Explicit credit is due for the analytical validation of the NSFD scheme and the reproducible comparison metrics.

major comments (1)
  1. [Urban flood simulation experiment] The manuscript states that a sufficient condition is provided so that solutions align with macroscopic physical conditions, yet the urban flood simulation experiment section does not confirm or demonstrate that the specific water motion and distribution rules chosen for that experiment satisfy the condition. This verification is load-bearing for the physical-consistency interpretation of the HEC-RAS comparison results (dp', d2', Rp, R2).
minor comments (1)
  1. The abstract and validation paragraphs would benefit from explicit cross-references to the section or equation defining the sufficient condition and to the precise rules employed in the urban experiment.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the thorough review and valuable comments on our manuscript. The major comment highlights an important aspect regarding the verification of the sufficient condition in the urban flood simulation experiment. We address this below and will incorporate the necessary clarification in the revised manuscript.

read point-by-point responses

  1. Referee: [Urban flood simulation experiment] The manuscript states that a sufficient condition is provided so that solutions align with macroscopic physical conditions, yet the urban flood simulation experiment section does not confirm or demonstrate that the specific water motion and distribution rules chosen for that experiment satisfy the condition. This verification is load-bearing for the physical-consistency interpretation of the HEC-RAS comparison results (dp', d2', Rp, R2).

    Authors: We appreciate the referee bringing this to our attention. While the sufficient condition is derived in the theoretical section to ensure that the dynamical system's solutions align with macroscopic physical conditions for water movement, we recognize that the experiment section would benefit from an explicit check that the chosen rules satisfy this condition. This will reinforce the validity of interpreting the HEC-RAS comparison results as physically consistent. In the revised manuscript, we will include a verification step for the specific water motion and distribution rules used in the urban flood simulation, demonstrating compliance with the sufficient condition. This addition will be placed in the relevant section to directly support the reported discrepancies of approximately 2 mm in dp' and 0.02 mm in d2', with relative distances of about 7.5% and 0.06%. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation remains independent with external validation

full rationale

The paper constructs the UFDSM by extending Cellular Automata concepts with customizable cell types and water rules, then supplies an explicit sufficient condition for alignment with macroscopic physics. The NSFD scheme is derived to enforce positivity, mass conservation, and fixed-point preservation as design goals, not as post-hoc fits. Validation proceeds via direct comparison to an independent analytical solution and to HEC-RAS outputs, yielding reported metric discrepancies. No quoted step reduces a claimed prediction or uniqueness result to a self-definition, a fitted parameter renamed as output, or a load-bearing self-citation chain. The central claims therefore rest on stated assumptions and external benchmarks rather than tautological closure.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim depends on the existence of a sufficient condition that enforces physical consistency for arbitrary simple rules and on the nonstandard finite difference scheme preserving positivity, conservation, and fixed points by construction; no free parameters or new entities are described in the abstract.

axioms (1)
  • domain assumption Water motion and distribution rules satisfy a sufficient condition ensuring alignment with macroscopic physical conditions governing water movement.
    Invoked in the abstract to justify that solutions remain physically meaningful despite simple rules.

pith-pipeline@v0.9.0 · 5826 in / 1363 out tokens · 57063 ms · 2026-05-20T23:44:51.050539+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

31 extracted references · 31 canonical work pages

  1. [1]

    CA TASTROPHIC NA TURAL DISASTERS AND ECONOMIC GROWTH

    Cavallo E, Galiani S, Noy I, Pantano J. CA TASTROPHIC NA TURAL DISASTERS AND ECONOMIC GROWTH. Review of Economics and Statistics. 2013;95(5):1549–1561. doi: 10.1162/REST_a_00413

    work page doi:10.1162/rest_a_00413 2013

  2. [2]

    Climate Change 2021: The Physical Science Basis

    Masson-Delmotte V , Zhai P , Pirani A, et al. Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change . Cambridge, United Kingdom and New Y ork, NY , USA: Cambridge University Press, 2021

    work page 2021

  3. [3]

    2020 Global Natural Disaster Assessment Report

    ADREM . 2020 Global Natural Disaster Assessment Report. Report; 2021

    work page 2020

  4. [4]

    A Brief review of flood forecasting techniques and their applications

    Jain SK, Mani P , Jain SK, et al. A Brief review of flood forecasting techniques and their applications. International Journal of River Basin Management. 2018;16(3):329–344. doi: 10.1080/15715124.2017.1411920

    work page doi:10.1080/15715124.2017.1411920 2018

  5. [5]

    Urban surface water flood modelling a comprehensive review of current models and future challenges

    Guo K, Guan M, Y u D. Urban surface water flood modelling a comprehensive review of current models and future challenges. Hydrol. Earth Syst. Sci.. 2021;25(5):2843–2860. doi: 10.5194/hess-25-2843-2021

    work page doi:10.5194/hess-25-2843-2021 2021

  6. [6]

    A review of flood modeling methods for urban pluvial flood application

    Bulti DT, Abebe BG. A review of flood modeling methods for urban pluvial flood application. Modeling Earth Systems and Environment. 2020;6(3):1293–1302. doi: 10.1007/s40808-020-00803-z

    work page doi:10.1007/s40808-020-00803-z 2020

  7. [7]

    Cellular Automata and Finite V olume solvers converge for 2D shallow flow modelling for hydrological modelling

    Caviedes-V oullième D, Fernández-Pato J, Hinz C. Cellular Automata and Finite V olume solvers converge for 2D shallow flow modelling for hydrological modelling. Journal of Hydrology. 2018;563:411–417. doi: 10.1016/j.jhydrol.2018.06.021

    work page doi:10.1016/j.jhydrol.2018.06.021 2018

  8. [8]

    A 2D FLOOD INUNDA TION MODEL BASED ON CELLULAR AUTOMA TA APPROACH

    Dottori F, Todini E. A 2D FLOOD INUNDA TION MODEL BASED ON CELLULAR AUTOMA TA APPROACH. In: International Conference on Computational Methods in Water Resources. 2010. 16 TA YLORET AL

    work page 2010

  9. [9]

    A weighted cellular automata 2D inundation model for rapid flood analysis

    Guidolin M, Chen AS, Ghimire B, Keedwell EC, Djordjevi S, Savi DA. A weighted cellular automata 2D inundation model for rapid flood analysis. Environmental Modelling & Software. 2016;84:378–394. doi: 10.1016/j.envsoft.2016.07.008

    work page doi:10.1016/j.envsoft.2016.07.008 2016

  10. [10]

    A Cellular Automata Fast Flood Evaluation (CAffé) Model

    Jamali B, Bach PM, Cunningham L, Deletic A. A Cellular Automata Fast Flood Evaluation (CAffé) Model. Water Resources Research. 2019;55(10):6605–6621. doi: 10.1029/2018WR023679

    work page doi:10.1029/2018wr023679 2019

  11. [11]

    A two scale cellular automaton for flow dynamics modeling (2CAFDYM)

    Kassogué H, Bernoussi A, Maâtouk M, Amharref M. A two scale cellular automaton for flow dynamics modeling (2CAFDYM). Applied Mathematical Modelling. 2017;43:61–77. doi: 10.1016/j.apm.2016.10.034

    work page doi:10.1016/j.apm.2016.10.034 2017

  12. [12]

    A Cellular Automata Based Rainfall-Runoff Model for Urban Inundation Analysis Under Different Land Uses

    Y ao S, Chen NC, Du WY , Wang C, Chen CZ. A Cellular Automata Based Rainfall-Runoff Model for Urban Inundation Analysis Under Different Land Uses. Water Resources Management. 2021;35(6):1991–2006. doi: 10.1007/s11269-021-02826-2

    work page doi:10.1007/s11269-021-02826-2 2021

  13. [13]

    A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling

    Bates PD, Horritt MS, Fewtrell TJ. A simple inertial formulation of the shallow water equations for efficient two-dimensional flood inundation modelling. Journal of Hydrology. 2010;387(1–2):33–45. doi: 10.1016/j.jhydrol.2010.03.027

    work page doi:10.1016/j.jhydrol.2010.03.027 2010

  14. [14]

    A rapid flood inundation model for urban flood analyses

    Wijaya OT, Y ang TH, Hsu HM, Gourbesville P . A rapid flood inundation model for urban flood analyses. Methodsx. 2023;10:102202. doi: 10.1016/j.mex.2023.102202

    work page doi:10.1016/j.mex.2023.102202 2023

  15. [15]

    Hydrological modelling of urbanized catchments: A review and future directions

    Salvadore E, Bronders J, Batelaan O. Hydrological modelling of urbanized catchments: A review and future directions. Journal of Hydrology. 2015;529:62–81. doi: 10.1016/j.jhydrol.2015.06.028

    work page doi:10.1016/j.jhydrol.2015.06.028 2015

  16. [16]

    Dynamic-wave cellular automata framework for shallow water flow modeling

    Chang TJ, Y u HL, Wang CH, Chen AS. Dynamic-wave cellular automata framework for shallow water flow modeling. Journal of Hydrology. 2022;613:128449. doi: 10.1016/j.jhydrol.2022.128449

    work page doi:10.1016/j.jhydrol.2022.128449 2022

  17. [17]

    Introduction to Applied Nonlinear Dynamical Systems and Chaos

    Wiggins S. Introduction to Applied Nonlinear Dynamical Systems and Chaos . Springer New Y ork, NY . 2 ed., 2010

    work page 2010

  18. [18]

    Urban flood numerical simulation: Research, methods and future perspectives

    Luo P , Luo M, Li F, et al. Urban flood numerical simulation: Research, methods and future perspectives. Environmental Modelling & Software. 2022;156:105478. doi: 10.1016/j.envsoft.2022.105478

    work page doi:10.1016/j.envsoft.2022.105478 2022

  19. [19]

    A Cellular Automaton Model for Freeway Traffic

    Nagel K, Schreckenberg M. A Cellular Automaton Model for Freeway Traffic. Journal De Physique I. 1992;2(12):2221–2229. doi: 10.1142/S0217987892900382

    work page doi:10.1142/s0217987892900382 1992

  20. [20]

    Universality and Complexity in Cellular Automata

    Wolfram S. Universality and Complexity in Cellular Automata. Physica D. 1984;10(1–2):1–17. doi: 10.1016/0167-2789(84)90245-8

    work page doi:10.1016/0167-2789(84)90245-8 1984

  21. [21]

    A class of nonstandard numerical methods for autonomous dynamical systems

    Wood DT, Kojouharov HV . A class of nonstandard numerical methods for autonomous dynamical systems. Applied Mathematics Letters. 2015;50:78–82. doi: 10.1016/j.aml.2015.06.008

    work page doi:10.1016/j.aml.2015.06.008 2015

  22. [22]

    A nonstandard finite difference method for n-dimensional productivedestructive systems

    Wood DT, Dimitrov DT, Kojouharov HV . A nonstandard finite difference method for n-dimensional productivedestructive systems. Journal of Difference Equations and Applications. 2015;21(3):240–254. doi: 10.1080/10236198.2014.997228

    work page doi:10.1080/10236198.2014.997228 2015

  23. [23]

    Contributions to the mathematics of the nonstandard finite difference method and applications

    Anguelov R, Lubuma JMS. Contributions to the mathematics of the nonstandard finite difference method and applications. Numerical Methods for Partial Differential Equations. 2001;17(5):518–543. doi: 10.1002/num.1025

    work page doi:10.1002/num.1025 2001

  24. [24]

    CONTRIBUTIONS TO THE THEORY OF NON-STANDARD FINITE DIFFERENCE METHODS AND AP- PLICA TIONS TO SINGULAR PERTURBA TION PROBLEMS

    Lubuma JMS, Patidar KC. CONTRIBUTIONS TO THE THEORY OF NON-STANDARD FINITE DIFFERENCE METHODS AND AP- PLICA TIONS TO SINGULAR PERTURBA TION PROBLEMS. https://doi.org/10.1142/9789812703316_0012. 2005. Online; accessed on 2023-10-01

    work page doi:10.1142/9789812703316_0012 2005

  25. [25]

    Nonstandard finite difference models of differential equations

    Mickens RE. Nonstandard finite difference models of differential equations . World Scientific, 1994

    work page 1994

  26. [26]

    Impact of urbanization on flooding: The Thirusoolam sub watershed - A case study

    Suriya S, Mudgal BV . Impact of urbanization on flooding: The Thirusoolam sub watershed - A case study. Journal of Hydrology. 2012;412:210–

    work page 2012

  27. [27]

    doi: 10.1016/j.jhydrol.2011.05.008

    work page doi:10.1016/j.jhydrol.2011.05.008 2011

  28. [28]

    The MA TLAB ODE suite

    Shampine LF, Reichelt MW. The MA TLAB ODE suite. Siam Journal on Scientific Computing. 1997;18(1):1–22. doi: 10.1137/S1064827594276424

    work page doi:10.1137/s1064827594276424 1997

  29. [29]

    A spatially distributed hydrodynamic model framework for urban flood hydrological and hydraulic processes involving drainage flow quantification

    Guo K, Guan M, Y an H, Xia X. A spatially distributed hydrodynamic model framework for urban flood hydrological and hydraulic processes involving drainage flow quantification. Journal of Hydrology. 2023;625:130135. doi: 10.1016/j.jhydrol.2023.130135

    work page doi:10.1016/j.jhydrol.2023.130135 2023

  30. [30]

    Harris, K

    Harris CR, Millman KJ, Walt v. dSJ, et al. Array programming with NumPy. Nature. 2020;585(7825):357–362. doi: 10.1038/s41586-020-2649-2

    work page doi:10.1038/s41586-020-2649-2 2020

  31. [31]

    D., 2007, @doi [Computing in Science and Engineering] 10.1109/MCSE.2007.55 , 9, 90

    Hunter JD. Matplotlib: A 2D Graphics Environment. Computing in Science & Engineering. 2007;9(3):90–95. doi: 10.1109/MCSE.2007.55

    work page doi:10.1109/mcse.2007.55 2007