In this paper, we consider a bivariate extension of blending type approximation by Lupaş-Durrmeyer type operators involving Pólya Distribution. We illustrate the convergence rate of these type operators using Peetre's K-functional, modulus of smoothness and for functions in a Lipschitz type space.