The present article is in the continuation of our previous work [26], where we have improved the order of approximation of ??Bernstein P?lt?nea operators. In the given note, we study the bivariate extension of first order modification of these operators and their approximation properties such as convergence, error of approximation in terms of complete and partial modulus of continuity and their asymptotic formula. We present numerical examples to show the convergence of functions of two variables with the help of MATLsoftware. Also, we construct the GBS operators associated to the bivariate extension and present their approximation behavior.