We introduce and study the sequence of bivariate Generalized Bernstein operators {Bm,s}m,s, m, s ∈ N, Bm,s = I − (I − Bm) s , B i m = Bm(B i−1 m), where Bm is the bivariate Bernstein operator. These operators generalize the ones introduced and studied independently in the univariate case by Mastroianni and Occorsio [Rend. Accad. Sci. Fis. Mat. Napoli 44 (4) (1977), 151-169] and by Micchelli [J. Approx. Theory 8 (1973), 1-18] (see also Felbecker [Manuscripta Math. 29 (1979), 229-246]). As well as in the one-dimesional case, for m fixed the sequence {Bm,s(f)}s can be successfully employed in order to approximate "very smooth" functions f by reusing the same data points