The aim of this paper is to extend the existing theory of second-order self-similar processes as defined by Cox (1984) from the univariate case to higher dimensions. Multivariate self-similar processes defined in terms of second-order theory for stationary time series can be used as models for long-range dependent observations when the marginal observations are long-range dependent. An interesting question concerns the correlation structure within the processes when the marginal processes are correlated. We show that the self-similarity requirement, as defined in this article, implies a cross-correlation structure similar to that for the marginal processes. This occurs both in the time domain and in the frequency domain. This fact can be used to obtain generalized least squares estimates for the long-range dependence parameters. We discuss some difficulties concerning estimation based on simulations.