PUBLICATIONS Journal of Advances in Modeling Earth Systems RESEARCH ARTICLE 10.1002/2014MS000415 Key Points: Climate models has cold SST biases in Northern basins The cold biases in the Atlantic bias have remote impacts on the Pacific The atmospheric bridge is the main pathway for connecting biases in both basins
Correspondence to: C. Zhao,
[email protected] Citation: Zhang, L., and C. Zhao (2015), Processes and mechanisms for the model SST biases in the North Atlantic and North Pacific: A link with the Atlantic meridional overturning circulation, J. Adv. Model. Earth Syst., 7, 739–758, doi:10.1002/2014MS000415. Received 9 DEC 2014 Accepted 6 APR 2015 Accepted article online 6 APR 2015 Published online 25 MAY 2015
Processes and mechanisms for the model SST biases in the North Atlantic and North Pacific: A link with the Atlantic meridional overturning circulation Liping Zhang1,2 and Chuanhu Zhao3 1
Atmospheric and Oceanic Sciences Program, Princeton University, Princeton, New Jersey, USA, 2NOAA Geophysical Fluids Dynamics Laboratory, Princeton, New Jersey, USA, 3Physical Oceanography Laboratory, Key Laboratory of OceanAtmosphere Interaction and Climate in Universities of Shandong, Ocean University of China, Qingdao, China
Abstract Almost all of CMIP5 climate models show cold SST biases in the extratropical North Atlantic (ENA) and tropical North Atlantic (TNA) as well as in the North Pacific which are commonly linked with the weak simulated Atlantic meridional overturning circulation (AMOC). A weak AMOC and its associated reduced northward oceanic heat transport are associated with a cooling of the ENA Ocean, whereas the TNA cooling is attributable to both weak AMOC and surface heat flux. The cold biases in the ENA and TNA have remote impacts on the SST bias in the North Pacific. Here we use coupled ocean-atmosphere model experiments to show the mechanisms and pathways by which the ENA and TNA affect the North Pacific. The model simulations demonstrate that the cooling SST bias in the North Pacific is largely due to the remote effect of the cooling SST bias in the ENA, while the remote impact of the TNA cooling SST bias is of secondary importance. The ENA cooling bias triggers the circumglobal teleconnection via the Northern Hemisphere annular mode, producing a strengthening of the Aleutian low, an enhancement of the southward Ekman and Oyashio cold advection, and thus a cooling SST in the North Pacific. In contrast, the TNA cooling produces a surface high extending to the eastern tropical North Pacific, inducing the northeasterly wind anomalies north, northerly cross-equatorial wind anomalies, and northwesterly wind anomalies south of the equator. This C-shape wind anomaly pattern generates an SST warming in the tropical southeastern Pacific, which eventually leads to an SST warming in the tropical central and western Pacific by the windevaporation-SST feedback. The tropical Pacific warming in turn leads to an SST cooling in the North Pacific by the Pacific North American teleconnection pattern.
1. Introduction Large-scale climate variations in the North Pacific and North Atlantic, such as the Pacific decadal oscillation (PDO) and Atlantic multidecadal oscillation (AMO), have been observed [e.g., Mantua et al., 1997; Enfield et al., 2001]. The PDO and AMO are associated with changes of climate and extreme weather events such as drought and hurricanes [e.g., Mantua and Hare, 2002; Hessl et al., 2004; Enfield et al., 2001; McCabe et al., 2004; Goldenberg et al., 2001; Wang et al., 2013; Zhang and Wang, 2012], and thus play an important role in economic, ecology, society, and people’s life.
C 2015. The Authors. V
This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.
ZHANG AND ZHAO
The observed large-scale AMO pattern resembles the simulated SST anomaly pattern induced by fluctuations of the Atlantic meridional overturning circulation (AMOC) [e.g., Delworth and Mann, 2000]. Both paleo data and coupled model simulations suggest that the AMO can influence the North Pacific variability such as the PDO. For example, higher oxygen levels off the California coast (indicating reduced upwelling and reduced California current) were synchronous with Greenland cooling stadials [Behl and Kennett, 1996], which are in turn hypothesized to link to the AMOC weakening. Using a hybrid version of the GFDL-CM2.1 climate model, Zhang and Delworth [2007] demonstrated that the AMO provides a source of multidecadal variability to the North Pacific. These research studies suggest that there is a link between the AMOC, Atlantic, and Pacific climate variations. The linkage between the AMOC, Atlantic, and Pacific climate variations can be also seen from the North Atlantic waterhosing experiments [e.g., Stouffer et al., 2006]. In response to a 1 Sv freshwater forcing in the subarctic North Atlantic, the AMOC slows down rapidly in all of the models. The associated climatic
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
739
Journal of Advances in Modeling Earth Systems
10.1002/2014MS000415
anomalies exhibit a large degree of similarity, especially in the Atlantic basin: a cooling of surface temperature and southward displacement of the ITCZ. The oceanic heat transport is fundamental for this response [e.g., Yang, 1999; Knutti et al., 2004; Wu et al., 2008]. Beyond the Atlantic basin, coupled models also exhibit a robust response over the North Pacific [e.g., Timmermann et al., 2007]. When the AMOC nearly shuts down, the North Pacific cools throughout the basin with maximum amplitude of 3–5 C along the oceanic frontal zone, and the westerly winds intensify, associated with deepening of the wintertime Aleutian low. Wang et al. [2014] analyzed 22 CMIP5 climate models and found a common pattern of global SST biases (defined as the SST difference between the model and observed SSTs). In the Northern Hemisphere, almost all of climate models show cold SST biases in the extratropical North Atlantic (ENA) and tropical North Atlantic (TNA) as well as in the North Pacific. These cold SST biases are linked with a weak simulated AMOC. However, the detailed mechanisms and processes for causing these cold biases need to be studied. In particular, it is not clear how the cold biases in the ENA and TNA affect the cold bias in the North Pacific. In this paper, we use both CMIP5 model outputs and coupled model experiments to show the potential cause of the SST biases in the Atlantic and their impacts on the state of the coupled system bias in the North Pacific. The results of these models suggest that the cause of the large systematic error in the Atlantic appears to be the systematic error in the Atlantic AMOC and its associated oceanic heat transport in the region, and these large cold biases can have a remote impact on the SST bias in the North Pacific. The cold bias in the ENA plays a more important role in the cold bias in the North Pacific than the cold bias in the TNA. In section 2, we describe the CMIP5 model outputs, dataset and coupled ocean-atmosphere model used in this study. Section 3 briefly summarizes and examines the SST biases in the Atlantic and North Pacific, their relationship with the AMOC in CMIP5 models. Section 4 investigates processes and sources of these SST biases. Section 5 uses coupled model experiments to show the possible pathways through which the ENA and TNA SST biases affect the North Pacific SST bias. Finally, section 6 discusses the implication of this study to coupled model bias in general and gives a brief summary.
2. CMIP5 Model Outputs, SST Data, and Coupled Model This study uses 22 coupled GCMs output data of the ‘‘historical’’ simulations provided to the Intergovernmental Panel on Climate Change Fifth Assessment Report (IPCC-AR5). The model data can be downloaded from the website of the Coupled Model Intercomparison Project phase 5 (CMIP5; Taylor et al., 2012) (http:// cmip-pcmdi.llnl.gov/cmip5/). The historical run is forced by observed atmospheric composition changes which reflect both anthropogenic (greenhouse gases, aerosols) and natural sources (volcanic influences, solar forcing and emissions of short-lived species and their precursors) and, for the first time, including time-evolving land cover. These historical runs cover much of the industrial period from the mid-19th century to the near present and are sometimes referred to as ‘‘twentieth century’’ simulations. The CMIP5 model sponsor, country, name, and letter denotation are shown in Table 1. Observational data set is used to validate the variability of coupled GCM simulations. SST data are from the NOAA Extended Reconstruction Sea Surface Temperature version 3 (ERSST v3) [Smith and Reynolds, 2004]. The temporal coverage is from January 1854 to the present and it has a spatial resolution on a 2 32 grid. The data can be obtained from http://www.ncdc.noaa.gov/oa/climate/research/sst/ersstv3.php. Because we are interested in only large-scale features, unless otherwise specified, all model outputs and ERSST data are interpolated to a 1o31o grid. To investigate the physical processes and mechanisms connecting the North Atlantic and North Pacific biases, we conduct three groups of partial coupling experiments (see section 5 for details) by using the fully coupled NCAR Community Earth System Model (CESM1.0.4, B_1850 case) [Meehl et al., 2012]. The atmospheric component in this version of CESM is the Community Atmospheric Model version 4 (CAM4) having 96 latitudes and 144 longitudes (1.9o32.5o) on a finite volume grid and 26 hybrid levels in the vertical. The ocean model is a version of the Parallel Ocean Program (POP) developed at Los Alamos National Lab with 1 horizontal resolution and enhanced meridional resolution (1/3 ) in the equatorial tropics and the North Atlantic and with 60 vertical levels. The fully coupled control simulation has been integrated for 1500 years without apparent climate shifts.
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
740
Journal of Advances in Modeling Earth Systems
10.1002/2014MS000415
Table 1. The 22 Climate Models Used in This Study and Their Sponsor, Country, Name, and Letter Denotation Sponsor, Country Commonwealth Scientific and Industrial Research Organisation (CSIRO), Australia Canadian Center for Climate Modeling and Analysis, Canada National Center for Atmospheric Research (NCAR), USA M et eo-France/Centre National de Recherches M et eorologiques, France Commonwealth Scientific and Industrial Research Organisation (CSIRO), Australia European Earth System Model, EU U.S. Department of Commerce/National Oceanic and Atmospheric Administration (NOAA)/Geophysical Fluid Dynamics Laboratory (GFDL), USA National Aeronautics and Space Administration (NASA)/Goddard Institute for Space Studies (GISS), USA Met office Hadley Centre, UK
Institute Pierre Simon Laplace, France
Center for Climate System Research (University of Tokyo), National Institute for Environmental Studies, and Frontier Research Center for Global Change (JAMSTEC), Japan Max Planck Institute for Meteorology, Germany Meteorological Research Institute, Japan Norwegian Climate Centre, Norway
Model Name
Letter Denotation
ACCESS1.0
M1
CanESM2 CCSM4 CNRM-CM5 CSIRO-Mk3.6.0
M2 M3 M4 M5
EC-EARTH GFDL-CM3 GFDL-ESM2G GFDL-ESM2M GISS-E2-H GISS-E2-R HadCM3 HadGEM2-CC HadGEM2-ES IPSL-CM5A-LR IPSL-CM5A-MR IPSL-CM5B-LR MIROC5
M6 M7 M8 M9 M10 M11 M12 M13 M14 M15 M16 M17 M18
MPI-ESM-LR MPI-ESM-P MRI-CGCM3 NorESM1-M
M19 M20 M21 M22
3. North Atlantic and Pacific SST Biases and the AMOC in CMIP5 Models Wang et al. [2014] found common patterns of global model SST biases in CMIP5 climate models. They also showed that regional model SST biases of the global patterns are commonly linked with AMOC simulations. Here we briefly summarize these results and perform some further analyses. As shown in Figures 1a–1d, both the North Atlantic and North Pacific show the cold SST biases with a magnitude of several degrees Celsius in all seasons. The largest model SST bias occurs in the ENA, with the cooling pattern resembling the cold phase of the AMO and featuring by a comma-like distribution with the largest amplitude over the subpolar region and an extension along the east side of the basin and into the tropics. If an AMOC index is defined by the maximum value of its streamfunction in the latitude band 20 –60 N, CMIP5 models show that the cold SST bias in the ENA is stronger when the AMOC is weaker, and vice versa, with an intermodel correlation of 0.85 between the SST bias and the AMOC index (Figure 2a). The relationship is consistent with the model SST response to the North Atlantic freshwater forcing [e.g., Zhang and Delworth, 2005; Wu et al., 2008], suggesting that the cold SST biases are due to the weak simulated AMOC. The North Pacific SST bias is also positively correlated with the ENA SST bias, with an intermodel correlation coefficient up to 0.61 (Figure 2b). This approximately linear relationship has a strong seasonality, with maximum (0.70) and minimum (0.48) correlations in the winter and summer, respectively (Figures 3a–3d). Additionally, the cold SST biases in the TNA positively covary with those in the North Pacific and ENA (Figures 2c and 2d). These correlations imply that the model SST biases in the North Atlantic and the North Pacific bias are connected with each other. To further verify the relationship of the AMOC with the model SST biases in the North Atlantic and North Pacific, we perform an intermodel singular value decomposition (SVD) analysis for the SST bias and the AMOC of 22 CMIP5 models over the Northern Hemisphere oceans. The first intermodel SVD mode (SVD1) explains 45% of total variance, and is characterized by a cold bias in the North Pacific and Atlantic Oceans associated with a weakened AMOC as manifested by the negative value of the AMOC stream function (Figures 4a and 4b). The two associated principal components have a positive correlation of 0.80 (Figure 4c). The maximum SST amplitudes occur in the North Atlantic subpolar region, which highly resembles the ensemble mean SST bias displayed in Figure 1. All of these results indicate that the North Pacific and Atlantic SST biases are connected with each other, both of which are associated with the AMOC. That is, if the
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
741
Journal of Advances in Modeling Earth Systems 60N
10.1002/2014MS000415
(a) MAM
40N 20N 0 20S 100E 60N
140E 180W 140W 100W
60W
20W
20E
60W
20W
20E
60W
20W
20E
60W
20W
20E
(b) JJA
40N 20N 0 20S 100E 60N
140E 180W 140W 100W
(c) SON
40N 20N 0 20S 100E 60N
140E 180W 140W 100W
(d) DJF
40N 20N 0 20S 100E
140E 180W 140W 100W
−3 −2.5 −2 −1.5 −1 −0.5
0
0.5
1
1.5
2
2.5
3
Figure 1. The seasonal variation of SST biases (oC) in 22 climate models. Shown are the SST biases during the boreal (a) spring (MAM), (b) summer (JJA), (c) fall (SON), and winter (DJF). The bias is calculated as the SST difference between the simulated SST and observed ERSST. The dots denote where at least 18 of 22 models have the same sign of the SST bias. The boxes in Figure 1a represent the focused regions: the North Atlantic (55 W–15 W, 30 N–55 N), the North Pacific (150 E–140 W, 15 N–45 N), and the tropical North Atlantic (80 W–25 W, 5 N–25 N).
model cannot reasonably simulate the AMOC mean state, they will generate at least partially SST bias in the North Atlantic. Then, the North Atlantic bias can induce the model biases in the TNA and North Pacific by the ocean and atmosphere bridges. In the next sections, we will examine the possible sources of biases in the Atlantic and North Pacific Oceans.
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
742
Journal of Advances in Modeling Earth Systems
R = 0.85
0 SST bias in the ENA
(b) ENA versus North Pacific
(a) ENA versus AMOC SST bias in the North Pacific
0.5 −0.5 −1.0 −1.5 −2.0 −2.5 −3.0 −3.5
0.5 0
ACCESS1.0
R = 0.61
CanESM2
−0.5
CCSM4
−1.0
CNRM−CM5 CSIRO−Mk3.6.0
−1.5
EC−EARTH
−2.0
GFDL−CM3
−2.5
GFDL−ESM2G
−3.0
GFDL−ESM2M
−4.0
0.4
8 10 12 14 16 18 20 22 24 26 AMOC strength
−4 −3.5 −3 −2.5 −2 −1.5 −1 −0.5 SST bias in the ENA
(c) TNA versus North Pacific
0
0
(d) TNA versus ENA 1
R = 0.4 SST bias in the ENA
SST bias in the North Pacific
6
−0.4 −0.8 −1.2 −1.6 −2 −2.4
10.1002/2014MS000415
R = 0.3 0
GISS−E2−H GISS−E2−R HadCM3 HadGEM2−CC HadGEM2−ES IPSL−CM5A−LR IPSL_CM5A−MR
−1 −2 −3
IPSL−CM5B−LR MIROC5 MPI−ESM−LR MPI−ESM−P MRI−CGCM3
−2.8 −3.2 −2.5 −2.0 −1.5 −1.0 −0.5 0 0.5 SST bias in the tropical North Atlantic
−4 −2.5 −2.0 −1.5 −1.0 −0.5 0 0.5 SST bias in the tropical North Atlantic
NorESM1−M
Figure 2. Relationship of SST bias in the North Atlantic with the AMOC and SST bias in the North Pacific. (a) Scatterplot of the annual mean SST bias (oC) in the extratropical North Atlantic (ENA) versus the AMOC strength (Sv). (b) Scatterplot of the annual mean SST bias in the ENA versus the annual mean SST bias in the North Pacific. (c) Scatterplot of the annual mean SST bias in the TNA versus the annual SST bias in the North Pacific. (d) Scatterplot of the annual mean SST bias in the TNA versus the annual mean SST bias in the ENA. The intermodel correlation R is shown in the left top side of each plot.
4. Processes Related to the SST Biases in the North Atlantic and Pacific To further examine the processes or sources of the model SST biases, we choose three models (M15, M17, and M21) with the coldest SST biases in the North Atlantic and Pacific (CNPAC-CATL) and three models (M3, M10, and M22) with the relatively warm SST biases in the North Pacific and Atlantic (WNPAC-WATL). This selection can be easily done from the SST principal components in Figure 4c or from the scatter plot of North Atlantic and North Pacific SST biases (Figure 2b). As expected, the annual mean SST difference between the CNPAC-CATL and WNPAC-WATL models is characterized by cooling SST anomalies in the North Atlantic and North Pacific (Figure 5a). For long-term mean SST, the mixed-layer heat budget can be expressed as follows: Q 1 Do 5 0, where Q denotes the net heat flux (positive downward) and Do represents ocean dynamical effects due to the three-dimensional advection and mixing. Figure 5b shows that the cooling SST over the ENA Ocean is damped by the surface net heat flux, implying the important role of ocean dynamics. As expected, the AMOC and its associated northward heat transport are significantly reduced in the CNPACCATL models compared to the WNPAC-WATL models (Figures 6a and 6b), which in turn produces the cooling SST bias in the ENA region (Figure 5a). The subpolar gyre over the North Atlantic also shows a strengthening as indicated by the negative sea surface height (SSH) anomaly (Figure 5b), which is a typical fingerprint of the weakened AMOC [Zhang, 2008]. This spin up of subpolar gyre further favors the cooling SST anomaly over the ENA due to the anomalous cold advection from the high latitude and the anomalous vertical cold advection and mixing. In contrast, the TNA cooling SST bias is not only attributed to the decrease of the northward heat transport, but also contributed by the heat flux as a result of the wind changes (Figures 5a and 5b).
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
743
Journal of Advances in Modeling Earth Systems
10.1002/2014MS000415
SST Bias in the Extratropical North Atlantic Versus North Pacific SST bias in the North Pacific
(a) MAM
(b) JJA
0.5
ACCESS1.0
0
0
R = 0.65
−0.5
−0.5
CanESM2
R = 0.48
CCSM4
−1.0
−1.0
CNRM−CM5
−1.5
−1.5 −2.0 −2.5 −3.0
CSIRO−Mk3.6.0
−2.0
EC−EARTH
−2.5
GFDL−CM3
−3.0
GFDL−ESM2G
−3.5
GFDL−ESM2M
−3.5 −4
−3
−2
−1
0
1
−3.5 −3 −2.5 −2 −1.5 −1 −0.5
0
SST bias in the Extratropical North Atlantic SST bias in the Extratropical North Atlantic SST bias in the North Pacific
(c) SON 0.5 0
0.5
−0.5 −1.0
HadGEM2−CC HadGEM2−ES
R = 0.70
IPSL−CM5A−LR
−0.5
IPSL_CM5A−MR
−1.0
IPSL−CM5B−LR
−1.5
MIROC5
−1.5
−2.0
GISS−E2−R HadCM3
(d) DJF
0
R = 0.55
GISS−E2−H
MPI−ESM−LR
−2.5
−2.0
MPI−ESM−P
−3.0
−2.5
MRI−CGCM3
−3.5 −3.5 −3 −2.5 −2 −1.5 −1 −0.5
0
−3.0 −4
NorESM1−M
−3
−2
−1
0
1
SST bias in the Extratropical North Atlantic SST bias in the Extratropical North Atlantic Figure 3. Relationship between the SST bias (oC) in the extratropical North Atlantic and Pacific. Shown are scatterplots of the SST bias in the North Pacific versus extratropical North Atlantic during the boreal (a) spring (MAM), (b) summer (JJA), (c) fall (SON), and winter (DJF). The intermodel correlation R is shown in the left top side of each plot.
As presented in Figure 5c, we perform a detailed heat budget analysis in the western and central North Pacific Ocean. The cooling SST bias is largely associated with the temperature advection by the anomalous meridional current, while the contribution from the surface heat flux is negative. This suggests that there is an anomalous southward current in the western and central midlatitudes, which can be also implied from the surface wind and SSH anomaly (Figures 5a and 5b). The surface wind in the North Pacific is characterized by a westerly anomaly in the midlatitude (Figure 5a), which is in favor of generating a southward cold Ekman transport. Moreover, the subpolar gyre over the North Pacific is significantly strengthened and shifts southward in response to the cyclonic wind (Figures 5a and 5b), which leads to a strengthening and southward movement of the Oyashio and thus generates a cooling SST anomaly in the western North Pacific. In addition, the temperature advections by the anomalous zonal and vertical currents also play a positive role. This is because the cyclonic wind strengthens the subpolar gyre and the associated eastward current in the southern branch, which in turn brings more western subpolar cold water to the central North Pacific. One the other hand, the cyclonic wind strengthens the Ekman pumping, which induces an anomalous upward current and therefore a cooling SST over the North Pacific. Overall, we can find that the cooling SST bias in the ENA is attributed to the reduced northward heat transport associated with the weak simulated AMOC. While in the TNA, both the northward heat transport and the reduced surface heat flux contribute to the cooling SST bias. Over the North Pacific, the cooling SST bias is primarily determined by the temperature advection by the anomalous southward current as a result of the surface cyclonic wind. However, the potential pathway by which the North Atlantic SST bias affects the North Pacific is still unknown. To investigate the possible mechanisms, we calculate the geopotential height differences at low (850 mb), middle (500 mb), and high (250 mb) heights, as exhibited in Figures 7a–7c. The anomalous atmosphere circulation displays an equivalent barotropic structure, with a low pressure belt in
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
744
Journal of Advances in Modeling Earth Systems 60N
10.1002/2014MS000415
(a) The first SVD mode of SST bias
50N 40N 30N 20N 10N 0 120E −2
0m
150E −1.6
180W −1.2
150W
−0.8
120W
−0.4
90W
0
60W
0.4
0.8
30W 1.2
1.6
0 2
(b) The first SVD mode of MOC
1000m 2000m 3000m 4000m 5000m 30S −3
20S −2.5
10S −2
−1.5
EQ −1
10N −0.5
20N 0
0.5
30N 1
40N 1.5
50N 2
2.5
60N 3
(c) SVD time series (correlation=0.8) 3 2
3
SST MOC
2
1
1
0
0
−1
−1
−2
−2 M1 M2 M3 M4 M5 M6 M7 M8 M9 M10 M11 M12 M13 M14 M15 M16 M17 M18 M19 M20 M21 M22
Figure 4. Intermodel SVD analysis of SST bias in the Northern Hemisphere and the Atlantic meridional overturning circulation (AMOC). Shown are (a) the first SVD mode spatial pattern of the SST bias in the Northern Hemisphere, (b) the first SVD mode spatial pattern of the AMOC, and (c) the corresponding SST and AMOC principal components. The x axis in Figure 4c represents different models (Table 1). The intermodel SVD analysis is performed by using different climate models as the time dimension. The first SVD mode accounts for 45% of the total variance.
the midlatitudes to the subtropics and a high pressure cap over the entire Arctic, resembling the negative phase of the Northern Hemisphere annular mode [Thompson and Wallace, 1998]. The low pressure belt in the North Atlantic and North Pacific is also seen in the surface wind (Figure 5a), which shows a reinforcement of the Aleutian low and Iceland low. This suggests that the North Atlantic cooling SST bias can be transmitted to the North Pacific through the annular mode and vice versa. Therefore, it is expected that if the model fails in simulating the AMOC, it will generate at least partially cooling biases in the ENA and TNA Oceans. These North Atlantic cooling SST biases can be further propagated to the North Pacific Ocean and then can be strengthened by the local air-sea coupling, which ultimately may positively feedback to the Atlantic Ocean. In other words, the Northern Hemisphere SST biases in coupled models are largely due to the AMOC simulation. Previous studies argue that the physical processes involved in the impacts of the North Atlantic SST on the North Pacific SST can be through two possible paths: the midlatitude atmosphere [e.g., Zhang and Delworth, 2007; Wu et al., 2008; Kang et al., 2014] and the tropical atmosphere and oceans [e.g., Timmermann et al., 2005; Dong et al., 2006; Xie et al., 2008; Wu et al., 2005, 2007; Zhang et al., 2011b; Okumura et al., 2009; Ham et al., 2013; Zhang et al., 2014]. Zhang and Delworth [2007] suggested that a warm (cold) phase of the AMO reduces (enhances) the meridional SST gradients in the midlatitude, which further reduces (increases) the surface atmosphere eddy heat transport and upper level vorticity flux. These cause northward (southward) shifts of midlatitude westerly winds both over the North Atlantic and North Pacific. The northward (southward) shift of the westerly wind over the North Pacific leads to a
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
745
Journal of Advances in Modeling Earth Systems
10.1002/2014MS000415
(a) SST and surface wind stress differences
o
60 N o
50 N o
40 N
0.06N/m2
o
30 N o
20 N o 10 N
0o o 120 E
−4
o
o
150 E
−3.5
−3
o
180 W
−2.5
−2
o
150 W
−1.5
−1
o
120 W
−0.5
0
o
90 W
0.5
1
o
60 W
1.5
2
o
30 W
2.5
3
0
3.5
4
(b) Net heat flux and Sea surface height differences
o
60 N 50oN 40oN o 30 N o
20 N o
10 N o
0
o
o
120 E
o
150 E
−50
−40
180 W
−30
−20
o
150 W
−10
o
120 W
0
o
o
90 W
10
o
60 W
20
o
30 W
30
0
40
50
(c) Heat budget over the North Pacific Ocean 3 1 −1 −3 −5 Ta
−uTx
−vTy
−Utx
−Vty
−wTz
−Wtz
HFLX
Mix
Figure 5. Shown are the (a) SST ( C, shading) and surface wind stress (N/m2, vectors), (b) net surface heat flux (positive downward, W/m2, shading) and sea surface height (SSH, contour interval is 0.01 m) differences between the CATL-CNPAC and WATL-WNPAC models, and (c) heat budget analysis over the North Pacific region (see box in Figure 5a). Ta denotes the mixed layer temperature. Here, the mixed layer depth is defined as the depth where the potential density is larger than that in the surface by 0.125 kg/m3. Mix denotes the mixing which is estimated as the residual. HFLX denotes the surface net heat flux. The rest of the terms describe advection, with capital (U, V, W, T) and lower case (u, v, w, t) letters denoting mean and anomaly, respectively. -uTx, -vTy, -Utx, -Vty, -wTz, and -Wtz represent term 2U’ @T =@x, 2V ’ @T =@y, 2U@T ’ =@x, 2V @T ’ =@y, 2W ’ @T =@z, and 2W @T ’ =@z, respectively. Units for the temperature Ta and the heat budget terms are C and 10 W/m2, respectively.
weakness (strengthening) of the Aleutian low, which further causes the warm (cold) SST anomalies in the North Pacific. Wu et al. [2008] argued that the strong North Atlantic cooling (warming) accelerates (decelerates) the westerly winds over the North Pacific by triggering an annular mode, and thus a cooling (warming) SST in the North Pacific. Both of these arguments are referred to the midlatitude pathway, although the former focuses on the effect of transient eddies on the Northern Hemisphere westerly, while the latter focuses on the atmosphere teleconnection. This midlatitude pathway is similar to our composite analysis from CMIP5 models to some degree. Wu et al. [2005, 2007] and Zhang et al. [2011b] further demonstrated that a warming of the TNA could lead to a La Nina-like SST response in the tropical Pacific, which in turn produces a North Pacific warming through the Pacific North American (PNA) teleconnection as suggested by Alexander et al. [2002]. Okumura et al. [2009] also showed a direct atmospheric pathway for the TNA to affect the Aleutian low and North Pacific SST. Although our analysis suggests that the North Pacific cooling SST bias can be remotely or at least partially affected by the North Atlantic SST bias through the annular mode, it is still not clear which path is more important and which path favors triggering the annular mode: the midlatitude path or the tropical path? Therefore, in the next section, we attempt to conduct a series of numerical experiments to examine which path is dominant to the North Pacific cooling SST anomaly.
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
746
Journal of Advances in Modeling Earth Systems 0m
10.1002/2014MS000415
(a) AMOC difference
1000m 2000m 3000m 4000m 5000m 30S −12
15S −9
EQ −6
−3
15N 0
30N 3
45N 6
60N 9
12
(b) Northward heat transport 1
WNPAC−WATL CNPAC−CATL
0.8 0.6 0.4 0.2 0 30S
15S
EQ
15N
30N
45N
60N
Figure 6. Shown are the (a) AMOC difference (Sv) between the CATL-CNPAC and WATL-WNPAC models, and (b) northward heat transport (PW) in WATL-WNPAC and CATL-CNPAC models, respectively.
5. Pathways Connecting the North Atlantic and Pacific 5.1. Midlatitude Pathway To examine the potential pathways connecting the model SST biases in the North Atlantic and North Pacific, we conduct two groups of partial coupling experiments by using the fully coupled NCAR Community Earth System Model (CESM1.0.4), with each group consisting of two parallel runs. The first group of partial coupling experiments is designed to test the midlatitude path from the ENA to the North Pacific, which is named as the PC_ENA. The first group includes a control restoring run and an ENA SST bias run. In the control restoring run, the entire tropical (20 S–20 N, 0oE–0oW) and ENA (20 N–60 N, east coast to west coast) SSTs are restored to their climatological SST seasonal cycle at every time integration step. The ENA SST bias run is configured as the same as the control restoring run but with the ENA region restoring to the CESM1.0.4 climatological SST seasonal cycle added with the monthly SST bias in CMIP5 models (Figure 1). Both the two experiments are integrated for 100 years. The difference averaged in the last 30 years is taken as the response. Overall, in the first group experiments, the tropical regions are decoupled, and thus the forced SST bias in the ENA can only transmit to the North Pacific region through the midlatitude path. Figure 8a shows the SST response in the PC_ENA run. The cooling over the North Pacific Ocean is a prominent feature in this experiment. The area-averaged cooling anomaly over the North Pacific is about 1.12 C, which accounts for 62% of the cooling bias in the coupled model (Figure 1). The cooling over the western North Pacific also exhibits a seasonal dependence with the maximum and minimum in late summer and winter, respectively (not shown). This implies that a large amount of the cooling bias in the North Pacific Ocean can be remotely impacted by the ENA bias. Over the North Pacific Ocean, the
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
747
Journal of Advances in Modeling Earth Systems (a) 850mb E
0
30
30
30 20
W
60
60
E
W
12
12 0W
90E
90W
10
0E
15
0E
180W
0 −10 −20
0W
−30
W
40
15
(b) 500mb 0
30
20 0 −20
12
0E
12
0W
90E
90W
W
60
60
E
E 30
15
0E
180W
0W
−40
W
60
15
(c) 250mb 30
40 20 90W
60
0
W 60
E
E 30
10.1002/2014MS000415
atmospheric response tends to be dominated by the equivalent barotropic low (Figures 9a–9c). This circumglobal teleconnection is a reminiscent of the annular mode, the leading mode of Northern Hemispheric atmosphere variability in the winter [Thompson and Wallace, 1998]. In response to the barotropic low, the surface wind anomaly is featured by a cyclone wind over the North Pacific Ocean and an intensified westerly in the midlatitude (Figure 10a). The midlatitude westerly anomaly first cools the ocean by reinforcing the southward Ekman cold advection and the surface turbulent heat flux cooling (not shown). The anomalous positive wind stress curl (cyclonic wind) then drives a strengthening and southward shift of the subpolar gyre after the adjustment of westward Rossby waves as indicated by the SSH anomaly (contours in Figure 10a), which in turn generates a SST cooling in the western North Pacific by the anomalous cold advection from the high latitude to the low latitude. The strengthening and southward shift of the subpolar gyre in the North Pacific is further seen from the barotropic streamfunction response (Figure 10b). The barotropic component of the ocean circulation response is predominantly due to changes of wind stress in response to the cooling bias. To demonstrate that, we calculate the Sverdrup transport based on the wind stress anomalies (Figure 10c). The differences between the Sverdrup transport and the simulated barotropic transport are less significant, except in the western boundary regions as expected.
12
0E 12
ZHANG AND ZHAO
0W
90E
In the North Pacific Ocean, the SST and geopotential height responses suggest that there is a local positive feedback, which is also pro−20 posed by the previous studies [Zhang and Delworth, 2005; Wu et al., 2008; Zhang et al., −40 2011a, 2011b; Zhang and Wu, 2012]. On one 15 hand, the extratropical atmospheric teleconW 0E −60 50 1 nection can transmit the effects of the North 180W Atlantic cooling to the midlatitude North Figure 7. Shown are the (a) 850 mb, (b) 500 mb, and (c) 250 mb geopoPacific by the global nature of the annular tential height (m) differences between the CATL-CNPAC and WATLmode. On the other hand, the cold SST in the WNPAC models. Kuroshio-Oyashio extension (KOE) region can feedback to the atmosphere, leading to an even stronger strengthening of the Aleutian Low. Therefore, once the North Pacific cold SST is initiated by the North Atlantic bias, then it can be amplified by the local positive air-sea feedback. In consequence, we conclude that the North Pacific bias can be remotely impacted by the North Atlantic bias by the midlatitude path, which can explain nearly 62% of the SST cooling bias in the North Pacific.
0
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
748
Journal of Advances in Modeling Earth Systems 60oN
10.1002/2014MS000415
(a) SST response in PC_ENA
40oN o
20 N 0o 20oS 40oS
0o
−2
30oE
60oE
−1.5
90oE
−1
120oE
150oE 180oW 150oW 120oW
−0.5
0
0.5
90oW
60oW
1
30oW
1.5
0o
2
(b) SST response in PC_TNA
o
60 N 40oN 20oN o
0
20oS 40oS
o
0
−2
o
30 E
−1.5
o
60 E
o
90 E
−1
o
120 E
o
o
o
o
150 E 180 W 150 W 120 W
−0.5
0
0.5
o
90 W
1
o
60 W
o
30 W
1.5
o
0
2
Figure 8. Annual mean SST response ( C) in the (a) PC_ENA and the (b) PC_TNA experiments.
5.2. Tropical Pathway We conduct the second group of partial coupling experiments to examine the relationship between the cooling biases in the TNA and North Pacific, which is named as the PC_TNA. The first experiment is also the control restoring run. The ocean and atmosphere are fully coupled except in the TNA region where we restore the model-produced SST to the climatological SST seasonal cycle at every integration step. The second experiment is configured as the same as the first experiment but with the TNA region restoring to the CESM1.0.4 climatological SST added with the monthly CMIP5 model SST cooling bias in the region of 0 – 20 N and from east to west coast (Figure 1). The difference averaged in the last 30 years is taken as the response. As shown in Figure 8b, the SST response in the North Pacific is characterized by a horseshoe-like SST pattern, with a cooling in the northwestern and central North Pacific and a warming in the east extending northwest into the subpolar ocean and southwest into the tropics. The magnitude of the area averaged cooling SST over the North Pacific Ocean is about 0.6 C, which accounts for 33% of the North Pacific cooling bias in CMIP5 models presented in Figure 1. This suggests that the SST bias in the North Pacific Ocean is partially attributed to the remote effect of the SST biases in the TNA Ocean. The North Pacific cooling anomaly triggered by the TNA bias is mainly through the tropical Pacific region. As seen in Figure 8b, the tropical Pacific is characterized by a warm SST extending from the southeastern Pacific (SEP) region to the western tropical Pacific. According to Alexander et al. [2002], the tropical Pacific warming can produce a cooling of the North Pacific Ocean by triggering the PNA teleconnection pattern [Horel and Wallace, 1981]. This can be seen from the geopotential response (Figures 9d– 9f). Figures 9d–9f show that the structure of the middle latitude geopotential response is similar to the PNA pattern, with a low pressure in the North Pacific, a high pressure in North America, and a low pressure in the south. The tropical Pacific warming is generated as follows: The cold TNA bias suppresses convection and rainfall over the AWP region (Figure 11b), producing a surface high extending to the eastern tropical North Pacific (Figure 11c). Over the eastern tropical Pacific, this anomalous high induces northeasterly wind anomalies north, northerly cross-equatorial winds, and northwesterly winds south of the equator (Figure 11a). This
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
749
Journal of Advances in Modeling Earth Systems
30 o W
o
30
20
90oW
E
15
0 oE
−30
15
40
o
30
90 W
o
15
0 oE
−40
60
0
30 o W
60
90 W
o
20
o
E
−60
o
−40
40
0 oW
−20 0
0 oW
0
12
90 W
20
12
12
30
E
90oE
E 180 W
o 0W
15
o
W
W
o
0
12 o
−40
o
o
o
0 oE
180 W
o 0W
15
60
40
60
90 E
15
−20
o
30 o W
60 o E
30
E
o
60 o E
o
0
(f) 250mb
(c) 250mb o 0
20
0 oW
90 W
o
0 oW 12
E
180 W
o 0W
15
40
W
o
0
12
o
30 o W
o
oE
0
12
−20
0o
90 E
0
E
o
20
W
o
0 oE
−20 −30
60
o
60
90 E
15
−10
12
30 o W
60 o E
30
E
180 W
o 0W
15
60 o E
o
o
0
(e) 500mb
(b) 500mb o 0
10
0 oW
90oW
0 oW
o
0
12
12
o
180 W
o 0W
−10 −20
20
W
10 0
30
90 E
o
30 o W o
W
oE
0
12
0 oE
E
0
60
o
o
90 E
15
PC_TNA
o
30 60
60 o E
E
(d) 850mb
PC_ENA
12
o
30
0
60 o E
(a) 850mb o
10.1002/2014MS000415
15
0 oE
o
180 W
o 0W
15
0 −20 −40 −60
Figure 9. Winter (DJF) mean (a) 850 mb, (b) 500 mb, and (c) 250 mb geopotential height responses (m) in the PC_ENA run. (d–e) Same as (a–c) but for the PC_TNA experiment.
C-shape wind anomaly generates a dipole SST anomaly through changes in evaporation. This coupled WES feedback, which has been extensively studied in relation to tropical Atlantic variability [e.g., Xie and Carton, 2004], acts to amplify the SST dipole. The formation of crossing interbasin wind is primarily due to the weakening of the Hadley-type circulation from the AWP region to the SEP as a result of the cooling bias imposed over the TNA region. As exhibited in Figures 12a and 12b, the TNA cooling is associated with a divergent
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
750
Journal of Advances in Modeling Earth Systems 60oN
10.1002/2014MS000415
(a) SST, SSH and surface wind stress responses in PC_ENA
50oN 0.03N/m2 o
40 N 30oN 20oN o 130 E
−2
o
o
150 E
−1.5
o
170 E
−1
o
170 W
150 W
0
0.5
−0.5
o
o
130 W
1
110 W
1.5
2
(b) Barotropical streamfunction response in PC_ENA
o
60 N o
50 N 40oN 30oN 20oN o 130 E 60oN
o
150 E
o
o
170 E
170 W
o
o
150 W
130 W
o
110 W
(c) Sverdrup streamfunction response in PC_ENA
o
50 N 40oN o
30 N o
20 N o 130 E
−5
o
150 E
−4
−3
o
o
170 E
−2
170 W
−1
0
o
o
150 W
1
130 W
2
3
o
110 W
4
5
Figure 10. Shown are the (a) SST (oC, shading), surface wind stress (N/m2, vectors), and sea surface height (SSH, contour interval is 0.01 m), (b) barotropic streamfunction (Sv), and (c) Sverdrup streamfunction (Sv) responses in the PC_ENA run.
circulation in the low troposphere that crosses the equator into the South Pacific and vice versa for the upper troposphere. This meridional circulation reduces the South Pacific subtropical anticyclone (Figure 11c) and the associated subsidence, which in turn leads to a reduction of the low clouds, the easterly trade wind, and thus an increase of the SST. The warming of the SEP Ocean is further transmitted to the western tropical Pacific Ocean by the WES mechanism [e.g., Xie, 1996]. As exhibited in Figure. 13, the SEP warming induces a northwesterly wind on its west side. These northwesterly wind anomalies decelerate the mean southeasterly winds, reducing oceanic evaporative heat loss to induce warming on the northwestward flank of the original SST warming and causing the coupled SST-wind pattern to propagate equatorward. Eventually, the tropical Pacific Ocean is occupied by the warming SST in the boreal winter, which triggers the PNA teleconnection. The PNA teleconnection can be also seen from the barotropic streamfunction response (Figure 12c). Here, the barotropic and baroclinic streamfunctions are calculated as wt 5ðw850mb 1w250mb Þ=2 and wc 5ðw850mb 2w250mb Þ=2, respectively. In response to the PNA teleconnection, the surface wind in the North ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
751
Journal of Advances in Modeling Earth Systems
10.1002/2014MS000415
(a) Surface windstress response 60N 2
40N
0.02N/m
20N 0 20S 40S 0
30E
60E
90E
120E 150E 180W 150W 120W
90W
60W
30W
0
90W
60W
30W
0
(b) Precipitation response 60N 40N 20N 0 20S 40S 0
−2
30E
60E
−1.5
90E
−1
120E 150E 180W 150W 120W
−0.5
0
0.5
1
1.5
2
(c) Sea level pressure response 60N 40N 20N 0 20S 40S 0
−120
30E
−90
60E
90E
−60
120E 150E 180W 150W 120W
−30
0
30
90W
60
60W
90
30W
0
120
Figure 11. Annual mean (a) surface wind stress (N/m2), (b) precipitation (mm/day), and (c) sevel level pressure (mb) responses in the PC_TNA run.
Pacific Ocean is cyclonic, which generates a horseshoe-like SST pattern, with a cooling in the northwestern and central North Pacific and a warming in the east extending northwest into the subpolar ocean and the southwest into the subtropics (Figure 8b). This SST response is mainly due to the surface heat flux and temperature advection by the anomalous meridional current (not shown). The eastern warming is also propagated to the western tropical Pacific by the positive WES feedback [e.g., Vimont et al., 2009]. To demonstrate this process, we examine the seasonal development of SST, surface wind, and latent heat flux (Figure 14). In March and April, the warming due to the warm northward advection is primarily limited to the north of 10 N and east of 180 W (Figures 14a and 14b). Associated with the warm anomalies, anomalous cyclonic winds develop in the east, which decelerate the northeast trades and thus reduce the oceanic latent heat loss (Figures 14f and 14g). Note that both wind and downward heat flux anomalies have extended to the lower latitude and more westward, which favor an equatorward extension of the warm anomalies in the following months (Figure 14c). At the same time, westerly wind anomalies develop in the western tropical Pacific and reduce the latent heat loss (Figure 14h). This can further substantiate equatorward progression of the warm anomalies (Figures 14d and 14e). Therefore, the northeastern Pacific warming can be eventually transmitted to the tropical Pacific Ocean, which in turn reinforces the subsequent PNA teleconncetion. Therefore, the air-sea coupling in the tropical-extratropical Pacific acts as a positive feedback to amplify the North Pacific cooling. In addition to the barotropic response, the atmosphere response also has the baroclinic component. The baroclinic streamfunction response shows a pair of anticyclones: one in the TNA and northeastern Pacific and the other in the SEP (Figure 12d). This model response is largely consistent with Gill’s [1980] solution to
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
752
Journal of Advances in Modeling Earth Systems (a) 850mb velocity potential and divergent wind
(c) Barotropic Stream Function
70N
70N
50N
50N
30N
30N
10N
10N
10S
10S
30S
30S
50S
50S
70S 100E
150E
160W
110W
60W
10W
10.1002/2014MS000415
70S 100E
150E
160W
110W
60W
10W
60W
10W
1m/s
(b) 250mb velocity potential and divergent wind
(d) Baroclinic Stream Function
70N
70N
50N
50N
30N
30N
10N
10N
10S
10S
30S
30S
50S
50S
70S 100E
150E
160W
110W
60W
10W
70S 100E
150E
160W
110W
1m/s −1.5 −1.2 −0.9 −0.6 −0.3
0
0.3
0.6
0.9
1.2
1.5
−2.5 −2 −1.5 −1 −0.5
0
0.5
1
1.5
2
2.5
Figure 12. Annual mean (a) 850 mb velocity potential and divergent wind, (b) 250 mb velocity potential and divergent wind, (c) barotropic streamfunction, and (d) baroclinic streamfunction responses in the PC_TNA run.
a cooling anomaly slightly north of the equator [Heckley and Gill, 1984]. There is also a pair of cyclones in the western Pacific. This can be interpreted as follows: The northeastern Pacific warming anomalies can gradually propagate to the western tropical Pacific due to the positive WES feedbacks (Figure 14), which induces a large amount of precipitation there (Figure 11b). The heating anomaly over the western tropical Pacific eventually induces a baroclinic response represented by a pair of cyclones. Note that the TNA cooling bias may influence the North Pacific SST through the TNA-North Atlantic-North Pacific path by exciting the barotropic Rossby waves [e.g., Hoskins and Karoly, 1981; Okumura et al., 2009]. As presented in Figure 12c, the barotropic component also shows a pattern of alternating high and low centers from the TNA to high latitudes, which is consistent with the pure atmosphere model response [e.g., Wang et al., 2010; Okumura et al., 2009]. To examine the importance of this path, we perform the third group experiment named PC_TNAxTPAC run, which is configured as the same as the second group experiment, but with the tropical Pacific Ocean (20oS–20oN, west coast to east coast) decoupled by restoring the model-generated SST to the climatological SST cycle. Figure 15 shows the SST response in this third group experiment. Although the North Pacific Ocean features a cooling SST anomaly, its magnitude is too small, particularly in the western midlatitude of the North Pacific. Therefore, the North Pacific SST cooling response to the TNA cold bias is primarily through the tropical Pacific path.
6. Discussion and Summary In this paper, we report the results from 22 coupled model simulations provided from the Coupled Model Intercomparison Projection (CMIP) phase 5 multimodel ensembles. Our intercomparison study of
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
753
Journal of Advances in Modeling Earth Systems SST
10.1002/2014MS000415
Heat flux and wind stress
(a) July
(f) July
0
0
20S
20S
40S 40S 120E 150E 180W 150W 120W 90W 60W 120E 150E 180W 150W 120W 90W 60W (b) August
(g) August
0
0
20S
20S
40S 40S 120E 150E 180W 150W 120W 90W 60W 120E 150E 180W 150W 120W 90W 60W (c) September
(h) September
0
0
20S
20S
40S 40S 120E 150E 180W 150W 120W 90W 60W 120E 150E 180W 150W 120W 90W 60W (d) October
(i) October
0
0
20S
20S
40S 40S 120E 150E 180W 150W 120W 90W 60W 120E 150E 180W 150W 120W 90W 60W (e) November
(j) November
0
0
20S
20S
40S 40S 120E 150E 180W 150W 120W 90W 60W 120E 150E 180W 150W 120W 90W 60W 2
0.03N/m
−0.8 −0.6 −0.4 −0.2
0
0.2 0.4 0.6 0.8
−20 −15 −10 −5
0
5
10
15
20
Figure 13. (a, f) July, (b, g) August, (c, h) September, (d, i) October, and (e, j) November SST (left plot), turbulent heat flux and wind stress (right plot) responses in the PC_TNA run. The units for the SST, heat flux (positive downward) and wind stress are C, W/m2, and N/m2, respectively.
coupled climate models shows that virtually all the current generation of coupled models without flux correction have significant, but similar biases in the Northern Hemisphere Oceans. These biases manifest themselves as the SST in the North Pacific, TNA and ENA Oceans being too cold. The magnitude of these biases can be as large as 3oC or more, particularly in the ENA Ocean, resulting in a significant distortion of the semiannual cycle of the northern hemisphere SST in these models. We explored the potential cause of the large biases in the Atlantic and its impact on the state of the coupled system bias in the North Pacific. The results of these models suggest that the cause of the large systematic error in the Atlantic appears to be the systematic error in the AMOC and its associated oceanic heat transport in the region, and these large cold biases can have a remote impact on the SST bias in the North Pacific Ocean. They contribute positively to the cold bias in the North Pacific Ocean. The study suggests that effort of
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
754
Journal of Advances in Modeling Earth Systems
60N
SST
(a) March
Heat flux and wind stress
60N
40N
40N
20N
20N
0 120E 150E (b) April 60N
180W 150W 120W
90W
40N
20N
20N 180W 150W 120W
90W
0 120E 150E (h) May 60N
40N
40N
20N
20N
0 120E 150E (d) June 60N
180W 150W 120W
90W
0 120E 150E (i) June 60N
40N
40N
20N
20N
0 120E 150E (e) July 60N
180W 150W 120W
90W
0 120E 150E (j) July 60N
40N
40N
20N
20N
0 120E
150E
180W 150W 120W
90W
(f) March
0 120E 150E (g) April 60N
40N
0 120E 150E (c) May 60N
10.1002/2014MS000415
0 120E
150E
180W 150W 120W
90W
180W 150W 120W
90W
180W 150W 120W
90W
180W 150W 120W
90W
180W 150W 120W
90W
0.03N/m2
−1.5
−1
−0.5
0
0.5
1
1.5
−20 −15 −10 −5
0
5
10
15
20
Figure 14. Same as Figure 13 but for the March to July SST, heat flux and wind stress responses in the North Pacific.
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
755
Journal of Advances in Modeling Earth Systems
10.1002/2014MS000415
reducing coupled model biases in the North Hemisphere Oceans should take into the consideration not only the local processes, but also the remote influence.
(a) SST response in PC_TNAxTPAC 60N 40N 20N 0 20S 40S 0
30E
60E
90E
120E 150E 180W 150W 120W
90W
60W
30W
0
The detailed mechanism through which the −2 −1.5 −1 −0.5 0 0.5 1 1.5 2 Atlantic SST biases Figure 15. Annual mean SST ( C) response in the PC_TNAxTPAC run. affect the North Pacific is currently under investigation. A preliminary analysis of CMIP5 models indicates that the Arctic oscillation or annular mode plays an important role in transmitting the Atlantic bias to the North Pacific Ocean. The strong Atlantic cooling reinforces the Aleutian low over the North Pacific by triggering an annular mode. The surface wind in the North Pacific is anomalously cyclonic, which is in favor of generating a southward Ekman cold advection. In addition, the subpolar gyre is significantly strengthened in response to the cyclonic wind, which leads to a strengthening of the Oyashio and thus favors a cooling anomaly in the North Pacific midlatitudes. The CESM1.0.4 simulations further demonstrate that the North Pacific cooling bias is largely due to the remote effect of the cooling bias in the North Atlantic Ocean, while the remote impact of the TNA cooling bias is of secondary importance. Similar to the CMIP5 model analyses, the North Atlantic cooling bias triggers the circumglobal teleconnection (Arctic Oscillation), producing a strengthening of the Aleutian low, an enhancement of the southward Ekman and Oyashio cold advection, and thus a cooling SST in the North Pacific Ocean. In contrast, the North Pacific cooling anomaly triggered by the TNA bias is mainly through the tropical Pacific region. The TNA cooling bias suppresses local convection and rainfall, producing a surface high extending to the eastern tropical North Pacific. Over the eastern tropical Pacific, this anomalous high induces northeasterly wind anomalies north, northerly cross-equatorial winds, and northwesterly winds south of the equator. This C-shape wind anomaly generates a warming in the SEP Ocean, which eventually leads to a warming in the central and western Pacific by the positive WES feedbacks. The tropical Pacific warming favors a cooling of the North Pacific Ocean by triggering the PNA teleconnection pattern. While emphasizing the role of ocean heat transport in the Atlantic sector, we do not wish to down play the potential importance of local feedbacks in causing the large systematic error in oceans such as the tropical Atlantic. It is very likely that regional feedbacks between stratocumulus clouds, surface winds, upwelling, coastal currents and ocean eddies are very important contributors to the large warm SST bias off the coast of Angola and Namibia and Peru. Future investigations are clearly needed to understand the causes of the large biases in the tropical Atlantic and Pacific Oceans. Our study here has some implications for the classic double intertropical convergence zone (ITCZ) problem in the coupled models. Recent studies suggested that the ITCZ problem is tied to bias in the Southern Hemisphere (SH) extratropics [e.g., Hwang and Frierson, 2013]. They demonstrated that the warm SST bias in the Southern Ocean as a result of unreasonable cloud parameterization and radiative forcing, leads to a warmer temperature anomaly in the SH than in the Northern Hemisphere. This interhemispheric temperature contrast favors an anomalous Hadley cell with a deep convection center in the warmer SH due to energy balance and thus produces an excess of precipitation in the SH in climate models. In the current paper, we found that the North Atlantic and North Pacific exhibit significant cooling SST bias in the coupled model, which is suggested to be related to the weak simulated AMOC. This cooling SST bias in the Northern Hemisphere extratropics could also induce a temperature contrast between two hemispheres, which is then expected to drive a southward shift of the Hadley circulation and thus an excessive of precipitation in the SH. Thereby, the relative roles of Southern Ocean warm bias and Northern Hemisphere cooling bias on the classic double ITCZ problem need to be addressed in the future.
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
756
Journal of Advances in Modeling Earth Systems Acknowledgments The CMIP5 data for this paper are available through the Earth System Grid — Center for Enabling Technologies (ESG-CET), on the page http://pcmdi9.llnl.gov/. The ERSST data are available at NOAA’s Earth System Research Laboratory, on the page http://www.esrl.noaa.gov/psd/data/ gridded/. Dataset: NOAA Extended Reconstructed SST. Dataset name: sea surface temperature. The model output used in this paper is available upon request (email:
[email protected]). Acknowledgement is made to the Chinese science foundation (NSFC, 41106009).
ZHANG AND ZHAO
10.1002/2014MS000415
References Alexander, M. A., I. Blad e, M. Newman, J. R. Lanzante, N. Lau, and J. D. Scott (2002), The atmospheric bridge: The influence of ENSO teleconnections on air–sea interaction over the global oceans, J. Clim., 15, 2205–2231. Behl, R., and J. P. Kennett (1996), Brief interstadial events in the Santa Barbara basin, NE Pacific, during the past 60 kyr, Nature, 379, 243–246. Delworth, T. L., and M. E. Mann (2000), Observed and simulated multidecadal variability in the Northern Hemisphere, Clim. Dyn., 16, 661–676. Dong, B., R. T. Sutton, and A. A. Scaife (2006), Multidecadal modulation of El Ni~ no-Southern Oscillation (ENSO) variance by Atlantic Ocean sea surface temperatures, Geophys. Res. Lett., 33, L08705, doi:10.1029/2006GL025766. Enfield, D. B., A. M. Mestas-Nu~ nez, and P. J. Trimble (2001), The Atlantic multidecadal oscillation and its relation to rainfall and river flows in the continental U.S., Geophys. Res. Lett., 28, 2077–2080. Gill, A. E. (1980), Some simple solutions for heat-induced tropical circulation, Q. J. R. Meteorol. Soc., 106, 447–462. Goldenberg, S. B., C. W. Landsea, A. M. Mestas-Nu~ nez, and W. M. Gray (2001), The recent increase in Atlantic hurricane activity: Causes and implications, Science, 293, 474–479. Ham, Y. G., J. S. Kug, J. Y. Park, and F. F. Jin (2013), Sea surface temperature in the north tropical Atlantic as a trigger for El Ni~ no/Southern Oscillation events, Nat. Geosci., 6, 112–116. Heckley, W. A., and A. E. Gill (1984), Some simple analytic solutions to the problems of forced equatorial long waves, Q. J. R. Meteorol. Soc., 110, 203–217. Hessl, A., D. McKenzie, and R. Schellhaas (2004), Drought and Pacific decadal oscillation linked to fire occurrence in the inland Pacific Northwest, Ecol. Appl., 14, 425–442. Horel, J. D., and J. M. Wallace (1981), Planetary-scale atmospheric Phenomena associated with the southern oscillation, Mon. Wea. Rev., 109, 813–829. Hwang, Y.-T., and D. M. M. Frierson (2013), Link between the double-intertropical convergence zone problem and cloud biases over the Southern Ocean, Proc. Natl. Acad. Sci. U. S. A., 110, 4935–4940. Kang, S. M., I. M. Held, and S. P. Xie (2014), Contrasting the tropical responses to zonally asymmetric extratropical and tropical thermal forcing, Clim. Dyn., 42, 2033–2043. Knutti, R., J. Fluckiger, T. F. Stocker, and A. Timmermann (2004), Strong hemispheric coupling of glacial climate through freshwater discharge and ocean circulation, Nature, 430, 851–856. Mantua, N. J., and S. R. Hare (2002), The Pacific decadal oscillation, J. Oceanogr., 58, 35–44. Mantua, N. J., S. R. Hare, Y. Zhang, J. M. Wallace, and R. C. Francis (1997), A Pacific interdecadal climate oscillation with impacts on salmon production, Bull. Am. Meteorol. Soc., 78, 1069–1079. McCabe, G. J., M. A. Palecki, and J. L. Betancourt (2004), Pacific and Atlantic Ocean influences. on multidecadal drought frequency in the United States, Proc. Natl. Acad. Sci. U. S. A., 101, 4136–4141. Meehl, G. A., W. M. Washington, J. M. Arblaster, A. Hu, H. Teng, C. Tebaldi, W. G. Strand, and J. B. White III (2012), Climate system response to external forcings and climate change projections in CCSM4, J. Clim., 25, 3661–3683. Okumura, Y. M., C. Deser, A. Hu, A. Timmermann, and S.- P. Xie (2009), North Pacific climate response to freshwater forcing in the subarctic North Atlantic: oceanic and atmospheric pathways, J. Clim., 22, 1424–1445. Smith, T. M., and R. W. Reynolds (2004), Improved extended reconstruction of SST (1854–1997), J. Clim., 17, 2466–2477. Stouffer, R. J., and Coauthors (2006), Investigating the causes of the response of the thermohaline circulation to past and future climate changes, J. Clim., 19, 1365–1387. Taylor, K. E., R. J. Stouffer, and G. A. Meehl (2012), An overview of CMIP5 and the experiment design, Bull. Am. Meteorol. Soc., 93, 485–498. Thompson, D. W. J., and J. M. Wallace (1998), The Arctic Oscillation signature in the wintertime geopotential height and temperature fields, Geophys. Res. Lett., 25, 1297–1300. Timmermann, A., S. I. An, U. Krebs, and H. Goosse (2005), ENSO suppression due to weakening of the North Atlantic thermohaline circulation, J. Clim., 18, 3122–3139. Timmermann, A., et al. (2007), The influence of a weakening of the Atlantic meridional overturning circulation on ENSO, J. Clim., 20, 4899–4919. Vimont, D. J., M. Alexander, and A. Fontaine (2009), Midlatitude excitation of tropical variability in the Pacific: the role of thermodynamic coupling and seasonality*, J. Clim., 22, 518–534. Wang, C., S.-K. Lee, and C. R. Mechoso (2010), Interhemispheric influence of the Atlantic warm pool on the Southeastern Pacific, J. Clim., 23, 404–418. Wang, C., L. Zhang, and S.-K. Lee (2013), Response of freshwater flux and sea surface salinity to variability of the Atlantic warm pool, J. Clim., 26, 1249–1267. Wang, C., L. Zhang, S.-K. Lee, L. Wu, and C. R. Mechoso (2014), A global perspective on CMIP5 climate model biases, Nat. Clim. Change, 4, 2201–2205. Wu, L., F. He, and Z. Y. Liu (2005), Coupled ocean–atmosphere response to north tropical Atlantic SST: Tropical Atlantic dipole and ENSO, Geophys. Res. Lett., 32, L21712, doi:10.1029/2005GL024222. Wu, L., F. He, Z. Y. Liu, and C. Li (2007), Atmospheric teleconnections of tropical Atlantic variability: Interhemispheric, tropical–extratropical, and cross-basin interactions, J. Clim., 20, 856–870. Wu, L., C. Li, C. Yang, and S.-P. Xie (2008), Global teleconnections in response to a shutdown of the Atlantic meridional overturning circulation, J. Clim., 21, 3002–3019. Xie, S.-P. (1996), Westward propagation of latitudinal asymmetry in a coupled ocean–atmosphere model, J. Atmos. Sci., 53, 3236–3250. Xie, S.-P., and J. A. Carton (2004), Tropical Atlantic variability: Patterns, mechanism, and impacts. Earth climate: The ocean–atmosphere interaction, Geophys. Monogr., 147, 121–142. Xie, S.-P., Y. Okumura, T. Miyama, and A. Timmermann (2008), Influences of Atlantic climate change on the Tropical Pacific via the Central American Isthmus, J. Clim., 21, 3914–3928. Yang, J. Y. (1999), A linkage between decadal climate variations in the Labrador Sea and the tropical Atlantic Ocean, Geophys. Res. Lett., 26, 1023–1026. Zhang, L., and C. Wang (2012), Remote influences on freshwater flux variability in the Atlantic warm pool region, Geophys. Res. Lett., 39, L19714, doi:10.1029/2012GL053530. Zhang, L., and L. Wu (2012), Can oceanic freshwater Flux Amplify Global Warming?, J. Clim., 25, 3417–3430.
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
757
Journal of Advances in Modeling Earth Systems
10.1002/2014MS000415
Zhang, L., L. Wu, and J. Zhang (2011a), Coupled ocean–atmosphere responses to recent freshwater flux changes over the KuroshioOyashio extension region, J. Clim., 24, 1507–1524. Zhang, L., L. Wu, and J. Zhang (2011b), Simulated response to recent freshwater flux change over the Gulf Stream and its extension: coupled ocean–atmosphere adjustment and Atlantic Pacific teleconnection, J. Clim., 24, 3971–3988. Zhang, L., C. Wang, Z. Song and S.-K. Lee (2014), Remote effect of the model cold bias in the tropical North Atlantic on the warm bias in the tropical southeastern Pacific, J. Adv. Model. Earth Syst., 6, 1016–1026. Zhang, R. (2008), Coherent surface–subsurface fingerprint of the Atlantic meridional overturning circulation, Geophys. Res. Lett., 35, L20705, doi:10.1029/2008GL035463. Zhang, R., and T. L. Delworth (2005), Simulated tropical response to a substantial weakening of the Atlantic thermohaline circulation. J. Clim., 18, 1853–1860. Zhang, R., and T. L. Delworth (2007), Impact of the Atlantic multidecadal oscillation on North Pacific climate variability, Geophys. Res. Lett., 34, L23708, doi:10.1029/2007GL031601.
ZHANG AND ZHAO
PROCESSES AND MECHANISMS FOR THE MODEL SST BIASES
758