# An assessment of upper ocean salinity content from the Ocean Reanalyses Inter-comparison Project (ORA-IP)

L. Shi, O. Alves, R. Wedd, M. A. Balmaseda, Y. Chang, G. Chepurin, N.

Ferry, Y. Fujii, Fabienne Gaillard, S. A. Good, et al.

► To cite this version:

L. Shi, O. Alves, R. Wedd, M. A. Balmaseda, Y. Chang, et al. An assessment of upper ocean salinity content from the Ocean Reanalyses Inter-comparison Project (ORA-IP). Climate Dynamics, 2017, 49 (3), pp.1009-1029. 10.1007/s00382-015-2868-7. hal-04201922

HAL Id: hal-04201922

https://hal.science/hal-04201922v1

Submitted on 3 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# An assessment of upper ocean salinity content from the ocean reanalyses inter-comparison project (ORA-IP)

Article

Accepted Version

Shi, L., Alves, O., Wedd, R., Balmaseda, M. A., Chang, Y., Chepurin, G., Ferry, N., Fujyii, Y., Gaillard, F., Good, S. A., Guinehut, S., Haines, K. ORCID: https://orcid.org/0000-0003-2768-2374, Hernandez, F., Lee, T., Palmer, M., Peterson, K. A., Masuda, S., Storto, A., Toyoda, T., Valdivieso, M. ORCID: https://orcid.org/0000-0002-1738-7016, Vernieres, G., Wang, X. and Yin, Y. (2017) An assessment of upper ocean salinity content from the ocean reanalyses inter-comparison project (ORA-IP). Climate Dynamics, 49 (3). pp. 1009-1029. ISSN 0930-7575 doi: 10.1007/s00382-015-2868-7 Available at https://centaur.reading.ac.uk/52361/

It is advisable to refer to the publisher's version if you intend to cite from the work. See Guidance on citing.

To link to this article DOI: http://dx.doi.org/10.1007/s00382-015-2868-7

Publisher: Springer
including copyright law. Copyright and IPR is retained by the creators or other
copyright holders. Terms and conditions for use of this material are defined in
the [End User Agreement](#).

www.reading.ac.uk/centaur

CentAUR

Central Archive at the University of Reading

Reading's research outputs online
An Assessment of Upper Ocean Salinity Content from the
Ocean Reanalyses Inter-Comparison Project (ORA-IP)

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

L. Shi¹), O. Alves¹), R. Wedd¹), M. A. Balmaseda²), Y. Chang³),¹⁶), G. Chepurin⁴), N. Ferry¹⁰), Y. Fujii⁵), F.

Gaillard⁶), S. A. Good⁷), S. Guinehut⁸), K. Haines⁹),F. Hernandez¹⁰), ¹⁷), T. Lee¹¹), M. Palmer⁷), D. Peterson⁷), S.

Masuda¹²), A. Storto¹³), T. Toyoda⁵), M. Valdivieso⁹),G. Vernieres¹⁴), X. Wang¹⁵), Y. Yin¹

1) Research and Development Branch, Bureau of Meteorology, Australia

2) European Centre for Medium-Range Weather Forecasting (ECMWF)

3) Geophysical Fluid Dynamics Laboratory (GFDL), Princeton University, U. S. A.

4) University of Maryland, U.S.A.

5) Meteorological Research Institute (MRI), Japan Meteorological Agency

6) Laboratoire de Physique des Océans (LPO/IFREMER), France

7) Met Office Hadley Centre, United Kingdom

8) CLS / Space Oceanography Division, France

9) Reading University, United Kingdom

10) Mercator Océan, Ramonville Saint Agne, France

11) Jet Propulsion Laboratory (JPL), NASA, U. S. A.

12) Research and Development Center for Global Change, Japan Agency for Marine-Earth Science and Technology

13) Euro-Mediterranean Centre for Climate Change (CMCC), Italy

14) Global Ocean and Assimilation Office (GMAO), NASA, U. S. A.

15) Joint Institute for Regional Earth System Science and Engineering (JIFRESSE), UCLA

16) Department of Earth Science, Kongju National University, Korea

17) Institut de Recherche pour le Développement (IRD), Toulouse, France
Abstract

Many institutions worldwide have developed ocean reanalyses systems (ORAs) utilizing a variety of ocean models and assimilation techniques. However, the quality of salinity reanalyses arising from the various ORAs has not yet been comprehensively assessed. In this study, we assess the upper ocean salinity content (depth-averaged over 0-700m) from 14 ORAs and 3 objective ocean analysis systems (OOAs) as part of the Ocean Reanalyses Intercomparison Project (ORA-IP).

Our results show that the best agreement between estimates of salinity from different ORAs is obtained in the tropical Pacific, likely due to relatively abundant atmospheric and oceanic observations in this region. The largest disagreement in salinity reanalyses is in the Southern Ocean along the Antarctic Circumpolar Current (ACC) as a consequence of the sparseness of both atmospheric and oceanic observations in this region. The West Pacific Warm Pool (WPWP) is the largest region where the signal to noise ratio (SNR) of reanalysed salinity anomalies is greater than 1. Therefore, the current salinity reanalyses in the tropical Pacific Ocean may be more reliable than those in the Southern Ocean and regions along the western boundary currents. Moreover, we found that the assimilation of salinity in ocean regions with relatively strong ocean fronts is still a common problem as seen in most ORAs.

The impact of the Argo data on the salinity reanalyses is visible, especially within the upper 500m, where the interannual variability is large. The increasing trend in global-averaged salinity anomalies can only be found within the top 0-300m layer, but with quite large diversity among different ORAs. Beneath the 300m depth, the global-averaged salinity anomalies from most ORAs switch their trends from a slightly growing trend before 2002 to a decreasing trend after 2002. The rapid switch in the trend is most likely an artefact of the dramatic change in the observing system due to the implementation of Argo.
1. Introduction

Studies of the seasonal variability of salinity in the tropical oceans have revealed that the salinity changes, in particular in the upper ocean, are strongly impacted by river discharges, surface freshwater flux (i.e., evaporation and precipitation; referred as E-P hereafter), and advection etc.

(Cronin and McPhaden, 1998; Johnson et al. 2002; Foltz et al. 2004). On decadal timescales, changes in global ocean salinity can be mainly attributed to changes in the global hydrological cycle, in particular the E-P pattern changes, (Curry et al. 2003; Durack and Wijffels 2010; Durack et al.

2012) that is possibly linked to global warming (Held and Soden 2006). In the high latitude Atlantic, the long-term salinity changes are impacted by the wind-driven export of ice, river discharge from the Arctic and advection by the currents (Vinje 2001; Belkin 2004).

