Summary In this paper we are continuing our work [6], concerning a-priori error estimates for the velocity tracking of two-dimensional evolutionary Navier-Stokes flows. The controls are of distributed type, and subject to point-wise control constraints. The discretization scheme of the state and adjoint equations is based on a discontinuous timestepping scheme (in time) combined with conforming finite elements (in space) for the velocity and pressure. Provided that the time and space discretization parameters, τ and h respectively, satisfy τ ≤ Ch 2 , error estimates of order O(h 2 ) and O(h 3 2 − 2 p ) with p > 3 depending on the regularity of the target and the initial velocity, are proved for the difference between the locally optimal controls and their discrete approximations, when the controls are discretized by the variational discretization approach and by using piecewise-linear functions in space respectively. Both results are based on new duality arguments for the evolutionary Navier-Stokes equations.