There are many effects that take part in forming the signal distribution. The multitude of references are popped up when searching with keywords "signal phase uniform distribution", and many are about the communication channel effects -- for the most part, about the Rayleigh fading in the channel. I think the channel impulse response is the subject of a few discussions that you mention in your question. Because the phase distribution of the Gaussian process per se is analyzed in multiple textbooks and courses of instruction, I am going to consider the contribution of the channel impulse response -- in particular, the scenarios where it can result in a uniform phase distribution.
The impulse response of a wireless channel is formed by the transmission medium. The three key components of the channel response are path loss, shadowing, and multipath. Multipath propagation and shadowing lead to fading. Similar to how stochastic processes in general can be of the Wide Sense Stationary kind, the impulse response of channels can often be described with the Wide Sense Stationary Uncorrelated Scattering model.
In WSSUS, if there is no dominant contribution of a unique propagation path (typically, a line-of-sight path) over the other paths (this is the Rayleigh fading scenario), the channel impulse response (CIR) reduces to a time-variant complex channel coefficient which is a sum over all constituent paths
$$
c(t) = \sum_i{a_i \exp(j\phi _i (t))}
$$
For the quadrature components of the channel coefficient that are i.i.d. Gaussian distributed variables, the probability density functions can be computed in a closed form
$$
p(r) = {\frac {r} {\sigma^2}} \exp\left(-{\frac {r^2} {2\sigma^2}}\right) \\
p(\phi) = {\frac {1} {2\pi}}
$$
If a contribution of the line-of-sight propagation path is strong as compared to the complex Gaussian component from the other propagation paths (the Rician fading scenario), this contribution persists in CIR
$$
c(t) = a_0 + \sum_i{a_i \exp(j\phi _i (t))}
$$
After some tedious computation, for the Rician magnitude pdf one can arrive at
$$
p(r) = {\frac {r} {\sigma^2}} \exp{\left(-{\frac {r^2 + a_0^2} {2\sigma^2}}\right)} I_0\left({\frac {r·a_0} {\sigma^2}}\right)
$$
where $I_0(·)$ is the modified Bessel function of the first kind.
I found this formula in multiple sources and also myself traced the derivation to some degree. The formula for the Rician phase pdf $p(\phi)$ depends on $\phi$ (non-uniform distribution), is very complicated, and I am not sure that it exists in the closed form. With MATLAB' Communications Toolbox you can implement a Rician fading channel object for modeling real-world phenomena in wireless communications, though.
The CIR phase distribution is uniform for the Rayleigh fading scenario and non-uniform for the Rician fading scenario. I emphasize that these formulas per se do not describe the phase distribution of a signal received from a multipath communication channel: the received signal is a convolution of CIR and a source signal.
You learned that the phase distribution for a AWGN signal is uniform, and the derivation for the WSS signal is trivial. Now, the CIR analysis for WSSUS fading channels is a hot topic and there is a field for rewarding research. So you are on the right track in your studies.