AMERICAN METEOROLOGICAL SOCIETY - AMS Journals

7 downloads 1115 Views 3MB Size Report
EARLY ONLINE RELEASE ..... http://www.cesm.ucar.edu/models/cesm1.0/cam/docs/ug5_0/ug.html ). 213. 3. ...... http://www.esrl.noaa.gov/psd/data/gridded/.
15 DECEMBER 2016

FENG ET AL.

8949

Contrasting Responses of the Hadley Circulation to Equatorially Asymmetric and Symmetric Meridional Sea Surface Temperature Structures JUAN FENG AND JIANPING LI College of Global Change and Earth System Science, Beijing Normal University, and Joint Center for Global Change Studies, Beijing, China

FEIFEI JIN School of Ocean and Earth Science and Technology, University of Hawai‘i at Manoa, Honolulu, Hawaii

ZHENGYU LIU Department of Atmospheric and Oceanic Sciences, and Nelson Institute Center for Climatic Research, University of Wisconsin–Madison, Madison, Wisconsin, and Laboratory for Climate and Ocean–Atmosphere Studies, Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing, China

XING NAN Beijing Meteorological Service, Beijing, China

YIPENG GUO School of Atmospheric Sciences, Nanjing University, Nanjing, China (Manuscript received 24 February 2016, in final form 22 August 2016) ABSTRACT The impacts of different meridional structures of tropical sea surface temperature (SST) on the Hadley circulation (HC) in the annual mean are investigated during the period 1948–2013. By decomposing the variations in SST and the HC into two components—that is, the equatorially asymmetric (SEA for SST, and HEA for HC) and the equatorially symmetric (SES for SST, and HES for HC) parts—it is shown that the long-term variability in SEA and SES captures well the temporal variations in equatorially asymmetric and symmetric variations in SST. The variation in HEA is closely linked to that of SEA, and the variation in HES is connected with that of SES. However, the response of HEA to a given amplitude variation in SEA is stronger (by ;5 times) than that of HES to the same amplitude variation in SES. This point is further verified by theoretical and numerical models, indicating that the meridional structure of SST plays a crucial role in determining the anomalies in HC. This result may explain why the principal mode of HC is dominated by an equatorially asymmetric mode in its long-term variability.

1. Introduction Denotes Open Access content.

Corresponding author address: Prof. Jianping Li, College of Global Change and Earth System Science, Beijing Normal University, 19 XinJieKouWai St., HaiDian District, Beijing 100875, China. E-mail: [email protected] DOI: 10.1175/JCLI-D-16-0171.1 Ó 2016 American Meteorological Society

The tropical Hadley circulation (HC) is defined as the zonal mean meridional mass circulation in the atmosphere bounded roughly by 308S and 308N. It is one of the most important and largest atmospheric circulations. The HC is mainly a thermally driven circulation, with warmer air rising in the tropics because of the release of

8950

JOURNAL OF CLIMATE

latent heat and colder air sinking in the subtropics in both hemispheres, generating a closed circulation in each hemisphere (Held and Hou 1980). This circulation is characterized by equatorward mass transport by the prevailing trade wind flow in the lower troposphere and poleward mass transport in the upper troposphere. The HC plays an essential role in the global climate system. It transports not only heat from the tropics to the subtropics and to high latitudes through extratropical eddies, but also momentum to the subtropics (Lindzen 1994; Hou 1998). Both heat and momentum transports have significant influences on the subtropical jet streams. Therefore, the HC has an important impact on global climate (Diaz and Bradley 2004; Feng et al. 2013). In recent years there has been growing interest in the changes in the HC. Studies have focused on two key problems: one is the intensity of the HC, and the other is its meridional extent (Hu and Zhou 2009). Wielicki et al. (2002) and Chen et al. (2002) suggested that the annual mean HC intensified in the 1990s. Subsequently, the variations in the seasonal mean HC have been discussed in many studies. The variations in the HC during boreal winter [December–February (DJF)] have been extensively explored, and an obvious strengthening trend has been observed. For example, Quan et al. (2004) reported that the HC has strengthened since the 1950s. Similar results have been obtained in other studies for different datasets or periods (Hu et al. 2005; Mitas and Clement 2005; Ma and Li 2007, 2008). A similar intensifying trend is seen in boreal spring [March–May (MAM); Feng et al. (2013)]. In boreal summer [June– August (JJA)], the HC shows minor changes (e.g., Quan et al. 2004; Tanaka et al. 2004; Mitas and Clement 2006; Feng et al. 2011). With respect to the meridional extent of the HC, results from several recent studies (e.g., Fu et al. 2006; Hu and Fu 2007; Seidel and Randel 2007) suggest that the HC has expanded poleward since 1979, by 1.218–58 latitude in different seasons and in different datasets (Hudson et al. 2006; Lu et al. 2007; Johanson and Fu 2009). In addition to the strength and extent of the HC, its long-term variability has been investigated, as well as its relationship with sea surface temperature (SST) (e.g., Feng et al. 2011; Feng and Li 2013; Guo et al. 2016). In particular, the long-term variability in HC has been discussed using empirical orthogonal function (EOF) analysis to explore its spatiotemporal characteristics. For example, focusing on the seasonal cycle, Dima and Wallace (2003) found that the HC is dominated by two components with comparable mean-square amplitudes, one equatorially symmetric and the other equatorially asymmetric. The two modes vary sinusoidally with the seasonal cycle. As for the year-to-year variability in the

VOLUME 29

