Abstract: In this paper, a 3-dimensional fine-mesh nonhydrostatic PBL model has been developed and used for numerical simulation on the PBL structure and turbulence characteristics over complex underlying regions. Some results from this model were in agreement with those from observation. Energy closure scheme(E-∈) was used. The initial fields used in this model were from observed data. Some numerical experiments were also done and effects on PBL simulation,under the conditions of using nonhydrostatic model or not and different closure schemes,were obtained.The results indicated that abrupt terrain and irregular coasts had more effects on the local wind fields and turbulence fields, and it was necessary to use nonhydrostatic model for complex underlying regions, and energy closure scheme could simulate turbulence diffusoin better than K-closure scheme in the complex airflow.