Estimating ground surface PM 2.5 with fine spatiotemporal resolution is a critical technique for exposure assessments in epidemiological studies of its health risks. Previous studies have utilized monitoring, satellite remote sensing or air quality modeling data to evaluate the spatiotemporal variations of PM 2.5 concentrations, but such studies rarely combined these data simultaneously. Through assembling techniques, including linear mixed effect regressions with a spatial-varying coefficient, a maximum likelihood estimator and the spatiotemporal Kriging together, we develop a three-stage model to fuse PM 2.5 monitoring data, satellite-derived aerosol optical depth (AOD) and community multi-scale air quality (CMAQ) simulations together and apply it to estimate daily PM 2.5 at a spatial resolution of 0.1 • over China. Performance of the three-stage model is evaluated using a cross-validation (CV) method step by step. CV results show that the finally fused estimator of PM 2.5 is in good agreement with the observational data (RMSE = 23.0 µg/m 3 and R 2 = 0.72) and outperforms either AOD-derived PM 2.5 (R 2 = 0.62) or CMAQ simulations (R 2 = 0.51). According to step-specific CVs, in data fusion, AOD-derived PM 2.5 plays a key role to reduce mean bias, whereas CMAQ provides spatiotemporally complete predictions, which avoids sampling bias caused by non-random incompleteness in satellite-derived AOD. Our fused products are more capable than either CMAQ simulations or AOD-based estimates in characterizing the polluting procedure during haze episodes and thus can support both chronic and acute exposure assessments of ambient PM 2.5 . Based on the products, averaged concentration of annual exposure to PM 2.5 was 55.7 µg/m 3 , while averaged count of polluted days (PM 2.5 > 75 µg/m 3 ) was 81 across China during 2014. Fused estimates will be publicly available for future health-related studies.