A number of papers have indicated that salinity has a great impact on the global ocean dynamical and thermal circulation through density and dynamical height variations (Cooper 1988; Rahmstorf, 1996; Murtugudde and Busalacchi 1998; Vialard et al. 2002; Fedorov et al. 2004; Zhang and Vallis 2006; Huang et al. 2008). For instance, “Great Salinity Anomalies (GSAs)” events, that have occurred during the 1970s, 1980s and 1990s in the North Atlantic Ocean (Dickson et al. 1988; Belkin et al. 1998; Belkin 2004), may play an important role in the variations of the thermohaline circulation, deep western boundary current, northern recirculation gyre, and Gulf Stream, etc.
(Wadley and Bigg 2006; Zhang and Vallis 2006). The close link between salinity variability in the western Pacific and subsequent ENSO events has also been revealed from both observational investigations and dynamical model simulations (Maes et al., 2005, 2006; O’Kane et al. 2014; Zhu et al. 2014). Furthermore, some studies have pointed out that the assimilation of observed salinity can provide more accurate initial ocean states for dynamical models since the imbalance between temperature and salinity is reduced, thereby resulting in better ENSO prediction skills (Ballabrera-Poy et al., 2002; Vialard et al., 2002; Yang et al. 2010; Hackert et al. 2011; Zhao et al. 2013, 2014; Zhu et al. 2014).
Most research and operational centres around the world have established their own ocean reanalysis systems (ORAs) for the purpose of building up historical ocean datasets and providing initial conditions for a range of forecast systems. In the early stages, the products from ORAs (Balmaseda et al 2009) mainly focused on the assimilation of observed temperature, while salinity was not adjusted at all or was constructed from the local climatological temperature-salinity (T-S) relationship, mainly due to the paucity of salinity observations (Behringer and Xue 2004). Since the international Argo Project (http://argo.jcommops.org), which collects real-time temperature and salinity profiles in upper 2000 metres of the ocean, started to provide comprehensive global ocean coverage from around 2006 onwards, most state-of-the-art ORAs assimilate both observed sea temperature and salinity profiles by using a variety of assimilation methods (Table 1).

To date, there is no intercomparison of the performance of salinity reanalyses from the latest vintage of ORAs from around the world. However, there have been several assessments of salinity analyses from some individual ORAs (Hernandez et al. 2009; Xue et al. 2011; Chang et al. 2012; Fujii et al. 2012). The Ocean Reanalyses Intercomparison Project (ORA-IP) was proposed by the participants of the joint GODAE (Global Ocean Data Assimilation Experiment) OceanView/CLIVAR GSOP (Global Synthesis and Observation Panel) workshop in Santa Cruz in June 2011 for the purpose of real time ocean monitoring and operational seasonal forecast systems improvement (Balmaseda et al. 2015). The work presented here is a contribution to the ORA-IP project. Its primary objective is to quantify the ensemble spread and signal to noise ratio in the estimation of salinity from an ensemble of existing global ocean reanalyses. This is the first step to evaluate the maturity level of existing global products. By identifying current deficiencies, it is expected that ORA-IP can help with the future development of ocean data assimilation and observing systems. The work presented here is only an initial and broad quantification of the signal to noise ratio, and it will not deal with the representation of process or specific modes of variability. We expect that further studies can follow once the ORA-IP data is made publicly available.
This paper is organized into six sections. Section 2 presents a brief description of the ORAs included in this study. The mean state of the reanalysed salinity for the period from 1993 to 2010<sup>Δ</sup> is assessed in section 3. In section 4, the variability of the salinity reanalyses, such as standard deviation, signal-to-noise ratio and temporal correlation is evaluated. In this section, the temporal correlation between local salinity anomalies and temperature anomalies and the impact of Argo data on the salinity reanalyses is also discussed. Furthermore, we will also examine the trend in global average salinity anomalies over depth levels above 1500 meters in section 5. The final section contains a discussion and conclusions.

9

10

## 2. Reanalyses Systems

Critical information (e.g., referred names and associated institutions; ocean model resolutions; atmospheric forcing; main assimilation methods, assimilated observations and relaxation to climatology) of the products assessed in this study is summarized in Table 1. The total of 17 estimates can be roughly classified into two groups: the first one consists of fourteen ORAs (the first fourteen in Table 1) which all assimilated various ocean observations into a variety of dynamical ocean model systems (coupled or uncoupled); the second group, in contrast, consists of ARMOR3D, ISAS13 and EN3v2a in which their salinity reanalyses are obtained from ocean observations through a statistical method with no dynamical ocean model. Thus, the three products in the second group are referred to as Objective Ocean Analyses (OOAs) hereafter.

The main assimilation techniques used by the fourteen ORAs can be simply summarised as:

Optimal Interpolation (OI; e.g., SODA and ORAS3), Ensemble OI (EnOI; e.g., GMAO), 3-dimension variational method (3DVAR; e.g., GloSea5, ORAS4, CGLORS, MOVE-C, MOVE-G2 and G2V3), 4-dimension variational method (4DVAR; e.g., ECCOV4 and K7ODA^), Kalman

<footnote>Δ Exception is the ISAS13, which is only available from 2002-2010 in this study. Hereafter, all the calculation of ISAS13, thus, is based on the period from 2002-2010.</footnote>

<footnote>^ Changes in water volume in conjunction with a free surface model used by K7ODA are ignored in this study.</footnote>
1 Filter method (KF; e.g., G2V3), and Ensemble Kalman Filter (EnKF; e.g., ECDA, PEODAS*,

2 PECDAS). The atmospheric surface forcing for most of the 14 ORAs are obtained from

3 atmospheric reanalyses through a variety of methodologies (e.g., corrected fluxes, different bulk

4 formulations) except that three ORAs (e.g., PECDAS, ECDA and MOVE-C) are provided by the

5 atmospheric component of the corresponding coupled model.

6 In addition to observed Temperature/Salinity (T/S) profiles, various other ocean observations,

7 such as altimeter-derived Sea Level Anomalies (SLA); Sea Surface Height (SSH) from tide gauges;

8 Sea Surface Temperature (SST) and satellite derived Sea-Ice Concentration (SIC), are also

9 assimilated by some of the ORAs in Table 1. In order to prevent the model from drifting, most

10 ORAs relax their 3 dimensional (3D) T/S or Sea Surface Salinity (SSS) to climatology with

11 differing relaxation time intervals (Table 1). Although a few studies have suggested the potential

12 impact of SSS on SST variability in the tropical Pacific (Ballabrera-Poy et al. 2002; Wang and

13 Chao 2004; Hackert et al. 2011), the performance of SSS reanalyses will not be assessed in this

14 paper because SSS is defined differently in each ORA (e.g., different layer depths). In this study,

15 we will mainly focus on the depth-averaged salinity over the upper 0-700m ocean layer# (S700)

16 since the largest salinity changes are usually observed in the upper 500 meters (Curry et al. 2003;

17 Boyer et al. 2005; Durack and Wijffels 2010). For this study, all the monthly salinity fields from the

18 ORA-IP participants have been interpolated to a standard 1°×1° latitude-longitude grid.

19 The biggest problem in the assessment of salinity reanalyses from various ORAs is the absence

20 of “reality” or proper standard due to the paucity of salinity observations, especially prior to Argo.

21 We will focus on the ensemble spread (SPD; refer to the Eq. (A2-3) for details; noted as ensemble

22 standard deviation in Balmaseda et al. 2015) of the variable ‘X’ from 14 ORAs about their

23 corresponding ensemble mean (EMORA; refer to the Eq. (A1) for details). This can be used to

24 measure the diversity / agreement of the salinity between different ORAs. Following previous

<footnote>* In a strict sense, PEODAS is an approximate form of an ensemble Kalman filter system (Yin et al. 2011).</footnote>

<footnote># The S700 values in the ocean coast regions, where the deepest depth is less than 700 meters, are defined as missing value.</footnote>
studies (Lee et al. 2009; Zhu et al. 2012), we also consider the standard deviation (STD; refer to
$STD_{EMORA}^A$ in the Eq. (A8)) of EMORA variability as the ‘signal’ (or certainty) part of the variability
of salinity reanalyses from the 14 ORAs. Thus, the corresponding SPD (refer to $SPD_{EMORA}^X$ in the Eq.
(A2-3)) can be considered a quantitative measurement of the ‘noise’ (or uncertainty), which is caused
by different assimilation methods, atmospheric forcing and ocean model dynamics etc., from the 14
ORAs. Therefore, the ‘uncertainty’ mentioned hereafter in this study cannot simply be considered
as the true ‘error’ between the reanalyses and the observations. For instance, it is possible that all
the products have similar systematic error, which will not be captured by this ensemble method. As
discussed by Balmaseda et al (2015), the ensemble method also assumes that all the estimates have
similar quality, which may not always be applicable.

The ensemble mean of the 3 OOAs (EMOO; refer to $\overline{X}_{EMOO}$ or $\overline{X}_{EMOO}^A$ in the Eq. (A1)) will be
compared with EMORA in this paper. This comparison should illustrate the main differences
between statistical and dynamical data assimilation estimates, assuming the bulk of in-situ salinity
observations are likely to be similar in ORAs and OOAs. The salinity from OOAs is likely to be
close to climatology prior to Argo due to the lack of observations. Even during the Argo period, we
note that the salinity reanalyses from the 3 OOAs or the corresponding EMOO cannot simply be
considered as the proxy of ‘reality’ due to inhomogeneous temporal and spatial distribution of Argo.
In contrast, the salinity in ORAs is affected not only by the salinity observations, but also by model
dynamics and mixing, surface fluxes, and imposed multivariate relations (for instance, observations
of temperature and sea level can affect the salinity). This variety of information sources in ORAs
can contribute to the coherence of the signal, but also to the ensemble spread, since there is large
uncertainty in ocean models, surface fluxes and multivariate relationships.

3. Mean State
Compared to EMOO (Fig. 1b), the AMS700 of EMORA (Fig. 1a) is relatively saltier (≥ 0.1 *psu*) in the Southern Ocean along the ACC and in the Kuroshio. A striking dipole appears along the Gulf Stream/Labrador Current system, where the AMS700 of EMORA is fresher/saltier than that of EMOO by more than 0.1 *psu*, respectively. This difference is consistent with the well-known systematic error in ocean models (incorrect strength and path of the Gulf Stream and Labrador Current). It is likely that in this area EMOO are relatively well constrained by the existing salinity observations. The large differences along the northern edge of the ACC are also likely to have a dynamical origin, but in this case EMOO may have significant uncertainty, due the paucity of in-situ observations. There are also large scale differences in the meridional distribution of salinity, which varies across basins. Thus, in the Atlantic, the AMS700 of EMORA shows a saltier equatorial band and fresher sub-tropical gyres than that of EMOO, while the opposite pattern occurs in the Pacific. Differences associated with the Equatorial Pacific current system are also visible.

Generally, the regions of relatively large $SPD_{EMORA}^{AM}$ of AMS700 ($SPD_{EMORA}^{AM} \ge 0.1 \text{ psu}$) shown in Fig. 1c, such as in the Southern Ocean along the ACC, the Kuroshio and the Gulf Stream, correspond with the regions of relatively large differences of AMS700 shown in Fig. 1a. This indicates that the existing observations are not able to constrain the large diversity among different ORAs in these areas, at least with the current assimilation systems. The spatial pattern of $SPD_{EMORA}^{AM}$ in the Atlantic resembles the footprint of the wind driven circulation.

In order to demonstrate the individual performance of each ORA, the zonal distributions of the differences of AMS700 between each product and EMOO are shown in Fig. 2a (meridionally averaged over 30°N-60°N), 2b (15°S-15°N) and 2c (60°S-30°S), respectively. The definition of the shaded band in Fig. 2, which represents the uncertainty range (i.e., $UCR_{EMORA}^{AM}$) of AMS700 differences from the 14 ORAs about their EMORA, is detailed in the Eq. (A4). Thus, the ORAs outside the shaded band can be considered as outliers.
The AMS700 differences between OOAs and EMOO are generally smaller than that between each ORA and EMOO in most parts of oceans except for the relatively large AMS700 differences for ISAS13 shown in the northern Atlantic and the tropical Pacific. This feature can probably be attributed to the calculation period for the ISAS13 (i.e., 2002-2010) which is shorter than that for other products (i.e., 1993-2010).

In the northern band (averaged over 30°N-60°N; Fig. 2a), the *UCR* of AMS700 differences gradually increases eastward in the northern Pacific Ocean, but gradually decreases eastwards in the North Atlantic Ocean. Compared to the tropical Indian Ocean and Atlantic Ocean, the obviously smaller *UCR* of AMS700 differences in the tropical Pacific Ocean can be attributed to the relative abundance of observations there. The *UCR* of AMS700 differences in the southern band (averaged over 60°S-30°S; Fig. 2c) is relatively large in the Indian Ocean sector and Pacific Ocean sector, presumably due to sparsity of observations. However, it is relatively lower in the Atlantic sector of the Southern Ocean. K7ODA is seemingly an outlier. This is mainly from the fact that its S700 is calculated under the assumption of constant water volume of the first ocean layer despite a free surface ocean model applied in K7ODA.

The seasonal cycle (i.e., January-December monthly climatology) of S700 for a moored buoy located at 8°N, 156°E (referred as T8N156E) is selected to compare with all products in Table 1 at the same location (by linear interpolation) and same period (1999-2010). The differences in the observed seasonal cycle for each ORA are shown in Fig. 3. This buoy has observations for the longest available period (from February 1999) and covers the most depth layers (from 1.5m to 750m depth) among all buoy sites of the TAO/TRITON array (available at http://www.pmel.noaa.gov/tao/proj_over/triton.html). Generally, the average $UCR_{EMORA}^{SC}$ (≈ 0.02 *psu*; shaded band in Fig. 3; refer to the Eq. (A4)) along the seasonal cycle of S700 differences between all ORAs and the T8N156E buoy are significantly smaller than that (≈ 0.1 *psu*) in other regions (Fig. 2). The bias of EMOO is comparable to that of EMORA. This is not surprising, given the
1 availability of observations at this specific location. A few individual ORAs show differences

2 comparable with EMOO or EMORA. However, the spread among the ORAs is much larger than

3 the spread among OOAs. This is indicative that errors in ocean models, surface forcing and data

4 assimilation methods are still an issue for the precise estimation of salinity. For instance, over-

5 estimation of precipitation in the tropical band by most atmospheric reanalyses has been reported by

6 several studies (Janowiak et al. 2010; Kim and Alexander 2013).

7

8 **4. Temporal Variability**

9 **4.1 Standard Deviation**

10 The first assessment of salinity variability is the standard deviation (STD; refer to the Eq. (A7-

11 8)), which is an important indicator of the amplitude of S700 anomalies⁺ (seasonal cycle removed)

12 from all products for the period 1993-2010 (Fig.4). In the northern band (30°N-60°N average; Fig.

13 4a), the STD<sub>n</sub><sup>A</sup> (refer to the Eq. (A7)) of S700 anomalies from most ORAs in the central north Pacific

14 is generally smaller than that in the north-western and north-eastern Pacific, but the UCR<sub>EMORA</sub><sup>STD</sup> (refer

15 to the Eq. (A4); shaded band in Fig. 4) is very similar over the whole northern Pacific. In the North

16 Atlantic Ocean, the STD of S700 anomalies from all products, as well as the UCR<sub>EMORA</sub><sup>STD</sup> , significantly

17 decreases from west (i.e., the Gulf Stream) to east. The STD<sub>EMORA</sub><sup>A</sup> (refer to the Eq. (A8)) of S700

18 anomalies from EMORA agrees well with that from EMOO (i.e., STD<sub>EMOO</sub><sup>A</sup> ; refer to the Eq. (A8))

19 except for the Gulf Stream.

20 In the tropical oceans (15°S-15°N average; Fig. 4b), the largest UCR<sub>EMORA</sub><sup>STD</sup> occurs in the Atlantic

21 Ocean, especially in the eastern part of the basin. In the Indian Ocean, the UCR<sub>EMORA</sub><sup>STD</sup> and amplitude

22 of S700 anomalies is largest in the central-eastern Indian Ocean (around 90°E). In the Pacific, both

23 the UCR<sub>EMORA</sub><sup>STD</sup> and the amplitude of the S700 anomalies is largest in the western edge of the WPWP

<footnote>+ Hereafter, the ‘anomalies’ in this study are relative to the corresponding January-December monthly climatology (i.e., seasonal cycle).</footnote>
region (around 165°E), and then, decreases both eastwards and westwards. In contrast, in the
Atlantic Ocean the largest STD<sup>A</sup><sub>n</sub> from most ORAs, as well as the UCR<sup>STD</sup><sub>EMORA</sub>, is seen in both the east
and west (e.g., the Gulf of Guinea). The STD<sup>A</sup><sub>EMORA</sub> of S700 anomalies from EMORA is consistently
larger than the STD<sup>A</sup><sub>EMOO</sub> from EMOO, especially in the central tropical Indian Ocean, central-western
tropical Pacific and western tropical Atlantic Ocean. Most individual ORAs (with the exceptions of
ECCOV4 and K7ODA) exhibit higher variability than that of individual OOAs, which highlights
the contribution of models and surface forcing to the estimation of salinity variability. The
differences in variability of S700 anomalies among the 3 OOAs are not small, even in the tropical
Pacific Ocean where there is a relative abundance of observations. ARMOR3D seems to be the
outlier, showing very small STD of S700 anomalies.

The STD<sup>A</sup><sub>n</sub> of S700 anomalies from most ORAs, as well as the UCR<sup>STD</sup><sub>EMORA</sub>, in the Southern Ocean along the ACC region (Fig. 4c) is generally larger than that in the tropical oceans and the northern band. This relatively large UCR<sup>STD</sup><sub>EMORA</sub> is likely caused by the lack of observation that results in the STD of S700 anomalies becoming more dependent on the ocean model, assimilation method and atmospheric forcing etc. The STD<sup>A</sup><sub>EMORA</sub> of S700 anomalies from EMORA is smaller than that of most individual ORAs, and comparable to the STD<sup>A</sup><sub>EMOO</sub> from EMOO, suggesting a lack of coherence in the variability of individual ORAs. Exceptions are the convergence zone (around 45°W) of the Brazil Current and the South Atlantic Current; the convergence zone (around 15°E) of the Benguela Current and the Agulhas Current. In these dynamically active regions the ocean model is likely to be playing a significant role.

It is worth noting that most ORAs show relatively large STD of S700 anomalies than the 3
OOAs over most parts of oceans. Although the amplitude of salinity variability may be
overestimated by the model-based ORAs, it is also possible that the 3 OOAs underestimate the
variability because in regions of sparse observations they will be closer to climatology. The
UCR<sub>EMORA</sub><sup>STD</sup> (up to 0.05 *psu*) of the STD<sub>II</sub><sup>A</sup> of S700 anomalies in Fig. 4 is smaller than that of corresponding AMS700 differences (around 0.1 *psu*) shown in Fig. 2, except for the tropical Pacific Ocean and the Atlantic sector of the Southern Ocean. Additionally, we note that the STD<sub>EMORA</sub><sup>A</sup> of EMORA is smaller than that of individual ORA in most cases, except for the tropical Indian-Pacific Ocean. This feature implies that the variability of S700 anomalies from individual ORA is quite diverse, and different from EMORA in most cases. The phase agreement of S700 variability between all individual ORAs and EMOO will be assessed further in the next sub-section 4.2.

There is no specific ORA that is an overall outlier in Figs 2 and 4. It suggests that no specific ORA is the best or worst one among all 14 ORAs in this study. However, in specific regions there are specific outliers, for example, the K7ODA in the North Atlantic Ocean and the Southern Ocean; the GloSea5 and G2V3 in the tropical Pacific Ocean; the PEODAS and PECDAS in the western tropical Atlantic Ocean.

13

## 4.2 Signal-to-Noise Ratio

The SPD of S700 anomalies from each ORA about the corresponding EMORA (refer to SPD<sub>EMORA</sub><sup>A</sup> in the Eq. (A3)) is shown in Fig.5a. The geographical distribution of the largest SPD<sub>EMORA</sub><sup>A</sup> of S700 anomalies (≥ 0.1 *psu*) is associated with the largest SPD<sub>EMORA</sub><sup>AM</sup> of AMS700 (see Fig. 1c), particularly in the western boundary currents, such as the Kuroshio, Gulf Stream and Brazil Current. Other areas of relatively large SPD<sub>EMORA</sub><sup>A</sup> of S700 anomalies (≥ 0.06 *psu*) can be seen in the sub-tropical eastern Indian Ocean and central Pacific Ocean. As was discussed in Balmaseda et al. (2015), the areas of the relatively large uncertainty in salinity reanalyses tends to occur in regions associated with both strong temperature and salinity fronts. Of course, the effects of the ocean models and assimilation techniques on the uncertainty cannot be discarded.
As mentioned above, the STD<sup>A</sup><sub>EMORA</sub> of S700 anomalies from EMORA (Fig. 5b), can be considered as a quantitative estimate of the signal. The regions with the largest signal mainly occur in the WPWP region, central Indian Ocean, Gulf of Alaska along the Alaska Current and a narrow band in the Southern Ocean (around 40°S, 20°W-70°E). It is also high in areas of strong variability such as the western boundary currents in Atlantic Ocean. Strong variability occurs in the WPWP due to strong rainfall and current variability.

Following the approach used by Lee et al. (2009) and Zhu et al. (2012); the ratio of the STD<sup>A</sup><sub>EMORA</sub> (Fig.5b) to the SPD<sup>A</sup><sub>EMORA</sub> (Fig.5a) of S700 anomalies can be considered as the so-called signal to noise ratio (SNR) that gives a good quantitative estimate of the reliability of S700 variability among the different ORAs. As shown in Fig. 5c, the regions where the SNR is greater than 1 mainly appear in the WPWP region, central tropical Indian Ocean, the Gulf of Alaska and other small regions of the mid-latitude oceans. The relatively large SNR over the WPWP, which is also an area of large interannual variability, is likely related to the constraint provided by the salinity observations from the TAO/TRITON moorings. Overall, the SNR is less than 1 over most parts of oceans, indicating that there is relatively large SPD<sup>A</sup><sub>EMORA</sub>, and therefore, disagreement in the estimates of S700 anomalies among different ORAs.

## 4.3 Correlation

Fig.6a illustrates how well the S700 variability in the two ensemble means agree with each other. Correlations are relatively high (≥0.75) in the central and western equatorial Pacific, western sub-tropical Pacific along the Kuroshio and north-eastern mid-latitude Pacific. They are also high in the eastern equatorial Indian Ocean, and throughout parts of the sub-tropical and mid-latitude oceans. Correlations are relatively low (≤0.5) around the northern edge of the ACC, western Indian Ocean and parts of the sub-tropical Atlantic, particularly downstream of the Mediterranean outflow.
The S700 variability of each ORA can be correlated with that of EMOO. And then, the $SPD_{EMORA}^{COR}$ of the correlations from 14 ORAs about their corresponding ensemble average of all 14 correlations (refer to the Eq. (A2) for details) provide an indication of the disagreement in the estimate of variability between the different systems (Fig. 6b). There is some correspondence between areas with large $SPD_{EMORA}^{COR}$ and low correlation in Fig. 6a, such as the northern edge of the ACC in the Pacific sector and the northern part of the tropical Atlantic. Equally, the high correlation in the Tropical Pacific, Eastern Indian Ocean, North East Pacific and North East Atlantic, where the spread is low, is indicative of consistency between the different estimates. The Southern Ocean is an exception, showing relatively large values of the correlation and the $SPD_{EMORA}^{COR}$.

The WPWP region and central tropical Indian Ocean, the regions with the smallest $SPD_{EMORA}^{COR}$ (best agreement among all ORAs, Fig. 6b), also have the highest SNR values (Fig. 5c). It is worth noting that these regions are also the places where the largest precipitation variability occurs (Storto et al. 2015).

Fig. 7 shows the zonal distributions of the correlation of S700 anomalies between each ORA and EMOO for the period 1993-2010. EMORA, unsurprisingly, obtains the highest correlation among all ORAs over most parts of oceans due to averaging out the impacts of different ocean models, forcing fields and assimilation techniques etc. Generally, the areas of relatively high correlation, associated with relatively small $UCR_{EMORA}^{COR}$ (refer to the Eq. (4)), are in the eastern tropical Indian Ocean, central-western tropical Pacific and north-east Pacific. In contrast, the areas of relatively low correlations, associated with relatively large $UCR_{EMORA}^{COR}$, are in the Southern Ocean and the tropical and north-west Atlantic Ocean.

22

## 4.4 Local T-S Correlation

The close relationship between seawater temperature and salinity (T-S) is too complicated to be precisely described and measured in one simple way. However, in this study, we utilize the
correlation between the local S700 anomaly and the corresponding depth-averaged temperature (over upper 0-700m ocean layer, referred as **T700**) anomaly to investigate how the co-variability between T700 and S700 anomalies represented by the ORAs compares with that of EMOO. Figs. 8a and 8b show the temporal correlation between the T700 anomaly and S700 anomaly from EMORA (Fig. 8a) and EMOO (Fig. 8b) for the period 1993-2010. The distribution of T700-S700 correlations from EMORA is quite similar to those of EMOO, showing coherent large scale patterns. For instance, relatively high positive correlation in most parts of Atlantic Ocean, equatorial Pacific, north-eastern and southern sub-tropical Pacific, southern sub-tropical Indian Ocean; whereas, negative correlations mainly occur in the eastern equatorial Indian Ocean, the north-eastern boundary of the Pacific (i.e., off the east coast of Mexico and the Gulf of Alaska), and in particular the Southern Ocean. However, EMORA (Fig. 8a) produces more extreme positive and negative correlation than EMOO (Fig. 8b). Yet, the region of negative correlation for EMORA is smaller than that of EMOO, in particular the narrower negative correlation belt in the Southern Ocean in Fig. 8b. It seems that most ORAs exhibit a stronger relationship between the local salinity content and heat content compared to EMOO, in particular in the north-west Indian Ocean, the equatorial Atlantic Ocean, and the northern edge of the ACC in the Indian-Pacific Ocean sector. In this area, where there are few temperature and salinity observations, the associated relationship between local salinity and temperature in the ORAs may come from the ocean-model information, which is absent in the OOAs.

There is no unique explanation for the large scale patterns of correlation between T700 and S700. It is possible that changes associated with local vertical displacement of the water column related with, say, variations in Ekman pumping, would result in positive/negative correlation of T700-S700 wherever the temperature and salinity vertical stratification (above and below 700m) has the same/opposite sign. But changes associated with horizontal displacement of water masses, or changes in the water mass properties cannot be discarded either.
Fig.8c shows the $SPD_{EMORA}^{COR}$ (refer to the Eq. (A2) for details) of the T700-S700 correlation from all ORAs about the corresponding ensemble average correlation (refer to $\overline{X}_{EMORA}$ in the Eq. (A1)) for the period 1993-2010. The largest $SPD_{EMORA}^{COR}$ of the T700-S700 correlation among the ORAs occurs in the Southern Ocean. This disagreement can be attributed to the lack of observations in this region, and to the different T700-S700 correlation among ocean models.

The zonal distributions of T700-S700 correlation from each ORA are shown in Fig.9. The correlations from most ORAs agree quite well with that of the OOAs in some parts of ocean, for instance, in the mid-latitude North Pacific and the eastern tropical Pacific. In some other parts of oceans, however, the correlation of most ORAs is higher than that of OOAs, such as in the central-western tropical Pacific, central tropical Indian Ocean, the Southern Ocean in both Indian Ocean sector and east Pacific sector, and in particular in the tropical Atlantic Ocean where the correlation of OOAs is around 0.4 but the correlations of most ORAs are generally more than 0.6 (Fig.9b). In the tropical Indonesian Sea/eastern Indian Ocean area (90°E-130°E), most ORAs show negative local T700-S700 correlations (up to -0.6; Fig. 9b). We suspect this may be attributed to the Indonesian Throughflow (ITF) transporting relatively warmer and fresher sea water from the Pacific Ocean into the Indian Ocean (Vranes et al. 2002; Sprintall et al. 2009).

17

## 4.5 Impacts of Argo

Figure 10 shows the depth/time evolution (1993-2010) of the centred pattern correlation coefficients$^{^}$ (i.e., CPCOR; refer to the Eq. (A9); Santer et al. 1993; Storch and Navarra 1999) of salinity anomalies between EMORA and EMOO, averaged over 0-360°E; 30°N-60°N (Fig.10a); 15°S-15°N (Fig.10b); 60°S-30°S (Fig.10c), respectively. The vertical distribution of CPCOR in both the mid-latitude northern oceans (Fig.10a) and the tropical oceans (Fig.10b) generally decrease downwards from 100-1500m depth. Relatively high correlation (≥ 0.5) is obtained over the upper

<footnote>$^{^}$ A 7-month running mean has been applied on the computed correlation coefficients to remove the intra-seasonal variability.</footnote>
1 300-400m depth before 2001, and then, increases up to more than 0.9 and extended to 500-600m
2 depth after 2002. This increase corresponds quite well with the beginning of the Argo project. The
3 relatively higher CPCOR within the upper 500m ocean layer in both mid-latitude and tropical
4 oceans after 2002 indicates higher consistency in the estimated pattern of salinity anomalies by both
5 EMORA and EMOO due to Argo. The CPCOR decreases downwards with the depth increasing,
6 and is likely attributed to a lack of coherence between the anomaly patterns from EMORA and
7 EMOO.

8 In contrast, the vertical distribution of *CPCOR* in the Southern Ocean band (Fig. 10c),
9 increases downward from the surface to 1500m depth before 2002. The reason for this feature needs
10 to be investigated further. This may be related to the existence of slowly varying spatial salinity
11 patterns at these latitudes, which can be sample even with a limited set of observations. A large
12 value of the correlation is not synonymous of adequate sampling though. The influence of a few
13 deep observations may also persist for longer in the slowly varying deep ocean. Hence, the
14 relatively higher correlation in the deep Southern Ocean prior to Argo may be also an artefact of
15 using climatology in all the estimates, either as a prior in the OOAs, or as a nudging term in the
16 ORAs. There is a period of lower spatial correlations during 2003, probably associated with the
17 diversity of ways in which different systems adjust to the spin-up of Argo (including different
18 quality control decisions). Since 2003, the correlation in the Southern Ocean also increases over the
19 0-1000m depth, in particular over 0-400m depth after 2009, as Argo floats were deployed.

20 Generally, after Argo the *CPCOR* over the upper 0-500m of the ocean in both northern and
21 southern mid-latitude oceans are significantly smaller than that in the tropical oceans. This indicates
22 that there is still room for the salinity reanalyses in both northern and southern mid-latitude oceans
23 to be further improved and highlights the need for more Argo floats in this region.

25 **5. Trend**
A phenomenon, which has been noticed by previous studies (Levitus et al. 2009; Xue et al. 2012; Balmaseda et al. 2015; Palmer et al. 2015) and announced by the IPCC (2013) and operational or research centres (http://www.epa.gov/climatechange/science/indicators/oceans/ocean-heat.html), is that the global averaged ocean heat content anomalies from either ORAs or OOAs have a growing trend from the 1990s till now. A similar growing trend has also been found in the steric sea level change (Levitus et al. 2012; also http://www.nodc.noaa.gov/OC5/3M_HEAT Content). This warming trend of ocean heat content anomalies has been considered as strong evidence for global warming in recent decades (IPCC 2013). Therefore, an interesting question is if the corresponding global averaged salinity anomalies retain a similar growing trend since the 1990s? If so, since the total salinity is approximately conserved in the global ocean, are the salinity anomalies that show a growing trend within a certain ocean layer compensated with a decreasing trend at other depths?

Fig. 11(a-c) shows the temporal evolution of global averaged (0-360°E; 60°S-60°N) salinity anomalies (relative to climatology for the period 1993-2010), depth-averaged within 0-300m (i.e., S300; Fig.11a), 300-700m (i.e., S3-700; Fig.11b) and 700-1500m (i.e., S7-1500; Fig.11c) ocean layers, from all products for the period 1993-2010. It is worth noting that the depth-average in this study is calculated within 3 continuous vertical layers (i.e., 0-300m; 300-700m and 700-1500m) rather than the top to bottom vertical average (e.g., 0-300m, 0-700m and 0-1500m) approach used by Balmaseda et al. (2015). Therefore, we can show the different features of the trend in salinity anomalies within different vertical layers. In addition, the reference period for the climatology in this study (1993-2010) is different to the 1993-2007 period used by Balmaseda et al. (2015).

In the 0-300m upper ocean (Fig. 11a), the temporal evolution of global averaged S300 anomalies from all 3 OOAs, EMOO and EMORA show a growing trend similar to the corresponding global averaged temperature anomalies (not shown). However, the temporal evolutions of global averaged S300 anomalies among different ORAs are quite divergent and with a relatively large $UCR_{EMORA}^A$ (refer to the Eq. (A5)). For instance, SODA, GloSea5, ORAS3 and G2V3
1 show rapidly growing trends after the beginning of Argo (2001), in contrast, PECDAS shows a
2 decreasing trend since the end of the 1990s and most other ORAs show weak decreasing and
3 increasing trends. In contrast, it can be seen from Fig.11b that the temporal evolution of the S3-700
4 anomalies from most ORAs shows a trend turning from generally increasing to decreasing after
5 2003. For all OOAs and a few ORAs, such as the ORAS4, CGLORS, there is no clear trend before
6 2003 followed by a very weak decreasing trend after 2003. This feature is quite different from that
7 of the corresponding global average temperature anomalies (not shown). Within the 700-1500m
8 depth layer, for the S7-1500 anomalies (Fig. 11c), most ORAs and OOAs, except for the G2V3,
9 ORAS4 and K7ODA, show a decreasing trend in the reference period, in particular after the
10 beginning of Argo.

11 Generally, the global averaged S300 anomalies from most ORAs and OOAs show a similar
growing trend in the reference period as that shown in corresponding global averaged temperature
12 anomalies, even though there is an increasing discrepancy of S300 anomalies among different
13 ORAs, in particular after the beginning of Argo. As the ocean depth increases, the global averaged
14 salinity anomalies from most ORAs and OOAs show a decreasing trend, in particular after Argo.
15 This feature can probably be explained by the approximate conservation of salt in the global ocean.
16 Interestingly, it can be seen from Figure 11 that most ORAs show a rapid change in both salinity
17 and temperature anomalies (not shown) after the beginning of Argo (i.e., 2002 or 2003). We note it
18 is more likely caused by the changes before and after Argo because of the shortage of reliable
19 observations of both salinity and temperature (in particular in the subsurface ocean) prior to the
20 Argo project.
21

22

23

**6. Conclusions**

24

In this paper, the reanalysed S700 of 14 ORAs from different institutions is assessed to address
the major agreement/disagreement among different ORAs. All ORAs assimilate both temperature
and salinity observations using a variety of ocean models and assimilation methods. In addition,

25

26
three OOAs are also used in this paper as independent data and reference for the assessment. The ensemble spread about the multi-system ensemble mean is utilized to demonstrate the agreement/disagreement and measure the uncertainty range among different ORAs.

Generally, the largest agreement (or smallest uncertainty range) of reanalysed S700 properties, such as mean state, standard deviation and correlation, among different ORAs occurred in the tropical Pacific. The largest disagreement (or uncertainty range) was found in the Southern Ocean along the ACC, and along the western boundary currents, such as the Kuroshio, Gulf Stream and Brazil Current. The main cause for the disagreement in the Southern Ocean can be attributed to both the shortage of ocean and atmospheric observations. Assimilation in the regions along the western boundary currents, in particular the Gulf Stream, can be more difficult as noted by Balmaseda et al. (2015) because the relatively stronger ocean fronts in these regions are not well simulated by the ocean models.

It is shown that the variability of S700 anomalies (i.e., standard deviation) from most of the ORAs is usually stronger over most parts of oceans compared with that from the OOAs. Moreover, the standard deviation of EMORA is smaller than that of most individual ORAs over most parts of oceans, except for the tropical Indian-Pacific Ocean. This is because there is relatively large phase dispersion among different ORAs in these regions. Consequently, EMORA obtains the highest correlation of S700 anomalies with the corresponding EMOO when compared with that of each ORA. A SNR value larger than one is mainly restricted to the WPWP region (probably because of the TAO/TRITON salinity observations), the central tropical Indian Ocean and a few parts of the north Pacific and is associated with the regions with the smallest disagreement of S700 anomalies correlations among different ORAs.

Correlations between T700-S700 anomalies show coherent high values and consistent spatial patterns in both ORAs and OOAs products. The reason for this coherent large scale behaviour needs to be explored further. It may be caused by the temporal variations in availability of observations, which would affect ORAs and OOAs in similar manners. But it may be indicative of real dynamical
signals associated to the large scale ocean circulation. Having good salinity estimations can thus
help us with the understanding and attribution of ocean variability.

The impact of Argo floats on the ocean reanalyses has been shown by some previous studies
(Balmaseda et al. 2007). In this study, our results demonstrated that the tropical oceans/Southern
Ocean have the largest/smallest improvement of salinity reanalyses during the Argo period.

Interestingly, the relatively large improvements in the salinity reanalyses due to Argo are mainly
confined within the upper 500m. It's probably because the models have better physics in the upper
ocean and therefore fit the Argo data in the upper ocean better. We note that the reason for this
phenomenon need to be further investigated in the future.

Although the assimilated global heat content anomalies within upper 700m from most ORAs
and OOAs show an increasing trend during the reference period 1993-2010 (refer to Palmer et al
2015), the global averaged salinity anomalies from most ORAs and OOAs only show an increasing
trend within the top 0-300m layer, in particular in the Argo period. In contrast, in the other two
layers beneath 300m (i.e., 300-700m; 700-1500m), the global averaged salinity anomalies from
most ORAs and OOAs switch their trends from a slightly increasing trend prior to Argo to a
decreasing trend after Argo. We note that there is a rapid change in the trend in global averaged
salinity anomalies around 2002 likely due to Argo.

While there is some agreement regarding the spatial patterns of interannual variability of
salinity and its relation with temperature, large uncertainty remains regarding global averaged
salinity anomaly trends that will affect the estimation of global steric height (Zuo et al 2015). Since
conservation of salt content is considered to be a good approximation, for diagnostic and
attributions studies of global sea level it may be more pertinent to ignore the halo-steric component,
rather than using unreliable halo-steric trends from ORAs. However, the ORA estimation of the
thermo-steric component appears to be more robust (Storto et al, 2015).

Finally, despite the progresses in salinity reanalyses made by most state-of-the-art ORAs, we
note that the current performance of salinity reanalyses from most ORAs is still a long way from
1 being considered a satisfactory and reliable estimation. As mentioned above, the relatively large

2 disagreement/agreement in reanalysed salinity among the different ORAs offers a useful guidance

3 to potential users and scientists. These results highlight ocean regions where the salinity reanalyses

4 may be more reliable (e.g., the tropical Pacific Ocean) and which regions the salinity reanalyses

5 need to be improved (e.g., the Southern Ocean and regions along the western boundary currents).

6 Sustaining and enhancing oceanic measurements of salinity such as those derived from Argo and

7 satellites (e.g., European Space Agency's Soil Moisture and Ocean Salinity and NASA's Aquarius

8 missions) and improving evaporation and precipitation products (e.g., from atmospheric reanalysis)

9 are important to improving the representation of salinity by ORAs in the future. It is also worth

10 noting that the impacts of the ocean models and assimilation techniques on the improvement of

11 salinity reanalyses are also important.
# 1 APPENDIX

2

3 a. Multi-system ensemble mean

In this study, the $X_n^A$ represents the anomaly (seasonal cycle removed) of corresponding total variable $X_n$ for individual n ORA. Thus, the multi-system ensemble mean (i.e., EMORA) of $X_n$ or $X_n^A$ from the 14 ORAs can be given by:

$$ \overline{X}_{EMORA} = \frac{1}{nsys} \sum_{n=1}^{nsys} X_n \quad \text{or} \quad \overline{X}_{EMORA}^A = \frac{1}{nsys} \sum_{n=1}^{nsys} X_n^A \qquad (A1) $$

The $nsys$ represents the total number of all ORAs for calculating EMORA ($nsys = 14$). The corresponding $\overline{X}_{EMOO}$ or $\overline{X}_{EMOO}^A$ can be similarly calculated by the Eq. (A1) except for the $X_n$ or $X_n^A$ of individual n OOA and $nsys = 3$.

11 b. Ensemble spread (SPD)

The ensemble spread of different variables $X$ from 14 ORAs about their corresponding EMORA shown in Fig. 1c, 6b and 8c is given by:

$$ SPD_{EMORA}^X(i, j) = \sqrt{\frac{1}{nsys} \sum_{n=1}^{nsys} (X_n(i, j) - \overline{X}_{EMORA}(i, j))^2} \qquad (A2) $$

Here, the $X$ represents the annual mean (AM) of S700 in Fig. 1c (i.e., $SPD_{EMORA}^{AM}$), the correlation of S700 anomalies in Fig. 6b and the correlation of T700-S700 anomalies in Fig. 8c (i.e., $SPD_{EMORA}^{COR}$), respectively. The $i/j$ represents the longitude/latitude, respectively.

Similarly, the $SPD_{EMORA}^A(i, j)$, that is shown in Fig. 5a, can be calculated as:

$$ SPD_{EMORA}^A(i, j) = \sqrt{\frac{1}{mons} \sum_{t=1}^{mons} \frac{1}{nsys} \sum_{n=1}^{nsys} (X_n^A(i, j, t) - \overline{X}_{EMORA}^A(i, j, t))^2} \qquad (A3) $$
Here, the $X_n^A(i, j, t)$ denotes the S700 anomalies for individual $n$ ORA. The *mons* is the total number of months for the variable $X$ (i.e., *mons* = 216 for the period 1993-2010).

### c. Uncertainty range (UCR)

The uncertainty range (i.e., the shaded band shown in Fig. 2) of the meridionally-averaged AMS700 (i.e., $X$) from 14 ORAs about their corresponding EMORA (i.e., $\overline{X_{EMORA}}$) is defined as:

$$UCR_{EMORA}^{AM}(i) = \overline{X_{EMORA}}(i) \pm SPD_{EMORA}^{AM}(i) \qquad (A4)$$

Here, the $SPD_{EMORA}^{AM}(i)$ can be calculated by the Eq. (2) but without the dimension $j$. The *UCR* shown in Fig. 4, Fig. 7 and Fig. 9 can be similarly calculated by the Eq. (A4) except that the variable $X$ should be replaced by standard deviation for Fig. 4, the correlation coefficients for Fig. 7 and Fig. 9, respectively. In addition, the *UCR* shown in Fig. 3, where the $X$ represents the seasonal cycle of S700, can be also calculated by the Eq. (A4) except for replacing the dimension $i$ by the dimension $t$.

The $UCR_{EMORA}^A(z, t)$ shown in Fig. 11 is given by:

$$UCR_{EMORA}^{A}(z,t) = \overline{X_{EMORA}^{A}(z,t)} \pm SPD_{EMORA}^{A}(z,t) \qquad (A5)$$

Here, the $\overline{X_{EMORA}^A(z,t)}$ denotes the global averaged salinity anomaly in different ocean layers $z$ for

EMORA. And, the $SPD_{EMORA}^A(z, t)$ is calculated as:

$$SPD_{EMORA}^{A}(z,t) = \sqrt{\frac{1}{n_{sys}} \sum_{n=1}^{n_{sys}} \left( X_{n}^{A}(z,t) - \overline{X_{EMORA}^{A}(z,t)} \right)^{2}} \qquad (A6)$$

Here, the $X_n^A(z, t)$ denotes the global averaged salinity anomaly in different ocean layers $z$ for the individual $n$ ORA.

### c. Standard deviation (STD)
1 The STD of the meridionally-averaged S700 anomalies (i.e., $X_n^A(i,t)$) for individual n ORA (i.e.,

2 STD_n^A(i)) and the corresponding EMORA (i.e., STD_{EMORA}^A(i)), which is shown in Fig. 4, is respectively

3 given by:

$$4 \qquad STD_n^A(i) = \sqrt{\frac{1}{mons} \sum_{t=1}^{mons} (X_n^A(i,t))^2} \qquad (A7)$$

5 and

$$6 \qquad STD_{EMORA}^A(i) = \sqrt{\frac{1}{mons} \sum_{t=1}^{mons} (\overline{X_{EMORA}^A}(i,t))^2} \qquad (A8)$$

7 The corresponding STD_{EMOO}^A(i) in Fig. 4 can be similarly calculated by the Eq. (A8) except for
8 replacing the $\overline{X_{EMORA}^A}(i,t)$ by the $\overline{X_{EMOO}^A}(i,t)$. Additionally, the STD of S700 anomaly for EMORA (i.e.,

9 STD_{EMORA}^A(i, j)), which is shown in Fig. 5b, can be also calculated by Eq. (A8) except for adding the
10 dimension j.

11

12 d. Centred pattern correlation (CPCOR)

13 The centred pattern correlation (i.e., CPCOR(z, t)) of salinity anomalies (seasonal cycle removed)

14 between EMORA (i.e., $\overline{X_{EMORA}^A}(i, j, z, t)$) and EMOO (i.e., $\overline{X_{EMOO}^A}(i, j, z, t)$) as a function of depth (0-

15 1500m) and time (1993-2010), which is shown in Fig. 10, is defined as:

$$16 \qquad CPCOR(z,t) = \left[ \frac{1}{m \times n} \sum_{j=1}^{n} \sum_{i=1}^{m} \left( (\overline{X_{EMORA}^A}(i,j,z,t) - M_{EMORA}) (\overline{X_{EMOO}^A}(i,j,z,t) - M_{EMOO}) \right) \right] / (S_{EMORA} S_{EMOO}) \qquad (A9)$$

17 Here, the m/n denotes total longitude/latitude grids of the calculated ocean band. The $M_{EMORA}$ /

18 $M_{EMOO}$ denotes the total mean of the $\overline{X_{EMORA}^A}(i, j, z, t) / \overline{X_{EMOO}^A}(i, j, z, t)$ over the calculated ocean band,

19 respectively. Thus, the $M_{EMORA}$ can be given by:
$$1 \qquad M_{EMORA} = \frac{1}{m \times n} \sum_{j=1}^{n} \sum_{i=1}^{m} \overline{X_{EMORA}^{A}}(i, j, z, t) \qquad (A10)$$

2 Thus, the $M_{EMOO}$ can be similarly calculated by the Eq. (A10) except for replacing the $\overline{X_{EMORA}^{A}}$ by

3 the $\overline{X_{EMOO}^{A}}$. The $S_{EMORA} / S_{EMOO}$ in Eq. (A10) denotes the spatial standard deviations of the $\overline{X_{EMORA}^{A}}$ /

4 $\overline{X_{EMOO}^{A}}$, respectively. The $S_{EMORA}$ can be obtained by:

$$5 \qquad S_{EMORA}(z,t) = \sqrt{\frac{1}{m \times n} \sum_{j=1}^{n} \sum_{i=1}^{m} \left( \overline{X_{EMORA}^{A}}(i, j, z, t) - M_{EMORA} \right)^2} \qquad (A11)$$

6 The corresponding $S_{EMOO}$ can be also obtained by the Eq. (A11) except for replacing the $\overline{X_{EMORA}^{A}}$

7 by the $\overline{X_{EMOO}^{A}}$.

8
1 Reference

2 Ballabrera-Poy, J., R. Murtugudde, and A. J. Busalacchi, 2002: On the potential impact of sea

3 surface salinity observations on ENSO predictions, *J. Geophys. Res.*, **107**(C12), 8007,

4 doi:10.1029/2001JC000834

5 Balmaseda, M.A., D. Anderson, and A. Vidard, 2007: Impact of Argo on analyses of the global

6 ocean, *Geophys. Res. Lett.*, **34**, L16605, doi:10.1029/2007GL030452.

7 Balmaseda, M. A., A. Vidard, and D. Anderson, 2008: The ECMWF ORA-S3 ocean analysis

8 system. *Mon. Wea. Rev.*, **136**, 3018-3034.

9 Balmaseda M. A, O. J. Alves, A. Arribas, T. Awaji, D.W. Behringer, N. Ferry, Y. Fujii, T. Lee, M.

10 Rienecker, T. Rosati, and D. Stammer, 2009: Ocean initialization for seasonal forecasts.

11 *Oceanography*, **22**, 154-159.

12 Balmaseda M.A., K. Mogensen, A. Weaver, 2013: Evaluation of the ECMWF Ocean Reanalysis

13 ORAS4. *Q. J. R. Meteor. Soc.*, **139**.1132–1161.

14 Balmaseda M.A., and Co-authors, 2015: The Ocean Reanalyses Intercomparison Project (ORA-IP).

15 *J. Operational Oceanography*, **7**, 81-99.

16 Behringer, D. W., and Y. Xue, 2004: Evaluation of the global ocean data assimilation system at

17 NCEP: The Pacific Ocean. *Eighth Symposium on Integrated Observing and Assimilation*

18 *Systems for Atmosphere, Oceans, and Land Surface*, AMS 84th Annual Meeting, Washington

19 State Convention and Trade Center, Seattle, Washington, 11-15 Jan.

20 Belkin, I. M., S. Levitus, J. Antonov, and S. A. Malmberg, 1998: “Great Salinity Anomalies” in the

21 North Atlantic. *Prog. Oceanogr.*, **41**, 1-68.

22 Belkin, I. M., 2004: Propagation of the “Great Salinity Anomaly” of the 1990s around the northern

23 North Atlantic. *Geophys. Res. Lett.*, **31**, L08306, doi:10.1029/2003GL019334.

24 Boyer, T. P., S. Levitus, I. Antonov, R. A. Locarnini, and H. E. Garcia, 2005: Linear trends in

25 salinity for the World Ocean, 1955-1998. *Geophys. Res. Lett.*, **32**, L01604,

26 doi:10.1029/2004GL021791.
1 Carton, J. A., and B. S. Giese, 2008: A reanalysis of ocean climate using Simple Ocean Data

2 Assimilation (SODA), *Mon. Wea. Rev.*, **136**, 2999-3017.

3 Chang, Y.-S., S. Zhang, A. Rosati, T. Delworth, and W. F. Stern, 2012: An assessment of oceanic

4 variability for 1960-2010 from the GFDL ensemble coupled data assimilation, *Clim. Dyn.*,

5 doi:10.1007/s00382-012-1412-2

6 Cooper, N. S., 1988: The effect of salinity on tropical ocean model. *J. Phys. Oceanography*. **18**,

7 697-707.

8 Cronin, M. F., and M. J. McPhaden, 1998: Upper ocean salinity balance in the western equatorial

9 Pacific, *J. Geophys. Res.*, **103**, 27567-27588.

10 Curry, R., B. Dickson, and I. Yashayaev, 2003: A change in the freshwater balance of the Atlantic

11 Ocean over the past four decades. *Nature*, **426**, 826–829.

12 Dickson, R. R., J. Meincke, S.A. Malmberg, and A. J. Lee, 1988: The “Great Salinity Anomaly” in

13 the northern North Atlantic, 1968–1982. *Prog. Oceanogr.*, **20**, 103–151.

14 Durack, P. J., and S. E. Wijffels, 2010: Fifty-year trends in global ocean salinities and their

15 relationship to broad-scale warming. *J. Climate*, **23**, 4342-4362.

16 Durack, P. J., S. E. Wijffels, and R. J. Matear, 2012: Ocean salinities reveal strong global water

17 cycle intensification during 1950 to 2000. *Science*, **336**, 455-458.

18 Fedorov, A. V., R. C. Pacanowski, S. G. Philander, and G. Boccaletti, 2004: The effect of salinity

19 on the wind-driven circulation and the thermal structure of the upper ocean. *J. Phys.*

20 *Oceanography*.34, 1949-1966.

21 Foltz, G. R., S. A. Grodsky, and J. A. Carton, 2004: Seasonal salt budget of the north western

22 tropical Atlantic Ocean along 38°W. *J. Geophys. Res.*, **109**, C03052,

23 doi:10.1029/2003JC002111.

24 Fujii, Y., T. Nakaegawa, S. Matsumoto, T. Yasuda, G. Yamanaka, and M. Kamachi, 2009: Coupled

25 Climate Simulation by Constraining Ocean Fields in a Coupled Model with Ocean Data. *J.*

26 *Climate*, **22**, 5541–5557.
1 Fujii, Y., M. Kamachi, S. Matsumoto, and S. Ishizaki, 2012: Barrier layer and relevant variability of the salinity field in the equatorial Pacific estimated in an ocean reanalysis experiment. *Pure Appl. Geophys.* **169** (3), 579-594.

4 Guinehut, S., A.-L. Dhomps, G. Larnicol and P.-Y. Le Traon, 2012: High resolution 3D temperature and salinity fields derived from in situ and satellite observations. *Ocean Sci.*, **8**, 845-857.

7 Hackert, E., J. Ballabrera-Poy, A. Busalacchi, R.H. Zhang, and R. Murtugudde, 2011: Impact of sea surface salinity assimilation on coupled forecasts in the tropical Pacific, *J. Geophys. Res.*, **116**, C05009, doi:10.1029/2010JC006708.

10 Held, I. M., and B. J. Soden, 2006: Robust responses of the hydrological cycle to global warming. *J. Climate*, **19**, 5686–5699.

12 Hernandez, F., L. Bertino, G. Brassington, E. Chassignet, J. Cummings, F. Davidson, M. Drévillon, G. Garric, M. Kamachi, J.-M. Lellouche, R. Mahdon, M.J. Martin, A. Ratsimandresy, and C. Regnier, 2009: Validation and intercomparison studies within GODAE. *Oceanography* **22**(3):128–143.

16 Huang, B., Y. Xue, and D.W. Behringer, 2008: Impacts of Argo salinity in NCEP global ocean data assimilation system: The tropical Indian Ocean. *J. Geophys. Res.*, **113**, C08002, doi:10.1029/2007JC004388

19 Ingleby, B., and M. Huddleston, 2007: Quality control of ocean temperature and salinity profiles - historical and real-time data. *J. Marine Sys.*, **65**, 158-175.

21 IPCC (Intergovernmental Panel on Climate Change) 2013, Climate change 2013: The physical science basis. Working Group I contribution to the IPCC Fifth Assessment Report. Cambridge, United Kingdom: Cambridge University Press (see www.ipcc.ch/report/ar5/wg1).

24 Janowiak, J. E., P. Bauer, W. Wang, P. A. Arkin, and J. Gottschalck, 2010: An Evaluation of Precipitation Forecasts from Operational Models and Reanalyses Including Precipitation Variations Associated with MJO Activity. *Mon. Wea. Rev.*, **138**, 4542–4560.
1 Johnson, E. S., G. S. E. Lagerloef, J. T. Gunn, and F. Bonjean, 2002: Salinity advection in the
2 tropical oceans compared to atmospheric forcing: A trial balance, *J. Geophys. Res.*, **107**, 8014.

3 Kim, J.-E., and M. J. Alexander, 2013: Tropical Precipitation Variability and Convectively Coupled
4 Equatorial Waves on Submonthly Time Scales in Reanalyses and TRMM. *J. Climate*, **26**,
5 3013–3030.

6 Lee, T., T. Awaji, M.A. Balmaseda, E. Greiner, and D. Stammer, 2009: Ocean state estimation for
7 climate research. *Oceanography*, **22**(3):160–167.

8 Levitus, S., J. I. Antonov, T. P. Boyer, R. A. Locarnini, H. E. Garcia, and A.V. Mishonov, 2009:
9 Global ocean heat content 1955–2008 in light of recently revealed instrumentation problems,
10 *Geophys. Res. Lett.*, **36**, L07608, doi:10.1029/2008GL037155

11 Levitus, S., J. I. Antonov, T. P. Boyer, O. K. Baranova, H. E. Garcia, R. A. Locarnini, A. V.
12 Mishonov, J. R. Reagan, D. Seidov, E. S. Yarosh, M. M. Zweng, 2012: World Ocean heat
13 content and thermosteric sea level change (0-2000 m) 1955-2010. *Geophys. Res. Lett.*, **39**,
14 L10603, doi:10.1029/2012GL051106.

15 O’Kane, T. J., R. J. Matear, M. A. Chamberlain, and P. R. Oke, 2014: ENSO regimes and the late
16 1970’s climate shift: The role of synoptic weather and South Pacific Ocean spiciness. *J. Comp.
17 *Phy.*, **271**, 19-38.

18 Maes, C., J. Picaut, and S. Belamari, 2005: Importance of salinity barrier layer for the buildup of El
19 Niño. *J. Climate*, **18**, 104–118.

20 Maes, C., K. Ando, T. Delcroix, W. S. Kessler, M. J. McPhaden, and D. Roemmich, 2006:
21 Observed correlation of surface salinity, temperature and barrier layer at the eastern edge of the
22 western Pacific warm pool. *Geophys. Res. Lett.*, **33**, L06601, doi: 10.1029/2005GL024772.

23 Marshall, J., A. Adcroft, C. Hill, L. Perelman, and C. Helsey, 1997: A finite-volume,
24 incompressible Navier–Stokes model for studies of the ocean on parallel computers. *J. Geophys.
25 *Res.*, **102**, 5753–5766
1 Masuda, S., and Co-authors, 2010: Simulated Rapid Warming of Abyssal North Pacific Waters.

2 *Science*, **329**, 319-322

3 Murtugudde, R., and A. J. Busalacchi 1998: Salinity effects in a tropical ocean model. *J. Geophys.*

4 *Res.*, **103**, 3282-3300.

5 Palmer, M., and Co-authors, 2015: Ocean heat content variability and change in an ensemble of
6 ocean reanalyses. Submitted to the Climate Dynamic

7 Rahmstorf, S., 1996: On the freshwater forcing and transport of the Atlantic thermohaline
8 circulation. *Clim. Dyn.*, **12**, 799-811

9 Santer, B.D., T.M.L. Wigley, and P.D. Jones, 1993: Correlation methods in fingerprint detection
10 studies. *Clim. Dyn.*, **8**, 265-276.

11 Sprintall, J., S. E. Wijffels, R. Molcard, and I. Jaya, 2009: Direct estimates of the Indonesian
12 Throughflow entering the Indian Ocean: 2004-2006. *J. Geophy. Res-Oce.*, **114**, C07001,
13 doi:10.1029/2008JC005257

14 Storto, A., S. Dobricic, S. Masina, and P. Di Pietro, 2011: Assimilating along-track altimetric
15 observations through local hydrostatic adjustments in a global ocean reanalysis system. *Mon.*
16 *Wea. Rev.*, **139**, 738-754.

17 Storto, A., and Co-authors, 2015: Steric sea level variability (1993-2010) in an ensemble of ocean
18 reanalyses and objective analyses. Submitted to the Climate Dynamic

19 Toyoda, T, Y. Fujii, T. Yasuda, N. Usui, T. Iwao, T. Kuragano, and M. Kamachi, 2013: Improved
20 analysis of the seasonal-interannual fields by a global ocean data assimilation system,
21 *Theoretical and Applied Mechanics Japan*, **61**, 31-48.

22 Vernieres, G., M. M. Rienecker, R. Kovach, and L. C. Keppenne, 2012: The GEOS-iODAS:
23 Description and evaluation. *Tech. Rep. TM-2012-104606*, NASA, National Aeronautics and
24 Space Administration Goddard Space Flight Center, Greenbelt, Md.
1 Vialard, J., P. Delecluse, and C. Menkes, 2002: A modelling study of salinity variability and its

2 effects in the tropical Pacific Ocean during the 1993-1999 period. *J. Geophys. Res.*, **107**, C12,

3 8005, doi:10.1029/2000JC000758

4 Vinje, T., 2001: Fram Strait ice fluxes and atmospheric circulation, 1950–2000. *J. Climate*. **14**,

5 3508–3517.

6 von Storch, H., and A. Navarra, 1999: Analysis of climate variability: applications of statistical

7 techniques. Published by Springer-Verlag Berlin.

8 Vranes, K., A. L. Gordon, and A. Ffield, 2002: The heat transport of the Indonesian throughflow

9 and implications for Indian Ocean heat budget. *Deep-Sea Research I II*, **49**, 1391-1410.

10 Wadley, M. R., G. R. Bigg, 2006: Are “Great Salinity Anomalies” advective? *J. Climate*, **19**, 1080-1088.

11

12 Wang, X., and Y. Chao, 2004: Simulated sea surface salinity variability in the tropical Pacific.

13 *Geophys. Res. Lett.*, **31**, L02302, doi:10.1029/2003GL01

14 Waters J, Martin M, While J, Lea D, Weaver A, and Mirouze I. 2014: Implementing a variational

15 data assimilation system in an operational 1/4 degree global ocean model. *Q. J. R. Meteorol. Soc.*

16 DOI: 10.1002/qj.2388

17 Xue Y., and Co-authors, 2011: An assessment of oceanic variability in the NCEP climate forecast

18 reanalysis. *Clim. Dyn.*, **37**, 2511-2539

19 Xue Y., and Coauthors, 2012: A Comparative Analysis of Upper-Ocean Heat Content Variability

20 from an Ensemble of Operational Ocean Reanalyses. *J. Climate*, **25**, 6905–6929.

21 Yang, S. C., M. Rienecker, and C. Keppenne, 2010: The impact of ocean data assimilation on

22 seasonal-to-interannual forecasts: A case study of the 2006 El Niño event. *J. Climate*, **23**, 4080–

23 4095.

24 Yin, Y., O. Alves, and P. R. Oke, 2011: An ensemble ocean data assimilation system for seasonal

25 prediction. *Mon. Wea. Rev.*, **139**, 786-808.
1 Zhao, M., H. H. Hendon, O. Alves, Y. Yin, and D. Anderson, 2013: Impact of salinity constraints on the simulated mean state and variability in a coupled seasonal forecast model. *Mon. Wea. Rev.*, **141**, 388-402.

4 Zhao, M., H. H. Hendon, O. Alves, and Y. Yin, 2014: Impact of improved assimilation of temperature and salinity for coupled model seasonal forecasts, *Clim. Dyn.*, **42**, 2565-2585

6 Zhang, R., and G. K. Vallis, 2006: Impact of great salinity anomalies on the low-frequency variability of the north Atlantic climate. *J. Climate*, **19**, 470-482.

8 Zhang S., M. J. Harrison, A. Rosati, and A. Wittenberg, 2007: System design and evaluation of coupled ensemble data assimilation for global oceanic studies. *Mon. Wea. Rev.*, **135**, 3541-3564.

10 Zhu, J., B. Huang, and M. Balmaseda, 2012: An ensemble estimation of the variability of upper-ocean heat content over the tropical Atlantic Ocean with multi-ocean reanalysis products. *Clim. Dyn.* **39**, 1001-1020.

13 Zhu, J., B. Huang, R.-H. Zhang, Z.-Z. Hu, A. Kumar, M. A. Balmaseda, L. Marx and J. L. Kinter III, 2014: Salinity anomaly as a trigger for ENSO events. *Nature, scientific report*, **4**:6821, doi: 10.1038/srep06821

16 Zuo, H., Balmaseda, M. A. and Mogensen, K., 2015. The new eddy-permitting ORAP5 ocean reanalysis: description, evaluation and uncertainties in climate signals. *Clim. Dyn.*, **10.1007/s00382-015-2675-1**

19
<table><caption>Table 1: List of the ocean data assimilation systems</caption><thead><tr><td>Products</td><td>Ocean Models</td><td>Forcing</td><td>Main Assim. Method</td><td>Data Assimilated</td><td>Relax to Clim.</td></tr></thead><tbody><tr><td>ECCOV4<sup>a)</sup><br>(MIT/AER/JPL))</td><td>MITGCM<br>1°×(0.4-1.0°); 50 vertical levels</td><td>ERAi + CORE Bulk</td><td>4DVAR Adjoint</td><td>T/S; SLA; SSH; SST</td><td>None</td></tr><tr><td>GMAO<sup>b)</sup><br>GSFC/NASA</td><td>MOM4p1<br>0.5°; 40 vertical levels</td><td>MERRA + Bulk</td><td>EnOI</td><td>T/S; SST;</td><td>SSS</td></tr><tr><td>SODA<sup>c)</sup><br>Uni. Maryland and TAMU</td><td>POP2.1<br>0.25°×0.4°; 40 vertical levels</td><td>ERA40 and ERAi (since 2002) + Bulk</td><td>OI</td><td>T/S; SST</td><td>SSS (3 months); 3D<br>T/S (10 years)</td></tr><tr><td>GloSea5<sup>d)</sup><br>Met Office, UK</td><td>NEMO3.2<br>1/4°; 75 vertical levels</td><td>ERAi + CORE Bulk</td><td>3DVAR</td><td>T/S; SST; SLA; SIC</td><td>Surface Haney<br>restoring + 3D T/S<br>(360 days)</td></tr><tr><td>ORAS3<sup>e)</sup><br>ECMWF</td><td>HOPE<br>1°×(0.3-1°); 29 vertical levels</td><td>ERA40 and ECMWF operational analysis<br>(since 2002)</td><td>OI</td><td>T/S; SST; SLA</td><td>SSS (3 years); 3D<br>T/S (10 years)</td></tr><tr><td>ORAS4<sup>f)</sup><br>ECMWF</td><td>NEMO3.0<br>1°×1°; 42 vertical levels</td><td>ERA40 to 1988, ERAi thereafter + Flux</td><td>3DVAR</td><td>T/S; SST; SLA</td><td>SSS (1 year); 3D<br>T/S (20 years)</td></tr><tr><td>CGLORS<sup>g)</sup><br>CMCC/IT</td><td>NEMO3.2.1-LIM2<br>0.5°×0.5°; 50 vertical levels</td><td>ERAi + CORE Bulk</td><td>3DVAR</td><td>T/S; SLA; SST; SIC</td><td>SSS (300 days); 3D<br>T/S</td></tr><tr><td>PEODAS<sup>h)</sup><br>BOM/AUS</td><td>MOM2<br>2°×(0.5-1.5°); 25 vertical levels</td><td>ERA40 and NCEP2 (since 2002)</td><td>EnKF</td><td>T/S; SST</td><td>SSS (1 year); 3D<br>T/S (2 years)</td></tr><tr><td>PECDAS<br>BOM/AUS</td><td>MOM2<br>2°×(0.5-1.5°); 25 vertical levels</td><td>Atmospheric component of the coupled model</td><td>EnKF</td><td>T/S; SST</td><td>SSS (1 year); 3D<br>T/S (2 years)</td></tr><tr><td>ECDA<sup>i)</sup><br>GFDL/NOAA</td><td>CM2.1/MOM4<br>1°×(0.3-1°); 50 vertical levels</td><td>Atmospheric component of the coupled model</td><td>EnKF</td><td>T/S; SST</td><td>None</td></tr><tr><td>MOVE-C<sup>j)</sup><br>MRI/JMA</td><td>MRI.COM V2<br>1°×(0.3-1.0°); 50 vertical levels</td><td>Atmospheric component of the coupled model</td><td>MOVE 3DVAR-IAU</td><td>T/S; SSH; SST</td><td>SSS (10 days); 3D<br>T/S (5 years)</td></tr><tr><td>MOVE-G2<sup>k)</sup><br>(MRI/JMA)</td><td>MRI.COM V3<br>1°×(0.3-0.5°); 52 vertical levels</td><td>JRA-55 CORE Bulk</td><td>MOVE 3DVAR-FGAT-IAU</td><td>T/S; SSH; SST</td><td>3D T/S (5 years)</td></tr><tr><td>K7ODA (ESTOC)<sup>l)</sup><br>(JAMSTEC/RCGC)</td><td>MOM3-based OGCM<br>1°×1°; 45 vertical levels</td><td>NCEP-R1 + Corr. Flux</td><td>4DVAR Adjoint</td><td>T/S; SLA; SST</td><td>BC (30 days by IAU)</td></tr><tr><td>G2V3<br>Mercator Océan/FR</td><td>NEMO3.1-LIM2evp<br>1/4°; 75 vertical levels</td><td>ERAi Corr + CORE Bulk</td><td>SAM2V1 - SEEK filter + 3DVAR</td><td>T/S; SLA; SST; SIC</td><td>3D T/S (360 days)</td></tr><tr><td>ARMOR3D<sup>m)</sup><br>CLS/FR</td><td>No Model<br>1/3°; 24 vertical levels</td><td>N/A</td><td>Objective Analysis</td><td>T/S; SLA; SST</td><td>None</td></tr></tbody></table>
<table><tr><td>ISAS13<br>LPO/Ifremer</td><td>No Model<br>0.5°; 152 vertical levels</td><td>N/A</td><td>Objective Analysis</td><td>T/S</td><td>None</td></tr><tr><td>EN3v2a<sup>0)</sup><br>Met Office, UK</td><td>No Model<br>1°×1°; 42 vertical levels</td><td>N/A</td><td>Objective Analysis</td><td>T/S</td><td>None</td></tr></table>

