The physics of seismic waves has been exploited to produce detailed shear wave velocity structure beneath highly seismogenic regions, Sumatra and Northeast India from modeling of Receiver Functions obtained from numerous earthquake events throughout the globe. A fine combination of inverse methods such as H-K stacking, Neighbourhood Algorithm, and Joint inversion was used to model the crustal structure as well as composition.