Lagrangian particle dispersion models (LPDMs) are frequently used for regional-scale inversions of greenhouse gas emissions. However, the turbulence parameterizations used in these models were developed for coarse resolution grids, hence, when moving to the kilometre-scale the validity of these descriptions should be questioned. Here, we analyze the influence of the turbulence parameterization employed in the LPDM FLEXPART-COSMO model. Comparisons of the turbulence kinetic energy between the turbulence schemes of FLEXPART-COSMO and the underlying Eulerian model COSMO suggest that the dispersion in FLEXPART-COSMO suffers from a double-counting of turbulent elements when run at a high resolution of $$1 \times 1 \,\hbox {km}^2$$
1
×
1
km
2
. Such turbulent elements are represented in both COSMO, by the resolved grid-scale winds, and FLEXPART, by its stochastic parameterizations. Therefore, we developed a new parametrization for the variations of the winds and the Lagrangian time scales in FLEXPART in order to harmonize the amount of turbulence present in both models. In a case study for a power plant plume, the new scheme results in improved plume representation when compared with in situ flight observations and with a tracer transported in COSMO. Further in-depth validation of the LPDM against methane observations at a tall tower site in Switzerland shows that the model’s ability to predict the observed tracer variability and concentration at different heights above ground is considerably enhanced using the updated turbulence description. The high-resolution simulations result in a more realistic and pronounced diurnal cycle of the tracer concentration peaks and overall improved correlation with observations when compared to previously used coarser resolution simulations (at 7 km $$\times $$
×
7 km). Our results indicate that the stochastic turbulence schemes of LPDMs, developed in the past for coarse resolution models, should be revisited to include a resolution dependency and resolve only the part of the turbulence spectrum that is a subgrid process at each different mesh size. Although our new scheme is specific to COSMO simulations at $$1 \times 1 \,\hbox {km}^2$$
1
×
1
km
2
resolution, the methodology for deriving the scheme can easily be applied to different resolutions and other regional models.