seasonal mean HC, Ma and Li (2007, 2008) reported that the principal modes of boreal winter HC include an equatorially asymmetric and an equatorially symmetric mode; the two modes are connected to SST variations over the Indo-Pacific warm pool and tropical Pacific, respectively. Similar results have been reported for the long-term variability in the boreal summer HC (Feng et al. 2011; Sun and Zhou 2014). The equatorially asymmetric modes in boreal winter and summer mirror the climatological structure of the HC, as it is equatorially asymmetric in both boreal winter and summer. Further research by Feng et al. (2013) and Guo et al. (2016) showed that the first two principal modes of boreal spring and autumn [September–November (SON)] HC are an equatorially asymmetric and an equatorially symmetric mode. Their results indicate that inhomogeneous warming of SST over the tropical oceans contributes to the formation of the equatorially asymmetric mode (Feng et al. 2013; Li et al. 2013). These studies highlight the role of the SST meridional gradient in the variations of the HC, as seen in observations and numerical experiments; however, they do not explain why the equatorially asymmetric mode dominates the longterm variability of the HC. In fact, theoretical and numerical studies have established that the meridional structure of tropical SST may change the strength and position of the ascending branch and convergence by adjusting the convergence process and thermal structure of the atmosphere (Schneider and Lindzen 1977; Rind and Rossow 1984). Theoretical work by Lindzen and Nigam (1987) indicates that meridionally inhomogeneous variations in SST lead to changes in the meridional gradient of SST, which in turn influence the boundary layer wind and vertical motion in the lower troposphere. Later work has also shown that convergence in the lower troposphere is sensitive to the meridional distribution of the heating profile (Lindzen and Hou 1988; Hou and Lindzen 1992). These studies indicate that the meridional distribution of tropical SST may have a considerable impact on the convergence in the lower troposphere. In previous studies, we also reported that the position of the anomalous ascending branch of the HC is linked to the meridional gradient of tropical SST (Feng et al. 2013; Feng and Li 2013). Thus, it is clear that the meridional structure of SST is a key influence on the spatial structure of the HC. Moreover, it is found that the equatorially asymmetric mode dominates the variability of the HC in both longterm seasonal variations and its seasonal cycle, even in the transition seasons (i.e., boreal spring and autumn), when the HC is equatorially symmetric in the climatology. This suggests that the HC may respond differently

15 DECEMBER 2016

to different SST meridional structures. Therefore, the response of the HC to equatorially asymmetric thermal forcing may be greater than to equatorially symmetric forcing, thus contributing to the formation of the equatorially asymmetric mode of the HC. Understanding this issue would further our understanding of the variability of the HC and help us to assess the responses of the HC to different SST meridional structures. Here we focus on this issue by separating the variations in SST into two components, the equatorially asymmetric and symmetric, and examining their respective influences on the HC. The remainder of this paper is organized as follows. The datasets, methodology, and models are described in section 2. Section 3 outlines the climatological features of the annual mean HC and zonal mean SST, as well as their principal modes. Section 4 demonstrates the responses of the HC to different SST meridional structures in reanalysis data, and section 5 shows the results from the theoretical and numerical models. Section 6 presents a simple theoretical deduction for the observation and model results. Finally, conclusions and a discussion are presented in section 7.

2. Datasets, methodology, and models a. Datasets

1 ›[y] ›[v] 1 5 0. a ›f ›p

(2)

Thus, we have [y] 5

g ›c 2pa cosf ›p

and

g ›c [v] 5 2 , 2pa2 cosf ›f

(3)

where [y] is the zonal mean meridional wind, a is the mean radius of Earth, f is latitude, g is gravitational acceleration, and p is pressure. The overbars and square brackets represent temporal and zonal averaging, respectively. Note that v is an indirect variable and is obtained from zonal and meridional winds. Here, to reduce errors we integrate Eq. (3) involving the meridional wind from the top of the atmosphere, where it is assumed that c 5 0 and p 5 0, yielding ð 2pa cosf p [y(f, p)] dp. (4) c(f, p) 5 g 0 To examine the influence of the meridional structure of SST on the HC, we separate the variations in SST and HC into two components: the equatorially asymmetric and the equatorially symmetric parts. For the variations in SST, the equatorially asymmetric (i.e., SEA) and symmetric (i.e., SES) parts are, respectively, defined as SST(j) 2 SST(2j) and 2 SST(j) 1 SST(2j) . SES(j) 5 2

SEA(j) 5

The atmospheric fields are from the National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR) reanalysis (Kalnay et al. 1996), with 2.58 3 2.58 horizontal resolution. The two global SST reanalysis datasets are from the improved Extended Reconstructed SST (ERSST) on a 28 3 28 grid (Smith and Reynolds 2004), and the Met Office Hadley Centre Sea Ice and SST dataset (HadISST) gridded at 18 3 18 resolution (Rayner et al. 2003). We focus on the period from 1948 to 2013 and consider the annual mean.

b. Methodology The HC is characterized by the mass streamfunction (MSF) c of the mean meridional circulation. Clockwise circulation (the northern cell) is defined as positive and anticlockwise circulation (the southern cell) is defined as negative. Using pressure as the vertical coordinate, the conservation of mass requires 1 ›u 1 ›(y cosf) ›v 1 1 5 0, a cosf ›l a cosf ›f ›p

8951

FENG ET AL.

(5)

For the HC, since the MSF values in the two hemispheres are opposite, the equatorially asymmetric (i.e., HEA) and symmetric (i.e., HES) components are defined as follows: MSF( j) 1 MSF(2j) and 2 MSF( j) 2 MSF(2j) , HES(j) 5 2

HEA( j) 5

(6)

where j and 2j are equatorially symmetric meridional positions. In addition, two SST indices defined by Feng and Li (2013), the meridional asymmetric SST index (MASI) and the meridional symmetric SST index (MSSI), are employed to further examine the meridional structure of tropical SST variations. The MASI and MSSI are defined as follows: MSSI 5 0:5 3 [SSTA(158S258S) 1 SSTA(58N2158N) ] and (7)

(1)

zonally averaged around the entire globe. The first term in Eq. (1) goes to zero, giving

MASI 5 SSTA(158S258S) 2 SSTA(158N258N) .

(8)

The SSTA in the Eqs. (7) and (8) is for the SST anomalies removing the annual cycle, and the subscripts are

8952

JOURNAL OF CLIMATE

VOLUME 29

for the areal mean over the specified region. Note that to enable a direct comparison between the indices above and the decompositions, the sign of the MASI here is reversed compared with Feng and Li (2013). And we further simplify the MSSI to facilitate a direct comparison with the MASI. The index of annual mean HC intensity (HCI) is defined as the maximum of the absolute value of the annual mean MSF between 308S and 308N (Oort and Yienger 1996). EOF analysis is employed to determine the principal mode of year-to-year variability of longterm annual mean MSF and the equatorially symmetric and asymmetric variabilities of the zonal mean SST and MSF. Regression and correlation methods are employed to detect the relationships between SST and the HC.