Note: **T/S**: Temperature/Salinity profiles; **SSH**: sea surface height from tide gauges; **SLA**: altimeter-derived sea level anomalies from satellite observation; **SIC**: satellite derived sea ice concentration. **BC**: lateral boundary conditions **IAU**: incremental analysis updates a) Marshall et al. 1997; b) Vernieres et al. 2012; c) Carton and Giese 2008; d) Water et al. 2014; e) Balmaseda et al. 2008; f) Balmaseda et al. 2013; g) Storto et al. 2011; h) Yin et al. 2011; i) Zhang et al. 2007; j) Fujii et al. 2009; k) Toyoda et al. 2013; l) Masuda et al. 2010; m) Guinehut et al. 2012; o) Ingleby and Huddleston 2007
<caption>Fig. 1a</caption>

Diff AnnualMean S700 EMORA-EMOO

<figure><img src="image_2.png" /><figcaption>Fig. 1b</figcaption></figure>

AnnualMean S700 EMOO

<figure><img src="image_3.png" /><figcaption>Fig. 1c</figcaption></figure>

Ensemble Spread AnnualMean S700

<figure><img src="image_4.png" /></figure>

Fig. 1: (a) Differences of the annual mean S700 (AMS700) between the ensemble mean of ORAs (EMORA) and the ensemble mean of 3 OOAs (EMOO) for the period 1993-2010. (b) The distribution of AMS700 from EMOO for the period 1993-2010. (c) The ensemble spread (i.e. SPD<sup>AM</sup><sub>EMORA</sub>) of AMS700 from individual ORAs about the corresponding EMORA. The unit of colour bar is psu.
<figure><img src="image_5.png" /><figcaption>**Fig. 2:** (a) Zonal distributions of differences of meridionally-averaged (over 30°N-60°N) AMS700 between individual reanalysis systems and EMOO for the period 1993-2010. The definition of shaded band (i.e. *UCR<sub>EMORA</sub><sup>AM</sup>*) can be referred to the Eq. (A4). b) Same as in (a), except for meridionally-averaged over 15°S-15°N. (c) Same as in (a), except for meridionally-averaged over 60°S-30°S. The unit of ordinate is *psu*.</figcaption></figure>
<caption>Fig. 3</caption>

