Seismic Velocity Assessment In The Kachchh Region, India, From Multiple Waveform Functionals

Abstract:

The primary goal of this study is to estimate well constrained crust and upper mantle seismic velocity structure in the Kachchh region of Gujarat, India - an area of active interest for earthquake monitoring purposes. Several models based on ‘stand-alone’ surface wave dispersion and receiver function modeling exist in this area. Here we jointly model the receiver function, surface wave dispersion and, S and shear-coupled PL wavetrains using broadband seismograms of deep (150–700 km), moderate to-large magnitude (5.5–6.8) earthquakes recorded teleseismically at semi-permanent seismograph stations in the Kachchh region, Gujarat, India. While surface wave dispersion and receiver function modeling is computationally fast, full waveform modeling makes use of reflectivity synthetic seismograms. An objective function that measures misfit between all three data is minimized using a very fast simulated annealing (VFSA) approach. Surface wave and receiver function data help reduce the model search space which is explored extensively for detailed waveform fitting. Our estimated crustal and lithospheric thicknesses in this region vary from 32 to 41 km and 70 to 80 km, respectively, while crustal P and S velocities from surface to Moho discontinuity vary from 4.7 to 7.0 km/s and 2.7 to 4.1 km/s, respectively. Our modeling clearly reveals a zone of crustal as well as an asthenospheric upwarping underlying the Kachchh rift zone relative to the surrounding unrifted area. We believe that this feature plays a key role in the seismogenesis of lower crustal earthquakes occurring in the region through the emanation of volatile CO₂ into the hypocentral zones liberating from the crystallization of carbonatite melts in the asthenosphere. Such a crust-mantle structure might be related to the plume–lithosphere interaction during the Deccan/Reunion plume episode (~65 Ma).