The purpose of this paper is to define normal Gaussian variables in the setting of hyperbolic probabilities, and introduce an associated Brownian motion, when both the index and the values of the process lie in the real algebra $\mathbb{H}$ of hyperbolic numbers. In Hida's white noise space, we construct two probability measures (say $P_1$ and $P_2$), and associate to them two families of $N(0,1)$ variables $(Z_n)_{n\in\mathbb N_0}$ (independent with respect to $P_1$) and $(W_n)_{n\in\mathbb N_0}$ (independent with respect to $P_2$). An important feature is that the $Z_n$ and $W_m$ need not be mutually independent either with respect to $P_1$ or $P_2$. An hyperbolic normal Gaussian variable is constructed (in non-degenerate cases) from two classical Gaussian variables and the hyperbolic Brownian motion is, in general, composed from two copies of the classical Brownian motion. Using the associated Gelfand triples we also compute the derivative of the hyperbolic Brownian motion as a stochastic distribution. The argument extends to the $\mathbb{H}$-valued fractional Brownian motion, and more generally to a wide family of $\mathbb{H}$-valued stationary-increment second order processes.