Numerical approximations to the solution of a singularly perturbed elliptic convection-diffusion problem in two space dimensions are generated using a monotone finite difference operator on a tensor product of piecewise-uniform Shishkin meshes. The bilinear interpolants of these numerical approximations are parameter-uniformly convergent to the solution of the continuous problem, in the pointwise maximum norm. In this article, discrete approximations to the first derivatives of the solution are shown to be globally first-order (up to logarithmic factors) uniformly convergent, when the errors are scaled within the analytical layers of the continuous problem. Numerical results are presented to illustrate the theoretical error bounds established in an appropriated weighted C 1 -norm.