In this paper results from the application of an ocean data assimilation (ODA) system, combining a multivariate reduced-order optimal interpolator (OI) scheme with a global ocean general circulation model (OGCM), are described. The present ODA system, designed to assimilate in situ temperature and salinity observations, has been used to produce ocean reanalyses for the 1962–2001 period. The impact of assimilating observed hydrographic data on the ocean mean state and temporal variability is evaluated. A special focus of this work is on the ODA system skill in reproducing a realistic ocean salinity state. Results from a hierarchy of different salinity reanalyses, using varying combinations of assimilated data and background error covariance structures, are described. The impact of the space and time resolution of the background error covariance parameterization on salinity is addressed.