Bivariate polynomial interpolation on the square at new nodal sets