Outlier detection for PIV velocity fields is still nowadays an active field of research. In the last decades, several image pre-processing and processing algorithms have been developed aiming at increasing the dynamic velocity range of PIV measurements and reducing the measurement uncertainty. Nevertheless, PIV velocity fields are still often characterised by the presence of outliers, which potentially hamper the correct interpretation of the flow physics and negatively affect the evaluation of the flow statistics. The outlier detection strategies presented in literature are mainly based on the statistical analysis of the velocity vector with its immediate neighbour. Most of these algorithms have been demonstrated to be effective for instantaneous flow fields, where the errors associated with the outliers are order of magnitude larger than the expected measurement uncertainty. However, these approaches are not as effective for the flow statistics, where the outliers yield errors of the same order of the measurement uncertainty. To overcome this limitation, this paper proposed an outlier detection approach based on the agreement of the flow statistics to the constitutive equations, more specifically to the turbulent kinetic energy (TKE) transport equation. The focus is posed on the ratio between the local advection terms of TKE and a robust estimation of the TKE production along the local streamline. It is demonstrated that, in presence of outliers, the proposed principle yields a clear separation between the correct and the erroneous vectors. In order to assess the performance of the proposed principle, three different test cases are considered. For all of them, the results are compared with a reference outlier detection methodology, namely the universal outlier detection method proposed by Westerweel and Scarano (2005). The proposed turbulence transport-based approach exhibits higher performance in terms of percentage of outliers correctly identified in the flow statistics.