Identifying and Estimating Biases between XBT and Argo Observations Using Satellite Altimetry

A methodology is developed to identify and estimate systematic biases between expendable bathythermograph (XBT) and Argo observations using satellite altimetry. Pseudoclimatological fields of isotherm depth are computed by least squares adjustment of in situ XBT and Argo data to altimetry-derived sea height anomaly (SHA) data. In regions where the correlations between isotherm depth and SHA are high, this method reduces sampling biases in the in situ observations by taking advantage of the high temporal and spatial resolution of satellite observations. Temperature profiles from deep XBTs corrected for a bias identified and adopted during the 1990s are considered in this study. The analysis shows that the pseudoclimatological isotherm depths derived from these corrected XBTs are predominantly deeper than the Argo-derived estimates during the 2000-07 period. The XBT-minus-Argo differences increase with depth consistent with hypothesized problems in the XBT fall-rate equation. The depth-dependent XBT-minus-Argo differences suggest a global positive bias of 3% of the XBT depths. The fact that this 3% error is robust among the different ocean basins provides evidence for changes in the instrumentation, such as changes in the terminal velocity of the XBTs. The value of this error is about the inverse of the correction to the XBT fall-rate equation (FRE) implemented in 1995, suggesting that this correction, while adequate during the 1990s, is no longer appropriate and could be the source of the 3% error. This result suggests that for 2000-07, the XBT dataset can be brought to consistency with Argo by using the original FRE coefficients without the 1995 correction.