One of the most challenging endeavors of contemporary research is to describe and analyze the dynamic behavior of time series arising from real-world systems. To address the need for analyzing long-range correlations and multifractal properties of multivariate time series, we generalize the multifractal detrended moving average algorithm (MFDMA) to the multivariate case and propose a multivariate MFDMA algorithm (MV-MFDMA). The validity and performance of the proposed algorithm are tested by conducting numerical simulations on synthetic multivariate monofractal and multifractal time series. The MV-MFDMA algorithm is then utilized to analyze raw, seasonally adjusted, and remainder components of five air pollutant time series. Results from all three cases reveal multifractal properties with persistent long-range correlations.