# Lagged Responses of the Tropical Pacific to the 11-yr Solar Cycle Forcing and Possible Mechanisms

## 热带太平洋对太阳准11年周期强迫的滞后响应及其可能机制

• Corresponding author: Ziniu XIAO, xiaozn@lasg.iap.ac.cn
• Funds:

Supported by the National Key Basic Research and Development (973) Program of China (2012CB957804); Project from State Key Laboratory of Tropical Oceanography, South China Sea Institute of Oceanology, Chinese Academy of Sciences (LTO1916); National Natural Science Foundation of China (42075040); and Science and Technology Project of State Grid Corporation of China (SGCC; NY71-19-013)

• doi: 10.1007/s13351-021-0137-8
• This paper uses two subsets of ensemble historical-Nat simulations and pi-Control simulations from CMIP5 as well as observational/reanalysis datasets to investigate responses of the tropical Pacific to the 11-yr solar cycle. A statistically significant 11-yr solar signal is found in the upper-ocean layers above the thermocline and tropospheric circulations. A warming response initially appears in the upper layers of the central equatorial Pacific in the solar maximum years in observations, then increases and shifts into the eastern Pacific at lagged 1–3 yr. Meanwhile, an anomalous updraft arises over the western equatorial Pacific and shifts eastwards in the following years with anomalous subsidence over the Maritime Continent. These lagged responses are confirmed by the historical-Nat simulations, except that the initial signal is located more to the west and all the responses are weaker than the observed. A simplified mixed-layer heat budget analysis based on the historical-Nat simulations suggests that the atmospheric forcing, especially the shortwave radiation, is the major contributor to the initial warming response, and the ocean heat transport effect is responsible for the eastward displacement of the lagged warming responses. In the solar maximum years, the zonal ocean temperature gradient in the western–central Pacific is reduced by the initial warming, and anomalous westerly winds appear over the western equatorial Pacific and extend into the eastern Pacific during the lagged years. These anomalous westerly winds reduce the wind-driven ocean dynamical transport, resulting in the initial warming in the central equatorial Pacific being amplified and the surface warming shifting eastward during the lagged 1–3 yr.

本文采用来自CMIP5自然强迫历史气候模拟和工业革命前控制实验的两个子集，以及观测和再分析数据集，分析了热带太平洋对太阳准11年周期强迫的响应，为提高年代际气候预测技巧提供参考。本文证实了在热带太平洋海洋温度中存在独立于ENSO循环的准11年太阳周期信号，对太阳最大值的显著增暖响应首先出现在赤道中太平洋，并随滞后时间向东移动。而自然强迫历史气候模拟中的响应强度比观测弱，且初始信号位置偏西。基于简化的混合层热量收支诊断可以发现，增强的太阳短波辐射（在太阳活动最大值年）是赤道太平洋初始增暖响应的主要贡献者，在滞后太阳周期峰值的1–3年里，由于海气耦合过程，增暖响应逐渐增强并引起热带太平洋的海洋次表层和对流层环流出现响应。

• Fig. 1.  The reconstructed TSI (solid gray line at top) used in the CMIP5 historical-Nat simulations and its 3-yr running mean (solid black line at top), as well as the observed annual mean SSN (dotted gray line at bottom) and 3-yr running mean SSN (dotted black line at bottom). The black dots (circles) indicate the solar maximum (minimum) years used for composite analysis in this study.

Fig. 2.  (a–d) Lagged correlation maps between TSI and the observational annual SST anomaly (ERSST; only 1950–2018 is used here). (e–h) As in (a–d), but between TSI and the MME SST anomaly from CMIP5 historical-Nat simulations (1861–2005). The dotted (solid) lines indicate a negative (positive) correlation. Black dots indicate the 95% confidence level (two-tailed Student’s t-test and controlling FDR method).

