Surface incident solar radiation (R s ) is a key parameter in many climatic and ecological processes. The data from satellites and reanalysis have been widely used. However, for reanalysis, R s data has been shown to have substantial spatial bias, and the time span of reliable satellite R s is too short for climatic and ecological studies. Combining reanalysis and satellite data would be an effective method for generating long-term and consistent R s datasets. Here, we apply a cumulative probability density function-based (CPDF) method to merge eight reanalyses with the latest available satellite R s data from Clouds and Earth's Radiant Energy System Energy Balanced and Filled (CERES EBAF) surface retrievals. The CPDF method not only reduces the spatial bias of the reanalysis R s data, but also makes the R s datasets in a global, long-term and consistent way. The observed R s data collected at 54 Baseline Surface Radiation Network (BSRN) stations from 1992 to 2016 are used to evaluate the method. Results show that the CPDF method could reduce the mean absolute biases (MAB) of the reanalysis R s effectively by 21.24-64.36%. The European Centre for Medium-Range Weather Forecasts Re-Analysis interim (ERA-interim) reanalysis R s data, which are available for 1979 onward, perform the best before MAB = 13.20 W·m −2 and after MAB = 10.40 W·m −2 merging. This small post-merging MAB of the ERA-interim reanalysis is caused by the MAB of 9.90 W·m −2 in the satellite R s retrievals. The Japanese 55-year reanalysis provides R s values back to 1958, and CPDF can reduce its MAB by 32.87%, to 11.17 W·m −2 . The National Oceanic and Atmospheric Administration (NOAA)-CIRES twentieth-century reanalysis (CIRES) and the ECMWF twentieth-century reanalysis (ERA20CM) provide century-long R s estimates. CIRES performs better after merging. The MAB of CIRES can be reduced by 32.10%, to 12.99 W·m −2 , while ERA20CM's can be reduced by 12.51%, to 16.40 W·m −2 .