TY - JOUR
T1 - Probabilistic, high-resolution tsunami predictions in northern Cascadia by exploiting sequential design for efficient emulation
AU - Salmanidou, Dimitra M.
AU - Beck, Joakim
AU - Pazak, Peter
AU - Guillas, Serge
N1 - KAUST Repository Item: Exported on 2022-01-26
Acknowledgements: Financial support. Serge Guillas was supported by the Alan Turing Institute project “Uncertainty Quantification of Complex Computer Models. Applications to Tsunami and Climate” under the Engineering and Physical Sciences Research Council (EPSRC) (grant no. EP/N510129/1).
PY - 2021/12/17
Y1 - 2021/12/17
N2 - The potential of a full-margin rupture along the Cascadia subduction zone poses a significant threat over a populous region of North America. Previous probabilistic tsunami hazard assessment studies produced hazard curves based on simulated predictions of tsunami waves, either at low resolution or at high resolution for a local area or under limited ranges of scenarios or at a high computational cost to generate hundreds of scenarios at high resolution. We use the graphics processing unit (GPU)-accelerated tsunami simulator VOLNA-OP2 with a detailed representation of topographic and bathymetric features. We replace the simulator by a Gaussian process emulator at each output location to overcome the large computational burden. The emulators are statistical approximations of the simulator's behaviour. We train the emulators on a set of input–output pairs and use them to generate approximate output values over a six-dimensional scenario parameter space, e.g. uplift/subsidence ratio and maximum uplift, that represent the seabed deformation. We implement an advanced sequential design algorithm for the optimal selection of only 60 simulations. The low cost of emulation provides for additional flexibility in the shape of the deformation, which we illustrate here considering two families – buried rupture and splay-faulting – of 2000 potential scenarios. This approach allows for the first emulation-accelerated computation of probabilistic tsunami hazard in the region of the city of Victoria, British Columbia.
AB - The potential of a full-margin rupture along the Cascadia subduction zone poses a significant threat over a populous region of North America. Previous probabilistic tsunami hazard assessment studies produced hazard curves based on simulated predictions of tsunami waves, either at low resolution or at high resolution for a local area or under limited ranges of scenarios or at a high computational cost to generate hundreds of scenarios at high resolution. We use the graphics processing unit (GPU)-accelerated tsunami simulator VOLNA-OP2 with a detailed representation of topographic and bathymetric features. We replace the simulator by a Gaussian process emulator at each output location to overcome the large computational burden. The emulators are statistical approximations of the simulator's behaviour. We train the emulators on a set of input–output pairs and use them to generate approximate output values over a six-dimensional scenario parameter space, e.g. uplift/subsidence ratio and maximum uplift, that represent the seabed deformation. We implement an advanced sequential design algorithm for the optimal selection of only 60 simulations. The low cost of emulation provides for additional flexibility in the shape of the deformation, which we illustrate here considering two families – buried rupture and splay-faulting – of 2000 potential scenarios. This approach allows for the first emulation-accelerated computation of probabilistic tsunami hazard in the region of the city of Victoria, British Columbia.
UR - http://hdl.handle.net/10754/668411
UR - https://nhess.copernicus.org/articles/21/3789/2021/
UR - http://www.scopus.com/inward/record.url?scp=85122008193&partnerID=8YFLogxK
U2 - 10.5194/nhess-21-3789-2021
DO - 10.5194/nhess-21-3789-2021
M3 - Article
SN - 1684-9981
VL - 21
SP - 3789
EP - 3807
JO - Natural Hazards and Earth System Sciences
JF - Natural Hazards and Earth System Sciences
IS - 12
ER -