De Haan and Karandikar (1989) [7] introduced generalized Ornstein-Uhlenbeck processes as one-dimensional processes (V t) t<0 which are basically characterized by the fact that for each h>0 the equidistantly sampled process (V nh) n∈N0 satisfies the random recurrence equation V nh=A (n-1)h,nhV (n-1)h+B( n-1)h,nh, n∈N, where ( A( n-1)h,nh,B( n-1)h,nh) n∈N is an i.i.d. sequence with positive A 0,h for each h>0. We generalize this concept to a multivariate setting and use it to define multivariate generalized Ornstein-Uhlenbeck (MGOU) processes which occur to be characterized by a starting random variable and some Lévy process (X,Y) in R m× m×R m. The stochastic differential equation an MGOU process satisfies is also derived. We further study invariant subspaces and irreducibility of the models generated by MGOU processes and use this to give necessary and sufficient conditions for the existence of strictly stationary MGOU processes under some extra conditions.