The spinless Falicov-Kimball model on a two-dimensional square lattice is studied using the method of restricted phase diagrams constructed in the grand canonical ensemble. The results are compared with the one-dimensional model. Although the two-dimensional phase diagrams are more complex, with several distinct families of ion configurations occurring as ground states, there are surprising similarities with the one-dimensional case. Within each family of configurations, the ground states form a devil`s staircase structure and the configurations are constructed according to a composition rule identical to that in one dimension. It is also found that, as in one dimension, segregation occurs in the non-neutral model for large ion-electron interaction strength. Some features of the phase diagrams are understood by examining the effective two-body ion interaction.