In this paper we introduce new bivariate Bernstein type operators BnM,i(f; x, y), i = 1, 2, 3. The rates of approximation by these operators are calculated and it is shown that the errors are significantly smaller than those of ordinary bivariate Bernstein operators for sufficiently smooth functions.