 Solution

Take $e^{-i\omega t}$ time dependence and measure angles from the upward normal. Write $p_i=k\sin {\theta }_i$, $q_i=k\cos {\theta }_i>0$ and choose incident amplitude $A_i$. Separate the flat Dirichlet reflection from the rough correction:Taylor expansion of the boundary condition about $z=0$ givesSince ${\partial }_z{\psi }^{[0]}(x,0)=-2i{q}_i{A}_i{e}^{i{p}_{i}x}$, the induced trace isUse $\hat{h}(Q)=\int h(x)e^{-iQx}dx$ and $\beta (\xi )=\sqrt{k^2-{\xi }^2}$ with the outgoing branch: $\beta >0$ for $|\xi |<k$ and $\mathrm{Im}\beta >0$ outside. The outgoing angular spectrum solves the Helmholtz equation with that trace, giving the first-order rough-surface scattered fieldIf “scattered field” includes the specular reflection, add $-A_i{e}^{i{p}_{i}x+i{q}_{i}z}$ to this rough correction. Fourier transforms of an infinite random surface are generalized objects; finite-window transforms give an ordinary integral and are also needed for its far-field intensity. The expansion assumes a fixed regular surface shape scaled to small height. Highly oscillatory profiles can require additional control of slopes and spectral moments; small height alone is not a uniform error estimate for every profile.