Drift in ocean and atmosphere state variables and boundary fluxes related to energy, mass and salt conservation. Each marker represents the linear trend over the full length of the pre-industrial control experiment, with CMIP5 and CMIP6 models designated with open and solid shapes, respectively. The colors represent models from the same institution (Table \ref{table:models}). Drift in ocean salinity was calculated by first converting to an equivalent change in ocean mass (as per Equation \ref{eq:salinity}; see panel d) and the time-integrated moisture flux into the atmosphere (panel e) has not been plotted against the change in atmospheric water vapor because the water vapor trends are negligible in comparison. The thick dashed lines indicate a 1-to-1 relationship (i.e. conservation) and estimates of the magnitude of the current planetary energy imbalance (estimates range from 0.4-1.0 Wm-2; shading in panels a and b), barystatic sea level rise (1.8 mm/year; thin dashed lines in panels c and d) and trend in the global mass of atmospheric water vapor (thin dashed lines in panel e) are shown.