PROBABILISTIC TSUNAMI HAZARD ASSESSMENT (GUIDELINE)

Introduction

Collection of scripts used for building a probabilistic tsunami hazard assessment (PTHA) for the southern coast of Java, Indonesia (link to the manuscript will be added upon acceptance). Steps taken are below:

1. Generate synthetic ETAS catalogue.
2. Generate stochastic heterogeneous earthquake slip models.
3. Stage 1: Offshore tsunami simulations using JAGURS.
4. Intermediate step before Stage 2
5. Stage 2: Onshore tsunami simulations using SFINCS.
6. Post-processing and analysis.

Note

Some scripts mentioned below will need to use utility scripts from ryanign Some I/O files are available at upon acceptance

Future questions:
Ryan Pranantyo
ignatius.pranantyo@ntu.edu.sg
ignatiusryanpranantyo@gmail.com

ETAS Catalogue

Stochastic heterogeneous slip model

I utilsed RPTHA code (from Geoscience Australia) to generate a set of heterogeneous earthquake slip models.Below were steps taken. I executed the scripts under rptha_py312 conda environment.

1. Generate unit sources displacement.

• Prepare a slab contour shapefile following Slab2.0 subduction model.
• Execute Rscript to generate a downdip shapefile.
• Visually checked and modified as required (inside upon acceptance/input_files__SouthernJava).
• Execute Rscript. I executed through pbs. If it works fine, you will get a set of .tif and .RDS files. Samples files are inside upon acceptance/OUTPUTS__Slab2__Jawa_samples.

2. Generate heterogeneous slip models.

• Execute Rscript.
  • lon and lat are the closest finite fault centroid to the lon and lat from ETAS catalogue.
  • Mw is from the ETAS catalogue.
  • Repeat for all lon, lat, and Mw. Here, I generated 20 samples per-one execution.
  • I executed through pbs.
  • If it works fine, your will get a set of .csv files. After this, you can skip the next step, and directly go to step 3. (saved inside upon acceptance/sffm_combinations/stochastic_slips__SLAB2__Jawa/)
• Execute py script -- which utilise py script -- through pbs. This step to eliminate unrealistic slip models.
• If it works fine, you will get a set of .csv files. (saved inside upon acceptance/SFFM_files_used/)

3. Combine unit sources displacement file and slip models to generate input files for JAGURS.

• Use py script execute using pbs. Results are inside upon acceptance/01__jagurs_displacement__mwX-X_to_Y-Y.zip files

4. Map which displacement model to use in the catalogue.

• Execute using py script

Stage 1: Offshore tsunami simulations using JAGURS

Please refer to JAGURS for the details. Below are files and scripts used.

• Files inside upon acceptance/STAGE-1__JAGURS_OFFSHORE_TEMPLATE.zip/ are template files used in this study. Update gridfile.dat to change between NLSWE and LSWE computation.
• Use one file from upon acceptance/01__jagurs_displacement__XXX.zip as an input. Update the filename or gridfile.dat first.
• Run all events in the catalogue!
• You may use py script submitted through pbs.

Intermediate steps before SFINCS simulations

This stage is to bridge between Stage 1 and Stage 2. Two main focuses are (1) select SFINCS domains to be simulated, and (2) extract tsunami elevation timeseries to be used in SFINCS simulations.

Use py script and execute using pbs.

Stage 2: Onshore tsunami simulations using SFINCS

Please refer to SFINCS for the details. In this study, I used v2.1.1_Dollerup_release./ There are two mains steps:

1. Generate templates for SFINCS simulations.

• Use sfincs template script, can be executed through pbs.

2. Launch the simulations.

• Use sfincs launch through pbs.

Post-processing

Generating offshore and onshore hazard database

Use scripts

Offshore tsunami hazard database

You can execute offshore database through pbs after the completion of intermediate steps. If it works fine, you will get upon acceptance/STAGE-1__RESULTS__OFFSHORE_DATABASE.zip/OFFSHORE_HAZARDS_DF__CATALOGUE_10000-YEARS__.csv.

Onshore tsunami hazard database

Execute onshore database through pbs after the completion of Stage 2. If it works fine, you will get upon acceptance/STAGE-2_RESULTS__ONSHORE_DATABASE.zip/flood_depth__tile_CATALOGUE__10000-YEARS.nc.

PTHA analysis

Below are two main scripts are used to analyse offshore and onshore hazard databases generated.

1. Use offshore ptha for offshore PTHA.

2. Use onshore ptha for onshore PTHA.

• execute onshore exceedance rates to generate annual exceedance rates first.