TY - JOUR

T1 - A Study of Implementation Features of the r-Solution Method for Tsunami Source Recovery in the Case of the Illapel Tsunami 2015: Implementation Features of the r-Solution Method

AU - Voronina, Tatyana A.

AU - Voronin, Vladislav V.

N1 - Funding Information: We would like to thank Artem Loskutov, the Institute of Marine Geology and Geophysics, the Far East Branch, Russian Academy of Sciences, for his assistance in the preparation of the observational data and to thank Vasily Titov, the Pacific Marine Environmental Laboratory in the NOAA Center for Tsunami Research, for permission and adaptation of the drawing from the paper Tang et.al. (2016 ). The research proposed was carried out under the State Budget Program with ICMMG SB RAS (0315-2020-0005). Publisher Copyright: © 2021, The Author(s), under exclusive licence to Springer Nature Switzerland AG.

PY - 2021/12

Y1 - 2021/12

N2 - A methodology is presented for tuning the main parameters of the r-solution method as applied to the problem of reconstructing the initial tsunami waveform. Based on singular value decomposition, the r-solution method provides an effective means for controlling the instability inherently present in the inverse problem. However, the quality of the inversion depends heavily on the main parameters of the method, which are the number of the right singular vectors which form the basis of a solution space and the number of spatial harmonics used in the tsunami source representation. Unlike previous studies, where these factors were analyzed separately, in this paper they are considered jointly, which allows one to discover the mutual connections between them. The latter makes it possible to improve the inversion results. An optimal number of spatial harmonics decreases the time for matrix calculation and reduces the number of inversion artifacts. The approach for choosing the optimal parameters is illustrated by a case study of the Illapel Tsunami of 2015.

AB - A methodology is presented for tuning the main parameters of the r-solution method as applied to the problem of reconstructing the initial tsunami waveform. Based on singular value decomposition, the r-solution method provides an effective means for controlling the instability inherently present in the inverse problem. However, the quality of the inversion depends heavily on the main parameters of the method, which are the number of the right singular vectors which form the basis of a solution space and the number of spatial harmonics used in the tsunami source representation. Unlike previous studies, where these factors were analyzed separately, in this paper they are considered jointly, which allows one to discover the mutual connections between them. The latter makes it possible to improve the inversion results. An optimal number of spatial harmonics decreases the time for matrix calculation and reduces the number of inversion artifacts. The approach for choosing the optimal parameters is illustrated by a case study of the Illapel Tsunami of 2015.

KW - ill-posed problem

KW - numerical modeling

KW - Tsunamis

KW - waveform inversion

UR - http://www.scopus.com/inward/record.url?scp=85114340618&partnerID=8YFLogxK

U2 - 10.1007/s00024-021-02843-7

DO - 10.1007/s00024-021-02843-7

M3 - Article

AN - SCOPUS:85114340618

VL - 178

SP - 4853

EP - 4863

JO - Pure and Applied Geophysics

JF - Pure and Applied Geophysics

SN - 0033-4553

IS - 12

ER -