Crustal scale heterogeneity across Bhutan Himalayas: insights from local earthquake shear-wave splitting
- 1Department of Applied Geophysics, Indian Institute of Technology (Indian School of Mines) Dhanbad
- 2Department of Physics, University of Naples Federico II
Abstract
In this study, we determine the fast polarization direction (FPD) and delay time ( δ t) of shear wave splitting from 313 Bhutan Himalayan crustal earthquakes recorded by 44 stations of the GANSSER (Geodynamics ANd Seismic Structure of the Eastern Himalaya Region) network deployed and operational in Bhutan from January 2013 to December 2014. We acquired 1686 high-quality measurements by using the Multiple Filter Automatic Splitting Technique (MFAST), which enables rapid processing of millions of events. The delay-time was found to vary between 0.125-0.6 s, indicating the presence of strong crustal anisotropy with an increase in δ t with depth. We also observe that the crustal anisotropy varies from ∼ 3% at 0-20 km depth to ∼ 5% in the 20-35 km depth interval, and more than 8% below 35 km hypocentral depth. The fast directions across the region is generally aligning with the regional stress directions, which is significantly different from the predominant trend of SW-NE to WSW-ENE as observed from the upper mantle splitting. An increase in the number of anisotropic measurements in the western and central Bhutan may be linked to Indian plate subduction. Bhutan Himalaya displays a noticeable along-strike variation in crustal thickness and geometry of Main Himalayan Thrust across a short lateral distance ( ∼ 90 km), which is significantly reflected in the anisotropic patterns, due to the uneven advancement of the Indian plate below western and eastern Bhutan. Integrated patterns of crustal and mantle splitting parameters provide new insights into the deformation processes beneath the Bhutan Himalaya and their implications in regional geodynamics.
Related
- Sustained crustal helium enrichment associated with postseismic viscoelastic relaxation following the 2025 Myanmar earthquake
- Was the mantle warmer when Pangea broke up? insights from initial oceanic crustal thickness alongside the rifted margins of the Atlantic and Indian Oceans
- The Brittle-Plastic Transition in Quartz-Albite Mixtures: New Insights From Shear Deformation Experiments at Mid-to-Lower Crustal Depth Conditions
- The 2016 Mw 7.0 Kumamoto Earthquake Sequence, Japan revisited: Insights from Spatio-Temporal Analysis of Seismicity Parameters
- Dependence of seismic hazard assessment on the observation time interval: insights from a synthetic earthquake catalogue in southeastern Spain
BibTeX
@article{jana2026crustal,
title = {{Crustal scale heterogeneity across Bhutan Himalayas: insights from local earthquake shear-wave splitting}},
author = {Niptika Jana and Mita Uthaman},
journal = {Scientific Reports},
year = {2026},
month = jul,
doi = {10.1038/s41598-026-61863-2},
url = {https://doi.org/10.1038/s41598-026-61863-2},
note = {CC BY},
}
RIS
TY - JOUR
TI - Crustal scale heterogeneity across Bhutan Himalayas: insights from local earthquake shear-wave splitting
AU - Niptika Jana
AU - Mita Uthaman
PY - 2026
DA - 2026/07/01
JO - Scientific Reports
DO - 10.1038/s41598-026-61863-2
UR - https://doi.org/10.1038/s41598-026-61863-2
AB - In this study, we determine the fast polarization direction (FPD) and delay time ( δ t) of shear wave splitting from 313 Bhutan Himalayan crustal earthquakes recorded by 44 stations of the GANSSER (Geodynamics ANd Seismic Structure of the Eastern Himalaya Region) network deployed and operational in Bhutan from January 2013 to December 2014. We acquired 1686 high-quality measurements by using the Multiple Filter Automatic Splitting Technique (MFAST), which enables rapid processing of millions of events. The delay-time was found to vary between 0.125-0.6 s, indicating the presence of strong crustal anisotropy with an increase in δ t with depth. We also observe that the crustal anisotropy varies from ∼ 3% at 0-20 km depth to ∼ 5% in the 20-35 km depth interval, and more than 8% below 35 km hypocentral depth. The fast directions across the region is generally aligning with the regional stress directions, which is significantly different from the predominant trend of SW-NE to WSW-ENE as observed from the upper mantle splitting. An increase in the number of anisotropic measurements in the western and central Bhutan may be linked to Indian plate subduction. Bhutan Himalaya displays a noticeable along-strike variation in crustal thickness and geometry of Main Himalayan Thrust across a short lateral distance ( ∼ 90 km), which is significantly reflected in the anisotropic patterns, due to the uneven advancement of the Indian plate below western and eastern Bhutan. Integrated patterns of crustal and mantle splitting parameters provide new insights into the deformation processes beneath the Bhutan Himalaya and their implications in regional geodynamics.
SN - 2045-2322
ER -