Moment tensor analysis and uncertainty quantification of local earthquake events: tectonic implication in the northwestern Himalayan region

Regular monitoring of small to moderate sources of continuous earthquake events in the complex tectonics of Himalayan region helps in clearly defining the ongoing seismotectonic process. The study of moment tensor inversion to decipher the fault planes responsible for current seismic activity in the Kishtwar region of Northwest part of Himalaya has been undertaken by establishing a six-station network in 2022 and among them 15 events of shallow origin with magnitude ranging from M_L ~ 3.0 to 4.0 occurred in the local region of seismic network are used for the moment tensor inversion. A few number of studies didn’t able to clearly demarcate the actual scenario of seismotectonics in the northwest part of Himalaya due to its difficult terrain and complex geology. This area has been studied for fault plane solution by a software package ISOLA based on MATLAB programming environment. The source inversion is performed via iterative deconvolution method and synthetic seismogram is generated through green’s function computation via discrete wavenumber method using the regional crustal velocity model. However, the inversion is performed at several trial source position and at various frequency bands based on the epicenter distance and the magnitude of earthquake to find the best solution resulting from the maximum correlation between the recorded and synthetically generated waveforms. A 2D space-time grid search is performed for determining the optimal time and positon of earthquake generation. Perhaps calculating source parameters such as moment magnitude, centroid depth and fault parameters equally with describing uncertainty quantities such as variance reduction factor and condition number will deliver the reliability and stability to the solution. A strong follow-up uncertainty quantification can justify the best estimated fault plane solution. Quality of earthquake event can be calculated through their DC and CLVD percentage and maximum & minimum compression stress direction. Focal mechanism solution of these events following thrust with strike-slip focal mechanism and represents the compressional regime in north-northeastern direction. The centroid depth obtained by moment tensor inversion of all events falls within the depth zone of Main Himalayan Thrust (MHT) suggesting seismicity is concentrated along the major detachment in the region.