In ITER, the helium (He) impurity produced by the deuterium-tritium reaction will bombard the tungsten (W) divertor armor at the strike points. Consequently, strong interaction occurs therein that both impact the performance of the plasmas and the lifetime of the divertor. Despite an ever-increasing understanding of this interaction, some experimental phenomena remain mysterious, especially the formation of orientation-dependent surface morphologies. Here, we combine multiscale experimental characterization and theoretical models to shed new light on this problem. After low-energy He plasma exposure in a linear plasma generator, the polycrystalline W surface developed various morphologies. Through electron backscatter diffraction analysis, we found that the {111} grains developed cube-corner structures, the {110} grains developed ripple structures, whereas the {100} grains remained smooth. Then, electron-transparent lamellae were extracted from such grains to observe the subsurface He bubbles by transmission electron microscopy. The volume density, size distribution, and depth range of the He bubbles weakly depend on the crystallographic orientation, suggesting that the migration of W atoms causes the morphology variety. Accordingly, we proposed a two-stage formation mechanism. First, W atoms generated by over-pressurized He bubbles glide on the slip plane and in the slip direction to reach the surface, forming characteristic patterns that are enclosed by the slip traces. Second, morphological instability drives the evolution of the surface patterns, in which the initial surface structure and surface self-diffusion kinetics mediate. The proposed mechanism has been incorporated into a topographical instability model to enable asemi-quantitative analysis. The obtained new insights are valuable to the impurity control of the core plasmas and the lifetime analysis of the divertor for ITER.