Accurate estimation of the planetary boundary layer (PBL) top is essential for air quality prediction, weather forecast, and assessment of regional and global climate models. In this article, the long-term climatology of seasonal, global distribution of PBL is presented by using global positioning system radio occultation (GPSRO) based payloads such as Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC), Communication/Navigation Outage Forecast System (C/NOFS), TerraSAR-X, and The Gravity Recovery and Climate Experiment (GRACE) from the year 2006–2015. We used Wavelet Covariance Transform (WCT) technique for precise PBL top identification. The derived PBL top from GPSRO data is rigorously evaluated with GPS radiosonde data over Gadanki. Significant seasonal variation is noticed in both radiosonde and GPSRO observations. Further, we compared the PBL obtained GPS RO with global radiosonde network and observed very good correlation. The number of occultations reaching down to 500 m and retrieval rate of PBL top from WCT method is very high in mid-latitudes compared to tropical latitudes. The global distribution of PBL top shows significant seasonal variation with higher during summer followed by spring, fall, and minimum in winter. In the vicinity of Inter Tropical Convergence Zone (ITCZ), the PBL top is high over eastern Pacific compared to other regions. The ERA-Interim reanalysis data underestimate the PBL top compared to GPS RO observations due to different measurement techniques. The seasonal variation of global averaged PBL top over land and ocean shows contrasting features at different latitude bands.