CaMa-Flood-GPU: a GPU-based hydrodynamic model implementation for scalable global simulations
- 1State Key Laboratory of Water Resources Engineering and Management, Wuhan University, Wuhan, P.R. China
- 2Institute of Industrial Science, The University of Tokyo, Tokyo, Japan
Abstract
Floods are among the costliest natural hazards, demanding scalable models to simulate river and floodplain dynamics at a global scale. The Catchment-based Macro-scale Floodplain (CaMa-Flood) model is a leading system for this purpose, but its CPU-based implementation is computationally demanding. This paper introduces CaMa-Flood-GPU, a fundamental refactoring of the model optimized for Graphics Processing Unit (GPU) architectures. We systematically reinterpreted its core algorithms – including river routing on irregular networks, runoff interpolation, and water depth diagnosis – into highly parallel, GPU-native operations. Key challenges were addressed by implementing scatter-add for flux updates, sparse matrix multiplication for runoff mapping, and branchless kernels for floodplain dynamics, all while preserving the original model's physical fidelity. Implemented in Python with Triton kernels and PyTorch, CaMa-Flood-GPU achieves multi-GPU scalability through optimized communication patterns that minimize synchronization overhead. The software adopts a modular structure with optional components (e.g., bifurcation routing, adaptive time stepping) and flexible data interfaces. Benchmarks demonstrate an order-of-magnitude speedup over a 192-core CPU baseline and near-linear scaling on multiple GPUs, with negligible numerical differences from the original model. This performance leap reduces simulation times for high-resolution global runs from days to hours, enabling larger ensembles and rapid scenario analysis. By providing a reproducible and efficient tool, CaMa-Flood-GPU lowers the barrier for adopting GPU acceleration in large-scale hydrology. The released implementation provides a reproducible reference for future method development.
Related
- CMIP6‐Driven CaMa‐Flood Simulations Reveal Anthropogenic Control on Global Flood Extremes via Precipitation and Snowmelt Drivers
- Implementation of MEGAN v3.2 in SURFEX v9.0: A new framework for global biogenic emission simulations
- Satellite altimetry reveals heterogeneity in global flood peak timing and its shifts
- Global estimation of bankfull river discharge reveals distinct flood recurrences across climate zones
- KNOU-Flood: Advancing global streamflow forecasting with lead-aware probabilistic machine learning
BibTeX
@article{kang2026camafloodgpu,
title = {{CaMa-Flood-GPU: a GPU-based hydrodynamic model implementation for scalable global simulations}},
author = {Shengyu Kang and Jiabo Yin and Dai Yamazaki},
journal = {Geoscientific Model Development},
year = {2026},
month = jun,
doi = {10.5194/gmd-19-5623-2026},
url = {https://doi.org/10.5194/gmd-19-5623-2026},
note = {CC BY},
}
RIS
TY - JOUR
TI - CaMa-Flood-GPU: a GPU-based hydrodynamic model implementation for scalable global simulations
AU - Shengyu Kang
AU - Jiabo Yin
AU - Dai Yamazaki
PY - 2026
DA - 2026/06/01
JO - Geoscientific Model Development
DO - 10.5194/gmd-19-5623-2026
UR - https://doi.org/10.5194/gmd-19-5623-2026
AB - Floods are among the costliest natural hazards, demanding scalable models to simulate river and floodplain dynamics at a global scale. The Catchment-based Macro-scale Floodplain (CaMa-Flood) model is a leading system for this purpose, but its CPU-based implementation is computationally demanding. This paper introduces CaMa-Flood-GPU, a fundamental refactoring of the model optimized for Graphics Processing Unit (GPU) architectures. We systematically reinterpreted its core algorithms – including river routing on irregular networks, runoff interpolation, and water depth diagnosis – into highly parallel, GPU-native operations. Key challenges were addressed by implementing scatter-add for flux updates, sparse matrix multiplication for runoff mapping, and branchless kernels for floodplain dynamics, all while preserving the original model's physical fidelity. Implemented in Python with Triton kernels and PyTorch, CaMa-Flood-GPU achieves multi-GPU scalability through optimized communication patterns that minimize synchronization overhead. The software adopts a modular structure with optional components (e.g., bifurcation routing, adaptive time stepping) and flexible data interfaces. Benchmarks demonstrate an order-of-magnitude speedup over a 192-core CPU baseline and near-linear scaling on multiple GPUs, with negligible numerical differences from the original model. This performance leap reduces simulation times for high-resolution global runs from days to hours, enabling larger ensembles and rapid scenario analysis. By providing a reproducible and efficient tool, CaMa-Flood-GPU lowers the barrier for adopting GPU acceleration in large-scale hydrology. The released implementation provides a reproducible reference for future method development.
SN - 1991-959X
SN - 1991-9603
ER -