For applications such as minehunting, communications, high frequency time‐reversal, there has been a need for reliable and fast propagation models able to correctly treat scattering at the rough waveguide boundaries. Data and scattering models have shown that as signal frequency or surface roughness increases, the scattering pattern from the ocean bottom and ocean surface evolves from specular to quasi omni‐directional. Ray propagation models are fast but, a priori, they would treat surface bounces as perfectly specular which would not account for the angular spread observed on multipath structures. The idea followed by the authors over the past few years has been to introduce in the Bellhop ray tracing model a reflection kernel including deterministic properties of the sea surface and bottom. Simulations using this kernel display a significant effect of the multipath structure similar to the angular spread observed in in‐situ data. Simulation results are compared with a recent experimental study at 10‐20...