This study deals with the development of statistical modeling for water wave surface elevation by using a method that combines a dynamic solution with random process statistics. Ocean wave data taken from four NOAA (National Oceanic and Atmospheric Administration) buoys moored in the northeast Pacific were used to validate the model. The results indicated that the nonlinear probability density distribution of ocean wave surface elevation derived from the model described the measurements much better than Gaussian distribution and Longuet-Higgins distribution.