Journal of Astronomy and Space Sciences
The Korean Space Science Society
Research Paper

Association between Solar Variability and Teleconnection Index

Jung-Hee Kim1, Heon-Young Chang1,2,http://orcid.org/0000-0003-2015-2725
1Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Daegu 41566, Korea
2Research and Training Team for Future Creative Astrophysicists and Cosmologists (BK21 Plus Program), Daegu 41566, Korea
Corresponding Author Tel: +82-53-950-6367, E-mail: hyc@knu.ac.kr

© The Korean Space Science Society. All rights reserved. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Aug 7, 2019; Revised: Aug 15, 2019; Accepted: Aug 16, 2019

Abstract

In this study, we investigate the associations between the solar variability and teleconnection indices, which influence atmospheric circulation and subsequently, the spatial distribution of the global pressure system. A study of the link between the Sun and a large-scale mode of climate variability, which may indirectly affect the Earth’s climate and weather, is crucial because the feedbacks of solar variability to an autogenic or internal process should be considered with due care. We have calculated the normalized cross-correlations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with teleconnection indices. We have found that the Southern Oscillation Index (SOI) index is anti-correlated with both solar activity and the solar North–South asymmetry, with a ∼3-year lag. This finding not only agrees with the fact that El Niño episodes are likely to occur around the solar maximum, but also explains why tropical cyclones occurring in the solar maximum periods and in El Niño periods appear similar. Conversely, other teleconnection indices, such as the Arctic Oscillation (AO) index, the Antarctic Oscillation (AAO) index, and the Pacific-North American (PNA) index, are weakly or only slightly correlated with solar activity, which emphasizes that response of terrestrial climate and weather to solar variability are local in space. It is also found that correlations between teleconnection indices and solar activity are as good as correlations resulting from the teleconnection indices themselves.

Keywords: solar activity; terrestrial climate; data analysis

1. INTRODUCTION

One of the most intensively studied current multidisciplinary scientific issues is estimating the anthropogenic contribution to global warming. This has led to investigating the possible role of the Sun with respect to a better understanding of natural climate change. As a result, it has been shown that the terrestrial climate system responds to the solar variability via changes in the Sun’s total irradiance (Frohlich 2006), the solar extreme ultraviolet (EUV) irradiance (Emmert & Picone 2010;Gray et al. 2017), the Earth’s magnetic/electric systems (Burns et al. 2007, 2008), galactic cosmic ray (GCR) influx (Tinsley & Deen 1991;Pudovkin et al. 1997;Marsh & Svensmark 2000;Roldugin & Tinsley 2004;Pudovkin 2004;Haigh 2007;Bazilevskaya et al. 2008;Artamonova & Veretenenko 2011), solar energetic particle events (Veretenenko & Thejll 2004;Cho et al. 2012;Mironova et al. 2012;Mironova & Usoskin 2013, 2014), highspeed solar winds (Zhou et al. 2016), etc. In addition, the solar North–South asymmetry has recently been considered as one of the parameters affecting climate conditions (e.g., Cho et al. 2012;Kim et al. 2017;Chang 2018;Kim et al. 2018).

Even though the manifestation of the response to solar activity appears to be more local in both time and space than global, there is indeed phenomenological evidence that solar activity is linked to changes in terrestrial weather and climate on the timescales of days to several tens of year (Svensmark & Friis-Christensen 1997;Marsh & Svensmark 2000;Tinsley 2000;Krivova & Solanki 2004;Pudovkin 2004;Roldugin & Tinsley 2004;Scafetta & West 2006;Haigh 2007;Burns et al. 2007, 2008;Cho & Chang 2008;Kniveton et al. 2008;Meehl et al. 2009;Gray et al. 2010;Cho et al. 2012;Gray et al. 2013;Lee & Yi 2018;Muraki 2018). For example, the solar variability is associated with the El Niño-Southern Oscillation (ENSO) variability (e.g., van Loon et al. 2007;van Loon & Meehl 2008;Meehl et al. 2008, 2009;Roy & Haigh, 2010, 2012;Zhou & Tung 2010;Haam & Tung 2012), the Quasi Biennial Oscillation (QBO) (Labitzke 1987;Labitzke & van Loon 1988;Sagir et al. 2015), and the North Atlantic Oscillation (NAO) (Gray et al. 2013). Furthermore, solar activity is related to various types of extreme weather, such as tropical cyclones (Kim et al. 2017;Kim et al. 2018), severe droughts (Bhalme et al. 1983;Meshcherskaya & Blazhevich 1997;Park & Chang 2013;Park et al. 2014), changes in the precipitation patterns of rainy seasons (Mazzarella & Palumbo 1992;Kavlakov 2005;Pérez-Peraza et al. 2008;Ho et al. 2016;Choi et al. 2017), and heatwaves (Coumou & Rahmstorf 2012).

