We take account of the effect of Earth’s surface topography in quasi-dynamic earthquake cycle simulations using a boundary integral equation method. While we have so far assumed a homogeneous elastic half-space medium with a flat free surface, Earth’s actual surface topography is complicated. Here, we constructed new slip response functions in half-space with an arbitrarily shaped surface t…

Hydrocarbon reservoir pressure depletion leads to stress changes inside the reservoir and ground deformation which is registered at the surface as subsidence. As reservoirs are often overlain by layers of rocksalt (or other evaporites), which are materials that flow so as to relax stresses inside them, there is the potential for time-varying surface subsidence. This work focuses on understandin…

Heterogeneous velocity structures are expected to affect fault rupture dynamics. To quantitatively evaluate some of these effects, we examine a model of dynamic rupture on a frictional fault embedded in an elastic full space, governed by plane strain elasticity, with a pair of off-fault inclusions that have a lower rigidity than the background medium.We solve the elastodynamic problem using the…

We present a new Sv-velocity model of the upper mantle under East and Southeast Asia constrained by the inversion of seismic waveforms recorded by broad-band stations. Seismograms from earthquakes occurred between 1977 and 2012 are collected from about 4786 permanent and temporary stations in the region whenever and wherever available. Automated Multimode Inversion of surface and multiple-S wav…

We present a new technique for modelling the global lithospheric magnetic field at Earth’s surface based on the estimation of equivalent potential field sources. As a demonstration we show an application to magnetic field measurements made by the CHAMP satellite during the period 2009–2010 when it was at its lowest altitude and solar activity was quiet. All three components of the vector fi…

Archaeological baked clay remains provide valuable information about the geomagnetic field in historical past, but determination of the geomagnetic field characteristics, especially intensity, is often a difficult task. This study was undertaken to elucidate the reasons for unsuccessful intensity determination experiments obtained from two different Bulgarian archaeological sites (Nessebar – …

We produce a S-wave velocity model of the crust and upper mantle around the Sorgenfrei–Tornquist Zone, southern Scandinavia, by analysing ambient seismic noise and earthquake recordings on temporary and permanent regional network stations. In a first step, we perform tomographical inversion of surface wave dispersion data from seismic noise to obtain Rayleigh and Love wave phase-velocity maps…

Data assimilation is a technique that optimizes the parameters used in a numerical model with a constraint of model dynamics achieving the better fit to observations. Optimized parameters can be utilized for the subsequent prediction with a numerical model and predicted physical variables are presumably closer to observations that will be available in the future, at least, comparing to those ob…

In this study, we propose a method to determine dislocation Love numbers using co-seismic gravity changes from GRACE measurements. First, we present an observation equation to model GRACE observations taking into account the effect of ocean water mass redistribution. The L-curve method was used to determine the regulation parameter in the inversion of the geopotential dislocation Love numbers c…

Ongoing works on full waveform inversion (FWI) are yielding an increasing number of objective functions as alternative to the traditional L2-waveform. These studies aim at designing more robust functions and inversion strategies to reduce the intrinsic dependence of the FWI results on (1) the initial model and (2) the lowest frequency present in field data. In this work, we perform a comparativ…