Diff AnnualMean S700 ReSys-TAO at 8N156E

<figure><img src="image_6.png" /></figure>

Fig. 3: Differences of the seasonal cycle (i.e. monthly climatology for the period 1999-2010) S700 between all products and the TAO/TRITON data at 8°N156°E. The definition of shaded band (i.e. UCR<sub>EMORA</sub><sup>SC</sup> ) can be referred to the Eq. (A4).

The unit of ordinate is *psu*.
<figure><img src="image_7.png" /><figcaption>Fig. 4a</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 4b</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 4c</figcaption></figure>

**Fig. 4:** (a) Zonal distribution of meridionally-averaged (over 30°N-60°N) temporal standard deviation of S700 anomalies (i.e. $STD_n^A$) from individual products for the period 1993-2010. The definition of shaded band (i.e. $UCR_{EMORA}^{STD}$) can be referred to the Eq. (A4). (b) Same as in (a), except for meridionally-averaged over 15°S-15°N. (c) Same as in (a), except for meridionally-averaged over 60°S-30°S. The unit of ordinate is *psu*.
<caption>Fig. 5a</caption>

Ensemble Spread of Anom S700 1993-2010

<figure><img src="image_10.png" /><figcaption>Fig. 5b</figcaption></figure>

Standard Deviation of Anom S700 EMORA 1993-2010

