SPGM: A Scalable PaleoGeomorphology Model

Abstract Numerical models of landscape evolution are playing an increasingly important role in providing an improved understanding of geomorphic transport processes shaping Earth’s surface topography. Improving theoretical underpinnings coupled with increasing computational capacity has led to the development of several open source codes written in low-level languages. However, adapting these codes to new functionality or introducing greater flexibility often requires significant recoding. Here we present a multi-process, scalable, numerical model of geomorphological evolution, built with a modular structure and geared toward seamless extensibility. We implement recent algorithmic advances that reduce the computational cost of flow routing – a problem that typically scales quadratically with the number of unknowns – to linear in time while allowing for parallel implementations of geomorphic transport processes. Our scalability tests demonstrate that such parallelizations can achieve an order of magnitude speedup on a typical desktop computer, making large-scale simulations more tractable.

[1]  M. Sambridge,et al.  Modelling landscape evolution on geological time scales: a new method based on irregular spatial discretization , 1997 .

[2]  F. Schlunegger,et al.  Topographic evolution and morphology of surfaces evolving in response to coupled fluvial and hillslope sediment transport , 2003 .

[3]  D. García-Castellanos.  Interplay between lithospheric flexure and river transport in foreland basins , 2002 .

[4]  Bjarne Stroustrup,et al.  C++ Programming Language , 1986, IEEE Softw..

[5]  Ken Perlin,et al.  [Computer Graphics]: Three-Dimensional Graphics and Realism , 2022 .

[6]  T. G. Freeman,et al.  Calculating catchment area with divergent flow based on a regular grid , 1991 .

[7]  J. Braun.  The many surface expressions of mantle dynamics , 2010, Nature Geoscience.

[8]  P. Roache.  Code Verification by the Method of Manufactured Solutions , 2002 .

[9]  Nicole M. Gasparini,et al.  Predictions of steady state and transient landscape morphology using sediment‐flux‐dependent river incision models , 2007 .

[10]  Sébastien Limet,et al.  Parallel Computing Flow Accumulation in Large Digital Elevation Models , 2011, ICCS.

[11]  G. Tucker,et al.  Erosional dynamics, flexural isostasy, and long-lived escarpments: A numerical modeling study , 1994 .

[12]  G. Tucker,et al.  Implications of sediment‐flux‐dependent river incision models for landscape evolution , 2002 .

[13]  M. Simoes,et al.  Continental‐scale erosion and transport laws: A new approach to quantitatively investigate macroscale landscapes and associated sediment fluxes over the geological past , 2010 .

[14]  H. Kooi,et al.  Escarpment evolution on high‐elevation rifted margins: Insights derived from a surface processes model that combines diffusion, advection, and reaction , 1994 .

[15]  T. Simon‐Labric,et al.  Eroding dynamic topography , 2013 .

[16]  Sean D. Willett,et al.  Orogeny and orography: The effects of erosion on the structure of mountain belts , 1999 .

[17]  Dani Lischinski.  Incremental Delaunay Triangulation , 1994, Graphics Gems.

[18]  Chengzhi Qin,et al.  Parallelizing flow-accumulation calculations on graphics processing units - From iterative DEM preprocessing algorithm to recursive multiple-flow-direction algorithm , 2012, Comput. Geosci..

[19]  David G. Tarboton,et al.  Parallel Flow-Direction and Contributing Area Calculation for Hydrology Analysis in Digital Elevation Models , 2009, PDPTA.

[20]  Jean Braun,et al.  A very efficient O(n), implicit and parallel method to solve the stream power equation governing fluvial incision and landscape evolution , 2013 .

[21]  M. Gurnis,et al.  Holographic measurements of surface topography in laboratory models of mantle hotspots , 1989 .

[22]  Philippe Fullsack,et al.  Erosional control of active compressional orogens , 1992 .

[23]  Leonidas J. Guibas,et al.  Primitives for the manipulation of general subdivisions and the computation of Voronoi diagrams , 1983, STOC.

[24]  G. Simpson.  A dynamic model to investigate coupling between erosion, deposition, and three‐dimensional (thin‐plate) deformation , 2004 .

[25]  Patrick Knupp,et al.  Verification of Computer Codes in Computational Science and Engineering , 2002 .

[26]  A. Stephen McGough,et al.  Massively parallel landscape-evolution modelling using general purpose graphical processing units , 2012, 2012 19th International Conference on High Performance Computing.

[27]  J. Z. Zhu,et al.  The finite element method , 1977 .