Origin of Warm SST Bias over the Atlantic Cold Tongue in the Coupled Climate Model FGOALS-g2

: Most of the coupled models contain a strong warm bias in sea surface temperature (SST) over the Atlantic Cold Tongue (ACT) region (10 ◦ S–3 ◦ N, 20 ◦ W–10 ◦ E) during June–August (JJA) and September–November (SON). In this study, the origins of the ACT SST bias and their relative contributions to the bias are explored by conducting a set of sensitivity experiments, which are based on an ocean-ice model, and by ignoring the nonlinear effects of each origin. The origins for the warm bias over the ACT in the coupled climate model during JJA are estimated as follows: westerly wind bias along the equator (5 ◦ S–5 ◦ N) during March–May (MAM; contributes approximately 32.6% of the warm bias), northerly bias over the southern tropical Atlantic (25 ◦ S–3 ◦ N, 40 ◦ W–20 ◦ E) during MAM and JJA (21.4%), bias in the surface speciﬁc humidity and surface air temperature (11.9%), and downward shortwave radiation bias (6.5%). The origins of the ACT bias during SON are as follows: northerly bias over the southern tropical Atlantic during SON (31.2%), bias in the surface speciﬁc humidity and surface air temperature (27.9%), downward shortwave radiation bias (17.4%), and zonal wind bias (13.4%). Note that these contribution ratios of these origins may be model-dependent. In addition, the local and non-local effects of the zonal wind bias are explored explicitly, while those of all the other biases are examined implicitly. Therefore, a better-performing atmospheric component is crucial when simulating zonal winds during MAM along the equator (5 ◦ S–5 ◦ N) and meridional winds during MAM, JJA, and SON over the southern tropical Atlantic, which will alleviate the warm bias over the ACT region in the coupled climate model.


