Respuesta corta:
Imagina la situación en un sistema de coordenadas esféricas, con el punto de reflexión en la superficie del agua en el origen y el eje $z$ paralelo a la superficie del agua quieta y apuntando en la dirección azimutal del sol.
A medida que el sol desciende, tanto los rayos de luz entrantes como los reflejados se acercan más a la paralelidad con el eje $z$. Imagina rotar la superficie del agua en pequeñas cantidades en las direcciones $\theta$ y $\varphi$ respectivamente y piensa en los efectos que eso tiene en la dirección de la luz reflejada.
A medida que los rayos de luz se acercan a la paralelidad con el eje $z$, los efectos de las rotaciones en la dirección de $\varphi$ tienen efectos cada vez más pequeños en la dirección de la luz reflejada. En el caso extremo del sol en el horizonte, las rotaciones en $\varphi$ no tienen ningún efecto en absoluto. Esto se debe a que todos los valores de $\varphi$ representan el mismo punto para $\theta=0$ y $\theta=\pi$.
Versión más larga de la misma respuesta (más matemática):
Encuentro útil pensar en esto en un sistema de coordenadas esféricas:
El sistema de coordenadas se elige de manera que el eje $x$ sea ortogonal a la superficie del agua quieta, y el eje $y$ sea ortogonal a la luz solar entrante. Imagina la superficie de agua reflectante en el origen, y 3 vectores en la esfera unitaria. El primer vector es la normal de la superficie reflectante, que escribiremos en esta forma:
$$ \boldsymbol{r}_n = (1,\, \pi/2 + \delta\theta_n,\, \delta\varphi_n) $$
Basicamente está apuntando en la dirección $x$, excepto por pequeñas variaciones debido a las olas representadas por $\delta\theta_n$ y $\delta\varphi_n$. El segundo vector apunta hacia la luz solar entrante:
$$ \boldsymbol{r}_i = (1,\, \theta_i,\, 0) $$
Y el tercero apunta en la dirección de la luz reflejada:
$$ \boldsymbol{r}_r = (1,\, \theta_r,\, \varphi_r) $$
Debido a la ley de la reflexión, obtenemos las siguientes relaciones entre los ángulos:
$$ \begin{align} (\pi/2 + \delta\theta_n) - \theta_i & = \theta_r - (\pi/2 + \delta\theta_n)\\ \Rightarrow\,\,\,\, \theta_r & = \pi - \theta_i + 2\delta\theta_n \end{align} $$
$$ \begin{align} \delta\varphi_n - 0 &= \varphi_r - \delta\varphi_n\\ \Rightarrow\,\,\,\, \varphi_r &= 2\delta\varphi_n \end{align} $$
Para ver la magnitud de la influencia de los pequeños cambios en la orientación de la superficie reflectante ($\delta\theta_n$ y $\delta\varphi_n$) en la dirección de la luz reflejada, podemos usar la fórmula para desplazamientos infinitesimales en sistemas de coordenadas esféricas:
$$ \mathrm{d}\boldsymbol{r} = \mathrm{d}r\,\boldsymbol{\hat{r}} + r\,\mathrm{d}\theta\,\boldsymbol{\hat{\theta}} + r \sin{\theta}\,\mathrm{d}\varphi\,\boldsymbol{\hat{\varphi}} $$
Ppara $\boldsymbol{r}_r$ obtenemos:
$$ \begin{align} \mathrm{d}\boldsymbol{r}_r &= \mathrm{d}\theta_r\,\boldsymbol{\hat{\theta}} + \sin{\theta_r}\,\mathrm{d}\varphi_r\,\boldsymbol{\hat{\varphi}}\\ &= 2\delta\theta_n\,\boldsymbol{\hat{\theta}} + \sin{(\pi - \theta_i + 2\delta\theta_n)}\,2\delta\varphi_n\,\boldsymbol{\hat{\varphi}}\\ &= 2\delta\theta_n\,\boldsymbol{\hat{\theta}} + 2\sin{(\theta_i - 2\delta\theta_n)}\,\delta\varphi_n\,\boldsymbol{\hat{\varphi}}\\ &\approx 2\delta\theta_n\,\boldsymbol{\hat{\theta}} + 2\sin{\theta_i}\,\delta\varphi_n\,\boldsymbol{\hat{\varphi}}\\ \end{align} $$
En el lado derecho, vemos que las variaciones de $\delta\theta_n$ se amplifican en la dirección de la luz reflejada por un factor de $2$, mientras que las variaciones de $\delta\varphi_n$ se amplifican por un factor de $2\sin{\theta_i}$.
Dado que $\theta_i$ es básicamente el ángulo de elevación del sol, esto significa que a medida que el sol se acerca al horizonte ($\theta_i = 0$), las variaciones de la superficie en una dirección comienzan a ser fuertemente atenuadas, mientras que en la otra dirección el factor de amplificación de $2$ permanece constante.