-
On the performance of sequential Bayesian update for database of diverse tsunami scenarios
Authors:
Reika Nomura,
Louise A. Hirao Vermare,
Saneiki Fujita,
Donsub Rim,
Shuji Moriguchi,
Randall J. LeVeque,
Kenjiro Terada
Abstract:
Although the sequential tsunami scenario detection framework was validated in our previous work, several tasks remain to be resolved from a practical point of view. This study aims to evaluate the performance of the previous tsunami scenario detection framework using a diverse database consisting of complex fault rupture patterns with heterogeneous slip distributions. Specifically, we compare the…
▽ More
Although the sequential tsunami scenario detection framework was validated in our previous work, several tasks remain to be resolved from a practical point of view. This study aims to evaluate the performance of the previous tsunami scenario detection framework using a diverse database consisting of complex fault rupture patterns with heterogeneous slip distributions. Specifically, we compare the effectiveness of scenario superposition to that of the previous most likely scenario detection method. Additionally, how the length of the observation time window influences the accuracy of both methods is analyzed. We utilize an existing database comprising 1771 tsunami scenarios targeting the city Westport (WA, U.S.), which includes synthetic wave height records and inundation distributions as the result of fault rupture in the Cascadia subduction zone. The heterogeneous patterns of slips used in the database increase the diversity of the scenarios and thus make it a proper database for evaluating the performance of scenario superposition. To assess the performance, we consider various observation time windows shorter than 15 minutes and divide the database into five testing and learning sets. The evaluation accuracy of the maximum offshore wave, inundation depth, and its distribution is analyzed to examine the advantages of the scenario superposition method over the previous method. We introduce the dynamic time warping (DTW) method as an additional benchmark and compare its results to that of the Bayesian scenario detection method.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Accelerating wave-propagation algorithms with adaptive mesh refinement using the Graphics Processing Unit (GPU)
Authors:
Xinsheng Qin,
Randall J. LeVeque,
Michael R. Motley
Abstract:
Clawpack is a library for solving nonlinear hyperbolic partial differential equations using high-resolution finite volume methods based on Riemann solvers and limiters. It supports Adaptive Mesh Refinement (AMR), which is essential in solving multi-scale problems. Recently, we added capabilities to accelerate the code by using the Graphics Process Unit (GPU). Routines that manage CPU and GPU AMR d…
▽ More
Clawpack is a library for solving nonlinear hyperbolic partial differential equations using high-resolution finite volume methods based on Riemann solvers and limiters. It supports Adaptive Mesh Refinement (AMR), which is essential in solving multi-scale problems. Recently, we added capabilities to accelerate the code by using the Graphics Process Unit (GPU). Routines that manage CPU and GPU AMR data and facilitate the execution of GPU kernels are added. Customized and CPU thread-safe memory managers are designed to manage GPU and CPU memory pools, which is essential in eliminating the overhead of memory allocation and de-allocation. A global reduction is conducted every time step for dynamically adjusting the time step based on Courant number restrictions. Some small GPU kernels are merged into bigger kernels, which greatly reduces kernel launching overhead. A speed-up between $2$ and $3$ for the total running time is observed in an acoustics benchmark problem.
△ Less
Submitted 8 August, 2018;
originally announced August 2018.
-
A High-Resolution Finite Volume Seismic Model to Generate Seafloor Deformation for Tsunami Modeling
Authors:
Christopher J. Vogl,
Randall J. LeVeque
Abstract:
A high-resolution finite volume method approach to incorporating time-dependent slip across rectangular subfaults when modeling general fault geometry is presented. The fault slip is induced by a modification of the Riemann problem to the linear elasticity equations across cell interfaces aligned with the subfaults. This is illustrated in the context of the high-resolution wave-propagation algorit…
▽ More
A high-resolution finite volume method approach to incorporating time-dependent slip across rectangular subfaults when modeling general fault geometry is presented. The fault slip is induced by a modification of the Riemann problem to the linear elasticity equations across cell interfaces aligned with the subfaults. This is illustrated in the context of the high-resolution wave-propagation algorithms that are implemented in the open source Clawpack software (www.clawpack.org), but this approach could be easily incorporated into other Riemann solver based numerical methods. Surface deformation results are obtained in both two and three dimensions and compared to those given by the steady-state, homogeneous half-space Okada solution.
△ Less
Submitted 13 May, 2017; v1 submitted 4 January, 2017;
originally announced January 2017.