c. Models Gill (1980) elucidated the responses of the tropical atmosphere to isolated equatorially symmetric heating in a resting atmosphere using linear shallow-water equations. Recently, Xing et al. (2014) explored the response of the atmosphere by adopting Gill’s model with detailed influences, including the location, strength, and extent of heating. Gill’s model is known to be a simple and useful theoretical model for exploring the response of the tropical atmosphere to a given distribution of heating in the tropics (e.g., Chao and Wang 1991; Ratnam et al. 2012). The steady-state version of the one-layer nondimensional equations following Gill (1980) may be written as follows: 1 ›p «u 2 yy 5 2 , 2 ›x 1 ›p «y 1 yu 5 2 , 2 ›y ›u ›y «p 1 1 5 2Q, and ›x ›y w 5 «p 1 Q ,

(9) (10) (11) (12)

where (x, y) is nondimensional distance, with x and y representing eastward and northward distances from the equator, respectively; (u, y) is proportional to horizontal velocity; p is proportional to the pressure perturbation; Q is proportional to the heating rate; w represents vertical velocity; and « is a decay factor. The atmospheric general circulation model employed is the NCAR Community Atmospheric Model, version 5 (CAM5; Collins et al. 2006). This model is able to reproduce observed climate features (e.g., Feng and Li 2013) and serves as the atmospheric component of the Community Climate System Model (CCSM). The horizontal resolution is spectral T42 (approximately 2.88 3 2.88), with 26 hybrid vertical levels (a complete

FIG. 1. (a) Climatological annual mean MSF over the period 1948–2013. The contour interval is 0.015 3 1010 kg s21. (b) Principal mode of the annual mean HC during 1948–2013. Positive (negative) contours are shown as solid (dotted) lines. The contour interval is 0.015 3 1010 kg s21. (c) Time series of the principal mode of the annual mean HC (31010 kg s21).

description of this model version is available online at http://www.cesm.ucar.edu/models/cesm1.0/cam/docs/ug5_0/ ug.html).

3. Climatological features of the Hadley circulation and zonal mean SST a. Hadley circulation Figure 1 displays the climatological annual mean MSF and its first principal mode during the period 1948–2013. The northern and southern components of the annual mean HC have equivalent magnitude and extent, with the rising branch around the equator and descending branches around 308 latitude in each hemisphere (Fig. 1a). The first principal mode is an equatorially asymmetric mode that dominates the variability of annual mean HC, with the ascending branch centered about 108S and the two descending branches located south of 308S and north of 308N. This mode accounts for 63.5% of the total variance. The stronger component of this mode is centered to the north of the equator, with the combined rising branch in the Southern Hemisphere

15 DECEMBER 2016

FENG ET AL.

8953

FIG. 2. (a) Climatological zonal mean tropical SST (8C) over the period 1948–2013. (b) Principal mode of zonal mean tropical SST during 1948–2013. (c) Time series of the principal mode of zonal mean SST. Red and blue lines are based on the ERSST and HadISST datasets, respectively.

and a descending branch in the Northern Hemisphere. In contrast, the counterpart in the Southern Hemisphere is weak, with a descending branch positioned south of 308S. The first principal mode here is similar to those observed in the seasonal long-term variations in the HC (Ma and Li 2008; Feng et al. 2011, 2013; Guo et al. 2016). Moreover, this mode has exhibited a significantly strengthened trend during the past six decades. The time series of this mode underwent a decadal transformation from negative to positive values during the late 1970s, which coincides with the timing of an interdecadal shift in tropical climate (e.g., Mantua et al. 1997; Xiao and Li 2007). The intensification of this mode is similar to the observed trend in the HCI, suggesting that intensification of the annual mean HC may be due to the variations in this mode. The strong correlation between these latter two factors (r 5 0.77) supports this view. The above result indicates that the year-to-year variability in long-term annual mean HC is dominated by an equatorially asymmetric mode that intensified during the period 1948–2013, and this variability contributes to the strengthening of the annual mean HC.

b. Zonal mean SST Figure 2 shows the distribution of the climatological annual zonal mean SST within the tropics and its first principal mode determined from the ERSST and HadISST data during 1948–2013. For the climatological mean, the zonal mean tropical SST shows two peaks

symmetrically positioned about the equator. However, the peak value in the Northern Hemisphere is significantly larger than that in the Southern Hemisphere. This is in agreement with the average position of the annual mean intertropical convergence zone (An et al. 2015), to the north of the equator, as determined by the larger ocean area in the Southern Hemisphere. The first mode explains 88.6% (86.3%) of the total variance of the zonal mean SST variation based on ERSST (HadISST), representing the fundamental variability of zonal mean SST. This pattern shows a stronger value around the equator that decreases with latitude, resembling the proposed enhanced equatorial response (i.e., a stronger warming near the equator than in the subtropical Pacific) reported in Liu et al. (2005) and Xie et al. (2010). Although there is some inconsistency in amplitude between the two SST reanalyses, both show a consistent equatorially asymmetric structure, with larger values in the Southern Hemisphere; that is, the warming rate in the tropical Southern Hemisphere is higher than that in the Northern Hemisphere. The tendency here is different from that presented in Xie et al. (2010), which indicated a tendency for a greater warming in the northern subtropics than in the southern subtropics. Here, we focused on the global mean; however, only the Pacific is focused on in Xie et al. (2010). The principal component (PC) of this mode shows a clear upward trend, indicating significant warming in the tropical oceans; it also implies that the warming difference

8954

JOURNAL OF CLIMATE

VOLUME 29

FIG. 3. (a) Climatological distribution of the annual mean HEA, and (b) its corresponding standard deviation. (c),(d) As in (a),(b), but for the annual mean HES. Positive (negative) contours are shown as solid (dotted) lines. The contour interval in (a),(c) is 0.01 3 1010 kg s2 and 0.003 3 1010 kg s2 in (c),(d).

between the two tropical hemispheres has increased over time. The HC is a thermally driven circulation, and its variations are closely linked to the underlying thermal structure. We therefore examine the possible linkage between the long-term variability of the annual mean HC and zonal mean SST. The correlation coefficient between the PCs of their first principal modes is 0.69, indicating that the tropical ocean warming is closely related to the long-term enhancement of the equatorially asymmetric modes.

4. Responses of the Hadley circulation to different meridional SST structures We decompose the HC into two components, HEA and HES, as described in section 2b. Figure 3 shows their climatological spatial distribution and corresponding standard deviations. The climatological HEA (upper) extends from 308S to 308N, centered at the equator; its largest variability is around the equator. In contrast, the climatological HES contain two mirrored cells (below), with a combined ascending branch located at the equator and two descending branches positioned in the subtropics in each hemisphere. The HES is stronger than the HEA, but they have equivalent extent, while the variability of the HEA is much larger than that of the HES. In addition, the most variable parts of both the HEA and HES are located in the tropics, further indicating that the tropical Hadley circulation is the most prominent meridional circulation. These results show that the HEA and HES are distinct in both intensity and variability: the intensity of HES is larger, while the variability of HEA is larger. This point

is consistent with the equatorial asymmetry of the first principal mode of the HC. Figure 4 displays the meridional climatological distributions of zonal mean SST for the two components and their corresponding standard deviations. Although there is some disagreement between the two SST reanalyses (i.e., ERSST and HadISST), the spatial structures show consistent features. The amplitudes of SEA and SES are equivalent in the tropics. The variability is zero at the equator for the SEA and increases with increasing latitude. In contrast, the maximum variability of the SES is located around the equator and decreases with latitude. We further examined the evolution of long-term variability in the HEA and HES. The principal modes as well as their corresponding time series are shown in Fig. 5. The first principal mode of SEA is similar to its climatological structure, with equatorial asymmetry, and explains 76.3% of the variance. The time series of this mode exhibits strong interdecadal variations and undergoes a decadal transformation from negative to positive values during the late 1970s. That is, the phase of the principal mode of the HEA changed at the same time as an interdecadal shift in tropical climate (e.g., Mantua et al. 1997; Xiao and Li 2007). Meanwhile, the first principal mode of HES is equatorially symmetric, with an explained variance of 63.5%, and the time series of this mode displays a long-term upward trend. The correlation coefficients of the PCs of the two modes and the climatology mean HC (HEA and HC, HES and HC) are 0.98 and 0.83, respectively. However, the correlation drops to 20.08 for the HES and HC when the linear trend is removed, while remaining at 0.98 for the HEA and HC, indicating that HES mainly represents

15 DECEMBER 2016

FENG ET AL.

8955

FIG. 4. (a) Climatological distribution of the SEA, and (b) its corresponding standard deviation (8C). (c),(d) As in (a),(b), but for the SES. Red and blue lines are based on the ERSST and HadISST datasets, respectively.

the long-term linear trend of the HC, while HEA refers to the interdecadal variations of the HC. As for the zonal mean SST, the first principal mode of SEA is equatorially asymmetric (Fig. 6a), explaining 89.3% (88.2%) of the variance based on the ERSST (HadISST). The primary mode of SES explains 93.3% (92.1%) of the variance based on the ERSST (HadISST; Fig. 6b). The time series of the SEA primary mode shows strong interdecadal variability with no obvious trend (Fig. 6b), while the time series of the SES principal mode has an obvious intensifying trend (Fig. 6d). The correlation coefficients of the corresponding PCs are

0.71 (0.63) and 0.77 (0.69) for those of the HEA and SEA and HES and SES based on the ERSST (HadISST), respectively. The high correlations indicate that the variations in HEA are closely linked with the SEA, and the variations in HES are connected with the SES. The time series of the SEA’s principal mode is also highly correlated with the MASI (correlation coefficient of 0.96; Fig. 7a), supporting the view that the MASI mainly represents the equatorially asymmetric variations in tropical SST. On the other hand, the time series of SES’s principal mode is associated with the MSSI (Fig. 7b), with a correlation coefficient of 0.99. These

FIG. 5. As in Fig. 1b, but for the (a) principal mode of the HEA during 1948–2013, and (b) its time series. (c),(d) As in (a),(b), but for the HES.

8956

JOURNAL OF CLIMATE

VOLUME 29

FIG. 6. As in Fig. 2b, but for the (a) principal mode of the SEA during 1948–2013, and (b) its time series. (c),(d) As in (a),(b), but for the SES. Red and blue lines are based on the ERSST and HadISST datasets, respectively.

results indicate that the MSSI and MASI perform well in capturing the equatorially symmetric and asymmetric variations in tropical SST. What then are the corresponding vertical circulations associated with the variations in tropical SST? Fig. 8 shows the regression patterns of MASI and MSSI onto the meridional circulation. It is evident that the equatorially symmetric and asymmetric meridional circulation anomalies are connected with the MSSI and MASI, respectively. The equatorially asymmetric SST variation (i.e., MASI) is accompanied by an equatorially asymmetric meridional circulation (Fig. 8a), with the anomalous ascending branch to the south of the equator (Fig. 8c). Meanwhile, the symmetric SST variation (i.e., MSSI) is associated with the anomalous equatorially symmetric meridional circulation (Fig. 8b), with the anomalous ascending branch around the equator (Fig. 8d). On the other hand, although the variations of MSSI are generally associated with equatorial symmetric meridional circulation anomalies (Figs. 8b,d,f), the extent and intensity of the anomalous circulation accompanied with the MSSI in the two hemispheres show certain differences. The possible reasons are as follows: 1) The decompose method presented in this study is linear: that is, the MSSI and MASI are not orthogonal. This point implies that there may be some nonlinear process involving the relationship between the SST and HC. 2) Although the MSSI represents the equatorially symmetric variations of tropical SST, the factors influencing the MSSI may not be equatorially symmetric, which also may play a role in impacting the distribution of meridional circulation. Nevertheless, the result here further verifies that the equatorially asymmetric SST variations contribute to variations in the equatorially asymmetric HC,

