Wind turbine wakes in the neutral atmospheric surface layer are simulated with Reynolds-averaged Navier–Stokes (RANS) using an explicit algebraic Reynolds stress model. Contrary to standard two-equation turbulence models, it can predict turbulence anisotropy and complex physical phenomena like secondary motions. For the cases considered, it improves Reynolds stress, turbulence intensity, and velocity deficit predictions, although a more top-hat-shaped profile is observed for the latter.