The mapped gaussian process (MGP) force-field used to elucidate the surface alloying of Cu-Zn. The force-field is made based on first-principles data by using machine-learning technique called Gaussian Process as implemented in FLARE package (https://github.com/mir-group/flare). Active and on-the-fly learning were employed to build the database efficiently. The simulation reveals atomistic details of the alloying process, i.e., the incorporation of deposited Zn adatoms to the Cu substrate. The surface alloying is found to start at upper and lower terraces near the step edge, which emphasize the role of steps and kinks in the alloying. The incorporation of Zn at the middle terrace was found at the later stage of the simulation.