Yuanxin LIU, Lijing CHENG, Yuying PAN, Zhetao TAN, John ABRAHAM,Bin ZHANG, Jiang ZHU, and Junqiang SONG

1College of Meteorology and Oceanography, College of Computer Science and Technology,National University of Defense Technology, Changsha 410073, China

2Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

3Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China

4University of Chinese Academy of Sciences, Beijing 100049, China

5School of Engineering, University of St. Thomas, 2115 Summit Ave., St Paul, MN 55105, USA

6Marine Science Data Center, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China

Fig. S1. Global climatology of SSS for individual CMIP5 models. The black countour line shows the 35 g kg-1 salinity.

Fig. S2. Global climatology of SSS for individual CMIP6 models. The black countour line shows the 35 g kg-1 salinity.

Fig. S3. As in Fig. 1 but using a subset of models (20 models) from the same institutes (Table S1).

Fig. S4. As in Fig. 1 but using 1980-2000 climatology.

Fig. S5. Global climatology of S1000 for all individual CMIP5 models. The black contour shows the 35 g kg-1 salinity.

Fig. S6. Global climatology of S1000 for all individual CMIP6 models. The black contour shows the 35 g kg-1 salinity.

Fig. S7. As in Fig. 2 but using a subset of models (20 models) from the same institutes (Table S1).

Fig. S8. As in Fig. 2 but using 1980-2000 climatology.

Fig. S9. Same as Fig. 11 but for SSS change within 60°S-60°N.

Table S1. Correspondence of CMIP5 models with CMIP6 models from the same institutes.