Global 20-year XCO2 mapping through synergy of multi-satellite observations

• Combine local LASSO and de-trending methods to reconstruct global XCO 2 , 2003–2022 • Correct discrepancies in multi-satellite XCO 2 data to achieve high consistency. • Reconstructed global long-term XCO 2 data align well with ground data (R 2 > 0.98). • Reveal global differences in the correlation between XCO 2 anomalies and NEE. Global-scale, long-term, high-consistency, and high-coverage carbon dioxide (CO 2 ) products are crucial for understanding the dynamics of CO 2 worldwide, which are often generated by integrating multi-source satellite and reanalysis data. However, existing research generally faces several challenges, including inconsistencies among different satellites, limited accuracy of multi-source data fusion modeling, and difficulties in extrapolating beyond the modeling period due to the interannual growth trend of CO 2 . To fill this gap, our study proposes a novel approach for reconstructing global column-averaged dry-air mole fraction of CO 2 (XCO 2 ) products of long time series (2003–2022) and high accuracy, with the combination of local Least Absolute Shrinkage and Selection Operator (LASSO) regression and de-trending methods. The proposed method corrects the differences between multi-source satellites to enhance the consistency of reconstruction results. Furthermore, it accounts for spatio-temporal heterogeneity and enhances extrapolation. Two long-term XCO 2 products are available: multi-satellite XCO 2 (MS-XCO 2 , 2003–2022), which integrates data from five satellites but still exhibits spatial gaps, and MS-CAMS-XCO 2 (2003–2020), which fuses CAMS data and MS-XCO 2 to provide spatially continuous coverage. Validation results against ground stations show high accuracy for both MS-XCO 2 and MS-CAMS-XCO 2 , with R 2 values of 0.985 and 0.989, and RMSE values of 1.08 ppm and 0.997 ppm, respectively. The reconstructed datasets reveal that the growth rate of global XCO 2 is 2.260 ppm/year, with higher CO 2 levels in the mid-latitude northern hemisphere and lower CO 2 levels in the southern hemisphere. Further analysis of the correlation between XCO 2 anomalies and net ecosystem exchange (NEE) indicates a significant negative correlation when the ecosystem is a carbon source. However, the correlation varies when the ecosystem is a carbon sink, with a significant negative correlation observed when vegetation is in good condition. This study generates two datasets of global long-term high-precision XCO 2 , providing valuable data for understanding the global carbon cycle.