Employing molecular dynamics simulation to study the thermodynamic properties of liquid metal, Nickel, we find the existence of surface-induced layering at 1220K (this temperature is very close to the phase transition temperature between solid and supercooled liquid of this system, 1190K). From the density profile at 1220K, the width of the outermost layer is found to be about 1.89 Å (this is a typical value for many models). On the other hand, the close-packed structure in these layers is recognized by analysing radial distribution function g(r). In addition, based on the distribution of coordinated numbers, we find out that the particle density in layering structure space is higher than that of the remaining part of the system. Finally, the 3D visualizations of roughness surface is taken at 1220 K in order to highlight the situation.