Low solar activity would allow a greater GCR influx into the Earth’s atmosphere. Consequently, more cloud condensation nuclei resulting from GCR via either the ionaerosol clear-sky mechanism or ion-aerosol near-cloud mechanism form to increase cloudiness and in turn affect the tropical cyclone activity. Another possible mechanism is that high solar activity results in warming of the lower stratosphere and the upper troposphere, decreasing the convective available potential energy and the intensity of tropical cyclones (Elsner & Jagger 2008;Hung 2013). Particularly, typhoons in the western North Pacific influencing the Korean peninsula are of great interest because typhoons occur here much more frequently than in any other tropical cyclone genesis basin. Thus, tropical cyclone activity over the Northwestern Pacific basin was examined to determine whether there is a long-term trend, and if so, whether it is modulated by solar activity (Kim et al. 2017;Kim et al. 2018). Kim et al. (2017, 2018) attempted to address questions regarding solar influences on tropical cyclone genesis in the Northwestern Pacific basin by analyzing locations where tropical cyclones form, reach their lowest central pressure, and downgrade back to a tropical depression. They found that tropical cyclones occurring in the solar maximum periods and in the solar descending periods have properties similar to those in El Niño periods. Unfortunately, however, their subsequent investigations show that that conclusion turns out to be inconclusive. Based on the linear correlation coefficients r between the locations of tropical cyclones and observed data of sunspot areas, they concluded that tropical cyclone properties appear weakly correlated or anti-correlated with the solar cycle, although some long-term trends can be observed. In spite of the fact that tropical cyclone activity seems to share periodicities with solar activity, they have also warned, with the wavelet transform method, that there is no evidence showing that tropical cyclone activity is causally related to solar activity, including positions of occurrence/development and other characteristics. This result essentially occurs because the wavelet transform method is inappropriate for searching for causality.

In this paper, we investigate associations between the solar variability and teleconnection patterns, which are spatial and/or temporal climate anomalies connected to each other over a large distance, typically thousands of kilometers. Such a large-scale anomaly surely influences the atmospheric circulation and subsequent spatial distribution of the global pressure system, so that climate factors governing tropical cyclone genesis, movement, and recurvature locations, as well as lifetime, are changed (Choi & Moon 2012; Choi & Cha 2017; Choi et al. 2017). That is, unlike earlier approaches, we study a link between the Sun and a large-scale mode of climate variability, which may indirectly affect the Earth’s climate and weather, instead of a direct link between the Sun and terrestrial climate and weather. This is crucial in the sense that feedbacks of solar variability to an autogenic or internal process, such as a long-term variation, should receive attention because solar forcing as an ingredient of the climate system is nonlinear in nature. To find a resemblance in shape of two physical quantities, we have calculated the normalized crosscorrelations, defined by

z ( t ) = u ( t )  v ( t ) E u E v ,
(1)

Where

w ( t ) = u ( t ) v ( t ) = + u * ( τ t ) v ( τ ) d τ ,
(2)

and

E u = + | u ( t ) | 2 d t ,
(3)

where u(τ) is the complex conjugate of u(τ). For all t, |z(t)| ≤ 1. Note that the argument of w(t) is called the lag if it is positive, or the lead if it is negative. That is, if the argument of w(t) is positive, the first quantity follows the second quantity by the amount of its value, and vice versa. The correlation of a signal with itself is called the autocorrelation. The autocorrelation at zero lag corresponds to the linear correlation coefficient. This paper is organized as follows. We begin with a description of the data analysis in Section 2. We present and discuss the results in Section 3. Finally, we briefly summarize and conclude in Section 4.