Fig. 3.  (a) The spatial pattern of the EOF1 of the annual mean MME SST anomaly over the tropical Pacific region in the historical-Nat simulations, (b) the corresponding PC1, and (c) the spectrum of PC1 (black thick solid line). The black thin (dotted) lines represent the 90% (95%) confidence level of the Markov “red noise” spectrum, and the red lines are the same as the black ones but for the spectrum of TSI. (d–f) As in (a–c), but for the pi-Control simulations.

Fig. 4.  Composite differences of the annual mean SST anomaly (ERSST; shaded; only values above the 90% confidence level are shown) and wind anomaly at 850 hPa (NCEP/NCAR Reanalysis 1; vector) between the (a) solar maximum and minimum years and (b–d) lagged 3 yr for 1950–2018. Green vectors indicate that it is above the 90% confidence level for zonal wind (bootstrapping difference).

Fig. 5.  As in Fig. 4, but for the annual mean MME SST anomaly (shaded and only values above the 90% confidence level are shown) and MME wind anomaly at 850 hPa (vector) from the CMIP5 historical-Nat simulations.

Fig. 6.  Composite differences of the annual mean MME vertical velocity anomaly averaged over the equatorial Pacific (5°S–5°N) in the CMIP5 historical-Nat simulations for (a) solar maximum minus minimum years and (b–d) lagged 3 yr. Black dots indicate the 95% confidence level (bootstrapping difference). Downward is positive (solid lines).

Fig. 7.  As in Fig. 6, but for the annual mean subsurface temperature anomaly from the objective analyses EN4 dataset (1950–2018). Black dots indicate the 95% confidence level (bootstrapping difference) and the thick black line represents the 20°C-isotherm depth (m).

Fig. 8.  As in Fig. 7, but for the annual mean MME subsurface temperature anomaly in CMIP5 historical-Nat simulations.

Fig. 9.  Lagged composite solar maximum minus minimum differences for the (a) annual mean SST change and associated (b) atmospheric forcing (${Q_{\rm{a}}}$), and (c) ocean transport effect (${D_{\rm{o}}}$) of the histori-cal-Nat simulations, in which all variables are averaged over 10°S–10°N. Black dots in (a) indicate exceedance of the 95% confidence level (bootstrapping difference test).

Fig. 10.  Lagged composite solar maximum minus minimum differences (W m−2) for the (a–d) surface net shortwave radiation (${Q_{\rm{S}}}$), (e–h) surface net longwave radiation (${Q_{\rm{L}}}$), (e–h) net sensible heat flux (${Q_{\rm{H}}}$), and (m–p) latent heat flux from atmospheric forcing ($Q_{\rm{E}}^{\rm{a}}$) after subtracting the ocean Newtonian cooling response. Downward is positive. Black dots indicate exceedance of the 95% confidence level (bootstrapping difference test).

Fig. 11.  Lagged composite solar maximum minus minimum differences of the annual mean MME total cloud fraction anomaly in historical-Nat simulations (color shading). Black dots indicate exceedance of the 95% confidence level (bootstrapping difference test) and the white contours in (a) are the climatological total cloud fraction averaged over 1861–2005.

