Crustal anisotropy from shear-wave splitting of local earthquakes in the Garhwal Lesser Himalaya

Abstract

SUMMARYCrustal anisotropy of the Garhwal Lesser Himalaya has been studied using local earthquake data from the Tehri seismic network. Earthquakes with magnitude (mL) up to 3, which occurred between January 2008 to December 2010, have been used for the shear wave splitting (SWS) analysis. SWS measurements have been done for steeply incident ray paths (ic ≤ 45°) to estimate the anisotropy fast axis orientation (ϕ) and the delay time (∂t). A total of 241 waveforms have been analysed, which yielded 209 splitting measurements, and 32 null results. The analysis reveals spatial and depth variation of ϕ and ∂t, suggesting complex anisotropic structure beneath the Garhwal Lesser Himalaya. The mean ∂t is estimated to be 0.07 ± 0.065 s with a mean depth normalized ∂t of 0.005 s km–1. We present the ϕ and Vs per cent anisotropy results by segregating these as a function of depth, for earthquakes originating above and below the Main Himalayan Thrust (MHT); and spatially, for stations located in the Outer Lesser Himalaya (OLH) and the Inner Lesser Himalaya (ILH). Earthquakes above the MHT sample only the Himalayan wedge, while those below the MHT sample both the underthrust Indian crust and the Himalayan wedge. Within the Himalayan wedge, for both OLH and ILH, the mean ϕ is oriented NE–SW, in the direction of maximum horizontal compressive stress axis (SHmax). This anisotropy is possibly due to stress-aligned microcracks controlled by the local stress pattern within the Himalayan wedge. The mean of normalized ∂t for all events originating within the Himalaya is 0.006 s km–1, which yields a Vs per cent anisotropy of ∼2.28 per cent. Assuming a homogeneous distribution of stress-aligned microcracks we compute a crack density of ∼0.0228 for the Garhwal Lesser Himalaya. At stations close to the regional fault systems, the mean ϕ is subparallel to the strike of the faults, and the anisotropy, locally, appears to be structure-related. For earthquakes originating below the MHT, in OLH, the mean ϕ orientation matches those from the Himalayan wedge and the normalized ∂t decreases with depth. This suggests depth localization of the anisotropy, primarily present within the Himalayan wedge. In the ILH, we observe large variations in the mean ϕ orientation and larger values of ∂t close to the regional fault/thrust systems. This is possibly a composite effect of the structure-related shallow crustal anisotropy and the frozen anisotropy of the underthrusting Indian crust. However, these cannot be segregated in this study.