In recent years, freak wave/rouge wave has become an important problem in science and engineering. Modulational instability is considered to be an important factor leading to freak wave in the wave evolution of deep water, and Janssen (2003) defined Benjamin-Feir index (BFI) to reflect it. Mori and Janssen (2006) gave the occurrence probability of freak waves based on a weakly non-Gaussian theory, and distribution of wave height is determined by skewness and kurtosis of surface elevation to a considerable extent in deep water. According to observational record, freak wave has not only been found in deep water in the ocean, but also been observed in shallow water and coastal areas. In the process of water wave entering continental shelf, water depth is changing with mild slope after a long distance propagation. This study focus on investigating how water depth affect skewness and kurtosis in the high order nonlinear wave evolution from deep water to finite water depth in two-dimension.Recorded Presentation from the vICCE (YouTube Link): https://youtu.be/a8LiJvXWRrw