In this paper, a three-dimensional (3D) cascaded lattice Boltzmann method (CLBM) is implemented to simulate the liquid-vapor phasechange process. The multiphase flow field is solved by incorporating the pseudopotential multiphase model into an improved CLBM, the temperature field is solved by the finite difference method, and the two fields are coupled via a non-ideal equation of state. Through numerical simulations of several canonical problems, it is verified that the proposed phase-change CLBM is applicable for both the isothermal multiphase flow and the liquid-vapor phase-change process. Using the developed method, a complete 3D pool boiling process with up to hundreds of spontaneously generated bubbles is simulated, faithfully reproducing the nucleate boiling, transition boiling, and film boiling regimes. It is shown that the critical heat flux predicted by the 3D simulations agrees better with the established theories and correlation equations than that obtained by two-dimensional simulations. Furthermore, it is found that with the increase in the wall superheats, the bubble footprint area distribution changes from an exponential distribution to a power-law distribution, in agreement with experimental observations. In addition, insights into the instantaneous and time-averaged characteristics of the first two largest bubble footprints are obtained.