while the equatorially symmetric SST variations accompany the equatorially symmetric HC. What, therefore, are the possible differences in meridional circulation resulting from the differences in the meridional structure of the underlying SST? To explore this issue, Fig. 9 shows scatterplots based on time series of the principal mode of HEA against SEA, as well as HES against SES. Besides the high linear correlations between the SEA and HEA, and between the SES and HES, we see that the responses of HEA to SEA and HES to SES are distinct. With ERSST data, the responses of HEA and HES to SST variations of equal amplitude in SEA and SES are significantly different, with coefficients of 24.84 and 4.30, respectively, a factor of roughly 5.8

FIG. 7. (a) Standardized time series of the principal mode of SEA (red) and the MASI (blue). (b) As in (a), but for the SES (red) and the MSSI (blue).

15 DECEMBER 2016

FENG ET AL.

8957

FIG. 8. (left) Regression pattern of (a) MSF, (c) zonal mean vertical velocity, and (e) zonal mean meridional wind with respect to the MASI. (right) As in (left), but with respect to the MSSI. Shading indicates statistical significance at the 0.05 level.

between the response of the HEA to SEA and that of the HES to SES. Similar results are observed for the HadISST data (a factor of about 5.4). Although the exact values are slightly different, it is a consistent result that the response of HEA to SEA is ;5 times larger than that of HES to SES, even though the thermal forcings have the same magnitude. This may explain why the first principal mode of the HC is equatorially asymmetric: the response of the equatorially asymmetric circulation to underlying SST forces is much larger than that of the equatorially symmetric circulation, resulting in the equatorially asymmetric circulation dominating the variability of the HC.

