We investigate the phase diagrams of two-dimensional lattice dipole systems with variable geometry. For bipartite square and triangular lattices with tunable vertical sublattice separation, we find rich phase diagrams featuring a sequence of easy-plane magnetically ordered phases separated by incommensurate spin-wave states. A recent breakthrough in the cooling of dipolar gases in optical lattices [1], following a decade of intensive research [2][3][4][5][6][7][8][9], has opened a door into the earlier inaccessible manybody physics of lattice systems with anisotropic long-range interaction. Although bulk crystalline dipolar systems are abundant in nature [10,11], their experimental investigation has been hindered by the extremely low temperatures required for the observation of ordering transitions [12] and the absence of continuously variable parameters. Recently, artificial twodimensional dipolar systems, such as lithographically created nanomagnet arrays, have been realized [13,14]. However, it is the advent of optical lattices with tunable lattice structure containing ultracold dipolar gases that has created numerous possibilities for studies of previously unexplored phase transitions-both classical and quantum-between ordered and disordered phases of this fundamental many-body system.In this paper, we analyze a series of magnetic phase transitions in a classical dipolar gas in deep optical lattices (square, Fig. 1 and triangular, Fig. 2) obtained from bipartite monolayer lattices by vertical separation, z, of the two sublattices. One way to realize such systems would be loading the ferromagnetic spinor Bose-Einstein minicondensates in the nodes [15] of a deep bilayer optical lattice created with the help of a painted potential technique [16], which would allow for a high degree of control over the shapes of optical lattices and interlayer separation.We find that, upon the variation of z, each system experiences a sequence of easy-plane magnetically ordered phases separated by incommensurate spin-wave states, which could be detected with the help of Bragg diffraction of light [17][18][19]. The phase diagram for the square lattice on the z-T plane is shown in Fig. 1. For sufficiently small separations z a, where a is lattice constant, we reproduce the earlier predicted [20,21] canted antiferromagnetic phase, AFM K , with the ordering vector K. For z > a, we find an antiferromagnetic phase, AFM M , with a larger unit cell and ordering wave vector at the M point of the Brillouin zone of the bipartite lattice. For intermediate interlayer distances, we find a stable ferromagnetic phase (FM), separated from the antiferromagnetic ones by incommensurate spin-wave states (ISW). At the critical temperature T c , all of the ordered phases feature a degeneracy in the orientation of magnetization, characterized in Fig. 1 by angle θ , or θ A and θ B for AFM M . Away from T c , such a degeneracy is lifted, and Fig. 1 shows the optimal orientation of the order parameter for the low-temperature states. The structure of the i...