The height of the planetary boundary layer (PBL) is an important parameter that relates to the various processes associated with the PBL. In this paper, we use Global Positioning System radio occultation (GPSRO) measurements to derive a global climatology of PBL heights. Utilizing the strength of GPSRO in capturing fine vertical structures, the top of the PBL is defined to be the height at which the vertical gradient of the refractivity or water vapor partial pressure is minimum, corresponding to the height where the refractivity or water vapor pressure changes most rapidly. A “sharpness parameter” is defined that quantifies the applicability of these definitions. The sharpness parameter is largest over the subtropical regions characterized by strong subsidence. When the sharpness parameter is large, the refractivity‐ and moisture‐based heights are shown to converge. We derived global PBL height climatology using three years (Dec. 2006–Nov. 2009) of COSMIC/FORMOSAT‐3 measurements and compared with values calculated from ECMWF Reanalysis Interim (ERA‐Int). We found that the mean PBL heights from GPSRO shared similar spatial and seasonal variations with ERA‐Int; however, GPSRO heights were higher by 500 m. The standard deviation was also higher from GPSRO, especially in the tropics. We present detailed comparisons between GPSRO and ERA‐Int over the Pacific Ocean and the Sahara desert and examine the PBL height distributions as well as its annual and diurnal variabilities. These results suggest that the underlying causes of the bias between GPSRO and ERA‐Int likely vary from region to region.

The global positioning system (GPS) radio occultation (RO) technique has demonstrated the ability to precisely probe earth’s atmosphere globally with high vertical resolution. However, the lowermost troposphere still presents some challenges for the technique. Over moist marine areas, especially in subtropical regions, a very large negative moisture gradient often exists across the thermal inversion capping the marine boundary layer (MBL), which frequently causes superrefraction (SR), or ducting. In the presence of SR, the reconstruction of refractivity from RO data becomes an ill-posed inverse problem. This study shows that one given RO bending angle profile is consistent with a continuum (an infinite number) of refractivity profiles. The standard Abel retrieval gives the minimum refractivity solution of the continuum and thus produces the largest negative bias, consistent with a negative bias often present in the retrieved refractivity profiles in the moist lower troposphere. By applying a simple linear parameterization of the refractivity structure within and just below the SR layer, an analytical relation between the Abel-retrieved refractivity and a continuum of solutions is derived. Combining the Abel retrieval and the analytical relation with some physical constraints, a novel approach is developed to reconstruct the vertical refractivity structure within and below the SR layer. Numerical simulation studies in this paper have demonstrated the great potential of the reconstruction method to provide a much-improved retrieval in the presence of SR, and the method should greatly enhance the ability to measure the MBL structure globally using the GPS RO technique.

Radio occultation (RO) can provide high-vertical-resolution thermodynamic soundings of the planetary boundary layer (PBL). However, sharp moisture gradients and strong temperature inversion lead to large gradients in refractivity N and often cause ducting. Ducting results in systematically negative RO N biases resulting from a nonunique Abel inversion problem. Using 8 years (2006–13) of Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) RO soundings and collocated European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-I) data, we confirm that the large lower-tropospheric negative N biases are mainly located in the subtropical eastern oceans and we quantify the contribution of ducting for the first time. The ducting-contributed N biases in the northeast Pacific Ocean (160°–110°W; 15°–45°N) are isolated from other sources of N biases using a two-step geometric-optics simulation. Negative bending angle biases in this region are also observed in COSMIC RO soundings. Both the negative refractivity and bending angle biases in COSMIC soundings mainly lie below ~2 km. Such bending angle biases introduce N biases that are in addition to those caused by ducting. Following the increasing PBL height from the southern California coast westward to Hawaii, centers of maxima bending angles and N biases tilt southwestward. In areas where ducting conditions prevail, ducting is the major cause of the RO N biases. Ducting-induced N biases with reference to ERA-I compose over 70% of the total negative N biases near the southern California coast, where strongest ducting conditions prevail, and decrease southwestward to less than 20% near Hawaii.