NCEP-FNL datasets are used as the initial and boundary fields of the WRF model. Six planetary boundary layer parameterization schemes (PBLPS) are applied in the model for Xinjiang region with 10-km horizontal resolution. The spatial distribution and temporal evolution of the meteorological elements are analyzed. The analysis results show the following aspects: 1) The WRF model with six PBLPS can simulate the seasonal circulation of monthly precipitation and the spatial pattern of annual and rainy season precipitation. 2) For the Xinjiang region, the deviation of rainy season precipitation between the simulation conducted with the Grenier-Bretherton-McCaa (GBM) scheme and the observations is within ±30%. For the Tianshan area, the deviation of annual precipitation between the simulation conducted with the Bougeault-Lacarrere (BouLac) scheme and the observations is －19.13%. The TS scores of moderate and heavy rains are 0.37 and 0.33, respectively, in the test results of daily precipitation simulated with the GBM scheme. For the different types of underlying surfaces in the Tianshan area, the day and night precipitation can be well simulated by the model with the GBM scheme with the deviation of precipitation within 5 mm during long precipitation days. 3) The WRF model with the BouLac scheme can simulate the annual spatial and temporal distribution characteristics of annual precipitation in the Tianshan area, and the rainy seasonal precipitation can be well simulated by the model with the GBM scheme in Xinjiang. Therefore, PLBS with the WRF model in the Xinjiang region should be considered.