Abstract:Ice-sheet flow models capable of accurately projecting their future mass balance constitute tools to improve flood risk assessment and assist sea-level rise mitigation associated with enhanced ice discharge. Some processes that need to be captured, such as grounding-line migration, require high spatial resolution (under the kilometer scale). Conventional ice flow models mainly execute on central processing units (CPUs), which feature limited parallel processing capabilities and peak memory bandwidth. This may hinder model scalability and result in long run times, requiring significant computational resources. As an alternative, graphics processing units (GPUs) are ideally suited for high spatial resolution, as the calculations can be performed concurrently by thousands of threads, processing most of the computational domain simultaneously. In this study, we combine a GPU-based approach with the pseudo-transient (PT) method, an accelerated iterative and matrix-free solution strategy, and investigate its performance for finite elements and unstructured meshes with application to two-dimensional (2-D) models of real glaciers at a regional scale. For both the Jakobshavn and Pine Island glacier models, the number of nonlinear PT iterations required to converge a given number of vertices ( N ) scales in the order of O( N 1.2 ) or better. We further compare the performance of the PT CUDA C implementation with a standard finite-element CPU-based implementation using the price-to-performance metric. The price of a single Tesla V100 GPU is 1.5 times that of two Intel Xeon Gold 6140 CPUs. We expect a minimum speedup of at least 1.5 times to justify the Tesla V100 GPU price to performance. Our developments result in a GPU-based implementation that achieves this goal with a speedup beyond 1.5 times. This study represents a first step toward leveraging GPU processing power, enabling more accurate polar ice discharge predictions. The insights gained will benefit efforts to diminish spatial resolution constraints at higher computing performance. The higher computing performance will allow for ensembles of ice-sheet flow simulations to be run at the continental scale and higher resolution, a previously challenging task. The advances will further enable the quantification of model sensitivity to changes in upcoming climate forcings. These findings will significantly benefit process-oriented sea-level-projection studies over the coming decades.

FastFlow: GPU Acceleration of Flow and Depression Routing for Landscape Simulation

Efficient Debris-flow Simulation for Steep Terrain Erosion

High-quality Texture-based Flow Visualization on Surfaces.

1D and 2D Flow Routing on a Terrain

Extraction of Drainage Network from Grid Terrain Datasets Using Parallel Computing

An Efficient GPU Implementation of a Coupled Overland-Sewer Hydraulic Model with Pollutant Transport

Parallel assignment of flow directions over flat surfaces in massive digital elevation models

GPU-based Dynamic Quad Stream for Forest Rendering

Efficient Modeling and Real-Time Rendering of Massive Urban Pipelines Based on GPU

Graphics-processing-unit-accelerated ice flow solver for unstructured meshes using the Shallow-Shelf Approximation (FastIceFlo v1.0.1)

Quad Mesh Based Dynamic Simulation of Large-Scale Forest on GPU

Parallel Transient Stability-Constrained Optimal Power Flow Using GPU as Coprocessor.

Multi-screen Tiled Displayed, Parallel Rendering System for a Large Terrain Dataset

Parallel Implementation for Phase-Field Simulation of Flow Effect on Dendritic Growth with GPU Acceleration

Implementation of a direct-addressing based lattice Boltzmann GPU solver for multiphase flow in porous media

FastGR : Global Routing on CPU-GPU with Heterogeneous Task Graph Scheduler

Simulating Dam-Break Flooding with Floating Objects Through Intricate City Layouts Using Gpu-Based Sph Method

Simple and fast terrain rendering using graphics hardware

Method for portable, scalable, and performant GPU-accelerated simulation of multiphase compressible flow

GPU‐Accelerated Urban Flood Modeling Using a Nonuniform Structured Grid and a Super Grid Scale River Channel

Method for scalable and performant GPU-accelerated simulation of multiphase compressible flow