Stress tensor inversion from focal mechanism data: variance component estimation approach for data weighting
摘要
Focal mechanism data are essential for estimating tectonic stress tensors, but heterogeneous data quality and unequal observation precision can bias inversions when uniform weights are assumed. We introduce a weighted least-squares inversion combined with variance component estimation (VCE; BIQUE in a Gauss–Helmert framework) that infers observation variances directly from residuals and assigns statistically consistent weights. In our implementation, one variance component is estimated per mechanism and applied jointly to its strike, dip, and slip angles, providing a parsimonious and catalog-consistent stochastic model. We validate the method on three datasets: a 40-event synthetic test, the 2008 West Bohemia swarm (167 mechanisms), and Agia Varvara in central Crete (31 mechanisms). Compared with uniform weighting, VCE produces sharper and more stable principal stress orientations, with azimuth and plunge uncertainties reduced by up to an order of magnitude. Variance–trace metrics decrease substantially, by roughly one order of magnitude in the synthetic test and two orders in West Bohemia. In Central Crete, the uncertainties shrink by a factor of 2–3, and the stress–shape ratio converges toward