<figure><img src="image_11.png" /><figcaption>Fig. 5c</figcaption></figure>

Ratio Signal-to-Noise Anom S700 1993-2010

<figure><img src="image_12.png" /></figure>

**Fig. 5:** (a) Distribution of ensemble spread of S700 anomalies (i.e. *SPD*<sup>A</sup><sub>EMORA</sub>) from 14 ORAs about the corresponding EMORA for the period 1993-2010. (b) The temporal standard deviation of S700 anomalies from EMORA (i.e. *STD*<sup>A</sup><sub>EMORA</sub>) for the period 1993-2010. The unit of colour bar in (a)-(b) is *psu*. (c) The distribution of the estimated Signal-to-Noise Ratio (SNR) of S700 anomalies from 14 ORAs for the period 1993-2010.
<figure><img src="image_13.png" /><figcaption>Fig. 6</figcaption></figure>

**Fig. 6:** (a) Temporal correlation coefficients of S700 anomalies between EMORA and EMOO for the period 1993-2010. (b) Ensemble spread (i.e. *SPD*<sub>EMORA</sub><sup>COR</sup>) of correlation coefficients of S700 anomalies between individual ORAs and EMOO about their mean correlation coefficients.
<figure><img src="image_14.png" /><figcaption>Fig. 7b</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Fig. 7c</figcaption></figure>

