Characteristics of turbulence anisotropy in flow over two-dimensional rigid dunes are analysed. The Reynolds stress anisotropy is envisaged from the perspective of the stress ellipsoid shape. The spatial evolutions of the anisotropic invariant map (AIM), anisotropic invariant function, eigenvalues of the scaled Reynolds stress tensor and eccentricities of the stress ellipsoid are investigated at various streamwise distances along the vertical. The data plots reveal that the oblate spheroid axisymmetric turbulence appears near the top of the crest, whereas the prolate spheroid axisymmetric turbulence dominates near the free surface. At the dune trough, the axisymmetric contraction to the oblate spheroid diminishes, as the vertical distance below the crest increases. At the reattachment point and one-third of the stoss-side, the oblate spheroid axisymmetric turbulence formed below the crest appears to be more contracted, as the vertical distance increases. The AIMs suggest that the turbulence anisotropy up to edge of the boundary layer follows a looping pattern. As the streamwise distance increases, the turbulence anisotropy at the edge of the boundary layer approaches the plane-strain limit up to two-thirds of the stoss-side, intersecting the plane-strain limit at the top of the crest and thereafter moving towards the oblate spheroid axisymmetric turbulence. © 2020 The Author(s).