We study the long time Schrödinger evolution of Lagrangian states $f_h$ on a compact Riemannian manifold $(X,g)$ of negative sectional curvature. We consider two models of semiclassical random Schrödinger operators $P_h^α=-h^2Δ_g +h^αQ_ω$, $0<α\leq 1$, where the semiclassical Laplace-Beltrami operator $-h^2Δ_g$ on $X$ is subject to a small random perturbation $h^αQ_ω$ given by either a random potential or a random pseudo-differential operator. Here, the potential or the symbol of $Q_ω$ is bounded, but oscillates and decorrelates at scale $h^β$, $0< β< \frac{1}{2}$. We prove a quantitative result that, under appropriate conditions on $α,β$, in probability with respect to $ω$ the long time propagation $$\mathrm{e}^{\frac{i}{h}t_h P_h^α} f_h, \quad o(|\log h|)=t_h\to\infty, ~~h\to 0,$$ rescaled to the local scale of $h$ around a uniformly at random chosen point $x_0$ on $X$, converges in law to an isotropic stationary monochromatic Gaussian field -- the Berry Gaussian field. We also provide and $ω$-almost sure version of this convergence along sufficiently fast decaying subsequences $h_j\to 0$.