<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "https://jats.nlm.nih.gov/nlm-dtd/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" specific-use="SMUR" dtd-version="3.0" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="publisher">EGUsphere</journal-id>
<journal-title-group>
<journal-title>EGUsphere</journal-title>
<abbrev-journal-title abbrev-type="publisher">EGUsphere</abbrev-journal-title>
<abbrev-journal-title abbrev-type="nlm-ta">EGUsphere</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub"></issn>
<publisher><publisher-name>Copernicus Publications</publisher-name>
<publisher-loc>Göttingen, Germany</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.5194/egusphere-2026-4124</article-id>
<title-group>
<article-title>Coupled hydrology, weathering geochemistry and surface-process evolution in a parallel global landscape evolution model</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Salles</surname>
<given-names>Tristan</given-names>
<ext-link>https://orcid.org/0000-0001-6095-7689</ext-link>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
</contrib-group><aff id="aff1">
<label>1</label>
<addr-line>School of Geosciences, The University of Sydney, Sydney, NSW 2006, Australia</addr-line>
</aff>
<pub-date pub-type="epub">
<day>22</day>
<month>07</month>
<year>2026</year>
</pub-date>
<volume>2026</volume>
<fpage>1</fpage>
<lpage>43</lpage>
<permissions>
<copyright-statement>Copyright: &#x000a9; 2026 Tristan Salles</copyright-statement>
<copyright-year>2026</copyright-year>
<license license-type="open-access">
<license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri"  xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p>
</license>
</permissions>
<self-uri xlink:href="https://egusphere.copernicus.org/preprints/2026/egusphere-2026-4124/">This article is available from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-4124/</self-uri>
<self-uri xlink:href="https://egusphere.copernicus.org/preprints/2026/egusphere-2026-4124/egusphere-2026-4124.pdf">The full text article is available as a PDF file from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-4124/egusphere-2026-4124.pdf</self-uri>
<abstract>
<p>&lt;p class=&quot;p1&quot;&gt;goSPL (Global Scalable Paleo Landscape Evolution) is an open-source, parallel, finite-volume landscape evolution model designed to simulate the co-evolution of topography, drainage and sedimentary basins from catchment to planetary scale over geological time. Here, I present the new capabilities which extend the model from a fluvial&amp;ndash;hillslope&amp;ndash;marine sediment-routing engine into a coupled surface-process, near-surface-hydrology and weathering-geochemistry framework. The principal additions are (i) an implicit Dupuit&amp;ndash;Boussinesq water table and a generic capillary-fringe duricrust that armours erodibility (forming either in situ or, with a groundwater-discharge gate, by lateral accumulation confined to valley discharge zones) and is archived in the stratigraphic record; (ii) a mass-conservative &lt;em&gt;solute geochemistry&lt;/em&gt; module in which one or more lumped weathering tracers are dissolved, transported along the groundwater flux, precipitated as crust and routed down the river network to the ocean, with per-species typing, source provenance, spatially variable (lithology-controlled) weatherability and optional in-transit and marine coupling; (iii) a soil/regolith model; (iv) a mesh-native orographic-rain solver; (v) a diagnostic glacial-erosion model; and (vi) a parallel finite-volume flexural isostasy solver for planar meshes. All components are opt-in and are formulated as implicit or steady solves on a distributed PETSc DMPlex mesh, with solver choices (Krylov method, preconditioner, iteration scheme) driven by the conditioning of each operator. I describe the governing equations and the numerical and design choices behind each component, the YAML-based experiment configuration, the analytical-benchmark and regression test suite, and strong-scaling measurements on global production meshes of up to 23.7 million nodes. Finally I summarise the online documentation, the conda/PyPI/Docker distributions, the application interface and post-processing tools, and a suite of fifteen generic examples spanning regional and global scales. goSPL is distributed under the GNU GPL v3 licence.</p>
</abstract>
<counts><page-count count="43"/></counts>
<funding-group>
<award-group id="gs1">
<funding-source>Australian Research Council</funding-source>
<award-id>DP260102926</award-id>
</award-group>
</funding-group>
</article-meta>
</front>
<body/>
<back>
</back>
</article>