Given a bivariate function f defined in a rectangular domain Ω, we approximate it by a C1 quadratic spline quasi-interpolant (QI) and we take partial derivatives of this QI as approximations to those of f. We give error estimates and asymptotic expansions for these approximations. We also...