<figure><img src="image_16.png" /></figure>

**Fig. 7:** (a) Zonal distribution of meridionally-averaged (over 30°N-60°N) correlation coefficients of S700 anomalies between individual ORAs and EMOO for the period 1993-2010. The definition of shaded band (i.e. *UCR<sub>EMORA</sub><sup>COR</sup>*) can be referred to the Eq. (A4). (b) Same as in (a), except for meridionally-averaged over 15°S-15°N. (c) Same as in (a), except for meridionally-averaged over 60°S-30°S.
<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

**Fig. 8:** (a) Temporal correlation coefficients between S700 anomalies and T700 anomalies from EMORA for the period 1993-2010. (b) Same as in (a), except for EMOO. (c) Ensemble spread (i.e. *SPD*<sub>EMORA</sub><sup>COR</sup>) of T700-S700 anomalies correlation coefficients from individual ORAs about their mean correlation coefficients.
<figure><img src="image_20.png" /><figcaption>Fig. 9a</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>Fig. 9b</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>Fig. 9c</figcaption></figure>

**Fig. 9:** Same as in Fig. 7, except for correlation coefficients between the T700-S700 anomalies.
<figure><img src="image_23.png" /></figure>

<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

**Fig. 10:** (a) Centred pattern correlation coefficients (i.e. CPCOR; calculated over the band area 0-360°E; 30°N-60°N) of salt anomalies between EMORA and EMOO as a function of depth (0-1500m) and time (1993-2010. Prior to plotting a 7-month running mean was applied on the computed correlation coefficients to remove the intra-seasonal variability. The ordinate has units meter (m). (b) Same as in (a), except for the band area 0-360°E; 15°S-15°N. (c) Same as in (a), except for the band area 0-360°E; 60°S-30°S.
<figure><img src="image_26.png" /><figcaption>Fig. 11</figcaption></figure>

**Fig.11:** (a) Evolution of global averaged (over 0-360°E; 60°S-60°N) and depth-averaged (within 0-300m ocean layer) salinity anomalies from all products for the period 1993-2010. The definition of shaded band (i.e. UCR<sub>EMORA</sub><sup>A</sup>) can be referred to the Eq. (A5). The unit of ordinate is psu. (b) - (c) Same as in (a), except for depth-averaged salinity anomalies within 300-700 meters layer and 700-1500 meters layer, respectively. Prior to plotting a 7-month running mean was applied to remove the intra-seasonal variability.