The article presents a study which investigates the global distribution of asthenosphere-lithosphere boundary. The study used global maps, numerical simulation method, and half-space cooling (HSC) and the Plate models the finite half-space (FHS) model used. Results show that there is an improvement on the theoretical values of the boundary that fit to the observational data for bathymetry and heat flow.