Storm Surge Modeling as an Application of Local Time-stepping in MPAS-Ocean

TL;DR

This paper demonstrates that local time-stepping schemes can be effectively applied in high-resolution global ocean models like MPAS-O, achieving comparable accuracy to traditional methods while significantly reducing computational costs during storm surge simulations.

Contribution

It introduces the first application of local time-stepping schemes in MPAS-Ocean, showing their efficiency and accuracy in storm surge modeling with high-resolution meshes.

Findings

LTS3 produces SSH solutions comparable to RK4.

LTS3 can be up to 35% faster than RK4.

Efficient LTS requires a coarse-to-fine cell ratio of at least 1:5.

Abstract

This paper presents the first scientific application of local time-stepping (LTS) schemes in the Model for Prediction Across Scales-Ocean (MPAS-O). We use LTS schemes in a single-layer, global ocean model that predicts the storm surge around the eastern coast of the United States during Hurricane Sandy. The variable-resolution meshes used are of unprecedentedly high resolution in MPAS-O, containing cells as small as 125 meters wide in Delaware Bay. It is shown that a particular, third-order LTS scheme (LTS3) produces sea-surface height (SSH) solutions that are of comparable quality to solutions produced by the classical four-stage, fourth-order Runge-Kutta method (RK4) with a uniform time step on the same meshes. Furthermore, LTS3 is up to 35% faster in the best cases, showing that LTS schemes are viable for use in MPAS-O with the added benefit of substantially less computational cost. The results of these performance experiments inform us of the requirements for efficient mesh design for LTS schemes. In particular, we see that for LTS to be efficient on a given mesh, it is important to have enough cells using the coarse time-step relative to those using the fine time-step, typically at least 1:5.