Viscoelastic orthorhombic full wavefield inversion: Development of multiparameter inversion methods

In this paper we compare two time-domain inversion methods that have been widely applied to the problem of modeling earthquake rupture using strong-motion seismograms. In the multi-window method, each point on the fault is allowed to rupture multiple times. This allows flexibility in the rupture time and hence the rupture velocity. Variations in the slip-velocity function are accommodated by variations in the slip amplitude in each time-window. The single-window method assumes that each point on the fault ruptures only once, when the rupture front passes. Variations in slip amplitude are allowed and variations in rupture velocity are accommodated by allowing the rupture time to vary. Because the multi-window method allows greater flexibility, it has the potential to describe a wider range of faulting behavior; however, with this increased flexibility comes an increase in the degrees of freedom and the solutions are comparatively less stable. We demonstrate this effect using synthetic data for a test model of the Mw 7.3 1992 Landers, California earthquake, and then apply both inversion methods to the actual recordings. The two approaches yield similar fits to the strong-motion data with different seismic moments indicating that the moment is not well constrained by strong-motion data alone. The slip amplitude distribution is similar using either approach, but important differences exist in the rupture propagation models. The single-window method does a better job of recovering the true seismic moment and the average rupture velocity. The multi-window method is preferable when rise time is strongly variable, but tends to overestimate the seismic moment. Both methods work well when the rise time is constant or short compared to the periods modeled. Neither approach can recover the temporal details of rupture propagation unless the distribution of slip amplitude is constrained by independent data.

Download Full-text

Tomotok: python package for tomography of tokamak plasma radiation

Journal of Instrumentation ◽

10.1088/1748-0221/16/12/c12015 ◽

2021 ◽

Vol 16 (12) ◽

pp. C12015

Author(s):

J. Svoboda ◽

J. Cavalier ◽

O. Ficker ◽

M. Imríšek ◽

J. Mlynář ◽

...

Keyword(s):

Diagnostic System ◽

Plasma Radiation ◽

Tokamak Plasma ◽

Tomographic Inversion ◽

Inversion Methods ◽

Building Geometry ◽

User Friendly ◽

Academy Of Sciences ◽

Near Future ◽

Python Package

Abstract A python package, called Tomotok, focused on performing tomographic inversion of tokamak plasma radiation is being developed at the Institute of Plasma Physics of the Czech Academy of Sciences. It aims at providing multiple inversion algorithms with an user friendly interface. In order to enable and ease performing tomographic inversion on different devices worldwide, it is planned to publish this software as open source in the near future. In this contribution, the package structure allowing an easy implementation of various tokamak and diagnostic geometries is described and an overview of the package contents is given. Apart from inversion methods, overview of Tomotok auxiliary content is given. The package provides tools for creating simple synthetic diagnostic system. These can be used for testing and benchmarking the code. This includes tools for building geometry matrices that describe the view of detectors using single line of sight approximation and artificial data generators capable of creating simple or hollow Gaussian profiles. The implemented inversion methods cover the minimum Fisher regularisation, biorthogonal decomposition and linear algebraic methods. The implementation of each method is explained, example results obtained by inverting phantom models are presented and discussed. The computation speed of implemented algorithms is benchmarked and compared.

Download Full-text