Mass and energy conservation after drift removal. For each model, linear regression coefficients were calculated between pairs of decadal-mean de-drifted time series of interest including the time-integrated netTOA (Qr),  time-integrated heat flux into ocean (Qh), time-integrated moisture flux into atmosphere (Qep), time-integrated freshwater flux into ocean (Qm), temperature component of globally integrated OHC (HT), ocean mass (M), ocean salinity (S) and mass of atmospheric water vapor (Ma). Each box shows the ensemble quartiles for the coefficients, while the whiskers extend to show the rest of the distribution, except for points determined to be outliers (values beyond 1.5 times the inter-quartile range). Values for each model (including the small number of outliers beyond the plot bounds) are given in Figures S1 and S2.