Solution (source code)

= Solution

Illuminate a surface window of length $D$, and define $\widehat h_D(Q)=\int_{-D/2}^{D/2}h(x)e^{-iQx}dx$. In the far field put $(x,z)=r(\sin\theta_s,\cos\theta_s)$, with $|\theta_s|<\pi/2$ away from grazing. The spectral phase $\xi\sin\theta_s+\beta(\xi)\cos\theta_s$ has stationary point $\xi_s=k\sin\theta_s$, value $k$, and second <derivative> $-1/(k\cos^2\theta_s)$. <Stationary phase> therefore gives
$$
\psi_s^{[1]}\sim iq_iA_i\sqrt{\frac{2k}{\pi r}}\cos\theta_s\,
\widehat h_D(Q)e^{ikr-i\pi/4},\qquad Q=k\sin\theta_s-p_i.
$$
Let $C_h(s)=\langle h(x+s)h(x)\rangle$ and $S_h(Q)=\int C_h(s)e^{-iQs}ds$. Statistical <stationarity> gives exactly
$$
\left\langle|\widehat h_D(Q)|^2\right\rangle
=\int_{-D}^D(D-|s|)C_h(s)e^{-iQs}ds.
$$
If $C_h$ is integrable, division by $D$ and <dominated convergence> yield $\langle|\widehat h_D(Q)|^2\rangle/D\to S_h(Q)$. Thus the leading <rough-surface far-field intensity per illuminated length> is
$$
\boxed{\frac{\langle|\psi_s^{[1]}(r,\theta_s)|^2\rangle}{D}
\sim\frac{2kq_i^2|A_i|^2\cos^2\theta_s}{\pi r}\,S_h(k\sin\theta_s-p_i).}
$$
This uses the field-intensity convention $I=|\psi|^2$; if $\psi$ is acoustic pressure, multiply by $1/(2\rho c)$ for physical far-field energy flux. The <variance> normalization is $\sigma^2=C_h(0)=(2\pi)^{-1}\int S_h(Q)dQ$. Consequently the diffuse intensity is of order $\sigma^2$ and samples the surface spectrum at the tangential momentum transfer. Its value is not determined by the r.m.s. height alone.

The finite window must first be in its far field, for example $r\gg kD^2$, before interpreting the large-window intensity density. An infinitely long stationary surface has infinite illuminated power and no finite total far-field intensity at a fixed point. For nonintegrable correlations the analogous result is a spectral measure, possibly with atoms. Zero mean height makes the linear rough mean field zero, but does not make its mean intensity zero; the flat specular reflection remains separate.