Current water resource recovery facility (WRRF) models only consider local concentration variations caused by inadequate mixing to a very limited extent, which often leads to a need for (rigorous) calibration. The main objective of this study is to visualize local impacts of mixing by developing an integrated hydrodynamic-biokinetic model for an aeration compartment of a full-scale WRRF. Such a model is able to predict local variations in concentrations and thus allows judging their importance at a process level. In order to achieve this, full-scale hydrodynamics have been simulated using computational fluid dynamics (CFD) through a detailed description of the gas and liquid phases and validated experimentally. In a second step, full ASM1 biokinetic model was integrated with the CFD model to account for the impact of mixing at the process level. The integrated model was subsequently used to evaluate effects of changing influent and aeration flows on process performance. Regions of poor mixing resulting in non-uniform substrate distributions were observed even in areas commonly assumed to be well-mixed. The concept of concentration distribution plots was introduced to quantify and clearly present spatial variations in local process concentrations. Moreover, the results of the CFD-biokinetic model were concisely compared with a conventional tanks-in-series (TIS) approach. It was found that TIS model needs calibration and a single parameter set does not suffice to describe the system under both dry and wet weather conditions. Finally, it was concluded that local mixing conditions have significant consequences in terms of optimal sensor location, control system design and process evaluation.