This study focuses on modeling air quality with an adaptive utilization of spatio-temporal information from multiple air quality monitoring stations under a multi-scale framework. To this end, it is necessary to consider different strategies to combine different methods to decompose the given series and to fuse multi-site information. Based on a systematic comparative analysis, we propose a novel multi-scale and multi-site modeling method named the multivariate empirical mode decomposition and spatial cosine-attention-based long short-term memory (MEMD-SCA-LSTM). The MEMD-SCA-LSTM first employs MEMD to decompose the multi-site air quality series into the scale-aligned components and then models the components at different scales. The high-frequency components are modeled by a novel SCA-LSTM, which employs LSTM with residual blocks to extract the temporal information and utilizes an attention mechanism based on the cosine similarity to adaptively extract interactions among different sites. Other relatively regular components are modeled by the LSTM. Empirical study on PM2.5 in Hong Kong has demonstrated the effectiveness of fusing multi-site information using the spatial attention (SA) mechanism under the multi-scale framework with MEMD. The proposed MEMD-SCA-LSTM can improve the one-day ahead modeling performance with the mean absolute error and the root mean square error reduced over 10%, compared to the baseline modeling methods. For the two-day and three-day ahead performance, the MEMD-SCA-LSTM is still the best one. Furthermore, by visualizing the attention weights, we illustrate that our proposed SCA-LSTM can overcome some limitations of the traditional attention mechanisms and that the attention weights exhibit more informative patterns which could be used to analysis the transport of air pollutant between sites. The proposed modeling method is a general method, which is feasible and applicable to other pollutants in other cities or regions.