Despite several reports on the surface phase transformations from a layered to a disordered spinel and a rock-salt structure at the surface of the Ni-rich cathodes, the precise structures and compositions of these surface phases are unknown. The phenomenon, in itself, is complex and involves the participation of several contributing factors. Of these factors, transition metal (TM) ion migration toward the interior of the particle and hence formation of TM-densified surface layers, triggered by oxygen loss, is thermodynamically probable. Here, we simulate the thermodynamic phase equilibria as a function of TM ion content in the cathode material in the context of lithium nickel oxides, using a combined approach of first-principles density functional calculations, the cluster expansion method, and grand canonical Monte Carlo simulations. We developed a unified lattice Hamiltonian that accommodates not only rock-salt like structures but also topologically different spinel-like structures. Also, our model provides a foundation to investigate metastable cation compositions and kinetics of the phase transformations. Our investigations predict the existence of several Ni-rich phases that were, to date, unknown in the scientific literature. Our simulated phase diagrams at finite temperature show a very low solubility range of the prototype spinel phase. We find a partially disordered spinel-like phase with far greater solubility that is expected to show very different Li diffusivity compared to that of the prototype spinel structure.