I read somewhere that Brownian motion can be defined to be a stochastic process with continuous paths, independent and stationary increments such that the mean of the process is zero at all times and the variance at $t$ is equal to $t$. So Gaussanity is not mentioned at all.
The claim is that any process that fits this definition is a Gaussian process. I am trying to understand why this is true. I thought about using Levy's characterization of BM but for that I need to show that a process with the definition above has a certain quadratic variation process. This requires me to make some statement on the fourth moments of the process but I don't see how that would follow from the definition above. Also, I want to use a more elementary method to show normality. Any hints/suggestions?
Recall that a stochastic process $(X_t)_{t \geq 0}$ with càdlàg sample paths is called a Lévy process if it has stationary and independent increments. The result which you mention can be formulated as follows:
Let me remark that the assumptions on $(X_t)_{t \geq 0}$ can be weakened; in fact, the assumptions $\mathbb{E}(X_t^2)=t$ and $\mathbb{E}(X_t)=0$ can be droped. Any Lévy process with continuous sample paths is Gaussian.
There are several ways to prove the above result:
References:
(1) Schilling & Partzsch: Brownian Motion - An Introduction to Stochastic Processes. De Gruyter. (2nd edition)
(2) Schilling: An Introduciton to Lévy and Feller Proceses. Online Version