The low accuracy of satellite Cloud fraction (CF) over the Arctic seriously restricts accurate assessment of regional and global radiant energy balance under the changing climate. Previous studies have reported that not a single satellite CF product could satisfy the needs of accuracy and spatio-temporal coverage simultaneously for long-term applications over the Arctic. Merging multiple CF products with complementary properties is an effective way to produce more spatiotemporally complete and accurate CF data record. This study proposed a spatiotemporal statistical data fusion framework based on cumulative distribution function (CDF) matching and Bayesian maximum entropy (BME) method to produce a syncretic 1°×1° CF dataset in the Arctic during 2000-2020. The original datasets contain CF from MOD08/MYD08, CERES-SSF Terra/Aqua, CLARA-A2 AM/PM, PATMOS-x AM/PM, ISCCP-H AM/PM. The fused CF product is more consistent with the active satellite data GEWEX-CALIPSO and the ground-based observation data CRU TS4.05.