Structural, petrophysical, and geological constraints in potential field inversion using the Tomofast-x v1.0 open-source code

The quantitative integration of geophysical measurements with data and information from other disciplines is becoming increasingly important in answering the challenges of undercover imaging and of the modelling of complex areas. We propose a review of the different techniques for the utilisation of structural, petrophysical, and geological information in single physics and joint inversion as implemented in the Tomofast-x open-source inversion platform. We detail the range of constraints that can be applied to the inversion of potential field data. The inversion examples we show illustrate a selection of scenarios using a realistic synthetic data set inspired by real-world geological measurements and petrophysical data from the Hamersley region (Western Australia). Using Tomofast-x's flexibility, we investigate inversions combining the utilisation of petrophysical, structural, and/or geological constraints while illustrating the utilisation of the L-curve principle to determine regularisation weights. Our results suggest that the utilisation of geological information to derive disjoint interval bound constraints is the most effective method to recover the true model. It is followed by model smoothness and smallness conditioned by geological uncertainty and cross-gradient minimisation.

Method of Source Identification following an Accidental Release at an Unknown Location Using a Lagrangian Atmospheric Dispersion Model

A computationally efficient source inversion algorithm was developed and applied with the Lagrangian atmospheric dispersion model DIPCOT. In the process of source location estimation by minimizing a correlation-based cost function, the algorithm uses only the values of the time-integrated concentrations at the monitoring stations instead of all of the individual measurements in the full concentration-time series, resulting in a significant reduction in the number of integrations of the backward transport equations. Following the source location estimation the release start time, duration and emission rate are assessed. The developed algorithm was verified for the conditions of the ETEX-I (European Tracer Experiment—1st release). Using time-integrated measurements from all available stations, the distance between the estimated and true source location was 108 km. The estimated start time of the release was only about 1 h different from the true value, within the possible accuracy of estimate of this parameter. The estimated release duration was 21 h (the true value was 12 h). The estimated release rate was 4.28 g/s (the true value was 7.95 g/s). The estimated released mass almost perfectly fitted the true released mass (323.6 vs. 343.4 kg). It thus could be concluded that the developed algorithm is suitable for further integration in real-time decision support systems.

Assessing the role of selected constraints in Bayesian dynamic source inversion: Application to the 2017 Mw 6.3 Lesvos earthquake

Summary A dynamic finite-fault source inversion for stress and frictional parameters of the Mw 6.3 2017 Lesvos earthquake is carried out. The mainshock occurred on June 12, offshore the southeastern coast of the Greek island of Lesvos in the north Aegean Sea. It caused 1 fatality, 15 injuries, and extensive damage to the southern part of the island. Dynamic rupture evolution is modeled on an elliptic patch, using the linear slip-weakening friction law. The inversion is posed as a Bayesian problem and the Parallel Tempering Markov Chain Monte Carlo algorithm is used to obtain posterior probability distributions by updating the prior distribution with progressively more constraints. To calculate the first posterior distribution, only the constraint that the model should expand beyond the nucleation patch is used. Then, we add the constraint that the model should reach a moment magnitude similar to that obtained from our centroid moment tensor inversion. For the final posterior distribution, 15 acceleration records from Greek and Turkish strong motion networks at near regional distances ($\approx 30 - 150$ km) in the frequency range of 0.05–0.15 Hz are used. The three posterior distributions are compared to understand how much each constraint contributes to resolving different quantities. The most probable values and uncertainties of individual parameters are also calculated, along with their mutual trade-offs. The features best determined by seismograms in the final posterior distribution include the position of the nucleation region, the mean direction of rupture (towards WNW), the mean rupture speed (with 68 per cent of the distribution lying between 1.4–2.6 km/s), radiated energy (12–65 TJ), radiation efficiency (0.09–0.38), and the mean stress drop (2.2–6.5 MPa).

Multicycle Simulation of Strike-Slip Earthquake Rupture for Use in Near-Source Ground-Motion Simulations

ABSTRACT Realistic dynamic rupture modeling validated by observed earthquakes is necessary for estimating parameters that are poorly resolved by seismic source inversion, such as stress drop, rupture velocity, and slip rate function. Source inversions using forward dynamic modeling are increasingly used to obtain earthquake rupture models. In this study, to generate a large number of physically self-consistent rupture models, rupture process of which is consistent with the spatiotemporal heterogeneity of stress produced by previous earthquakes on the same fault, we use multicycle simulations under the rate and state (RS) friction law. We adopt a one-way coupling from multicycle simulations to dynamic rupture simulations; the quasidynamic solver QDYN is used to nucleate the seismic events and the spectral element dynamic solver SPECFEM3D to resolve their rupture process. To simulate realistic seismicity, with a wide range of magnitudes and irregular recurrence, several realizations of 2D-correlated heterogeneous random distributions of characteristic weakening distance (Dc) in RS friction are tested. Other important parameters are the normal stress, which controls the stress drop and rupture velocity during an earthquake, and the maximum value of Dc, which controls rupture velocity but not stress drop. We perform a parametric study on a vertical planar fault and generate a set of a hundred spontaneous rupture models in a wide magnitude range (Mw 5.5–7.4). We validate the rupture models by comparison of source scaling, ground motion (GM), and surface slip properties to observations. We compare the source-scaling relations between rupture area, average slip, and seismic moment of the modeled events with empirical ones derived from source inversions. Near-fault GMs are computed from the source models. Their peak ground velocities and peak ground accelerations agree well with the ground-motion prediction equation values. We also obtain good agreement of the surface fault displacements with observed values.

The Comparison of Performance between Contrast Source Inversion and Its Related Algorithms and Several Improved Methods

The contrast source inversion (CSI) is an effective method for solving microwave imaging problems which is widely utilized. The core of the CSI is to change the conventional inverse scattering problem into an optimization problem. The two items in the objective function describe the state error and data error, respectively. As it is all known, there is almost no complete performance comparison based on Fresnel data for the CSI and its related improved algorithms. In addition, the performance of the algorithm under different weights was not analyzed before and the convergence speed of original CSI is slow. Firstly, this paper compares the performance of traditional CSI and its improved algorithms from three aspects of qualitative imaging effect, convergence speed, and objective function value based on Fresnel data. Secondly, the influence of the state error and the data error under different weights on the convergence rate and the objective function value are studied. For the limitation of a slower convergence rate, the CSI with weights (W-CSI), the CSI with dynamic reduction factor (CSI-DRF), and its related algorithms, which can get better convergence rate compared with their relative original algorithms, are proposed. Eventually, the future research work is prospected.