A novel class of bivariate \(\alpha \) -fractal functions is constructed on rectangular grids in this paper, with sufficient conditions established for these \(\alpha \) -fractal functions to become \(\alpha \) -fractal interpolation functions ( \(\alpha \) -FIFs). The box-counting dimension of the resulting \(\alpha \) -FIFs is subsequently estimated. Three numerical examples are provided to illustrate the effectiveness and practical applicability of the proposed construction method. Finally, it is demonstrated that the Riemann-Liouville fractional integrals of the constructed bivariate \(\alpha \) -fractal functions remain \(\alpha \) -fractal functions.