•  [1] Bretherton, C. S., M. Widmann, V. P. Dymnikov, et al., 1999: The effective number of spatial degrees of freedom of a time-varying field. J. Climate, 12, 1990–2009.. [2] Brönnimann, S., 2007: Impact of El Niño–Southern Oscillation on European climate. Rev. Geophys., 45, RG3003.. [3] Cai, M., and K. K. Tung, 2012: Robustness of dynamical feedbacks from radiative forcing: 2% solar versus 2 × CO2 experiments in an idealized GCM. J. Atmos. Sci., 69, 2256–2271.. [4] Cionni, I., V. Eyring, J. F. Lamarque, et al., 2011: Ozone database in support of CMIP5 simulations: results and corresponding radiative forcing. Atmos. Chem. Phys., 11, 11,267–11,292.. [5] Dai, A. G., and T. M. L. Wigley, 2000: Global patterns of ENSO-induced precipitation. Geophys. Res. Lett., 27, 1283–1286.. [6] de Szoeke, S. P., S. P. Xie, T. Miyama, et al., 2007: What maintains the SST front north of the eastern Pacific equatorial cold tongue? J. Climate, 20, 2500–2514.. [7] Diaconis, P., and B. Efron, 1983: Computer-intensive methods in statistics. Sci. Amer., 248, 116–130.. [8] Dima, M., and M. Voiculescu, 2016: Global patterns of solar influence on high cloud cover. Climate Dyn., 47, 667–678.. [9] Du, Y., and S. P. Xie, 2008: Role of atmospheric adjustments in the tropical Indian Ocean warming during the 20th century in climate models. Geophys. Res. Lett., 35, L08712.. [10] Flato, G., J. Marotzke, B. Abiodun, et al., 2013: Evaluation of Climate Models. In: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Stocker, T. F., D. Qin, G. -K. Plattner, et al., Eds., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 741–866. [11] Frame, T. H. A., and L. J. Gray, 2010: The 11-yr solar cycle in ERA-40 data: An update to 2008. J. Climate, 23, 2213–2222.. [12] Gershunov, A., and T. P. Barnett, 1998: ENSO influence on intraseasonal extreme rainfall and temperature frequencies in the contiguous United States: Observations and model results. J. Climate., 11, 1575–1586.. [13] Good, S. A., M. J. Martin, and N. A. Rayner, 2013: EN4: Quality controlled ocean temperature and salinity profiles and monthly objective analyses with uncertainty estimates. J. Geophys. Res. Oceans, 118, 6704–6716.. [14] Gray, L. J., S. A. Crooks, M. A. Palmer, et al., 2006: A possible transfer mechanism for the 11-year solar cycle to the lower stratosphere. Space Sci. Rev., 125, 357–370.. [15] Gray, L. J., S. T. Rumbold, and K. P. Shine, 2009: Stratospheric temperature and radiative forcing response to 11-year solar cycle changes in irradiance and ozone. J. Atmos. Sci., 66, 2402–2417.. [16] Haam, E., and K. K. Tung, 2012: Statistics of solar cycle–La Niña connection: Correlation of two autocorrelated time series. J. Atmos. Sci., 69, 2934–2939.. [17] Huo, W. J. and Z. N. Xiao, 2017a: Anomalous pattern of ocean heat content during different phases of the solar cycle in the tropical Pacific. Atmos. Ocean. Sci. Lett., 10, 9–16.. [18] Huo, W. J., and Z. N., Xiao, 2017b: Modulations of solar activity on El Niño Modoki and possible mechanisms. J. Atmos. So-lar-Terr. Phys., 160, 34–47.. [19] Khodri, M., T. Izumo, J. Vialard, et al., 2017: Tropical explosive volcanic eruptions can trigger El Niño by cooling tropical Africa. Nat. Commun., 8, 778.. [20] Kodera, K., and Y. Kuroda, 2002: Dynamical response to the solar cycle. J. Geophys. Res. Atmos., 107, 4749.. [21] Kodera, K., R. Thiéblemont, S. Yukimoto, et al., 2016: How can we understand the global distribution of the solar cycle signal on the Earth’s surface? Atmos. Chem. Phys., 16, 12,925–12,944.. [22] Kopp, G., and J. L. Lean, 2011: A new, lower value of total solar irradiance: Evidence and climate significance. Geophys. Res. Lett., 38, L01706.. [23] Kosaka, Y., and S. P. Xie, 2013: Recent global-warming hiatus tied to equatorial Pacific surface cooling. Nature, 501, 403–407.. [24] Li, G., and S. P. Xie, 2014: Tropical biases in CMIP5 multimodel ensemble: The excessive equatorial Pacific cold tongue and double ITCZ problems. J. Climate, 27, 1765–1780.. [25] Matthes, K., Y. Kuroda, K. Kodera, et al., 2006: Transfer of the solar signal from the stratosphere to the troposphere: Northern winter. J. Geophys. Res. Atmos., 111, D06108.. [26] Meehl, G. A., and J. M. Arblaster, 2009: A lagged warm event-like response to peaks in solar forcing in the Pacific region. J. Climate, 22, 3647–3660.. [27] Meehl, G. A., J. M. Arblaster, G. Branstator, et al., 2008: A coupled air–sea response mechanism to solar forcing in the Pacific region. J. Climate, 21, 2883–2897.. [28] Meehl, G. A., M. Arblaster, J. K. Matthes, et al., 2009: Amplifying the Pacific climate system response to a small 11-year solar cycle forcing. Science, 325, 1114–1118.. [29] Misios, S., and H. Schmidt, 2012: Mechanisms involved in the amplification of the 11-yr solar cycle signal in the tropical Pacific Ocean. J. Climate, 25, 5102–5118.. [30] Misios, S., D. M. Mitchell, L. J. Gray, et al., 2016: Solar signals in CMIP-5 simulations: effects of atmosphere–ocean coupling. Quart. J. Roy. Meteor. Soc., 142, 928–941. doi: 10.1002/qj.2695. [31] Misios, S., L. J. Gray, M. F. Knudsen, et al., 2019: Slowdown of the Walker circulation at solar cycle maximum. Proc. Natl. Acad. Sci. USA, 116, 7186–7191.. [32] Philander, S. G. H., 1981: The response of equatorial oceans to a relaxation of the trade winds. J. Phys. Oceanogr., 11, 176–189.. [33] Pyper, B. J., and R. M. Peterman, 1998: Comparison of methods to account for autocorrelation in correlation analyses of fish data. Can. J. Fish. Aquat. Sci., 55, 2127–2140. doi: 10.1139/f98-104. [34] Roy, I., and J. D. Haigh, 2010: Solar cycle signals in sea level pressure and sea surface temperature. Atmos. Chem. Phys., 10, 3147–3153.. [35] Roy, I., and J. D. Haigh, 2012: Solar cycle signals in the Pacific and the issue of timings. J. Atmos. Sci., 69, 1446–1451.. [36] Schneider, E. K., and M. Z. Fan, 2012: Observed decadal North Atlantic tripole SST variability. Part II: Diagnosis of mechanisms. J. Atmos. Sci., 69, 51–64.. [37] Thiéblemont, R., K. Matthes, N. E. Omrani, et al., 2015: Solar forcing synchronizes decadal North Atlantic climate variability. Nat. Commun., 6, 8268.. [38] Tung, K. K., and J. S. Zhou, 2010: The Pacific’s response to surface heating in 130 yr of SST: La Niña–like or El Niño–like? J. Atmos. Sci., 67, 2649–2657.. [39] van Loon, H., G. A. Meehl, and D. J. Shea, 2007: Coupled air–sea response to solar forcing in the Pacific region during northern winter. J. Geophys. Res. Atmos., 112, D02108.. [40] Wang, W. K., K. Matthes, W. S. Tian, et al., 2019: Solar impacts on decadal variability of tropopause temperature and lower stratospheric (LS) water vapour: a mechanism through ocean-atmosphere coupling. Climate Dyn., 52, 5585–5604.. [41] Wang, Y. M., J. L. Lean, and N. R. Jr. Sheeley, 2005: Modeling the Sun’s magnetic field and irradiance since 1713. Astrophys. J., 625, 522–538. doi: 10.1086/429689. [42] Ward, P. J., B. Jongman, M. Kummu, et al., 2014: Strong influence of El Niño Southern Oscillation on flood risk around the world. Proc. Natl. Acad. Sci. USA, 111, 15,659–15,664.. [43] White, W. B., and Z. Y. Liu, 2008: Resonant excitation of the quasi-decadal oscillation by the 11-year signal in the Sun’s irradiance. J. Geophys. Res. Oceans, 113, C01002.. [44] White, W. B., J. Lean, D. R. Cayan, et al., 1997: Response of global upper ocean temperature to changing solar irradiance. J. Geophys. Res. Oceans, 102, 3255–3266.. [45] Wilks, D. S., 2016: “The stippling shows statistically significant grid points”: How research results are routinely overstated and overinterpreted, and what to do about it. Bull. Amer. Meteor. Soc., 97, 2263–2273.. [46] Xiao, Z. N., Y. C. Liao, and C. Y. Li, 2016: Possible impact of so-lar activity on the convection dipole over the tropical pacific ocean. J. Atmos. Solar-Terr. Phys., 140, 94–107.. [47] Xie, S. P., C. Deser, G. A. Vecchi, et al., 2010: Global warming pattern formation: Sea surface temperature and rainfall. J. Climate, 23, 966–986.. [48] Zhang, Y., J. M. Wallace, and D. S. Battisti, 1997: ENSO-like interdecadal variability: 1900–93. J. Climate, 10, 1004–1020.. [49] Zhou, T. J., B. Wu, and L. Dong, 2014: Advances in research of ENSO changes and the associated impacts on Asian-Pacific climate. Asia-Pac. J. Atmos. Sci., 50, 405–422..
###### 通讯作者: 陈斌, bchen63@163.com
• 1.

