Abstract. Three dealiasing schemes for large-eddy simulation of turbulent flows are inter-compared for the canonical case of pressure-drive atmospheric boundary-layer type flows. Aliasing errors arise in the multiplication of partial sums, such as those encountered when integrating the non-linear terms of the Navier-Stokes equations in spectral methods (Fourier or polynomial discrete series), and are detrimental to the accuracy of the numerical solution. This is of special relevance when using highorder schemes. In this work, a performance/cost analysis is developed for three well-accepted approaches: the exact 3/2 rule, 5 the Fourier truncation method, and a high order Fourier smoothing method. Tests are performed within a newly developed mixed pseudo-spectral collocation -finite differences large-eddy simulation code, parallelized using a two-dimensional pencil decomposition. The static Smagorinsky eddy-viscosity model with wall damping of the model coefficient is used. A series of simulations are performed at varying resolution and key flow statistics are inter-compared among the considered dealiasing schemes. The numerical results validate the numerical performance predicted by theory when using the Fourier truncation and 10 Fourier smoothing methods. In terms of turbulence statistics, the Fourier Truncation method proves to be over-dissipative when compared against the Fourier Smoothing method and the traditional 3/2-rule, leading to an enhanced horizontal integrated mass flux and to higher dispersive momentum fluxes. Its use in large-eddy simulation of atmospheric boundary-layer type flows is therefore not recommended. Conversely, the Fourier Smoothing method yields accurate flow statistics, comparable to those resulting from the application of the 3/2 rule, with a significant reduction in computational cost, which makes it a convenient 15 alternative for use in the studies related to the atmospheric boundary layer.