 Solution

For a partition of $[a,b]$, every Riemann sum $S_n={\sum }_j{B}_{t_j}\Delta t_j$ is a centered Gaussian random variable, withPath continuity gives $S_n\to {\int }_a^b{B}_{s}ds$ almost surely, and the covariance bound $\mathbb{E}|B_s-B_t{|}^2=|s-t|$ also gives convergence in $L^2$. Hence the limit is Gaussian, centered, and passage to the limit in the displayed sums gives variance