In this study, the interaction of horizontal axis wind turbine (HAWT) with neutrally stratified atmospheric boundary layer (ABL) and its wake characteristics are investigated. Important wake characteristics of wind turbine such as velocity deficit and turbulence level are analyzed. For this purpose, Unsteady Reynolds-Averaged Navier-Stokes (URANS) using k-e turbulence closure models are performed using commercial Computational Fluid Dynamics (CFD) software called ANSYS FLUENT. Full rotor CFD simulations of the NREL Phase VI wind turbine by virtually placing on a flat surface with different aerodynamic roughness lengths are performed. Discussions on effective modelling of horizontal homogeneity for the undisturbed ABL is included. The influence of inflow ABL's turbulence level in the wake velocity recovery and the ground effect on the wake turbulence intensity (TI) is analyzed. In addition, comparison of rotor aerodynamics of wind turbine in different terrains is performed using pressure coefficient distributions. Finally, the necessity of inclusion of TI recovery in addition to velocity recovery in the wake for the wind farm alignment is discussed. (C) 2018 Elsevier Ltd. All rights reserved.