Abstract: When considering the modelling of small particle dispersion in the lower part of the Atmospheric Boundary Layer (ABL) using Reynolds Averaged Navier Stokes simulations, the particle paths depend on the velocity profile and on the turbulence kinetic energy, from which the fluctuating velocity components are derived to predict turbulent dispersion. It is therefore important to correctly reproduce the ABL, both for the velocity profile and the turbulence kinetic energy profile. For RANS simulations with the standard kå model, Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the kå turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 4647, 145153.) proposed a set of boundary conditions which result in horizontally homogeneous profiles. The drawback of this method is that it assumes a constant profile of turbulence kinetic energy, which is not always consistent with field or wind tunnel measurements. Therefore, a method was developed which allows the modelling of a horizontally homogeneous turbulence kinetic energy profile that is varying with height. By comparing simulations performed with the proposed method to simulations performed with the boundary conditions described by Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the kå turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 4647, 145153.), the influence of the turbulence kinetic energy on the dispersion of small particles over flat terrain is quantified.