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 θ y φ 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 φ 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 φ no tienen ningún efecto en absoluto. Esto se debe a que todos los valores de φ representan el mismo punto para θ=0 y θ=π.
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:
![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.