沈阳化工大学材料科学与工程学院 沈阳 110142

## Lagged Responses of the Tropical Pacific to the 11-yr Solar Cycle Forcing and Possible Mechanisms

###### Corresponding author: Ziniu XIAO, xiaozn@lasg.iap.ac.cn;
• 1. State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
• 2. Research Division Ocean Circulation and Climate, GEOMAR Helmholtz Centre for Ocean Research, Kiel 24105, Germany
• 3. State Key Laboratory of Tropical Oceanography, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou 510301, China
• 4. Innovation Academy of South China Sea Ecology and Environmental Engineering,Chinese Academy of Sciences, Guangzhou 510301, China
• 5. Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511458, China
Funds: Supported by the National Key Basic Research and Development (973) Program of China (2012CB957804); Project from State Key Laboratory of Tropical Oceanography, South China Sea Institute of Oceanology, Chinese Academy of Sciences (LTO1916); National Natural Science Foundation of China (42075040); and Science and Technology Project of State Grid Corporation of China (SGCC; NY71-19-013)

Abstract:

This paper uses two subsets of ensemble historical-Nat simulations and pi-Control simulations from CMIP5 as well as observational/reanalysis datasets to investigate responses of the tropical Pacific to the 11-yr solar cycle. A statistically significant 11-yr solar signal is found in the upper-ocean layers above the thermocline and tropospheric circulations. A warming response initially appears in the upper layers of the central equatorial Pacific in the solar maximum years in observations, then increases and shifts into the eastern Pacific at lagged 1–3 yr. Meanwhile, an anomalous updraft arises over the western equatorial Pacific and shifts eastwards in the following years with anomalous subsidence over the Maritime Continent. These lagged responses are confirmed by the historical-Nat simulations, except that the initial signal is located more to the west and all the responses are weaker than the observed. A simplified mixed-layer heat budget analysis based on the historical-Nat simulations suggests that the atmospheric forcing, especially the shortwave radiation, is the major contributor to the initial warming response, and the ocean heat transport effect is responsible for the eastward displacement of the lagged warming responses. In the solar maximum years, the zonal ocean temperature gradient in the western–central Pacific is reduced by the initial warming, and anomalous westerly winds appear over the western equatorial Pacific and extend into the eastern Pacific during the lagged years. These anomalous westerly winds reduce the wind-driven ocean dynamical transport, resulting in the initial warming in the central equatorial Pacific being amplified and the surface warming shifting eastward during the lagged 1–3 yr.

Reference (49)

/