We performed a receiver function analysis to study crust and uppermost upper mantle structure beneath
We calculated radial receiver functions using extended time multitaper frequency domain cross-correlation receiver function estimation method. We selected stations for which receiver functions with high signal-to-noise ratios are obtained among the 17 stations. The obtained receiver functions show a significant azimuthal dependence, which suggests the complexity of the structures beneath these stations. This is consistent with the large tangential components mentioned above. We performed inversion of the obtained receiver functions to determine velocity structures of crust and upper mantle using a genetic algorithm inversion technique. The results suggest that the Moho beneath the stations in the Potwar region is thick (around a depth of 50 km), while it becomes thinner in the center of