Eddy-resolving Simulation of CAS-LICOM3 for Phase 2 of the Ocean Model Intercomparison Project

A 61-year (1958–2018) global eddy-resolving dataset for phase 2 of the Ocean Model Intercomparison Project has been produced by the version 3 of Chinese Academy of Science, the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics (LASG/IAP) Climate system Ocean Model (CAS-LICOM3). The monthly and a part of the surface daily data in this study can be accessed on the Earth System Grid Federation (ESGF) node. Besides the details of the model and experiments, the evolutions and spatial patterns of large-scale and mesoscale features are also presented. The mesoscale features are reproduced well in the high-resolution simulation, as the mesoscale activities can contribute up to 50% of the total SST variability in eddy-rich regions. Also, the large-scale circulations are remarkably improved compared with the low-resolution simulation, such as the climatological annual mean SST (the RMSE is reduced from 0.59°C to 0.47°C, globally) and the evolution of Atlantic Meridional Overturning Circulation. The preliminary evaluation also indicates that there are systematic biases in the salinity, the separation location of the western boundary currents, and the magnitude of eddy kinetic energy. All these biases are worthy of further investigation.


Background and Summary
To understand the origin of model biases in phase 6 of the Coupled Model Intercomparison Project (CMIP6), the Ocean Model Intercomparison Project (OMIP; Griffies et al., 2016) was proposed to trace the coupled model biases in ocean-sea-ice models. Recently, the protocol of OMIP has been updated and two kinds of OMIP experiments forced by two different datasets have been proposed. Phase 1 of OMIP (OMIP-1) is forced by the CORE-II dataset (Large and Yeager, 2009), mainly derived from the NCEP-NCAR atmospheric reanalysis (Kalnay et al., 1996;Kistler et al., 2001), and phase 2 of OMIP (OMIP-2) is forced by the JRA55-do dataset (Tsujino et al., 2018), based on the Japanese 55-year Reanalysis (JRA-55; Kobayashi et al., 2015) and river runoff.
The results of OMIP-1 and OMIP-2 from 11 state-ofthe-art global low-resolution ocean-sea-ice models were evaluated by Tsujino et al. (2020). Some improvements to key features were identified from the comparison, such as the attenuation of the warming biases off the eastern coast of the Pacific Ocean, as well the reproduction of the observed global warming and warming hiatus in the sea surface temperature (SST). However, some common model biases still persist, indicating errors in representing important processes or biases in the surface forcing. Among the possible causes, the horizontal resolution of the model is considered to be one of the important factors, particularly for the biases in eddy-rich regions, western boundary currents, and narrow straits. Therefore, additional OMIP-2 experiments using global eddy-resolving ocean-sea-ice models were proposed and compared by Chassignet et al. (2020) to identify the robust improvements of ocean-sea-ice simulation by refining the model resolution from 100 km to 10 km.
Under the framework of this resolution-related comparison, four ocean modeling groups have conducted and submitted a pair of experiments: 100 km and 10 km versions of the model forced by the JRA55-do dataset for only one cycle . The experiments of the coarse resolution are just the same as the first cycle of OMIP-2 in Tsujino et al. (2020). The Chinese Academy of Science, State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics (LASG/IAP) Climate system Ocean Model, version 3 (CAS-LICOM3) is one of the models that has participated in the comparisons, which was developed based on the previous version, CAS-LICOM2 . CAS-LICOM3 is also the ocean component of FGOALS3, participating in CMIP6. The eddy-resolving simulation of CAS-LICOM3 for OMIP-2 (hereafter called CAS-LICOM3_ High) was finished in November 2019 and the data submitted to the Earth System Grid Federation (ESGF) data server (https://esgf-node.llnl.gov/projects/cmip6/). Here, we document detailed information on the model configurations and experiments, as well as the results of diagnostic variables from CAS-LICOM3_High, to provide a general description of this dataset for a variety of users.
The organization of the paper is as follows. In section 2, we describe the ocean model, the experiment design, and the methods used here. Section 3 presents the basic technical validation of CAS-LICOM3_High. In this section, firstly, the temporal evolution of SST, sea surface salinity (SSS), Atlantic Meridional Overturning Circulation (AMOC), kinetic energy (KE), and upper-ocean temperature are presented for examining the trends and variabilities of the simulation. Secondly, the simulated spatial pattern and standard deviation (STD) of sea surface height (SSH) and SST are validated by splitting the large-scale signal and mesoscale eddies. Thirdly, the ocean surface currents are provided for showing the major currents, such as the western boundary currents and Antarctic Circumpolar Current (ACC). Finally, the North Equatorial Undercurrent (NEUC) is selected to evaluate this dataset, because it is a good indicator to examine the long-term and high-resolution simulation . In section 4, the data record is described. Section 5 provides usage notes.
2. Model, experiment and methods 2.1. Introduction to the model CAS-LICOM is a global ocean general circulation model developed by LASG/IAP of the Chinese Academy of Sciences (Zhang and Liang, 1989;Liu et al., 2004Liu et al., , 2012Yu et al., 2018). CAS-LICOM3 is an ocean general circulation model with a free sea surface and Arakawa B grid, and uses the primitive equations with Boussinesq and hydro-static approximations. It has been extensively upgraded since CAS-LICOM2. First, the coordinates of the dynamical core have been replaced by arbitrary orthogonal curvilinear coordinates (Yu et al., 2018). Therefore, the tripolar grid from Murray (1996) can be applied in CAS-LICOM3, with two North Poles, on the Eurasian (55°N, 95°E) and North American (55°N, 85°W) continents, respectively.
Second, the coupler was updated from NCAR flux coupler 6 to coupler 7 when the ocean component was coupled with the Community Ice Code, version 4 (CICE4). Also, it has been proved to be helpful for the high-resolution coupling (Lin et al., 2016).
Third, the central difference scheme is used in the momentum equations and the Leapfrog format with an Asselin filter is used for the time integration of the momentum equations. Additionally, the two-step preserving shape advection scheme (Yu, 1994;Xiao, 2006) and implicit vertical viscosity/diffusivity (Yu et al., 2018) were adopted for the tracer equations.
Fourth, with regard to the physical processes, the St. Laurent et al. (2002) internal tidal mixing scheme has been introduced into CAS-LICOM3 (Yu et al., 2017). The vertical viscosity and diffusion coefficients in the mixed layer are computed by the scheme of Canuto et al. (2001Canuto et al. ( , 2002, with background values of 2 × 10 −6 m 2 s −1 and upper limit of 2 × 10 −2 m 2 s −1 . Besides, the chlorophyll-a dependent solar penetration of Ohlmann (2003) introduced by Lin et al. (2007) was also implemented in CAS-LICOM3.

High-resolution experiment design
The eddy-resolving simulation of CAS-LICOM3 (CAS-LICOM3_High) is forced by JRA55-do data (Tsujino et al., 2018), which were developed based on the JRA-55 product (Kobayashi et al., 2015). The temporal resolution of JRA55-do is 3-h and the horizontal resolution is 0.5625°. The atmospheric state variables and radiative fluxes are employed to compute the net surface flux into the ocean-sea-ice system in JRA55-do data, which included the air temperature at 10 m, the air density at 10 m, the relative humidity at 10 m, the surface wind vectors at 10 m, the sea surface pressure, the surface downward and upward shortwave flux, and the surface downward longwave flux (here, flux into the ocean is called "downward"). The turbulence heat and momentum fluxes of JRA55-do data have been computed using the bulk formula of Large and Yeager (2004). In the formula, the relative 10 m winds are employed and they are obtained by subtracting full ocean currents from the 10 m winds. The freshwater flux includes the precipitation (snow and rainfall) from JRA-55, as well as the continental runoff, which incorporates the runoff of ice sheets and glaciers from Greenland and Antarctica. The simulation is forced every 6 hours, even though the JRA55-do data are three-hourly, and integrated for 61 years from 1958 to 2018.
The experiment follows the OMIP-2 protocol, with the initial condition of the temperature and salinity from observation and a state of rest. The initial values in this study are different from the standard OMIP-2 experiment; they are the temperature and salinity on 1 January 2014 from the Mercator Ocean analysis (Lellouche et al., 2018;Artana et al., 2019), not the World Ocean Atlas Locarnini et al, 2006Locarnini et al, , 2010. Meanwhile, the surface salinity is restored to the monthly Polar Science Center Hydrographic Climatology, version 3 (PHC3.0; Steele et al., 2001) over the entire domain with a salinity piston velocity of 50 m (4 yr) −1 [plus 50 m (30 d) −1 for the sea ice regions]. η The horizontal grid contains 3600 × 2302 points with a resolution of approximately 1/10°, which is 11 km zonally and varies from 11 km at the equator to 8 km in midlatitudes and 2.7 km around the Antarctic. The higher resolution of CAS-LICOM3_High than the forcing data may imply that eddy-resolving activities come from the internally generated variability. The -coordinate (Mesinger and Janjic, 1985) with 55 levels is applied in the vertical grid. The biharmonic scheme with a horizontal viscosity coefficient of −2.8 × 10 10 m 4 s −1 is chosen in CAS-LICOM3_ High. The full ocean surface velocity is used in the calculation of wind stress (relative wind stress), which can affect the strength of the North Equatorial Countercurrent (Sun et al., 2019(Sun et al., , 2020.

Methods
To avoid the interference of the long-term trend and seasonal cycle when analyzing the temporal variations, the least-squares linear trend is first removed. Then, the de-seasonalized method is applied. The de-seasonalized method here means the removal of the climatological monthly mean from the detrended time series.
The STD here is calculated after removing the leastsquares linear trend and the seasonal cycle. Also, the data from CAS-LICOM3_High are interpolated into a 25-km resolution when comparing with the observations or reanalysis data.
The spatial correlation coefficient calculated in this study is the Pearson product-moment coefficient of linear correlation between two datasets. For the Pearson correlation coefficient, the linear change of the two variables will not change the value of Pearson correlation. That is, high correlation does not mean two variables are exactly the same; rather, that they have the same spatial gradient.
The root-mean-square error (RMSE) of SST calculated in this study is also based on the model data interpolated to the observation grid. The formula is: where A i stands for the area of each grid cell; and SST model and SST obs are the simulated and observed SST, respectively.
A spatial filter method is used here, which is adopted from Bryan et al. (2010) and Lin et al. (2019), to extract the oceanic mesoscale signal. First, a filter box with 3° of longitude and 3° of latitude is used for monthly mean data to obtain the box-mean value as the lowpass-filtered value.
The value is considered as the large-scale signal. Second, the STDs of both daily SST and the lowpass-filtered monthly SST are calculated. Third, the difference between those two STDs is used to reflect the role of the mesoscale signal. In that way, the impact of the interannual variability of large-scale signals can be removed.
The eddy kinetic energy (EKE) computed here is the square of the temporal anomaly of the velocity ( and ), which is calculated based on the daily anomaly of SSH ( ). The formula is as follows: where is gravity, is the Coriolis parameter, and is the daily anomaly of SSH. is calculated by subtracting the climatological annual mean of the dataset.

Temporal evolution
The temporal evolutions of the annual mean globally averaged SST, SSS, AMOC, KE and temperature anomaly over the Niño3.4 index region from CAS-LICOM3_High are evaluated against the observation, reanalysis and the low-resolution experiment (CAS-LICOM3_Low; see details in Appendix A) (Fig. 1). The temporal evolution of the annual globally averaged SST is captured well by CAS-LICOM3_High (Fig. 1a). Here, three validated datasets have been applied, including the Extended Reconstructed SST, version 5 (ERSST.v5; resolution of 2°; Huang et al., 2017), Optimum Interpolation Sea Surface Temperature, version 2 (OISSTv2; resolution of 1/4°; Banzon et al., 2016;Reynolds et al., 2007), and SST from IAP (IAP_T; resolution of 1°; Cheng et al., 2015). An obvious warming trend of about 0.11°C (10 yr) −1 can be seen during 1958-2018, and 0.13°C (10 yr) −1 during 1982-2018, in CAS-LICOM3_ High. The trends during the two periods are larger than those from ERSST.v5 and slightly smaller than those from OISSTv2. Compared with OISSTv2, the trends in CAS-LICOM3_High are better than those in CAS-LICOM3_ Low, which may be related to the better simulation of both large-scale currents and mesoscale eddies in CAS-LICOM3_High. The STDs range from 0.06°C to 0.08°C for all simulations and observations, and the correlations are all higher than 0.9, as listed in Table 1. This suggests both the magnitude and phase of interannual-decadal variability in CAS-LICOM3_High match the observed one very well. Figure 1b shows the time series of the global mean annual SSS from CAS-LICOM3_High, Argo data (http://doi. org/10.17882/42182), and salinity data from IAP (IAP_S; Cheng et al., 2015). The values of the trends and STDs of the annual mean time series in CAS-LICOM3_High, Argo and IAP_S are all small (close to 0). The correlation coefficients of SSS between simulation and observation for both high-and low-resolution results are much lower than that of SST (Table 1), while the correlation coefficients between the two simulations are much higher. This indicates the large uncertainties in the surface freshwater fluxes. The mean values of SSS for the two experiments are smaller than those of observations, suggesting differences between the observational datasets and the model initial conditions.
The temporal evolution of AMOC is displayed in Fig. 1c. Compared with observation from RAPID (RAPID/MOCHA/ WBTS array; Cunningham et al., 2007), CAS-LICOM3_ High has a better simulation of AMOC evolution than CAS-LICOM3_Low. Although we cannot exactly figure out the reason for the improvement in CAS-LICOM3_High yet, the refined horizontal resolution does contribute to improving the simulation of the ocean large-scale circulations. The improvement of the mesoscale eddy, the location of the Gulf Stream and the deep convection area, and the heat and salinity budget in the upper mixed layer, are all possibly responsible for the improvement in the AMOC variation. Figure 1d shows the evolution of the global mean KE from CAS-LICOM3_High during the whole integration period. The global KE increases quickly during first two years, and then decreases slowly in the next five or ten years. It finally keeps a steady state at values of around 15 cm 2 s −2 in the following years. Therefore, the 61-year integration can reach a steady state for the global mean KE, but it may be not long enough for the deep circulation. Figure 1e displays the time-depth section and STD profiles of the annual mean upper 300-m temperature anomaly averaged over the Niño3.4 region (5°N-5°S, 170°-120°W) from CAS-LICOM3_High and EN4.2.1 objective analyses ( Good et al., 2013). The close match of the temperature anomaly and the STD profiles between those two datasets indicates a good simulation of the interannual variability by CAS-LICOM3_High. The STDs of the annual mean SST anomaly over the Niño3.4 region are 0.65°C for CAS-LICOM3_High and 0.62°C for EN4.2.1. The temporal correlation coefficient of the annual mean SST anomaly over the Niño3.4 region between EN4.2.1 and CAS-LICOM3_High is significantly high (0.93). All these results indicate CAS-LICOM3_High can reproduce El Niño and La Niña events well. In addition, the warm temperature anomalies appear more frequently below the depth of 100 m before 1980 than they do after 1980. This phenomenon might be deserving of further research.

SSH and STD
The simulated spatial pattern of large-scale SSH from CAS-LICOM3_High is almost identical to the Archiving, Validation and Interpretation of Satellite Oceanographic Data (AVISO; http://www.aviso.altimetry.fr/), with a spatial correlation coefficient of 0.99 (Figs. 2a and b). In both CAS-LICOM3_High and AVISO, the high SSHs locate around the western boundary, and the low ones close to the eastern boundary, in the low and middle latitudes of the Pacific and Atlantic Ocean. The west-east contrast is associated with the easterly trade winds and the westward ocean currents. In the Indian Ocean, high SSH values are located in the east and low ones are in the west. In the high latitudes, the relatively low SSHs are located in the subpolar gyre in the Northern Hemisphere and south to the ACC in the Southern Ocean. This high spatial correlation between CAS-LICOM3_High and AVISO is expected, as the spatial correlation coefficient between AVISO and CAS-LICOM3_Low is also high.
Apart from the large-scale circulation, mesoscale eddies can also be reproduced in CAS-LICOM3_High. The amplitudes of mesoscale eddies are shown through the STD of the daily SSH anomaly from both AVISO and CAS-LICOM3_High (Figs. 2c and d). The high values of STD are presented in the western boundary currents (e.g., Kuroshio extension, Gulf Stream) and ACC, both for CAS-LICOM3_High and AVISO. There are also large STDs in the western and central tropical Pacific, which may be caused by the tropical instability waves (TIWs). The spatial correlation coefficient of the STD between AVISO and CAS-LICOM3_High is 0.93, which may be related to the eddy-resolving horizontal grid spacing of CAS-LICOM3_ High.
Although the present resolution can resolve these western boundary currents well, it is still a challenge to simulate a correct separation owing to the sensitivity to the magnitude and formula of viscosity (Chassignet and Marshall, 2008). Thus far, we believe that the wrong separation location is related to the large viscosity in CAS-LICOM3_High. The biharmonic viscosity equation is employed with a viscosity coefficient of −2.8 × 10 10 m 4 s −1 in the present experiment. In addition to the explicit viscosity, the Euler backward differences in the two-step shape preservation advection will provide additional numerical viscosity. Generally, the simulated amplitudes of eddy activities are underestimated compared with the observed ones, which is also shown and discussed in the time series of EKE later (Fig. 2e).
The temporal evolution of simulated EKE from CAS-LICOM3_High is compared with the value from AVISO, which has a horizontal resolution of 25 km, in Fig. 1e. Because the value from AVISO is available after October 1992, the statistics for both 1993-2018 and 1958-2018 are chosen and computed. The global annual mean EKE from CAS-LICOM3_High is underestimated by about 15%-20% compared with AVISO, which may be caused by the relat- Table 1. Climatological annual mean (units: °C, psu, cm 2 s −2 ), STD (units: °C, psu, cm 2 s −2 ) and linear trend [units: °C (100 yr) −1 , psu (10 yr) −1 , cm 2 s −2 (10 yr) −1 ] of SST, SSS and EKE during 1958-2018 from CAS-LICOM3_High, CAS-LICOM3_Low, and multiple observations. The temporal correlation coefficients of these variables between models and observations during the full period are also shown. The values in brackets represent the values during different periods. SST 1982, SSS 2005and EKE 1993  ively high viscosity or high numerical damping effects in CAS-LICOM3_High. The linear trend is overestimated during 1993-2018, with 6.15 cm 2 s −2 (10 yr) −1 in CAS-LICOM3_High and 1.87 cm 2 s −2 (10 yr) −1 in AVISO. The variability of observed EKE can be captured by CAS-LICOM3_High, with a temporal correlation coefficient of 0.53 during the same period. The magnitude of the variability of the simulated annual EKE is underestimated in CAS-LICOM3_High, as the STD is about 3.83 cm 2 s −2 for AVISO and 2.42 cm 2 s −2 for CAS-LICOM3_High.

SST and STD
The large-scale spatial features of SST are also simulated well by CAS-LICOM3_High (Figs. 3a and b), with a spatial correlation coefficient close to 1.0 between CAS-LICOM3_High and OISSTv2. We find that the meridional contrast of warm SST in the low latitudes and cold SST in the high latitudes, as well as the zonal contrast of the warm pool [enclosed by the 28°C contour (thick blue line) in Fig.  3a] in the western Pacific and cold tongue in the eastern Pacific, are represented well in CAS-LICOM3_High. Also, the global RMSE is 0.47°C, which is smaller than that of CAS-LICOM3_low (0.59°C). The improvements of the SST simulation in these regions can mostly be attributed to the better simulation of the eddy transports in CAS-LICOM3_High than those in CAS-LICOM3_Low. For the Labrador Sea, the reduced bias may be related to convection processes, which have a correct location in CAS-LICOM3_High. However, in CAS-LICOM3_High, there are still larger biases in the regions over the western boundary, which may be caused by the bias of the simulation of the western boundary jet axes, as shown in Figs. 2 and 4.
The de-seasonalized STDs of the daily SST anomaly from OISSTv2 and CAS-LICOM3_High are shown in Figs. 3c and d. The spatial pattern of STDs is simulated well, with a global spatial correlation coefficient of 0.98 between CAS-LICOM3_High and OISSTv2. The large STD values (> 1°C) are located in the equatorial Pacific, western boundary currents and their extensions, the Argentine Basin, and the Southern Ocean between the Indian Ocean and Atlantic Ocean. In the equatorial Pacific, the large STD is mainly due to the interannual variability associated with ENSO, and partly due to the intraseasonal variability related to eddy-induced TIWs. However, there are some underestimations of the simulated STD (Fig. 3d) over the warm SST regions (> 26°C contour lines in Fig. 3a).
Figures 3e and f show the difference between the daily SST STD and 3° × 3° lowpass-filtered monthly SST STD (ΔSTD). The difference indicates the variability due to the highpass signal related to mesoscale activities. As Figs. 3c and d show in the observation and CAS-LICOM3_High, large daily SST STDs are located in the eddy-rich regions, such as the western boundary and ACC. In those places, ΔSTD can reach 1°C (Figs. 3e and f). These large ΔSTD values are attributable to mesoscale eddies, which means the mesoscale signal can contribute to 50% of the total STD of daily SST (Figs. 3e and f). Meanwhile, in the eastern part of the Pacific, mesoscale eddies can contribute to 20% of the STD in the observation and relatively less in CAS-LICOM3_ High.

Surface ocean currents
The global large-scale surface current velocity during 1993-2018 for both the Ocean Surface Current Analyses Real-time (OSCAR; http://podaac.jpl.nasa.gov) dataset and CAS-LICOM3_High are shown in Fig. 4, with two western boundary regions zoomed in. The high speeds of the western boundary currents (such as the Kuroshio, Gulf Stream, North Brazil Current, etc.), equatorial current systems (including the South Equatorial Current, North Equatorial Countercurrent) and ACC, can be better represented in CAS-LICOM3_High than in CAS-LICOM3_Low.
For the western boundary currents, the simulated Kuroshio path flows along the shelf break inside the East China Sea between Taiwan Island and Yonaguni-jima Island, which is also proven by observations and other numerical simulations (Johns et al., 2001;Yang et al., 2012). The Gulf Stream in CAS-LICOM3_High can extend to 40°W, but it does not extend as a coherent feature past the New England seamounts, which is also exhibited in the 1/25° simulation from Chassignet and Xu (2017). However, there is a bias about the western boundary currents, as the separation of both the Kuroshio and Gulf Stream shifts northward by about 2° of latitude, compared with OSCAR.

NEUC
The zonal belt structure of the zonal currents are important features for the global ocean, which also cannot be simulated by coarse-resolution models. As discussed in Li et al. (2018), the simulation of the NEUC requires both a high-resolution and long-term simulation, thus meaning the NEUC is a good metric to evaluate high-resolution and long-term simulations. As presented in Fig. 5, the meridional sections of zonal velocity averaged over 135°-140°E, 175°-180°E, and 150°-145°W, from 11-year mean Argo data (Roemmich and Gilson, 2009) and CAS-LICOM3_High, show three branches of NEUCs (Qiu et al., 2013), in all three meridional sections. The Argo zonal velocities are estimated from the density gradients calculated based on the observed temperature and salinity. The upper boundary of the three NEUC jets for CAS-LICOM3_High become deeper with increasing latitude and shallower from west to east, which is exactly the same as in Argo. Although the locations of the NEUC jets are reproduced very well by CAS-LICOM3_ High, their strength in CAS-LICOM3_High is weaker than in Argo, which may be related to the non-geostrophic velocity  or the relatively large viscosity in the model.
In summary, CAS-LICOM3_High can reproduce well the long-term trend and interannual-decadal variability of various variables. The refined resolution simulation can not only capture the mesoscale eddies, the western boundary currents, and the zonal jets, but also improve the representation of large-scale circulations. In addition, the mesoscale activities can affect the magnitude of the variabilities, as the mesoscale signal is responsible for up to 50% of the variability of daily SST in eddy-rich regions. The processes related to these variabilities can be further investigated using this dataset and the added value of the high resolution can also be evaluated by comparing with the low-resolution simulation, particularly where the biases have been significantly reduced. In addition, this preliminary evaluation also indicates some systematic biases in the global temperature and salinity, the EKE, and regional current patterns. We will further improve the model in the future.

Data records
An eddy-resolving dataset of CAS-LICOM3_High for OMIP-2 has been uploaded to the ESGF node and can be found at https://esgf-node.llnl.gov/projects/cmip6/. Although the model outputs are unformatted with big-endianness and double precision, we converted all the variables into NetCDF, version 4, and single precision for analysis. These data can be easily handled by common computer programming languages and professional software, such as Climate Data Operators (https://code.mpimet.mpg.de/projects/cdo/) or NetCDF Operator (http://nco.sourceforge.net).

Usage notes
The original model outputs are on a tripolar grid with two poles at (55°N, 95°E) and (55°N, 85°W) in the Northern Hemisphere continent. There are 3600 and 2302 grid numbers in the zonal and meridional directions of the horizontal domain, respectively. The dataset was interpolated to a 1/4°l atitude-longitude rectangular grid via the first-order conservation interpolation method for ease of analysis. The dataset has 55 vertical levels, which is the original vertical level, and has not changed on ESGF nodes. The detail of each vertical level can be found in Appendix B, Table B.
There are 23 variables and 4 parameters provided in this dataset, including daily and monthly SSS, SST and SSH. The monthly three-dimensional ocean temperature, salinity and ocean velocities are also included. The monthly global mean temperature, salinity and AMOC are also uploaded in the dataset. Information on the data grid, such as area, mask, cell thickness and volume for every ocean grid, is also provided. The detail of each variable can be found in Appendix C, Table C.

DESCRIPTION OF CAS-LICOM3_LOW
This appendix describes the experiment design of CAS-LICOM3_Low and compares it with CAS-LICOM3_High.
CAS-LICOM3_Low shares most settings with CAS-LICOM3_High, but there are some differences. First, CAS- LICOM3_Low uses the PHC3.0 temperature and salinity as the initial state. In CAS-LICOM3_Low, the Laplacian scheme was used, which is different from the biharmonic scheme in CAS-LICOM3_High. The biharmonic scheme has the so-called "scale selected" feature. It is more effective to damp the checkerboard noise than the Laplacian scheme. Therefore, the resolved scales are less damped. This is good for the simulation of the mesoscale eddies.
The coupled CICE4 in CAS-LICOM3_Low contains both dynamic and thermodynamic sea-ice processes, while only the thermodynamic process of CICE4 is coupled in CAS-LICOM3_High. The lack of dynamic sea ice can lead to bias in the Arctic, especially for sea-ice volume. A detailed description of the low-resolution dataset can be found in Lin et al. (2020). Table A compares the parameters from CAS-LICOM3_Low and CAS-LICOM3_High. Table A. Comparison between CAS-LICOM3_High and CAS-LICOM3_Low.