This study presents the coherently stacked P-wave receiver functions and bouguer anomaly mapping ofWestern Himalayas (longitude 71◦ -74◦ E and latitude 31◦ -34◦ N ) to estimate crustal thickness. Data used for P-wave receiver function is from a local seismic network of Pakistan whereas, gravity data is extracted from Topex and available in the public domain for research. The crust thickness and average crust Vp/Vs ratio at each station of the network are obtained by coherently stacking the Ps, PpPs, PpSs + PsPs phases of 15seismic stations. The data used in this study was collected from 2012 to 2019, events with magnitudes mb ≥ 6 and epicentral distances 30◦ to 95◦ were chosen. There is a significant difference in Moho depth beneath the broadband seismic stations used for the investigations. Moho depth is 36 km in the south, on average 46 km in the center, and 52 km at the northernmost seismic station of the study area. The crust is generally dipping from South to North for the study area. In order to support this interpreted argument of moho-depth variation from the P-wave receiver function, residual calculation of Bouguer anomaly data was carried out as well. The residuals showing a variation of anomalous data from -89 to 193 mGal in the study area have presented a good correlation and favored the argument of crustal dipping as suggested by the P-wave function from the seismic network. The trend confirms that the crustal thickening and shortening are caused by the collision of Indian and Asian plates.