5. Theoretical and numerical modeling Results from the observational datasets indicate that the meridional structure of the underlying thermal forcing can result in very different induced responses of the HC. To further examine the influences of the meridional structure of the underlying SST forcing on the vertical circulation, we performed theoretical and numerical experiments.

a. Theoretical model Gill’s model is used to investigate the responses of meridional circulation to the SST forcing. Considering

8958

JOURNAL OF CLIMATE

VOLUME 29

FIG. 9. (a) Scatterplot of the first PC of the SEA against the first PC of the HEA using the ERSST. The red dots are the PCs of HEA and SEA (SES), the black dots are the linear fit of the scatter for which the regression coefficient is shown. (b) As in (a), but for first PC of the HES. (c),(d) As in (a),(b), but based on the HadISST.

that the heating forcing source in the natural atmosphere often decays in strength from the center outwards and varies zonally (Jiang and Li 2011), the ideal heating source is calculated following Li et al. (2016): Q(x, y) 5 Q1 1 Q2 5 A1 g1 (x)e2(1/R1 )(y1d1 )

2

2

1 A2 g2 (x)e2(1/R2 )(y1d2 ) and ( cosk1 (x 1 l1 ) jxj # L1 and g1 (x) 5 0 jxj . L1 ( cosk2 (x 1 l2 ) jxj # L2 and g2 (x) 5 0 jxj . L2

k1 5 p/2L1 k2 5 p/2L2

.

Here, A1 and A2 represent bipolar heating intensities; R1 and R2 represent the meridional gradients of heating; d1 and d2 are meridional positions of the bipolar heating centers; d1 . 0 and d2 . 0 and d1 , 0 and d2 , 0 represent heating component centers located in the Southern and Northern Hemispheres, respectively; l1 and l2 are the zonal distances of the bipolar heating centers from the origin 0 (the intersection of the x and y axes); and 2L1 and 2L2 represent the zonal and meridional widths of heating

sources. For equatorially symmetric heating, A1 5 3.15, d1 5 0, R1 5 0.5, l1 5 0, 2L1 5 36, and A2 5 d2 5 R2 5 l2 5 2L2 5 0. For equatorially asymmetric heating, A1 5 3.15, A2 5 23.15, d1 5 21.5, d2 5 1.5, R1 5 R2 5 0.5, l1 5 l2 5 0, and 2L1 5 2L2 5 36. Two types of SST anomalies are constructed: one equatorially asymmetric and one equatorially symmetric (Fig. 10). Given that the intensity and meridional gradient of SST have important influences on the HC (e.g., Lindzen and Hou 1988; Chiang et al. 2001), the intensity of the underlying SST forcing and the meridional gradient centered on the equator are fixed so that the anomalies in meridional circulation are not contaminated by these parameters. The circulation anomalies are therefore mainly due to the differences in the meridional structure of the anomalous SST. The results show that equatorially asymmetric meridional circulation is associated with the equatorially asymmetric SST forcing, with anomalous ascent in the Southern Hemisphere and the circulation extending from 208S to 208N. In contrast, an equatorially symmetric meridional circulation is observed when the SST forcing is equatorially symmetric. The extent of this anomalous meridional circulation is smaller than in the equatorially asymmetric

15 DECEMBER 2016

FENG ET AL.

8959

FIG. 10. (a) Distribution of the equatorially asymmetric heating in the theoretical experiment (8C). (b) As in (a), but for the equatorially symmetric heating. (c) Corresponding MSF solutions for the equatorially asymmetric heating (0.01 3 1010 kg s2). (d) As in (c), but for the corresponding MSF solutions for the equatorially symmetric heating.

case. The intensities of the two anomalous circulations are very different (with center values of 20 vs 4): the response of the meridional circulation to asymmetric SST forcing is about 5 times the response to the symmetric SST forcing. Note that the response difference shown here is insensitive to the shape of heating; the results based on the rectangular shape of heating show a similar result (figures not shown). This point indicates the responses from the equatorially asymmetric and symmetric structure are steady. The result here agrees well with the observations, further establishing the sensitivity of the response of the meridional vertical circulation to the meridional structure of the thermal forcing.

b. Numerical model The CAM5 is employed to further examine the response of meridional circulation to the thermal meridional structure of SST. We use the aquaplanet version of CAM5. This is sufficient for our study, since the land area within the tropics (i.e., 208S–208N in this study) is small compared with the ocean. A control run was integrated for 8 yr and was used to derive a reference state. Two sets of sensitivity experiments were designed and integrated for 10 yr, with the last 8 yr of the integration used to construct an eight-member ensemble mean to reduce the uncertainties arising from varying initial conditions. To isolate and reproduce the impact of the

SST meridional gradient and intensity on the associated circulation, the only difference between the two sets of sensitivity experiments is the SST meridional structure in the tropical oceans (see Fig. 11), as in the theoretical experiment. Note that the SST anomalies are identical to those used in the theoretical experiment. The prescribed profiles in the experiments are equatorially asymmetric and symmetric (Figs. 11a,c). Figures 11b,d shows the ensemble-mean response of the MSF anomalies associated with the two sets of experiments. When forced with the equatorially asymmetric SST anomalies (see Fig. 11a), the observed anomalous equatorially asymmetric meridional circulation is well captured by the simulation (Fig. 11b). The convergence and rising branch lie to the south of the equator, with the descending branch in the Northern Hemisphere. With the equatorially symmetric SST forcing, we see anomalous convergence and a rising branch at the equator, with descent at 208 poleward in both hemispheres (Figs. 11b, d). These characteristics, in response to the different SST meridional structures, are similar to those seen in the observations. We also find that the anomalous meridional circulation induced by the asymmetric SST forcing is about 5 times stronger than that induced by the symmetric SST forcing (with center values of 20 vs 4). The consistency between the observations and model results further verifies the crucial impact of the meridional structure of SST on the meridional circulation, with the

8960

JOURNAL OF CLIMATE

VOLUME 29

FIG. 11. (left) Specified anomalous SST pattern used in the numerical simulations for (a) the equatorially asymmetric heating and (c) the equatorially symmetric heating (8C). (right) (b) Response of MSF under the forcing of (a) (0.01 3 1010 kg s2); (d) as in (b), but for the forcing of (c).

The fundamental mechanism of the stronger response of HEA to SEA can be further understood in the Gill model, which can be used to investigate the importance of rotation and the b plane. In the simplest form without rotation, the Gill model described in section 2 can be reduced to the following form: ›p and «y 5 2 ›y ›y «p 1 5 2Q . ›y

direction. We could get the solution p 5 2(«Ql /l 2 1 «2 )eily ; thus, the response function is p/Q 5 «/(l 2 1 «2 ). Namely, the response amplitude depends on the spatial scale and decreases with the spatial scale. Neglecting the l 5 0 response, which has a small projection into the symmetric mode of forcing, an asymmetric mode has a small wavenumber l 5 1, and the symmetric mode has l 5 2, so the ratio is about 4 if « is small, which is close to the observations and model results. On the other hand, the effect of the b plane is further examined based on the Gill model. We further calculated the extreme value of zonally integrated meridional wind for the two kinds of modes [Eqs. (4.15) and (5.7) in Gill (1980)]. For the equatorially symmetric heating,

That is (›2 p/›y2 ) 2 «2 p 5 2«Q. Given a waveform forcing Q 5 Ql eily , where l represents the wavenumber in y

y 5 4«q2 (x)ye2y

influence of the equatorially asymmetric SST a factor of 5 larger than that of the equatorially symmetric forcing.

6. Theoretical explanations

2 /4

1 yQ and

Q 5 coskx 3 e2y

2 /4

,

15 DECEMBER 2016

FENG ET AL.

Ðx 2 we get that 2x y dx 5 (2y/3)(4L/p)e2y /4 . For the equatorially asymmetric, y 5 yQ 1 6«q3 (x)(y2 2 1)e2y 2y2 /4

Q 5 coskx 3 ye

2 /4

and

,

Ðx

that is 2x y dx 5 [(6 2 y2 )/5](4L/p)e2y /4 . The definitions of q2 and q3 are the same as those described in Gill (1980). It is found that the extreme value for the equatorially asymmetric mode (when) is ;4 times that of the equatorially symmetric mode (when y 5 1.41). The quantity of the above deduction without rotation and considering the b plane indicates that the rotation may not be the main factor affecting the intensity of atmospheric response, but the effects of the b plane play roles in determining the different responses of HC to SST. 2

7. Discussion and conclusions Using observational datasets for the past 66 yr, we investigated the influence of the meridional structure of SST on the HC, focusing on the annual mean response. By decomposing the long-term variations in SST and HC into two components, one equatorially asymmetric and one equatorially symmetric, we found that the equatorially asymmetric components of SST and HC (i.e., HEA and SEA) mainly reflect the interdecadal variations; in contrast, the equatorially symmetric parts (i.e., HES and SES) mainly represent the long-term linear trends. These results indicate that the equatorially asymmetric variations in HC are linked to the equatorially asymmetric variations in SST, while their equatorially symmetric variations are closely correlated. However, the responses of the HEA and HES to SST forcing are very different: the response of the HEA to SEA is much bigger than that of the HES to SES. This result is further confirmed by both theoretical and numerical experiments. It has been reported that the principal mode of long-term variability of the HC is equatorially asymmetric, not only for boreal winter and summer (Ma and Li 2008; Feng et al. 2011), when the climatological structure of HC is equatorially asymmetric, but also for the boreal spring and autumn (Feng et al. 2013; Guo et al. 2016), when the climatological structure of the HC is equatorially symmetric. The present results provide a plausible explanation for this phenomenon, as disturbances of the same magnitude in SEA and SES are associated with different responses in the HC, and the response to SEA is about 5 times larger than that to SES, thereby contributing to the formation of the equatorially asymmetric mode. In addition, previous studies have reported that the meridional distribution of heating plays an important role in inducing atmospheric anomalies. For example,

8961

Lindzen (1994) found that a larger meridional gradient of SST induces stronger atmospheric anomalies. In previous studies we have shown that the meridional gradient profile of tropical SST determines the location of anomalous atmospheric convergence (Feng et al. 2013; Feng and Li 2013). These studies gave a qualitative description of the importance of the meridional structure of SST for atmospheric circulation that was useful for better understanding the role of SST in impacting the atmosphere. In the present work, by decomposing the variations of SST into two components, we have quantitatively assessed the important role of SST meridional structure on the HC. Moreover, the effects of rotation and the b plane in impacting the contrast responses of HC to different SST meridional structures is further explored based on the Gill model. The theoretical deduction indicates that the cases without rotation and including the b plane both show equivalent quantity for different heating sources, suggesting the b plane plays a role in determining the contrast responses of HC to different SST meridional structures, whereas the role of rotation is limited. The results are helpful for understanding the important impact of the meridional structure of SST on atmospheric anomalies. On the other hand, we found that the HC responds differently to different meridional structure in SST, and it is worth exploring further the reasons for the different responses. Are they due to internal ocean processes or to external processes involving air–sea interaction? In this study, we focused on the long-term variability in HC, but the HC is reported to be dominated by an equatorially asymmetric mode in the seasonal cycle (Dima and Wallace 2003). Thus, it would be interesting to further explore the responses of the HC to different SST meridional structures in the seasonal cycle, to discover whether a similar contrasting sensitivity is observed. Since the HC behaves differently during different seasons, such as greater interdecadal variability during boreal summer (Feng et al. 2011) and a close link to ENSO during boreal winter (Ma and Li 2007), it would also be interesting to further examine and compare the similarities and differences in the response of the HC to SST meridional structures in different seasons. Meanwhile, since the atmosphere–ocean system exhibits strong decadal variability in the late 1970s, it would be interesting to examine whether the response of the HC to different meridional structures of SST may also have some decadal variation. Do the background climate conditions play a role in determining the response of the atmosphere to SST? These questions remain unresolved and require further work. Finally, our results highlight the important role of the meridional structure of SST in determining the

8962

JOURNAL OF CLIMATE

corresponding atmospheric circulations, and we conclude that the atmospheric response to a meridional distribution of SST can vary by more than a factor of 5. Therefore, the differences in the response of the HC to the meridional structure of SST may provide a useful criterion for assessing the output of coupled ocean–atmosphere general circulation models (CGCMs), while the sensitivity of the responses to climate change forcing can be assessed using CGCMs. Further work will investigate the application of the responses of the HC to different SST forcing using an ensemble of CGCM results. This appears to provide a useful ‘‘fingerprinting’’ technique for this specific aspect of climate change in the tropics. Acknowledgments. This work was jointly supported by the National Natural Science Foundation of China (41475076), the 973 Program (2016YFA0601801), and the SOA Program on Global Change and Air–Sea interactions (GASI-IPOVAI-03). Prof. Zhengyu Liu is supported by the NSF. The NCEP–NCAR reanalysis and ERSST were obtained from the NOAA (available online at http://www.esrl.noaa.gov/psd/data/gridded/). The HadISST dataset was obtained from the Met Office Hadley Centre (available online at http://www.metoffice. gov.uk/hadobs/hadisst/data/download.html). REFERENCES An, Z. S., and Coauthors, 2015: Global monsoon dynamics and climate change. Annu. Rev. Earth Planet. Sci., 43, 29–77, doi:10.1146/annurev-earth-060313-054623. Chao, J. P., and Z. G. Wang, 1991: Influence of the ocean heating scale on the structure of the vertical circulation cell in the tropical atmosphere. J. Nanjing Inst. Meteor., 14, 10–17. Chen, J. Y., B. E. Carlson, and A. D. Genio, 2002: Evidence for strengthening of the tropical general circulation in the 1990s. Science, 295, 838–841, doi:10.1126/science.1065835. Chiang, J. C. H., S. E. Zebiak, and M. A. Cane, 2001: On the relative roles of elevated heating and surface temperature gradients in driving anomalous surface winds over tropical oceans. J. Atmos. Sci., 58, 1371–1394, doi:10.1175/1520-0469(2001)058,1371: RROEHA.2.0.CO;2. Collins, W. D., and Coauthors, 2006: The formulation and atmospheric simulation of the Community Atmosphere Model version 3 (CAM3). J. Climate, 19, 2144–2161, doi:10.1175/JCLI3760.1. Diaz, H. F., and B. Bradley, 2004: The Hadley Circulation: Present, Past and Future. Advances in Global Change Research, Vol. 21, Kluwer Academic Publishers, 511 pp. Dima, I. M., and J. M. Wallace, 2003: On the seasonality of the Hadley cell. J. Atmos. Sci., 60, 1522–1526, doi:10.1175/ 1520-0469(2003)060,1522:OTSOTH.2.0.CO;2. Feng, J., and J. P. Li, 2013: Contrasting impacts of two types of ENSO on the boreal spring Hadley circulation. J. Climate, 26, 4773–4789, doi:10.1175/JCLI-D-12-00298.1. ——, ——, and F. Xie, 2013: Long-term variation of the principal mode of boreal spring Hadley circulation linked to SST over the Indo-Pacific warm pool. J. Climate, 26, 532–544, doi:10.1175/ JCLI-D-12-00066.1.

VOLUME 29

Feng, R., J. P. Li, and J. C. Wang, 2011: Regime change of the boreal summer Hadley circulation and its connection with the tropical SST. J. Climate, 24, 3867–3877, doi:10.1175/2011JCLI3959.1. Fu, Q., C. M. Johanson, J. M. Wallace, and T. Reichler, 2006: Enhanced mid-latitude tropospheric warming in satellite measurements. Science, 312, 1179, doi:10.1126/science.1125566. Gill, A. E., 1980: Some simple solutions for heat-induced tropical circulation. Quart. J. Roy. Meteor. Soc., 106, 447–462, doi:10.1002/ qj.49710644905. Guo, Y. P., J. P. Li, J. Feng, F. Xie, C. Sun, and J. Zheng, 2016: The multidecadal variability of the asymmetric mode of the boreal autumn Hadley circulation and its link to the Atlantic multidecadal oscillation. J. Climate, 29, 5625–5641, doi:10.1175/ JCLI-D-15-0025.1. Held, I. M., and A. Y. Hou, 1980: Nonlinear axially symmetric circulations in a nearly inviscid atmosphere. J. Atmos. Res., 37, 515– 533, doi:10.1175/1520-0469(1980)037,0515:NASCIA.2.0.CO;2. Hou, A. Y., 1998: Hadley circulation as a modulator of the extratropical climate. J. Climate, 55, 2437–2457, doi:10.1175/ 1520-0469(1998)055,2437:HCAAMO.2.0.CO;2. ——, and R. S. Lindzen, 1992: The influence of concentrated heating on the Hadley circulation. J. Atmos. Sci., 49, 1233–1241, doi:10.1175/1520-0469(1992)049,1233:TIOCHO.2.0.CO;2. Hu, Y. Y., and Q. Fu, 2007: Observed poleward expansion of the Hadley circulation since 1979. Atmos. Chem. Phys., 7, 5229– 5236, doi:10.5194/acp-7-5229-2007. ——, and C. Zhou, 2009: Decadal changes in the Hadley circulation. Atmospheric Science, J. H. Oh, Ed., Advances in Geosciences, Vol. 16, World Scientific Publishing Company, 61–74. ——, K.-K. Tung, and J. Liu, 2005: A closer comparison of early and late-winter atmospheric trends in the Northern Hemisphere. J. Climate, 18, 3204–3216, doi:10.1175/JCLI3468.1. Hudson, R. D., M. F. Andrade, M. B. Follette, and A. D. Frolov, 2006: The total ozone field separated into meteorological regimes—Part II: Northern Hemisphere mid-latitude total ozone trends. Atmos. Chem. Phys., 6, 5183–5191, doi:10.5194/ acp-6-5183-2006. Jiang, X. W., and J. P. Li, 2011: Influence of the annual cycle of sea surface temperature on the monsoon onset. J. Geophys. Res., 116, D10105, doi:10.1029/2010JD015236. Johanson, C., and Q. Fu, 2009: Hadley cell widening: Model simulations versus observations. J. Climate, 22, 2713–2725, doi:10.1175/2008JCLI2620.1. Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437–472, doi:10.1175/1520-0477(1996)077,0437:TNYRP.2.0.CO;2. Li, J. P., and Coauthors, 2013: Progress in air–land–sea interactions in Asia and their role in global and Asian climate change (in Chinese). Chin. J. Atmos. Sci., 37, 518–538. ——, N. Xing, L. N. Wang, and Y. K. Li, 2016: Responses of the tropical atmosphere to bipolar heating and the explanation of anomalous atmospheric circulation associated with Indian Ocean dipole events. Climate Dyn., in press. Lindzen, R. S., 1994: Climate dynamics and global change. Annu. Rev. Fluid Mech., 26, 353–378, doi:10.1146/ annurev.fl.26.010194.002033. ——, and S. Nigam, 1987: On the role of sea surface temperature gradients in forcing low-level winds and convergence in the tropics. J. Atmos. Sci., 44, 2418–2436, doi:10.1175/ 1520-0469(1987)044,2418:OTROSS.2.0.CO;2. ——, and A. Y. Hou, 1988: Hadley circulations for zonally averaged heating centered off the equator. J. Atmos. Sci., 45, 2416–2427, doi:10.1175/1520-0469(1988)045,2416:HCFZAH.2.0.CO;2.

15 DECEMBER 2016

FENG ET AL.

Liu, Z. Y., S. Vavrus, F. He, N. Wen, and Y. F. Zhong, 2005: Rethinking tropical ocean response to global warming: The enhanced equatorial warming. J. Climate, 18, 4684–4700, doi:10.1175/JCLI3579.1. Lu, J., G. A. Vecchi, and T. Reichler, 2007: Expansion of the Hadley cell under global warming. Geophys. Res. Lett., 34, L06805, doi:10.1029/2006GL028443. Ma, J., and J. P. Li, 2007: Strengthening of the boreal winter Hadley circulation and its connection with ENSO. Prog. Nat. Sci., 17, 1327–1333. ——, and ——, 2008: The principal modes of variability of the boreal winter Hadley cell. Geophys. Res. Lett., 35, L01808, doi:10.1029/2007GL031883. 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. Amer. Meteor. Soc., 78, 1069–1079, doi:10.1175/1520-0477(1997)078,1069: APICOW.2.0.CO;2. Mitas, C. M., and A. Clement, 2005: Has the Hadley cell been strengthening in recent decades? Geophys. Res. Lett., 32, L03809, doi:10.1029/2004GL021765. ——, and ——, 2006: Recent behavior of the Hadley cell and tropical thermodynamics in climate models and reanalyses. Geophys. Res. Lett., 33, L01810, doi:10.1029/2005GL024406. Oort, A. H., and J. J. Yienger, 1996: Observed interannual variability in the Hadley circulation and its connection to ENSO. J. Climate, 9, 2751–2767, doi:10.1175/1520-0442(1996)009,2751: OIVITH.2.0.CO;2. Quan, X. W., H. F. Diaz, and M. P. Hoerling, 2004: Change in the tropical Hadley cell since 1950. Hadley Circulation: Present, Past and Future, H. F. Diaz and R. S. Bradley, Eds., Springer, 85–120. Ratnam, J. V., S. K. Behera, Y. Masumoto, K. Takahashi, and T. Yamagata, 2012: Anomalous climatic conditions associated with the El Niño Modoki during boreal winter of 2009. Climate Dyn., 39, 227–238, doi:10.1007/s00382-011-1108-z. Rayner, N. A., D. E. Parker, E. B. Horton, C. K. Folland, L. V. Alexander, D. P. Rowell, E. C. Kent, and A. Kaplan, 2003:

8963

Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century. J. Geophys. Res., 108, 4407, doi:10.1029/2002JD002670. Rind, D., and W. B. Rossow, 1984: The effects of physical processes on the Hadley circulation. J. Atmos. Sci., 41, 479–507, doi:10.1175/1520-0469(1984)041,0479:TEOPPO.2.0.CO;2. Schneider, E., and R. S. Lindzen, 1977: Axially symmetric steady-state models of the basic state for instability and climate studies. Part I: Linearized calculations. J. Atmos. Sci., 34, 253–279, doi:10.1175/1520-0469(1977)034,0263: ASSSMO.2.0.CO;2. Seidel, D. J., and R. J. Randel, 2007: Recent widening of the tropical belt: Evidence from tropopause observations. J. Geophys. Res., 112, D20113, doi:10.1029/2007JD008861. Smith, T. M., and R. W. Reynolds, 2004: Improved extended reconstruction of SST (1854–1997). J. Climate, 17, 2466–2477, doi:10.1175/1520-0442(2004)017,2466:IEROS.2.0.CO;2. Sun, Y., and T. J. Zhou, 2014: How does El Niño affect the interannual variability of the boreal summer Hadley circulation? J. Climate, 27, 2622–2642, doi:10.1175/JCLI-D-13-00277.1. Tanaka, H. L., N. Ishizaki, and A. Kitoh, 2004: Trend and interannual variability of Walker, monsoon and Hadley circulations defined by velocity potential in the upper troposphere. Tellus, 56A, 250–269, doi:10.1111/j.1600-0870.2004.00049.x. Wielicki, B. A., and Coauthors, 2002: Evidence for large decadal variability in the tropical mean radiative energy budget. Science, 295, 841–843, doi:10.1126/science.1065837. Xiao, D., and J. P. Li, 2007: Spatial and temporal characteristics of the decadal abrupt changes of global atmosphere–ocean system in the 1970s. J. Geophys. Res., 112, D24S22, doi:10.1029/ 2007JD008956. Xie, S. P., C. Deser, G. A. Vecchi, J. Ma, H. Y. Teng, and A. T. Wittenberg, 2010: Global warming pattern formation: Sea surface temperature and rainfall. J. Climate, 23, 966–986, doi:10.1175/2009JCLI3329.1. Xing, N., J. P. Li, and Y. K. Li, 2014: Response of the tropical atmosphere to isolated equatorial asymmetric heating (in Chinese). Chin. J. Atmos. Sci., 38, 1147–1158.