Fig. 1

Land sink efficiencies based on 5-yr moving averages from different datasets during 1959–2016: a annual values, b annual values (dot) with linear regression (line) from 1959 to 1988 and piecewise regression (line) from 1996–2016. B in the inversion datasets (CAMS and Jena) is calculated by BL−L where L from BLUE is used. The diamonds and error bars in (b) indicate the detected breakpoints with 95% confidence interval. The Pinatubo eruption period is marked in light yellow shade. Note that we used 5-yr moving average fluxes, and thus the analysis period is shown as 1959–2016 instead of 1957–2018 (the original annual values)