2. DATA

As a proxy of solar activity, we have taken the sunspot area and sunspot number data from the NASA Marshall Space Flight Center website (http://solarscience.msfc.nasa. gov/greenwch.shtml). We extracted the daily sunspot area (in units of millionths of a hemisphere) and daily sunspot number during the period from 1977 to 2016, approximately spanning from solar cycles 21 to 24. The sunspot dataset is available as ASCII text files for both the Northern Hemisphere and the Southern Hemisphere, separately. To calculate the solar North–South asymmetry, we have used sunspot area data, defined by (ASAN)/(AS + AN), where AS and AN are the sunspot area appearing in the solar Southern and Northern hemispheres, respectively. For the timeseries of the teleconnection indices, such as the Southern Oscillation Index (SOI) the NAO index, the Arctic Oscillation (AO) index, the Antarctic Oscillation (AAO) index, and the Pacific-North American (PNA) index, we have downloaded daily data from the NOAA website (https://www.ncdc. noaa.gov/teleconnections), which provides various indices representing large-scale anomalies in space and time that control the variability of atmospheric circulation. Having downloaded all the daily data, we averaged those data yearly so that we can compute the normalized cross-correlations by year. We are interested in yearly values because climate indices primarily vary by season.

3. RESULTS AND DISCUSSION

In Fig. 1, from left to right, we show the normalized crosscorrelations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with the SOI index, respectively. According to the definition, when the correlation values are positive (negative), there is a positive (negative) correlation. When the argument of the function at extremum is positive, the first observable lags (or follows) the SOI index by that amount. Conversely, when the argument of the function at maximum (minimum) is negative, the first observable leads (or goes ahead of ) the SOI index by that amount. The SOI index is anticorrelated with both solar activity and the solar North– South asymmetry with a lag of ∼−3 years. Because the SOI index is based on the standardized difference between observed sea-level pressures in Tahiti and Darwin, the negative SOI index corresponds to the El Niño period. This result therefore means that the El Niño episodes likely occur three years after the solar maximum. This not only agrees with the fact that El Niño episodes likely occur around the solar maximum, but also explains why tropical cyclones occurring in the solar maximum periods and in El Niño periods appear similar (e.g., Kim et al. 2017;Kim et al. 2018).

jass-36-3-149_F1
Fig. 1. Normalized cross-correlations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with the SOI index are shown from left to right, in order. The time series of data is taken during the period from 1977 to 2016. SOI, Southern Oscillation Index.
Download Original Figure

In Fig. 2, from left to right, we show the normalized crosscorrelations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with the NAO index, respectively. The NAO index is the sealevel pressure difference between the Azores High and the Subpolar Low. The positive and negative phases are associated with the North Atlantic jet stream and storm tracks. When the NAO phase is positive, temperatures in the eastern Americas and northern Europe become high. The NAO index is weakly correlated with solar activity with a lag of ∼−1 year. However, it is very weakly correlated with the solar North–South asymmetry. This might reflect the fact that the climate of regions under the influence of the NAO, including Europe, is less sensitive to solar variability, even though some European countries experience severe winters during grand solar minima, such as the Maunder minimum (e.g., Eddy 1977).

jass-36-3-149_F2
Fig. 2. Similar to Fig. 1, except the NAO index is used instead of the SOI index. NAO, North Atlantic Oscillation; SOI, Southern Oscillation Index.
Download Original Figure

In Fig. 3, from left to right, we show the normalized crosscorrelations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with the AO index, respectively. The AO is a wind pattern circulating counterclockwise around the North-polar region. It is sometimes referred to as the Northern Hemisphere annular mode (NAM). When it is positive, cold air across the polar region is confined. The AO index appears to be weakly correlated with solar activity, with a lag of ∼−1 year. It is anti-correlated with the solar North–South asymmetry with a lag of ∼−7 years. These results could mean that high solar activity may lead to a positive AO pattern, so that relatively warm winters in the Northern hemisphere occur one year after a solar maximum.

jass-36-3-149_F3
Fig. 3. Similar to Fig. 1, except the AO index is used instead of the SOI index. AO, Arctic Oscillation; SOI, Southern Oscillation Index.
Download Original Figure

In Fig. 4, from left to right, we show the normalized crosscorrelations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with the AAO index, respectively. Similar to the AO, the AAO is defined as a belt of westerly winds around Antarctica. It is also known as the Southern Annular Mode (SAM). This may lead to ice shelf basal melt. The AAO index appears to be anti-correlated with solar activity and the solar North– South asymmetry, with a lag of ∼−2 years. Note that this behavior is opposite to what is seen in Fig. 3 resulting from the AO index. It is thus interesting to investigate the cause of this difference in view of the fact that recent global warming apparently leads to warming in the Northern hemisphere and somewhat different effects in the Southern hemisphere.

jass-36-3-149_F4
Fig. 4. Similar to Fig. 1, except the AAO index is used instead of the SOI index. AAO, Antarctic Oscillation; SOI, Southern Oscillation Index.
Download Original Figure

In Fig. 5, from left to right, we show the normalized cross-correlations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with the PNA index, respectively. The PNA pattern is associated with the East Asian jet stream. When the PNA pattern is in a positive phase, the temperature in western Canada is high. The PNA index appears to be only slightly correlated with sunspot area/number, or the North–South asymmetry. Again, as in the earlier cases seen in Figs. 24, this finding emphasizes that responses of terrestrial climate and weather to solar variability are local in space and demonstrates that the solar effect on the terrestrial climate is not straightforward; rather, feedbacks play a critical role.

jass-36-3-149_F5
Fig. 5. Similar to Fig. 1, except the PNA index is used instead of the SOI index. PNA, Pacific-North American; SOI, Southern Oscillation Index.
Download Original Figure

In Fig. 6, for comparison, we show the normalized cross-correlations resulting from some combinations of teleconnection indices. On the upper right corner in each panel, the chosen computational indices are indicated. Because teleconnections influencing the weather pattern were first noticed in the late 19th century, the correlation between their time series has been probed, and spatial patterns have been compared. Associated teleconnections are now well known. For example, the PNA pattern is strongly affected by the El Niño phenomenon, and NAO phases are related to AO patterns. As seen in Figs. 16, however, the correlations between solar activity and teleconnection indices are as good as correlations resulting from the teleconnection indices themselves. Hence, although the normalized cross-correlations of the total sunspot area, the total sunspot number, and the solar North– South asymmetry with teleconnection indices are marginal, the normalized cross-correlations of solar variability with teleconnection indices cannot be ignored in comparison with those of the teleconnection indices themselves.

jass-36-3-149_F6
Fig. 6. For comparison, normalized cross-correlations resulting from combinations of teleconnection indices. Indices used in computations are indicated in the upper right corner in each panel.
Download Original Figure

4. SUMMARY AND CONCLUSIONS

Sunspot activity is associated with terrestrial weather and climate. It operates in a complex way, so that the resulting effects are nonlinear and local in nature. It may directly modulate the system of terrestrial climate via changes in the Earth’s thermal system, or Earth’s magnetic/ electric systems. More importantly, solar variability may be linked to a largescale mode of climate variability, which may indirectly affect the Earth’s climate and weather. We consider this a crucial matter because feedbacks of solar forcing on an autogenic or internal process such as longterm variation should be taken into account in the sense that a solar variability as an ingredient of the terrestrial climate system is nonlinear in nature. In this paper, we have examined associations between solar variability and teleconnection patterns, which are largescale spatial and/ or temporal climate anomalies connected to each other. We have calculated the normalized cross-correlations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with teleconnection indices including the SOI, the NAO index, the AO index, the AAO index, and the PNA index during the period from 1977 to 2016, approximately spanning from solar cycles 21 to 24.

As a result, we have found that

  • (1) the SOI index is anti-correlated with both solar activity and the solar North–South asymmetry, with a lag of ∼−3 years. This implies that El Niño episodes likely occur three years after a solar maximum. This fact both agrees with the fact that El Niño episodes likely occur around the solar maximum and explains why both tropical cyclones occurring in the solar maximum periods and in El Niño periods appear similar.

  • (2) Other teleconnection indices, such as the AO index, the AAO index, and the PNA index, are weakly or only slightly correlated with solar activity. This fact emphasizes that the responses of terrestrial climate and weather to solar variability are local in space and demonstrates that the solar effect on the terrestrial climate is not straightforward; rather, feedbacks play a crucial role.

  • (3) The correlations between the indices of solar activity and teleconnection are as good as the correlations resulting from the teleconnection indices themselves. Hence, although the normalized cross-correlations of the total sunspot area, the total sunspot number, and the solar North–South asymmetry with teleconnection indices are marginal, the normalized cross-correlations of solar variability with teleconnection indices cannot be ignored in comparison with the normalized cross-correlations of the teleconnection indices themselves.

ACKNOWLEDGMENTS

The authors thank the anonymous referees for their critical comments and helpful suggestions, which greatly improved the original version of the manuscript. This study was funded by the KMA/NMSC (Korea Meteorological Administration/National Meteorological Satellite Center)’s R&D Project, “Technical Development on Satellite Data Application for Operational Weather Service”, and by BK21 Plus of the National Research Foundation of Korea (22A20130000179). HYC was supported by a National Research Foundation of Korea grant funded by the Korean government (NRF-2018R1D1A3B070421880 and 2018R1A6A1A06024970).

References

1.

Artamonova I, Veretenenko S, Galactic cosmic ray variation influence on baric system dynamics at middle latitudes, J. Atmos. Sol.-Terr. Phys. 73, 366-370 (2011).

2.

Bazilevskaya GA, Usoskin IG, Fluckiger EO, Harrison RG, Desorgher L, et al., Cosmic ray induced ion production in the atmosphere, Space Sci. Rev. 137, 149-173 (2008). https://

3.

Bhalme HN, Mooley DA, Jadhav SK, Fluctuations in the drought/ flood area over India and relationships with the Southern Oscillation, Mon. Wea. Rev. 111, 86-94 (1983). https://doi. org/10.1175/1520-0493 (1983).

4.

Burns AG, Solomon SC, Wang W, Killeen TL, The ionospheric and thermospheric response to CMEs: challenges and successes, J. Atmos. Sol.-Terr. Phys. 69, 77-85 (2007).

5.

Burns AG, Zeng Z, Wang W, Lei J, Solomon SC, et al., Behavior of the F2 peak ionosphere over the South Pacific at dusk during quiet summer conditions from COSMIC data, J. Geophys. Res. 113, A12305 (2008).

6.

Chang HY, Latitudinal distribution of sunspot and North-South asymmetry revisited, J. Astron. Space Sci. 35, 55-66 (2018).

7.

Cho IH, Chang HY, Long term variability of the sun and climate change, J. Astron. Space Sci. 25, 395-404 (2008).

8.

Cho IH, Kwak YS, Chang HY, Cho KS, Kim YH, et al., The global temperature anomaly and solar North-South asymmetry, Asia-Pac. J. Atmos. Sci. 48, 253-257 (2012).

9.

Choi JW, Cha Y, Kim HD, Interdecadal variation of precipitation days in August in the Korean Peninsula, Dyn. Atmos. Oceans 77, 74-88 (2017).

10.

Choi KS, Moon IJ, Influence of the Western Pacific teleconnection pattern on Western North Pacific tropical cyclone activity, Dyn. Atmos. Oceans 57, 1-16 (2012).

11.

Coumou D, Rahmstorf S, A decade of weather extremes, Nat. Clim. Chang. 2, 491-496 (2012).

12.

Eddy JA, Climate and the changing sun, Clim. Chang. 1, 173-190 (1977).

13.

Elsner JB, Jagger TH, United States and Caribbean tropical cyclone activity related to the solar cycle, Geophys. Res. Lett. 35, L18705 (2008).

14.

Emmert JT, Picone JM, Climatology of globally averaged thermospheric mass density, J. Geophys. res. 115, A09326 (2010).

15.

Frohlich C, Solar Irradiance variability since 1978, Space Sci. Rev. 125, 53-65 (2006).

16.

Gray LJ, Ball W, Misios S, Solar influences on climate over the Atlantic/European sector, AIP Conf. Proc. 1810, 020002 (2017).

17.

Gray LJ, Beer J, Geller M, Haigh JD, Lockwood M, et al., Solar influences on climate, Rev. Geophys. 48, RG4001 (2010).

18.

Gray LJ, Scaife AA, Mitchell DM, Osprey S, Ineson S, A lagged response to the 11 year solar cycle in observed winter Atlantic/European weather patterns, J. Geophys. Res. Atmos. 118, 13405-13420 (2013).

19.

Haam E, Tung KK, Statistics of solar cycle-La Nina connection: correlation of two autocorrelated time series, J. Atmos. Sci. 69, 2934-2939 (2012).

20.

Haigh JD, The Sun and the Earth’s climate, Living. Rev. Sol. Phys. 4, 2 (2007).

21.

Ho CH, Choi W, Kim J, Kim MK, Yoo HD, Does El Nino-Southern Oscillation affect the precipitation in Korea on seasonal time scales? Asia-Pac. J. Atmos. Sci. 52(4), 395-403 (2016).

22.

Hung CW, A 300-year typhoon record in Taiwan and the relationship with Solar activity, Terr. Atmos. Ocean. Sci. 24, 737-743 (2013).

23.

Kavlakov SP, Global cosmic ray intensity changes, solar activity variations and geomagnetic disturbances as North Atlantic hurricane precursors, Int. J. Mod. Phys. A. 20, 6699 (2005).

24.

Kim JH, Kim KB, Chang HY, Solar influence on tropical cyclone in western north pacific ocean, J. Astron. Space Sci. 34, 257- 270 (2017).

25.

Kim KB, Kim JH, Chang HY, Do solar cycles share spectral properties with tropical cyclones that occur in the western north pacific ocean? J. Astron. Space Sci. 35, 151-161 (2018).

26.

Kniveton DR, Tinsley BA, Burns GB, Bering EA, Troshichev OA, Variations in global cloud cover and the fair-weather vertical electric field, J. Atmos. Sol.-Terr. Phys. 70, 1633-1642 (2008).

27.

Krivova NA, Solanki SK, Solar variability and global warming: A statistical comparison since 1850, Adv. Space Res, 34(2), 361- 364 (2004).

28.

Labitzke K, Sunspots, the QBO, and the stratospheric temperature in the north polar region, Geophys. Res. Lett. 14, 535-537 (1987).

29.

Labitzke K, van Loon H, Associations between the 11-year solar cycle, the QBO and the atmosphere. Part I: the troposphere and stratosphere in the northern hemisphere in winter, J. Atmos. Terr. Phys. 50, 197-206 (1988).

30.

Lee S, Yi Y, Pacific equatorial sea surface temperature variation during the 2015 El Nino period observed by advanced very-high-resolution radiometer of NOAA satellites, J. Astron. Space Sci. 35, 105-109 (2018).

31.

Marsh N, Svensmark H, Cosmic rays, clouds, and climate, Space Sci. Rev. 94, 215-230 (2000).

32.

Mazzarella A, Palumbo F, Rainfall fluctuations over Italy and their association with solar activity, Theor. Appl. Clim. 45, 201-207 (1992).

33.

Meehl GA, Arblaster JM, Branstator G, von Loon H, A coupled air-sea response mechanism to solar forcing in the Pacific region, J. Clim. 21, 2883-2897 (2008).

34.

Meehl GA, Arblaster JM, Matthes K, Sassi F, von Loon H, Amplifying the Pacific climate system response to a small 11 year solar cycle forcing, Science 325, 1114-1118 (2009).

35.

Meshcherskaya AV, Blazhevich VG, The drought and excessive moisture indices in a historical perspective in the principal grain-producing regions of the former Soviet Union, J. Clim. 10, 2670-2682 (1997).

36.

Mironova IA, Usoskin IG, Kovaltsov GA, Petelina SV, Possible effect of extreme solar energetic particle event of 20 January 2005 on polar stratospheric aerosols: direct observational evidence, Atmos. Chem. Phys. 12, 769-778 (2012).

37.

Mironova IA, Usoskin IG, Possible effect of extreme solar energetic particle events of September–October 1989 on polar stratospheric aerosols: a case study, Atmos. Chem. Phys. 13, 8543-8550 (2013).

38.

Mironova IA, Usoskin IG, Possible effect of strong solar energetic particle events on polar stratospheric aerosol: a summary of observational results, Environ. Res. Lett. 9, 015002 (2014).

39.

Muraki Y, Application of a coupled harmonic oscillator model to solar activity and El Nino phenomena, J. Astron. Space Sci. 35, 75-81 (2018).

40.

Park JH, Chang HY, Drought over Seoul and its association with solar cycles, J. Astron. Space Sci. 30, 241-246 (2013).

41.

Park JH, Kim KB, Chang HY, Statistical properties of effective drought index (EDI) for Seoul, Busan, Daegu, Mokpo in South Korea, Asia-Pac. J. Atmos. Sci. 50(4), 453-458 (2014).

42.

Perez-Peraza J, Kavlakov S, Velasco V, Gallegos-Cruz A, Azpra- Romero E, et al., Solar, geomagnetic and cosmic ray intensity changes, preceding the cyclone appearances around Mexico, Adv. Space Res. 42, 1601-1613 (2008).

43.

Pudovkin MI, Influence of solar activity on the lower atmosphere state, Int. J. Geomagn. Aeron. 5, GI2007 (2004).

44.

Pudovkin MI, Veretenenko SV, Pellinen R, Kyro E, Meteorological characteristic changes in the high-latitudinal atmosphere associated with Forbush decreases of the galactic cosmic rays, Adv. Space Res. 20, 1169-1172 (1997).

45.

Roldugin VC, Tinsley BA, Atmospheric transparency changes associated with solar wind-induced atmospheric electricity variations, J. Atmos. Sol.-Terr. Phys. 66, 1143-1149 (2004).

46.

Roy I, Haigh JD, Solar cycle signals in sea level pressure and sea surface temperature, Atmos. Chem. Phys. 10, 3147-3153 (2010).

47.

Roy I, Haigh JD, Solar cycle signals in the pacific and the issue of timings, J. Atmos. Sci. 69, 1446-1451 (2012).

48.

Sagir S, Karatay S, Atici R, Yesil A, Ozcan O, The relationship between the quasi biennial oscillation and sunspot number, Adv. Space Res. 55, 106-112 (2015).

49.

Scafetta N, West BJ, Phenomenological solar contribution to the 1900-2000 global surface warming, Geophys. Res. Lett. 33, L05708 (2006).

50.

Svensmark H, Friis-Christensen E, Variation of cosmic ray flux and global cloud coverage-a missing link in solar-climate relationships, J. Atmos. Sol.-Terr. Phys. 59, 1225-1232 (1997).

51.

Tinsley BA, Deen GW, Apparent tropospheric response to MeV-GeV particle flux variations: a connection via electrofreezing of supercooled water in high-level clouds? J. Geophys. Res. 96, 22283-22296 (1991).

52.

Tinsley BA, Influence of solar wind on the global electric circuit, and inferred effects on cloud microphysics, temperature, and dynamics in the troposphere, Space Sci. Rev. 94, 231-258 (2000).

53.

Van Loon H, Meehl GA, Shea DJ, Coupled air-sea response to solar forcing in the pacific region during northern winter, J. Geophys. Res. 112, D02108 (2007).

54.

Van Loon H, Meehl GA, The response in the pacific to the sun’s decadal peaks and contrasts to cold events in the southern oscillation, J. Atmos. Sol.-Terr. Phys. 70, 1046-1055 (2008).

55.

Veretenenko S, Thejll P, Effects of energetic solar proton events on the cyclone development in the North Atlantic, J. Atmos. Sol.-Terr. Phys. 66, 393-405 (2004).

56.

Zhou J, Tung KK, Solar cycles in 150 years of global sea surface temperature data, J. Clim. 23, 3234-3248 (2010).

57.

Zhou L, Tinsley B, Chu H, Xiao Z, Correlations of global sea surface temperatures with the solar wind speed, J. Atmos. Sol.-Terr. Phys. 149, 232-239 (2016).