Determining Atmospheric Boundary Layer Height with the Numerical Differentiation Method Using Bending Angle Data from COSMIC

This paper presents a new method to estimate the height of the atmospheric boundary layer (ABL) by using COSMIC radio occultation bending angle (BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data (January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops (large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e., the uncertainties are small; whereas, for non-sharp ABL tops (small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.