Open AccessEstimating source time function and moment tensor from moment tensor rate functions by constrained L 1 norm minimization
Author(s) -
Wéber Zoltán
Publication year - 2009
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2009.04202.x
Subject(s) - moment tensor , moment (physics) , tensor (intrinsic definition) , mathematics , waveform , inversion (geology) , inverse problem , focal mechanism , algorithm , computer science , mathematical analysis , physics , geology , geometry , classical mechanics , magnitude (astronomy) , radar , seismology , astronomy , telecommunications , fault (geology) , tectonics
SUMMARY Linear inversion of three‐component waveform data for the time‐varying moment tensor rate functions (MTRFs) is a powerful method for studying seismic sources. After finding the MTRFs, however, we should try to represent an earthquake by just one moment tensor and one source time function (STF), if possible. This approach is particularly justified when dealing with weak local events. Unfortunately, extraction of a moment tensor and STF from the MTRFs is essentially a non‐linear inverse problem. In this paper, we introduce an iterative L p norm minimization technique to retrieve the best moment tensor and STF from the MTRFs obtained by waveform inversion. To allow only forward slip during the rupture process, we impose a positivity constraint on the STF. The error analysis, carried out by using Monte Carlo simulation, allows us to estimate and display the uncertainties of the retrieved source parameters. On the basis of the resulting moment tensor uncertainties, the statistical significance of the double‐couple, compensated linear vector dipole and volumetric parts of the solution can be readily assessed. Tests on synthetic data indicate that the proposed algorithm gives good results for both simple and complex sources. Confidence zones for the retrieved STFs are usually fairly large. The mechanisms, on the other hand, are mostly well resolved. The scalar seismic moments are also determined with acceptable accuracy. If the MTRFs cannot resolve the complex nature of a source, the method yields the average source mechanism. If the subevents are well separated in time, their mechanisms can be estimated by appropriately splitting the MTRFs into subintervals. The method has also been applied to two local earthquakes that occurred in Hungary. The isotropic component of the moment tensor solutions is insignificant, implying the tectonic nature of the investigated events. The principal axes of the source mechanisms agree well with the main stress pattern published for the epicentral region.