Earthquake Moment Tensor Inversion Using 3D Velocity Model in the Himalayas

The Himalayan region, shaped by the ongoing collision of the Indian and Eurasian tectonic plates, is one of Earth’s most seismically active and geologically complex areas. The Indian plate moves northeastward at a rate of approximately 5 cm per year, driving tectonic activity in this region. Understanding earthquake source mechanisms in this region is crucial for seismic hazard assessment and geodynamic studies. Moment tensor (MT) inversion, a widely used technique for analysing earthquake faulting mechanisms, matches synthetic waveforms to observed data by minimising the misfit. However, conventional 1D velocity models often fail to capture the region’s complex lateral heterogeneities, leading to inaccuracies in source characterisation. Synthetic waveforms, generated via Green’s functions using frequency waveform (FK) methods and 1D velocity models, are critical for MT solutions, with time shifts playing a pivotal role in achieving optimal waveform correlations.

This study employs a 3D velocity model to improve MT inversion for a Mw 3.5 earthquake on 9 January 2021 (30.76°N, 78.54°E). Green’s functions were generated using the spectral element method for six simulations. Each simulation resulted in three-component waveforms, with a total of 18 synthetics per station. Observed data from 24 broadband stations were analysed, and results were compared to those obtained using 1D models. Slight variations in strike, dip, and rake values underscore the limitations of 1D models in capturing Earth’s heterogeneities.

The study reveals that 3D velocity models significantly enhance MT solution accuracy, particularly in determining focal depths, faulting mechanisms, and seismic moment magnitudes. A probabilistic approach was also applied to quantify the uncertainty associated with MT estimates, providing confidence measures. Extending this approach, MT inversion was performed for another earthquake in the Uttarakhand Himalaya using the same 3D velocity model, further demonstrating the advantages of 3D wavefield simulations in seismically active regions.

Keywords: Himalayas, Moment Tensor, Green’s Function, Spectral element method.