A statistical model based on the geometric optics approximation for predicting wave scattering from random rough surfaces is developed. This model uses statistical concepts. rather than bundle tracing employed in the geometric optics approximation. This model includes first-, second- and higher-order scatterings. Shadowing effects are included, and the model can be expressed in a closed form. For one-dimensional random rough surfaces, the statistical model is quantified by comparisons with the exact electromagnetic theory calculations, and the results are in good agreement for surfaces with a surface slope less than unity. For two-dimensional random rough surfaces, the approximate results are compared with the existing experimental findings and electromagnetic theory calculations. The statistical method is a computationally inexpensive wave scattering approach for surfaces with a surface slope less than 0.5, and the method can approximate the magnitude and the trend for surfaces with a surface slope greater than 0.5.