Introduction
Most of the coupled climate models (CCMs) show a severe warm sea surface temperature (SST) bias for the tropical Atlantic [1][2][3][4][5]. In the southeastern tropical Atlantic (SETA) region (25 • S-10 • S, 10 • W-15 • E), the warm SST bias is mainly concentrated in the Angola-Benguela Frontal Zone [3,[6][7][8][9][10][11][12] (ABFZ) with a magnitude that exceeds 5 • C (IPCC 2013). The origins of the warm bias over the ABFZ have already been attributed to both local sources [4,5,[9][10][11] and remote sources [4,13]. Local sources include underestimated low-level winds [1,11,[14][15][16][17][18] and insufficient stratocumulus cloud cover [19][20][21]. The wind bias results in an SST bias by weakening upwelling in the ABFZ, while the bias in the stratocumulus cloud cover induces an SST bias by enhancing the downward radiative fluxes [8,22]. Cabos et al. [9] and Koseki et al. [11] further linked these local sources to the 2.8 • between 65.58 • S and 65.58 • N and a weighted equal-area grid poleward of 65.58 • with 26 vertical σ-layers up to 2.194 hPa [48]. The horizontal resolution of the ocean model is 1 • with a meridional refinement to 0.5 • in the tropics. There are 30 vertical layers with a spacing of 10 m in the upper 150 m, which gradually increased to 700 m in the deep ocean. The horizontal resolution of the sea ice component is exactly the same as the ocean component [49,50], while that of the land surface component is the same as the atmospheric component. In the tropical Pacific, the model captures well the El Niño-Southern Oscillation (ENSO) [51,52] variability and cloud radiative forcing [53]. However, there are severe warm biases over the ACT region, as most of the CCMs have (detailed in Section 2.3) [4]. To quantify biases in different atmospheric variables to the warm bias over ACT, we adopt the ocean-ice model, i.ethe version 3.4 of the Nucleus for European Modelling of the Ocean (NEMO v3.4) framework [56] in its global configuration (i.e., ORCA2). This ocean-ice model includes the Océan PArallélisé (OPA) for ocean dynamics and thermodynamics and version 3 of the Louvain-la-neuve Ice Model (LIM3) for the sea ice dynamics and thermodynamics [57,58]. The horizontal resolution of this ocean-ice model is 2 • × 2 • with a meridional refinement of 0.5 • in the latitudinal bands 5 • S-5 • N. There are 31 vertical layers from the surface to 5000 m with a grid spacing of 10 m for each layer in the upper 150 m, which gradually increases to an interval of 500 m in the deep ocean. Input fields of the ocean-ice model include 6-h horizontal winds, air temperature, and humidity at 10 m height, daily radiation (shortwave and longwave), and monthly total precipitation (rain and snow) and runoff. The default external forcing of the ocean-ice model NEMO v3.4 is from Common Ocean-ice Reference Experiment version II (CORE-II) [59]. Among all the forcings, the atmospheric forcing is the interannual forcing during the period from 1979 to 2005 and the river runoff forcing is 12-month climatology. We interpolate the forcing data to the grid of NEMO v3.4. Note that the river runoff with a salinity value of 0 psu is discharged into the first layer of the grid cell closest to the river mouth.

Data
The CORE-II dataset has been used as observational data when computing the model bias. The CORE-II data was constructed based on National Centers for Environmental Prediction and the National Center for Atmospheric Research (NCEP-NCAR) reanalysis [60] and many high-quality observation datasets. Therefore, compared to NCEP-NCAR and many other reanalysis datasets, CORE-II has a higher reliability and has been used as the forcing data in many ocean-ice models [59].
Data from the historical experiment of FGOALS-g2 is used for the external forcing of the ocean-ice sensitivity experiments. Variables used include the monthly mean surface winds, downward shortwave radiation, surface specific humidity, surface air temperature, and precipitation. In particular, all variables are from the first member of the historical experiment (i.e., r1i1p1).
The data nudged to the climate model FGOALS-g2 is from NCEP-NCAR for the period from 1979 to 2005 (details about the nudging section experiments can be referred to the supplementary materials). Variables used include the 6-h horizontal winds from 1000 hPa to 10 hPa.
The observational SST data used is from the Hadley Center Sea Ice and Sea Surface Temperature (HadISST) dataset [61] for the period from 1979 to 2005. The thermocline depth is calculated as the 20 • C isotherm depth of the climatological mean temperature data from the National Oceanographic Data Center (NODC) World Ocean Atlas (WOA2009) [62].

Tropical Atlantic Biases of FGOALS-g2
In this section, we access the performance of FGOALS-g2 for simulating the observed seasonal mean tropical Atlantic climate and compare it with the other twenty CCMs from CMIP5 (details in Table S1). All biases shown in this study are estimated according to the following formula: Here, A is the variable under examination. The subscripts of m and o denote the model and observation, respectively, and the overbar is the time-mean of a variable.
In this study, the time period of all experiments is from 1979 to 2005. FGOALS-g2 has a severe warm SST bias in the ACT region, which is very similar to that of most CCMs from CMIP5 and to their ensemble means (Figures 1 and S1). The largest value of SST bias from FGOALS-g2 and their ensemble means are 4.63 • C and 4.43 • C, respectively, over the ACT region during JJA ( Table 2). The RMSE of SST from FGOALS-g2 and their ensemble means are 2.25 • C and 2.38 • C, respectively, over the ACT region during JJA. In addition, the spatial correlation of the SST bias between FGOALS-g2 and the ensemble mean for most of the CMIP5 models is 0.99. Moreover, the largest bias and RMSE of the SST over the ACT region during JJA is comparable to that of all CMIP5 models (Table 2). This indicates that the performance of FGOALS-g2 in the SST simulation is comparable to that of all CMIP5 models (and their ensemble means) in terms of their spatial patterns and magnitudes.

Tropical Atlantic Biases of FGOALS-g2
In this section, we access the performance of FGOALS-g2 for simulating the observed seasonal  mean tropical Atlantic climate and compare it with the other twenty CCMs from CMIP5 (details in  Table S1). All biases shown in this study are estimated according to the following formula: Here, A is the variable under examination. The subscripts of m and o denote the model and observation, respectively, and the overbar is the time-mean of a variable.
In this study, the time period of all experiments is from 1979 to 2005. FGOALS-g2 has a severe warm SST bias in the ACT region, which is very similar to that of most CCMs from CMIP5 and to their ensemble means ( Figure 1 and Figure S1). The largest value of SST bias from FGOALS-g2 and their ensemble means are 4.63 °C and 4.43 °C , respectively, over the ACT region during JJA ( Table  2). The RMSE of SST from FGOALS-g2 and their ensemble means are 2.25 °C and 2.38 °C , respectively, over the ACT region during JJA. In addition, the spatial correlation of the SST bias between FGOALS-g2 and the ensemble mean for most of the CMIP5 models is 0.99. Moreover, the largest bias and RMSE of the SST over the ACT region during JJA is comparable to that of all CMIP5 models (Table 2). This indicates that the performance of FGOALS-g2 in the SST simulation is comparable to that of all CMIP5 models (and their ensemble means) in terms of their spatial patterns and magnitudes. As for the zonal wind, in the observational data, there are strong prevailing easterly winds along the equator during MAM (Figure 2(a1)). In FGOALS-g2, the equatorial easterly wind is underestimated and associated with a westerly bias (Figure 2(a2,b2)). The westerly bias is distinct during MAM and is mainly concentrated in the western and central parts of the Atlantic equator ( Figure 2(a2)). This is similar to the bias in most of the CCMs of CMIP5 ( Figure S2) and their ensemble means (Figure 2(a3,b3)). There is also a serious northerly wind bias over the ACT region in the FGOALS-g2 model (Figure 2(c2)). In the observations, there are prevailing southerly winds As for the zonal wind, in the observational data, there are strong prevailing easterly winds along the equator during MAM (Figure 2(a1)). In FGOALS-g2, the equatorial easterly wind is underestimated and associated with a westerly bias (Figure 2(a2,b2)). The westerly bias is distinct during MAM and is mainly concentrated in the western and central parts of the Atlantic equator (Figure 2(a2)). This is similar to the bias in most of the CCMs of CMIP5 ( Figure S2) and their ensemble means (Figure 2(a3,b3)).
There is also a serious northerly wind bias over the ACT region in the FGOALS-g2 model (Figure 2(c2)). In the observations, there are prevailing southerly winds across the equator year-round ( Figure 2(c1,d1); not shown during SON and DJF). The northerly bias is more severe during MAM than during JJA (Figure 2(c2,d2)). This meridional wind bias in FGOALS-g2 is similar to that in most of the CCMs of CMIP5 ( Figure S3) and similar to their ensemble means (Figure 2(c3,d3)).
Atmosphere 2018, 9, x FOR PEER REVIEW 5 of 21 across the equator year-round ( Figure 2(c1,d1); not shown during SON and DJF). The northerly bias is more severe during MAM than during JJA (Figure 2(c2,d2)). This meridional wind bias in FGOALS-g2 is similar to that in most of the CCMs of CMIP5 ( Figure S3) and similar to their ensemble means (Figure 2(c3,d3)). The horizontal distribution of the biases in the surface air temperature and surface specific humidity during JJA over the ACT region mainly follow the SST bias distribution in that area ( Figures 1 and 3(a2,b2,a3,b3)). Therefore, similar to the SST bias, the bias in the surface air temperature and surface specific humidity of FGOALS-g2 can also be representative of that for all the twenty CMIP5 models over the ACT region (Figures 3(a2,b2,a3,b3) and 5, Figures S4 and S5).
Over the equatorial Atlantic, there is intense precipitation (Figure 3(c1)). FGOALS-g2 has a large wet bias over the southern equatorial Atlantic and a large dry bias over the northern equatorial Atlantic (Figure 3(c2)). These remarkably large biases are similar to those in other CCMs in the CMIP5 ( Figure S6) and their ensemble means (Figure 3(c3)).
FGOALS-g2 also has a large bias in the downward shortwave radiation over the tropical Atlantic (Figures 3(d2) and Figure S7). In the observations, the annual mean shortwave radiations over the ACT region and SETA region are lower than other regions in the tropical Atlantic ( Figure  3(d1)). FGOALS-g2 shows a larger positive radiation bias around the SETA region, a smaller positive radiation bias in the ACT region and a negative bias in the western and central equatorial region (Figure 3(d2)). This negative bias may be caused by the positive precipitation bias (corresponding to more clouds). The horizontal distribution of the biases in the surface air temperature and surface specific humidity during JJA over the ACT region mainly follow the SST bias distribution in that area ( Figures 1 and 3(a2,b2,a3,b3)). Therefore, similar to the SST bias, the bias in the surface air temperature and surface specific humidity of FGOALS-g2 can also be representative of that for all the twenty CMIP5 models over the ACT region (Figures 3(a2,b2,a3,b3) and 5, S4 and S5).
Over the equatorial Atlantic, there is intense precipitation (Figure 3(c1)). FGOALS-g2 has a large wet bias over the southern equatorial Atlantic and a large dry bias over the northern equatorial Atlantic (Figure 3(c2)). These remarkably large biases are similar to those in other CCMs in the CMIP5 ( Figure S6) and their ensemble means (Figure 3(c3)).
FGOALS-g2 also has a large bias in the downward shortwave radiation over the tropical Atlantic (Figures 3(d2) and S7). In the observations, the annual mean shortwave radiations over the ACT region and SETA region are lower than other regions in the tropical Atlantic (Figure 3(d1)). FGOALS-g2 shows a larger positive radiation bias around the SETA region, a smaller positive radiation bias in the ACT region and a negative bias in the western and central equatorial region (Figure 3(d2)). This negative bias may be caused by the positive precipitation bias (corresponding to more clouds). The largest biases and RMSEs of these variables are shown in Table 2. The bias of the forcing data is comparable to that of all twenty CMIP5 models (and their ensemble means) in terms of their spatial patterns and magnitudes. Therefore, it is reasonable to adopt the ocean model NEMO to quantify the contributions of different atmospheric forcing datasets from FGOALS-g2, which has a different ocean component, i.e., LICOM2. In addition, the ocean model NEMO is a widely used ocean general circulation model in the scientific community. Moreover, the ocean model NEMO has been used as the ocean components of a series of coupled climate models in CMIP5. Therefore, studying the SST bias with the NEMO model may benefit the future works using the same ocean model.  The largest biases and RMSEs of these variables are shown in Table 2. The bias of the forcing data is comparable to that of all twenty CMIP5 models (and their ensemble means) in terms of their spatial patterns and magnitudes. Therefore, it is reasonable to adopt the ocean model NEMO to quantify the contributions of different atmospheric forcing datasets from FGOALS-g2, which has a different ocean component, i.e., LICOM2. In addition, the ocean model NEMO is a widely used ocean general circulation model in the scientific community. Moreover, the ocean model NEMO has been used as the ocean components of a series of coupled climate models in CMIP5. Therefore, studying the SST bias with the NEMO model may benefit the future works using the same ocean model.

Transformation of the FGOALS-g2 Data into Forcing Data for the Ocean-Ice Model
The atmospheric variables for FGOALS-g2 have a monthly temporal resolution, which is coarser than the corresponding variables of CORE-II (6 h or daily). Following Koseki et al. [11], the variables from FGOALS-g2 are transformed to have the same temporal resolutions as the variables in CORE-II. Here, we take a variable in CORE-II with a daily resolution as an example. A daily variable in CORE-II could be represented as follows: C daily = C monthly + C daily (2) Here, an overbar denotes the monthly mean of a variable, a prime denotes the corresponding daily deviations, and 'C' denotes the CORE-II forcing. The transformed variable from FGOALS-g2 is computed as follows: Here, 'M' denotes the monthly variable of the FGOALS-g2. The transformed variable serves as input forcing data for the ocean-ice model. The method for computing the 6-h transformed variable is similar.
In some of the experiments, there are atmospheric forcing variables with two data sources (i.e., a sub-region A with forcing from FGOALS-g2 and the rest region B with forcing from CORE-II) in the spatial dimension. To blend these two types of forcing data spatially, five grid cells (within the rest region B) normal to the spatial boundary between A and B are used to build up a buffer region. Taking the daily forcing as an example, the variable of this buffer region is computed as: where α is linearly decreased from 1 to 0 in five cells within region B close to the outer boundary of sub-region A.

Experimental Design
To quantify the relative contributions of different atmospheric forcing biases to the warm SST bias in the ACT region, sensitivity experiments are conducted with an ocean-ice model (NEMO v3.4), which is forced by different mixtures of atmospheric forcings from historical FGOALS-g2 and CORE-II datasets. Table 3 gives a detailed description of these experiments.
There are three control experiments in this study: the first is H_CCM, which is the historical experiment of FGOALS-g2 [6]; the second is CTRL_OI, which is an ocean-ice (OI) experiment forced by the atmospheric forcing from the interannual CORE-II dataset; the third is H_All_OI, which is the same as CTRL_OI, except that the atmospheric forcing is from the H_CCM experiment. All three experiments ran from 1979 to 2005. The CTRL_OI experiment is used to evaluate the SST bias that exists in the ocean-ice model, while the H_ALL_OI experiment is conducted to verify the ability of the NEMO v3.4 model in reproducing the mean state and seasonal cycle of the warm SST bias in the southern tropical Atlantic region given by FGOAL S-g2. All experiments based on the NEMO v3.4 model are restarted from the ocean state after a 500-year spin-up experiment from a state of rest with the initial temperature and salinity fields from Levitus94 [63,64], which is forced by the Normal Year CORE-II forcing dataset.
In addition to the three control experiments, a set of sensitivity experiments have been performed using the NEMO v3.4 model to quantify the roles of the biases in different atmospheric variables. All ocean-ice experiments had the same model configuration, but each had different atmospheric forcing data.
These ocean-ice sensitivity experiments can be divided into two types. In the first type, there are five ocean-ice sensitivity experiments used to quantify the relative contributions from the biases in atmospheric variables (i.e., precipitation, downward shortwave radiation, surface specific humidity, surface air temperature, zonal wind, and meridional wind) to the warm SST bias in the ACT region. These five experiments include the following: (1) H_UVQTR_OI, which is the same as H_All_OI except that the precipitation forcing is from the CORE-II dataset and the aim is to examine the role of precipitation bias; (2) H_UVQT_OI, which is the same as H_UVQTR_OI, except that the downward shortwave radiation is from the CORE-II dataset and is used to quantify the role of shortwave radiation bias; (3) H_UV_OI, which is the same as H_UVQT_OI, except that the specific humidity and air temperature are from the CORE-II dataset and is used to assess the combined effects of biases in surface humidity and surface air temperature that may co-vary with each other [65][66][67][68]; (4) H_U_OI, which is the same as H_UV_OI except that the meridional wind is from the CORE-II dataset; and (5) H_V_OI, which is the same as H_UV_OI, except that the zonal wind is from the CORE-II dataset. The H_U_OI and H_V_OI experiments are carried out to separately quantify the roles of surface meridional wind and zonal wind. Thus, this type of experiment could compare the relative contribution of the bias in each atmospheric variable to the SST bias in the ACT region. In the names of these experiments, the letters after '_' denote different atmospheric forcings. Specifically, 'U' stands for zonal wind, 'V' for meridional forcing, 'Q' for specific humidity, 'T' for air temperature, and 'R' for downward shortwave radiation. In addition, the letter 'H' before '_' indicates that the atmospheric forcing after '_' is from the historical run of FGOALS-g2. 'OI' denotes that the experiment is an ocean-ice run. Based on the above rule, 'H_UV_OI' means that this experiment is forced by zonal (U) and meridional (V) wind forcings from the historical run of FGOALS-g2 and the other atmospheric forcing from the CORE-II dataset.
The second type of sensitivity experiment also contains five experiments, i.e., H_U(MAM_Eq)_OI, H_V(MAM_STA_STA)_OI, H_V(JJA_STA)_OI, H_V(SON_STA)_OI and H_V(DJF_STA)_OI. Based on the results of the first type of experiment, the bias in the surface wind (zonal and meridional) plays a leading role in the formation of the warm SST bias over the ACT region. The second type of experiment is carried out to quantify the impact of the surface wind bias during different seasons on the SST bias. For the zonal wind, we mainly focus on the bias along the equator (5 • S-5 • N) during MAM, since previous studies have pointed out its possible role in inducing the warm bias [8,24,37]. In the second type of experiment, the name 'H_U(MAM_Eq)_OI' means that zonal wind (U) forcing during MAM along the equator (Eq; 5

Computing the Contribution Ratio of a Variable Bias
The contribution ratio of a variable bias to the SST bias is estimated as follows: Here, R is the contribution ratio; RMSE denotes the RMSE of the SST over the ACT region; subscripts a and b denote the experiment after changing the atmospheric forcing and the experiment before changing the atmospheric forcing, respectively; and RMSE hist denotes the RMSE of the SST over the ACT region in the historical FGOALS-g2 run. In addition, the term bracketed by two vertical bars denotes the corresponding absolute value of this term.

Bias in the Control Runs
The seasonal means of the biases from 1979 to 2005 for the H_CCM experiment in terms of surface winds, downward shortwave radiation, surface specific humidity, surface air temperature, and precipitation over the ACT region are displayed in Figure 4. The ACT SST bias is smaller during DJF and MAM and then rapidly increases and reaches its maximum in JJA (Figure 5a). The westerly bias is mainly concentrated along the equator during MAM and north of the equator during JJA. During SON and DJF, the westerly bias remains around the equator, but with much smaller amplitudes ( Figure 4a). The northerly bias is over the ACT region year-round ( Figure 4b) and this bias is strongest during MAM and weakest during JJA (Figure 4b). The biases of the surface specific humidity and surface air temperature over the ACT region are positive during JJA and SON and negative during MAM and DJF (Figure 4c,d). Wet bias (i.e., positive precipitation bias) exists year-round over the ACT region and reaches its peak during MAM and JJA (Figure 4e). In particular, the maxima of this wet bias during MAM and JJA cover the western part of the ACT and the oceans to the west of the ACT. Results from the atmospheric component of the FGOALS-g2 model forced by SST perturbations over the ACT region indicate that the warm bias over the ACT region plays an indispensable role in generating this wet bias (not shown). Similarly, the bias in shortwave radiation over the ACT region is positive during JJA and SON and negative during MAM and DJF. The roles of these atmospheric biases on the ACT SST bias will be discussed in Sections 4.2 and 4.3. Next, we assess the SST bias of the control experiment CTRL_OI in the ocean-ice model NEMO v3.4. In the CTRL_OI experiment, we applied restoring conditions for both the SST and sea surface salt (SSS). As a consequence, the seasonal mean and annual mean of the SST biases over the ACT region and SETA region are much smaller (Figures 5c and 6c) relative to the previous studies [11,18]. The RMSE of the SST in the ACT region during JJA is 0.82 ℃ with the largest value exceeding 1 ℃ ( Figure 5(c2)) and the RMSE of the SST over the southern tropical Atlantic is 0.56 ℃ (Figure 6c). A similar experiment, except without any restoring conditions, is conducted to evaluate the impact of the restoring conditions (Figures 5d and 6d). In this experiment (No_restore_OI), the RMSE of the SST during JJA is 1.41 ℃, with the largest value exceeding 2 ℃ in the ACT region ( Figure 5(d2)) and the RMSE is 0.94 ℃ with the largest value exceeding 1.5 ℃ in the southern tropical Atlantic region (Figures 5(d3) and 6d), which is consistent with previous studies [11,18]. It is worth noting that the impacts of the restoring conditions on the SST RMSE during the four seasons are different.  Next, we assess the SST bias of the control experiment CTRL_OI in the ocean-ice model NEMO v3.4. In the CTRL_OI experiment, we applied restoring conditions for both the SST and sea surface salt (SSS). As a consequence, the seasonal mean and annual mean of the SST biases over the ACT region and SETA region are much smaller (Figures 5c and 6c) relative to the previous studies [11,18]. The RMSE of the SST in the ACT region during JJA is 0.82 • C with the largest value exceeding 1 • C ( Figure 5(c2)) and the RMSE of the SST over the southern tropical Atlantic is 0.56 • C (Figure 6c). A similar experiment, except without any restoring conditions, is conducted to evaluate the impact of the restoring conditions (Figures 5d and 6d). In this experiment (No_restore_OI), the RMSE of the SST during JJA is 1.41 • C, with the largest value exceeding 2 • C in the ACT region ( Figure 5(d2)) and the RMSE is 0.94 • C with the largest value exceeding 1.5 • C in the southern tropical Atlantic region ( Figures 5(d3) and 6d), which is consistent with previous studies [11,18]. It is worth noting that the impacts of the restoring conditions on the SST RMSE during the four seasons are different. The SST RMSEs during MAM, JJA, SON, and DJF in the experiment with restoring conditions (CTRL_OI experiment) are decreased by 0.51 • C, 0.59 • C, 0.41 • C, and 0.32 • C, respectively (Figure 5c,d), relative to those in the experiment without restoring conditions (No_restore_OI experiment). These differences in SST RMSEs between the two experiments correspond to about 48.6%, 15.9%, 14.9%, and 35.2% of the SST RMSEs over the ACT region in the experiment H_CCM during MAM, JJA, SON, and DJF. Therefore, the impacts of the restoring conditions on the SST bias over the ACT region are relatively smaller during JJA and SON compared to the other two seasons, indicating that the use of the restoring conditions in examining the SST bias over the ACT region during JJA and SON will not strongly affect the results. Moreover, the restoring conditions could alleviate the SST bias induced by the oceanic sources, which, in turn, could help us emphasize the impacts of the atmospheric sources. Notably, we applied restoring conditions on the SST and SSS to all the ocean-ice experiments with a restoring time scale of 50 days except the No_restore_OI experiment. Finally, we evaluate the ability of the NEMO v3.4 model to reproduce the mean state and seasonal cycle of the warm SST bias of FGOALS-g2 over the southern tropical Atlantic region through an ocean-ice experiment forced by the atmosphere forcing from FGOALS-g2 (H_All_OI). H_All_OI has a very similar SST bias as H_CCM in the ACT region (Figure 5a,b) during all four seasons. The RMSE of SST in H_All_OI is 0.09 ℃ larger than that in H_CCM over the ACT region (Figure 5a,b). This may be related to the fact that the underestimated shortwave radiation and underestimated surface specific humidity and surface air temperature during MAM and DJF ( Figure 4) could lead to a cold bias, which, in turn, increases these biases over the ACT region in H_CCM. By contrast, the above process is lacking in H_All_OI. However, in the SETA region, the SST bias is smaller (Figures 5a,b and 6a,b) in H_All_OI compared to that in H_CCM. This may be because that there is a stratocumulus-SST feedback in H_CCM. Specifically, the overestimated shortwave radiation [69] leads to a warm bias, which in turn results in less stratocumulus (corresponding to more shortwave radiation [69]) over the SETA region in the H_CCM experiment. By contrast, the above feedback process is lacking in H_ALL_OI. Finally, we evaluate the ability of the NEMO v3.4 model to reproduce the mean state and seasonal cycle of the warm SST bias of FGOALS-g2 over the southern tropical Atlantic region through an ocean-ice experiment forced by the atmosphere forcing from FGOALS-g2 (H_All_OI). H_All_OI has a very similar SST bias as H_CCM in the ACT region (Figure 5a,b) during all four seasons. The RMSE of SST in H_All_OI is 0.09 • C larger than that in H_CCM over the ACT region (Figure 5a,b). This may be related to the fact that the underestimated shortwave radiation and underestimated surface specific humidity and surface air temperature during MAM and DJF ( Figure 4) could lead to a cold bias, which, in turn, increases these biases over the ACT region in H_CCM. By contrast, the above process is lacking in H_All_OI. However, in the SETA region, the SST bias is smaller (Figues 5a,b and 6a,b) in H_All_OI compared to that in H_CCM. This may be because that there is a stratocumulus-SST feedback in H_CCM. Specifically, the overestimated shortwave radiation [69] leads to a warm bias, which in turn results in less stratocumulus (corresponding to more shortwave radiation [69]) over the SETA region in the H_CCM experiment. By contrast, the above feedback process is lacking in H_ALL_OI.

Role of Different Atmospheric Biases
In this section, we perform sensitivity experiments to quantify the roles of biases in different atmospheric variables in inducing the SST bias in the ACT region. We do this by removing the atmospheric variable bias from the atmospheric forcing for the ocean-ice experiment one by one. Figure 7 shows the seasonal mean SST bias from 1979 to 2005 in the second type of sensitivity experiments over the ACT region. The H_UVQTR_OI experiment is identical to the H_All_OI experiment except that the precipitation forcing is from CORE-II. The warm SST bias in H_UVQTR_OI (Figure 7a) is nearly the same as that in H_All_OI ( Figure 5) over the ACT region during all four seasons, with a mean difference of less than 0.1 ℃. This indicates that the precipitation bias has a relatively weak impact on the SST bias over the ACT region. This is consistent with Richter et al. [24] who has suggested that the precipitation bias generally induces an SST bias by setting up a westerly wind bias in the CCM. Thus, the nearly unchanged SST bias over the ACT region is reasonable in H_UVQTR_OI, where the zonal wind is from observational data.
As for the role of the bias in the shortwave radiation, we performed an ocean-ice experiment H_UVQT_OI, which is similar to H_UVQTR_OI except that the shortwave radiation is from CORE-II. H_UVQT_OI shows that the absence of a shortwave radiation bias enhances the RMSE of the SST over the ACT region by approximately 0.11 ℃ during MAM, approximately 0.03 ℃ during DJF and reduces it by 0.24 ℃ during JJA and by 0.48 ℃ during SON (Figure 7a,b). This indicates that the shortwave radiation bias contributes approximately 6.5% to the warm SST bias over the ACT region during JJA, approximately 17.4% during SON and it is not a source of the warm SST bias over the ACT region during MAM and DJF. The reduced SST bias during JJA and SON can be explained by the reduction in the shortwave radiation bias over the ACT region. Moreover, the more the shortwave radiation bias is reduced, the more the warm SST bias is reduced (Figure 4(f2,f3)). Similarly, the enhanced SST bias during MAM and DJF can be explained by the enhancement of the shortwave radiation bias. This also shows that the impact of the shortwave radiation bias on the SST bias over the ACT region is local.

Role of Different Atmospheric Biases
In this section, we perform sensitivity experiments to quantify the roles of biases in different atmospheric variables in inducing the SST bias in the ACT region. We do this by removing the atmospheric variable bias from the atmospheric forcing for the ocean-ice experiment one by one. Figure 7 shows the seasonal mean SST bias from 1979 to 2005 in the second type of sensitivity experiments over the ACT region. The H_UVQTR_OI experiment is identical to the H_All_OI experiment except that the precipitation forcing is from CORE-II. The warm SST bias in H_UVQTR_OI (Figure 7a) is nearly the same as that in H_All_OI ( Figure 5) over the ACT region during all four seasons, with a mean difference of less than 0.1 • C. This indicates that the precipitation bias has a relatively weak impact on the SST bias over the ACT region. This is consistent with Richter et al. [24] who has suggested that the precipitation bias generally induces an SST bias by setting up a westerly wind bias in the CCM. Thus, the nearly unchanged SST bias over the ACT region is reasonable in H_UVQTR_OI, where the zonal wind is from observational data.
As for the role of the bias in the shortwave radiation, we performed an ocean-ice experiment H_UVQT_OI, which is similar to H_UVQTR_OI except that the shortwave radiation is from CORE-II. H_UVQT_OI shows that the absence of a shortwave radiation bias enhances the RMSE of the SST over the ACT region by approximately 0.11 • C during MAM, approximately 0.03 • C during DJF and reduces it by 0.24 • C during JJA and by 0.48 • C during SON (Figure 7a,b). This indicates that the shortwave radiation bias contributes approximately 6.5% to the warm SST bias over the ACT region during JJA, approximately 17.4% during SON and it is not a source of the warm SST bias over the ACT region during MAM and DJF. The reduced SST bias during JJA and SON can be explained by the reduction in the shortwave radiation bias over the ACT region. Moreover, the more the shortwave radiation bias is reduced, the more the warm SST bias is reduced (Figure 4(f2,f3)). Similarly, the enhanced SST bias during MAM and DJF can be explained by the enhancement of the shortwave radiation bias. This also shows that the impact of the shortwave radiation bias on the SST bias over the ACT region is local. As for the role of bias in surface specific humidity and air temperature, we perform another ocean-ice experiment (H_UV_OI), which is similar to H_UVQT_OI except that the specific humidity and air temperature forcing is from CORE-II. We remove the biases in both the specific humidity and air temperature from the forcing, as these parameters co-vary with each other and have similar effects on the SST bias [67,68]. In H_UV_OI the RMSE of SST over the ACT is enhanced by approximately 0.28 ℃ during MAM, 0.09 ℃ during DJF, and is reduced by approximately 0.44 ℃ during JJA and 0.77 ℃ during SON (Figure 7b,c). This indicates that the bias in the specific humidity and air temperature mainly contributes to the warm SST bias over the ACT region during JJA (approximately 11.9%) and SON (27.9%), while these parameters have cooling effects on the warm SST bias during MAM and DJF. The reduced SST bias during JJA and SON can be explained by the reduction in the positive bias in the specific humidity and air temperature, which leads to more latent and sensible heat flux out from the sea surface. Moreover, the enhanced SST bias during MAM and DJF can be attributed to the reduction in the negative biases in the specific humidity and air temperature. Furthermore, the effect of the biases in the specific humidity and air temperature on the SST bias are local over the ACT region.
For the bias effects in the surface meridional wind, we perform an ocean-ice experiment (H_U_OI), which is similar to H_UV_OI, except that the meridional wind is from CORE-II. H_U_OI shows that the absence of the meridional wind bias can reduce the warm SST bias over the ACT region in all four seasons. In H_U_OI, the SST RMSE over the ACT region is reduced by approximately 1.05 ℃ during MAM, 0.62 ℃ during JJA, 0.86 ℃ during SON, and 0.97 ℃ during DJF (Figure 7c,d). This indicates that the original warm SST bias of FGOALS-g2 over the ACT region during MAM (1.05 ℃) and DJF (0.91 ℃) is mainly induced by the meridional wind bias, and 16.7% of that during JJA (3.71 ℃) and 31.2% of that during SON (2.76 ℃) is from the meridional wind bias.
For the role of the bias in the surface zonal wind, we perform an ocean-ice experiment (H_V_OI), which is similar to H_UV_OI except that the zonal forcing is from CORE-II. H_V_OI shows that the absence of the zonal wind bias can reduce the warm SST bias over the ACT region in As for the role of bias in surface specific humidity and air temperature, we perform another ocean-ice experiment (H_UV_OI), which is similar to H_UVQT_OI except that the specific humidity and air temperature forcing is from CORE-II. We remove the biases in both the specific humidity and air temperature from the forcing, as these parameters co-vary with each other and have similar effects on the SST bias [67,68]. In H_UV_OI the RMSE of SST over the ACT is enhanced by approximately 0.28 • C during MAM, 0.09 • C during DJF, and is reduced by approximately 0.44 • C during JJA and 0.77 • C during SON (Figure 7b,c). This indicates that the bias in the specific humidity and air temperature mainly contributes to the warm SST bias over the ACT region during JJA (approximately 11.9%) and SON (27.9%), while these parameters have cooling effects on the warm SST bias during MAM and DJF. The reduced SST bias during JJA and SON can be explained by the reduction in the positive bias in the specific humidity and air temperature, which leads to more latent and sensible heat flux out from the sea surface. Moreover, the enhanced SST bias during MAM and DJF can be attributed to the reduction in the negative biases in the specific humidity and air temperature. Furthermore, the effect of the biases in the specific humidity and air temperature on the SST bias are local over the ACT region.
For the bias effects in the surface meridional wind, we perform an ocean-ice experiment (H_U_OI), which is similar to H_UV_OI, except that the meridional wind is from CORE-II. H_U_OI shows that the absence of the meridional wind bias can reduce the warm SST bias over the ACT region in all four seasons. In H_U_OI, the SST RMSE over the ACT region is reduced by approximately 1.05 • C during MAM, 0.62 • C during JJA, 0.86 • C during SON, and 0.97 • C during DJF (Figure 7c,d). This indicates that the original warm SST bias of FGOALS-g2 over the ACT region during MAM (1.05 • C) and DJF (0.91 • C) is mainly induced by the meridional wind bias, and 16.7% of that during JJA (3.71 • C) and 31.2% of that during SON (2.76 • C) is from the meridional wind bias.
For the role of the bias in the surface zonal wind, we perform an ocean-ice experiment (H_V_OI), which is similar to H_UV_OI except that the zonal forcing is from CORE-II. H_V_OI shows that the absence of the zonal wind bias can reduce the warm SST bias over the ACT region in all four seasons.
In H_V_OI the absence of the zonal wind bias can reduce the SST RMSE by 0.54 • C during MAM, 1.36 • C during JJA, 0.37 • C during SON, and 0.23 • C during DJF over the ACT region (Figure 7c,e).
This indicates that the zonal wind bias is one of the sources of the SST bias over the ACT region, and it contributes approximately 51.4% to the SST bias during MAM, 36.7% during JJA, 13.4% during SON and 25.3% during DJF over the ACT region. This shows that the contribution ratio of the zonal wind bias is the largest among the six atmospheric variables during JJA, which is consistent with previous studies.
The sum of the contribution ratios of zonal wind bias and meridional wind bias to the SST bias over the ACT region during MAM and DJF is larger than 100%, which is because the precipitation bias, specific humidity and air temperature biases, and shortwave radiation bias have cooling effects on the SST bias over the ACT region during MAM and DJF.

Role of the Surface Wind Bias
The previous subsection shows that both the zonal wind bias and meridional wind bias contribute largely to the warm SST bias over the ACT region year-round. Previous studies have indicated that the leading role of the zonal wind bias during MAM along the equator induces the warm SST bias over the ACT region. In this section, we perform the second type of sensitivity experiment to verify the important role of the zonal wind bias during MAM and further investigate the role of the meridional wind bias during different seasons. Figure 8 shows the seasonal mean SST bias over the southern tropical Atlantic for all four seasons. To verify the role of the zonal wind bias during MAM along the equator, we perform the H_U(MAM_Eq)_OI experiment, which is similar to the CTRL_OI run, except where the zonal wind forcing during MAM along the equator is from the H_CCM experiment. H_U(MAM_Eq)_OI shows a significant warm SST bias over the ACT region during JJA. This indicates that the westerly bias during MAM along the equator is a main source of the warm bias during JJA over the ACT region and the contribution ratio is approximately 32.6% (Figure 8a). This zonal wind bias enhances the SST bias over the ACT region by 0.17 • C during MAM, 0.10 • C during SON, and 0.01 • C during DJF, which indicates that the zonal wind bias during MAM is not the main source of the warm bias during these three seasons, which is consistent with previous studies. The westerly bias during MAM along the equator enhances the warm bias during JJA and can be explained by the deepening of the thermocline in the ACT region during MAM and JJA (Figure 9a), which has been mentioned in previous studies by Richter et al. [24,37]. The deepening of the thermocline during MAM is accompanied by a warm bias in the subsurface temperature (Figure 10a). The warmer water will then upwell to the sea surface by the prevailing southerly winds, which inhibits the upwelling of the cooling water and results in a warm bias in the sea surface.
As for the meridional wind bias role over the southern tropical Atlantic in inducing the warm SST bias over the ACT region during different seasons, we perform four experiments that are similar to CTRL_OI, except that the meridional wind forcing during different seasons over the southern tropical Atlantic is from the H_CCM experiment. H_V(MAM_STA)_OI shows that the RMSE of SST over the ACT region is enhanced by approximately 0.55 • C during MAM, 0.34 • C during JJA, and 0.02 • C during SON and DJF compared with that of CTRL_OI. This indicates that the northerly bias during MAM over the southern tropical Atlantic mainly contributes to the warm SST bias over the ACT region during MAM (approximately 52.4%) and JJA (approximately 9.2%), and has a very small effect on the warm bias over the region during SON and DJF. This northerly bias during MAM, which enhanced the warm bias during JJA over the ACT region, can be explained by the deepening of the thermocline as the zonal wind bias does (Figure 9b). The deepening of the thermocline is accompanied by a warm bias in the subsurface temperature (Figure 10b). The warmer water is then upwelled to the sea surface and prevents development of the cold tongue during JJA over the ACT region. Similarly, H_V(JJA_STA)_OI shows an enhanced SST bias only during JJA compared with that of CTRL_OI in the ACT region (Figure 8c), which means that the northerly wind bias during JJA can only contribute to the warm bias during JJA (approximately 12.2%) in the region. This meridional wind bias during JJA also affects the SST bias in that area through deepening of the thermocline and warming of the subsurface temperature during JJA over the ACT region (Figures 9c and 10c). Moreover, the northerly wind bias during SON (approximately 29.0%) and DJF (approximately 81.3%) only enhances the SST bias during the same time period (Figure 8d,e). Therefore, only the northerly wind bias during MAM has a delayed effect on the SST bias over the ACT region during the next season (JJA), while that during each of the other seasons could affect the SST bias over the ACT region during the same season.
Furthermore, results from the nudging experiments (Table S2) in the FGOALS-g2 model support the above findings. In the nudging experiments, the zonal wind bias during MAM along the equator plays a leading role in inducing the warm SST bias over the ACT region during JJA with a contribution ratio of approximately 45.6%. Meanwhile, the northerly wind bias during MAM and JJA over the southern tropical Atlantic contributes approximately 32.3% to the warm SST bias over the ACT region during JJA. In the nudging experiments, the combined effects of the zonal wind bias and meridional wind bias contribute approximately 71.4% of the warm bias during JJA (Table S3 and Figure S8).
Atmosphere 2018, 9, x FOR PEER REVIEW 15 of 21 12.2%) in the region. This meridional wind bias during JJA also affects the SST bias in that area through deepening of the thermocline and warming of the subsurface temperature during JJA over the ACT region (Figures 9c and 10c). Moreover, the northerly wind bias during SON (approximately 29.0%) and DJF (approximately 81.3%) only enhances the SST bias during the same time period (Figure 8d,e). Therefore, only the northerly wind bias during MAM has a delayed effect on the SST bias over the ACT region during the next season (JJA), while that during each of the other seasons could affect the SST bias over the ACT region during the same season. Furthermore, results from the nudging experiments (Table S2) in the FGOALS-g2 model support the above findings. In the nudging experiments, the zonal wind bias during MAM along the equator plays a leading role in inducing the warm SST bias over the ACT region during JJA with a contribution ratio of approximately 45.6%. Meanwhile, the northerly wind bias during MAM and JJA over the southern tropical Atlantic contributes approximately 32.3% to the warm SST bias over the ACT region during JJA. In the nudging experiments, the combined effects of the zonal wind bias and meridional wind bias contribute approximately 71.4% of the warm bias during JJA (Table S3 and Figure S8).   12.2%) in the region. This meridional wind bias during JJA also affects the SST bias in that area through deepening of the thermocline and warming of the subsurface temperature during JJA over the ACT region (Figures 9c and 10c). Moreover, the northerly wind bias during SON (approximately 29.0%) and DJF (approximately 81.3%) only enhances the SST bias during the same time period (Figure 8d,e). Therefore, only the northerly wind bias during MAM has a delayed effect on the SST bias over the ACT region during the next season (JJA), while that during each of the other seasons could affect the SST bias over the ACT region during the same season. Furthermore, results from the nudging experiments (Table S2) in the FGOALS-g2 model support the above findings. In the nudging experiments, the zonal wind bias during MAM along the equator plays a leading role in inducing the warm SST bias over the ACT region during JJA with a contribution ratio of approximately 45.6%. Meanwhile, the northerly wind bias during MAM and JJA over the southern tropical Atlantic contributes approximately 32.3% to the warm SST bias over the ACT region during JJA. In the nudging experiments, the combined effects of the zonal wind bias and meridional wind bias contribute approximately 71.4% of the warm bias during JJA (Table S3 and Figure S8).

Summary and Conclusions
This study quantifies the sources of the warm SST bias over the ACT region in most of the CCMs of CMIP5. This is completed through two types of sensitivity experiments based on an ocean-ice model and the historical simulations of the CCM FGOALS-g2. The FGOALS-g2 model has an SST bias and atmospheric variable bias over the ACT region similar to that of most CMIP5 CCMs. The SST bias during JJA of FGOALS-g2 over the ACT region is the largest, followed by that during SON and the biases during MAM and DJF are much smaller. This is similar to that of most CCMs in CMIP5. The westerly surface wind bias of FGOALS-g2 is mainly concentrated in the central and western equatorial Atlantic during MAM and is relatively weak during the other seasons in that area. The year-round northerly surface wind bias of FGOALS-g2 is mainly concentrated over the southern tropical Atlantic. The bias distributions in the shortwave radiation, surface specific humidity, and surface air temperature are similar to each other over the southern tropical Atlantic, with a positive value during JJA and SON and a negative value during MAM and DJF. In addition, there is a year-round wet bias over the ACT region. All these biases in FGOALS-g2 are similar to that of most CMIP5 CCMs and their ensemble means.
The two types of sensitivity experiments are forced by different mixtures of atmospheric forcings from historical FGOALS-g2 and CORE-II experimental datasets. The analysis based on these experiments confirms that the SST biases over the ACT region are different during all four seasons and the sources of them are different. The largest warm SST bias appears during JJA, which has four main origins, i.e., westerly wind bias during MAM in the western and central equatorial Atlantic (5° S-5° N, 40° W-5° W), northerly wind bias over the southern tropical Atlantic during MAM and JJA, and positive biases in shortwave radiation, surface specific humidity, and air temperature above the southern tropical Atlantic during JJA. Among these origins, the westerly bias during MAM in the western and central equatorial Atlantic is the most important. It can explain approximately 32.6% of the ACT warm bias during JJA. This is consistent with the previous study of Richter et al. [16] who had pointed out the leading role of this zonal wind bias. The northerly wind bias during MAM and JJA over the southern tropical Atlantic contributes approximately 21.4% to the warm bias through both local and non-local effects. The local effect is exerted by the northerly wind bias over the ACT during MAM and JJA through deepening the thermocline and warming the subsurface water during MAM and JJA, which can prevent the development of the cold tongue during JJA over the ACT region and result in a warm SST bias in that area. The non-local effect is exerted by the northerly wind bias off the ACT region through weakening the Angola Current, which can prevent the heat loss of the ACT region and result in a warm SST bias in that area. Moreover, the positive biases in the surface specific humidity and air temperature during JJA over the ACT region contribute approximately 11.9% to the warm SST bias during JJA over the ACT region through local effects. The overestimated surface specific humidity and air temperature induces SST bias through preventing the latent and sensible heat fluxes from dissipating from the sea surface, which leads to more heat remaining in the sea surface and results in a warm SST bias. It should be noted that the biases in surface specific humidity and air temperature essentially stem from the warm SST bias under examination. Therefore, the impacts of the biases in surface specific humidity and air temperature on the warm bias can be viewed as a feedback. Finally, the bias in

Summary and Conclusions
This study quantifies the sources of the warm SST bias over the ACT region in most of the CCMs of CMIP5. This is completed through two types of sensitivity experiments based on an ocean-ice model and the historical simulations of the CCM FGOALS-g2. The FGOALS-g2 model has an SST bias and atmospheric variable bias over the ACT region similar to that of most CMIP5 CCMs. The SST bias during JJA of FGOALS-g2 over the ACT region is the largest, followed by that during SON and the biases during MAM and DJF are much smaller. This is similar to that of most CCMs in CMIP5. The westerly surface wind bias of FGOALS-g2 is mainly concentrated in the central and western equatorial Atlantic during MAM and is relatively weak during the other seasons in that area. The year-round northerly surface wind bias of FGOALS-g2 is mainly concentrated over the southern tropical Atlantic. The bias distributions in the shortwave radiation, surface specific humidity, and surface air temperature are similar to each other over the southern tropical Atlantic, with a positive value during JJA and SON and a negative value during MAM and DJF. In addition, there is a year-round wet bias over the ACT region. All these biases in FGOALS-g2 are similar to that of most CMIP5 CCMs and their ensemble means.
The two types of sensitivity experiments are forced by different mixtures of atmospheric forcings from historical FGOALS-g2 and CORE-II experimental datasets. The analysis based on these experiments confirms that the SST biases over the ACT region are different during all four seasons and the sources of them are different. The largest warm SST bias appears during JJA, which has four main origins, i.e., westerly wind bias during MAM in the western and central equatorial Atlantic (5 • S-5 • N, 40 • W-5 • W), northerly wind bias over the southern tropical Atlantic during MAM and JJA, and positive biases in shortwave radiation, surface specific humidity, and air temperature above the southern tropical Atlantic during JJA. Among these origins, the westerly bias during MAM in the western and central equatorial Atlantic is the most important. It can explain approximately 32.6% of the ACT warm bias during JJA. This is consistent with the previous study of Richter et al. [16] who had pointed out the leading role of this zonal wind bias. The northerly wind bias during MAM and JJA over the southern tropical Atlantic contributes approximately 21.4% to the warm bias through both local and non-local effects. The local effect is exerted by the northerly wind bias over the ACT during MAM and JJA through deepening the thermocline and warming the subsurface water during MAM and JJA, which can prevent the development of the cold tongue during JJA over the ACT region and result in a warm SST bias in that area. The non-local effect is exerted by the northerly wind bias off the ACT region through weakening the Angola Current, which can prevent the heat loss of the ACT region and result in a warm SST bias in that area. Moreover, the positive biases in the surface specific humidity and air temperature during JJA over the ACT region contribute approximately 11.9% to the warm SST bias during JJA over the ACT region through local effects. The overestimated surface specific humidity and air temperature induces SST bias through preventing the latent and sensible heat fluxes from dissipating from the sea surface, which leads to more heat remaining in the sea surface and results in a warm SST bias. It should be noted that the biases in surface specific humidity and air temperature essentially stem from the warm SST bias under examination. Therefore, the impacts of the biases in surface specific humidity and air temperature on the warm bias can be viewed as a feedback. Finally, the bias in shortwave radiation during JJA over the ACT region contributes approximately 6.5% of the warm SST bias during the same season through local effects. Specifically, the positive bias in shortwave radiation over the southeastern part of the ACT region could increase the warm bias, while the negative bias in shortwave radiation over the northwestern part of the ACT region could alleviate the warm bias. However, the combined effects of the biases in shortwave radiation over the two regions increase the warm bias. Following the warm SST bias over the ACT region during JJA, the mean warm SST bias during SON is approximately 2.76 • C over the ACT region. This bias is mainly from four origins, i.e., northerly wind bias, positive bias in shortwave radiation, surface specific humidity and surface air temperature over the southern tropical Atlantic during SON, and zonal wind bias over the central part of the equatorial Atlantic during SON. Among these origins, the northerly bias during SON over the ACT region plays a leading role with a contribution ratio of approximately 31.2%. This northerly bias induces a warm SST bias in the region during SON through both local effects (by weakening the upwelling) and non-local effects (by weakening the Angola Current). The positive biases in the surface humidity and surface air temperature during SON over the ACT region contribute approximately 27.9% of the warm bias during SON through local effects, while the positive bias in the shortwave radiation during the same season contributes approximately 17.4% through local effects. The zonal wind bias contributes approximately 13.4% of the warm bias during SON. Furthermore, the warm bias over the ACT region during MAM and DJF is approximately 1 • C. This warm bias is mainly from the northerly bias during the same time period, of which the contribution ratio is larger than 95.0%. It is worth noting that all the above contribution ratios are estimated based on the ocean-ice model NEMO forced by atmospheric forcing from the CCM FGOALS-g2, which means that these ratios may be model-dependent.
The results of our study support the findings by Richter et al. [24,37] that the westerly bias along the equator during MAM (with a contribution ratio of approximately 32.6% in our study) plays a leading role in inducing the warm bias over the ACT region during JJA. Compared to previous studies, there are new findings in our study. First, results in this study indicate the northerly wind biases over the southern tropical Atlantic during MAM and JJA (with a contribution ratio of approximately 21.4%) also play important roles in inducing the warm bias during JJA over the ACT region and the northerly wind bias during SON play a leading role (approximately 31.2%) in inducing the warm bias over the ACT region during SON. The northerly wind bias over the southern tropical Atlantic has both local and non-local effects on the warm SST bias over the ACT region. Second, the biases in surface specific humidity and surface air temperature mainly contribute to the warm SST bias over the ACT region during JJA (approximately 11.9%) and SON (27.9%) and have cooling effects on the warm SST bias during MAM and DJF. Moreover, there are feedback processes existing between biases in surface specific humidity and air temperature and the warm bias over the ACT region. Third, the bias in shortwave radiation mainly contributes to the warm SST bias over the ACT region during JJA (approximately 6.5%) and SON (approximately 17.4%) and is not a source of the warm SST bias over the ACT region during MAM and DJF. Note that the biases in surface specific humidity, surface air temperature, and shortwave radiation mainly have local effects on the warm SST bias over the ACT region. Fourth, the precipitation bias has a relatively weak impact on the SST bias over the ACT region. Therefore, an atmospheric component with better performance in simulating the zonal winds during MAM over the equatorial Atlantic and meridional winds during MAM, JJA, and SON over the southern tropical Atlantic may be helpful in alleviating the warm bias over the ACT region in the CCM.
There are two limitations in our study. One is that we use an ocean-ice model with a relatively coarse resolution. Small et al. [12], Seo et al. [70], and Zheng et al. [71] attributed the warm SST bias to the low horizontal resolution of the oceanic models. Further, Harlaß el al. [10] pointed out that higher resolution of the atmospheric component in the CCM is also crucial for reducing the warm SST bias over the ACT region. In the future, a high-resolution ocean model should be used to quantify the role of atmospheric variable biases in inducing SST biases over the ACT region. Another is that we have isolated the effect of each atmospheric variable bias in this study. However, there exists nonlinear interaction between two or more variables. In the future, the nonlinear interaction of the biases in several variables will be studied and the possible impacts of the nonlinear interaction will be examined.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4433/9/7/275/s1, Figure S1: The annual-mean sea surface temperature (SST; units: • C) bias during 1979-2005 of different coupled climate models (CCMs) in CMIP5 (a and c-v) and their ensemble mean (b) over the southern tropical Atlantic. The blue box denotes the Atlantic cold tongue (ACT; 10 • S-3 • N, 20 • W-10 • E) region, Figure S2: The same as Figure S1 but for the annual-mean surface zonal wind (units: m s −1 ) bias. The blue box denotes the central and western equatorial Atlantic (5 • S-5 • N, 40 • W-5 • W), Figure S3: The same as Figure S1 but for the annual-mean surface meridional wind (units: m s −1 ) bias. The blue box includes the ACT and the region south of it, Figure S4: The same as Figure S1 but for the surface air temperature (units: • C) bias. The blue box includes the ACT region and the southeastern tropical Atlantic (SETA; 25 • S-10 • S, 10 • W-15 • E), Figure S5: The same as Figure S1 but for the surface specific humidity (units: g kg −1 ) bias, Figure S6: The same as Figure S1 but for the precipitation (units: mm day −1 ) bias. Note that, the color bar is reversed compared with that of the figures above, Figure S7: The same as Figure S1 but for the downward shortwave radiation (units: W m −2 ) bias. The blue box includes the ACT and a part of the SETA region, Figure Table S1: CMIP5 models used in this study, Table S2: Nudging experiments based on the coupled model FGOALS-g2,