In this paper, it aims to model wind speed time series at multiple sites. The five-parameter Johnson distribution is deployed to relate the wind speed at each site to a Gaussian time series, and the resultant m- dimensional Gaussian stochastic vector process $boldsymbol{Z}(t)$