The planetary boundary layer (PBL) parameterization schemes play a critical role in the weather and climate models, while they describe physical processes associated with the momentum, heat and humidity exchange between land surface and atmosphere. The sensitivity of boundary layer variables to eight PBL parameterization schemes (three non-local and five local closure schemes), available in the Weather Research and Forecasting (WRF) model, is evaluated over the central Tibetan Plateau with field measurements of the Third Tibetan Plateau atmospheric scientific experiment (TIPEX III) in July. Model results showed acceptable behavior, but no particular scheme produced the best performances for all observation stations and meteorological parameters. All PBL schemes underestimated the surface temperature over the central Tibetan Plateau. The BouLac scheme showed the minimum cold bias of the surface temperature. For the surface energy budget components, it was found that the sensible heat flux and the downward longwave radiation were the main factors causing the lower surface temperature. The sub-grid scale gravity wave drag was added to reduce biases result from unresolved topography over the central Tibetan Plateau. It led to smaller cold bias, causing warmer lower-tropospheric temperature, smaller water vapor content and higher PBL height. The modified model results show more close to the observation.