Research Paper

Normalized Cross-Correlations of Solar Cycle and Physical Characteristics of Cloud

Heon-Young Chang1,2,http://orcid.org/0000-0003-2015-2725
Author Information & Copyright
1Department 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: Sep 23, 2019; Revised: Nov 25, 2019; Accepted: Nov 26, 2019

Abstract

We explore the associations between the total sunspot area, solar north-south asymmetry, and Southern Oscillation Index and the physical characteristics of clouds by calculating normalized cross-correlations, motivated by the idea that the galactic cosmic ray influx modulated by solar activity may cause changes in cloud coverage, and in turn the Earth’s climate. Unlike previous studies based on the relative difference, we have employed cloud data as a whole time-series without detrending. We found that the coverage of high-level and low-level cloud is at a maximum when the solar north-south asymmetry is close to the minimum, and one or two years after the solar north-south asymmetry is at a maximum, respectively. The global surface air temperature is at a maximum five years after the solar north-south asymmetry is at a maximum, and the optical depth is at a minimum when the solar north-south asymmetry is at a maximum. We also found that during the descending period of solar activity, the coverage of low-level cloud is at a maximum, and global surface air temperature and cloud optical depth are at a minimum, and that the total column water vapor is at a maximum one or two years after the solar maximum.

Keywords: solar activity; terrestrial climate; data analysis

1. INTRODUCTION

The global temperature rise of recent decades is believed to have been driven mainly by the release of carbon dioxide into the Earth’s atmosphere by the burning of fossil fuels (e.g., Solomon et al. 2007). Determining an exact level of warming due to man-made greenhouse gases requires an unequivocal understanding of the natural causes of climate change. Provided that the Sun is the fundamental energy source in the Earth’s thermal energy budget, and that the radiative equilibrium temperature of the Earth is determined by a balance between the solar irradiation and the Earth’s energy loss by outgoing long-wave radiation, sun-climate connections as possible forces of nature that cause the global warming have been of general interest (Ney 1959;Dickinson 1975;Tinsley 2000;Roldugin & Tinsley 2004;Scafetta & West 2006;Burns et al. 2007;Haigh 2007;Burns et al. 2008;Cho & Chang 2008;Gray et al. 2010;Cho et al. 2012;Lee & Yi 2018;Muraki 2018). For instance, Eddy (1976) who showed the practical absence of sunspot activity during the period from 1645 to 1715, which is now known as the Maunder minimum (Spörer 1887;Maunder 1904), has claimed that this period coincides with an era when extremely cold weather prevailed in Western Europe including Central England. Thus, any variability in solar insolation has the potential to affect the terrestrial climate (Haigh 1996, 2007;Solanki et al. 2002;Solanki et al. 2013). It has been shown, however, that space-born radiometers providing accurate and stable measurements suggest that the value of the total solar irradiance (TSI) on the top of the Earth’s atmosphere, defined as the total power from the Sun impinging on a unit area perpendicular to the Sun’s rays, is 1,360.8 ± 0.5 Wm–2 under solar minimum conditions (Kopp et al. 2012), and it changes by less than ~0.1 % over the time scale of the solar cycle, which is clearly insufficient to account for the rise in global temperature. This small amount, however, may be sufficient to cause a significant climate effect if it is present over a long period of time or if there are nonlinear responses that amplify feedback.

On the basis of the fact that such variability is insufficient to influence Earth’s climate, nonetheless, answering the question of solar contribution to variation in the terrestrial climate seems premature. This is because the magnitude of variations in solar radiation is a strong function of wavelength. Changes in the ultraviolet (UV) range of the solar spectrum with the highest amplitude influence the structure of the middle atmosphere by affecting stratospheric chemistry through the modification of the photochemical dissociation rates of ozone (Emmert & Picone 2010;Gray et al. 2017). The resulting warming of the lower stratosphere via absorption of solar UV radiation by atmospheric ozone and other species of gas molecules may propagate its effects dynamically downwards to the troposphere, and cause both poleward shifts in position of the subtropical westerly jets, and poleward shifts of the tropical cyclone tracks, as those in El Niño periods (Kim et al. 2017;Kim et al. 2018).

The Sun may also influence the terrestrial climate in other indirect ways (Tinsley & Deen 1991;Pudovkin et al. 1997;Pudovkin 2004;Roldugin & Tinsley 2004;Haigh 2007;Bazilevskaya et al. 2008;Svensmark et al. 2009). An eminent example is by affecting the Earth’s climate by modulating the magnitude of influx of high energy galactic cosmic rays (GCRs), which has been known to change according to the magnetic characteristics of the heliosphere (Forbush 1954;Kolomeets et al. 1973;Kane 2005;Usoskin 2013). That is, incoming GCRs interact with atmospheric molecules, and subsequently induce important results in atmospheric ionization and condensation processes. One of the most important implications is modifying the coverage and physical properties of clouds, and thus the Earth’s radiative balance (Svensmark & Friis-Christensen 1997;Marsh & Svensmark 2000;Pallé Bagó & Butler 2000;Todd & Kniveton 2004;Harrison & Stephenson 2006;Pierce & Adams 2009;Singh & Singh 2010;Rawal et al. 2013;Svensmark et al. 2017). Changes in the total albedo via a GCRs-cloud interaction, therefore, may help in explaining how a relatively small change in TSI can provoke changes in Earth’s climate. There is further evidence that GCR influx is phenomenologically related to tropical cyclone activity (Cohen & Sweetser 1975;Kavlakov 2005;Elsner & Jagger 2008;Pérez-Peraza et al. 2008;Hodges & Elsner 2011;Oey & Chou 2016). Increases in GCR influx penetrating the Earth’s atmosphere during solar minima contribute to an increase in cloudiness, and therefore tropical cyclone intensification, via increases in the convective available potential energy (CAPE).

It may be helpful to briefly mention a process associated with GCRs that causes a modification in cloud microphysics, because clouds play a fundamental role in the Earth’s radiative energy equilibrium. The occurrence of clouds can be predicted by meteorological parameters, such as, humidity and temperature. In condensing water vapor to a cloud droplet, physical characteristics of clouds, including reflectivity, lifetime, coverage, size distribution, are controlled by cloud condensation nuclei (CCN) of 50– 100 nm through a series of subtle and complex processes. Two ways by which GCRs may affect the formation of cloud droplets are by the ion-aerosol clear-air mechanism and the ion-aerosol near-cloud mechanism (see Carslaw et al. 2002 and reference therein). The ion-aerosol clearair mechanism is based on the idea that the presence of ions enhanced by electric attraction force the birth and early growth of ultra-fine condensation nuclei (UCN) in the upper troposphere, which then become condensation nuclei (CN) that may eventually become CCN. The ionaerosol near-cloud mechanism is less well understood than the former. This mechanism is based on the fact that the aerosol electrical charge is different near clouds than it is in clear air. Highly charged aerosols and the electric field at cloud boundaries, the development of which is modulated by GCRs, may migrate within cloud and enhance aerosol efficacy as ice-forming nuclei. The CLOUD (Cosmics Leaving Outdoor Droplets) experiment conceived at CERN in the late 1990s was created to systematically test the link between GCRs and CCN (Kirkby 2001;Dunne et al. 2016;Gordon et al. 2017). The CERN CLOUD experiments are the most comprehensive laboratory measurements of aerosol nucleation rates so far achieved. It should be pointed out, however, that according to current measurements and simulations, the hypothesized effect is too small to account for the observed climate change due to the slow growth of UCN, though unknown mechanism may be operating to elevate the growth of UCN (Pierce & Adams 2007, 2009;Snow-Kropla et al. 2011;Dunne et al. 2012;Kazil et al. 2012;Yu et al. 2012;Yu & Luo 2014).

In this study, we explore associations between the total sunspot area and characteristics of clouds by calculating normalized cross-correlations. We note that all the phenomenological studies on the GCR-cloud relation are focused on the relative difference because the correlation between changes in low cloud fraction and changes in GCR influx was first reported by Svensmark & Friis-Christensen (1997). The total amount of low-lying cloud must be considered to calculate the albedo. However, the relative change is useful only when a background variation is trivial or the trend of the relative change is the same as that of the background variation that is to be detrended for the observed relative change. Unlike as in previous studies, we have directly employed cloud coverage with respect to the altitude ranges as a whole without detrending along with some physical properties of clouds, such as, optical depth and top temperature. We also investigated associations with the solar north-south asymmetry, having been motivated by the fact that the mean global temperature anomaly is smaller when the solar northern hemisphere is active than when the southern hemisphere is active (Cho et al. 2012). It is well known that the magnitude of GCR influx depends on the solar north-south asymmetry (e.g., Cho et al. 2011). To compare the level of possible contributions, we repeat the same analysis with the Southern Oscillation Index (SOI), which is a measure of the El Niño-Southern Oscillation (ENSO) global climate regulating phenomena. This paper is organized as follows. We begin with descriptions of the data analyzed for the present study in Section 2. We present and discuss results in Section 3. Finally, we briefly summarize and present conclusions in Section 4.

2. DATA

For a proxy of solar variability, we used sunspot area data from the NASA Marshall Space Flight Center website (http://solarscience.msfc.nasa.gov/greenwch.shtml). We downloaded the daily sunspot area in units of millionths of a hemisphere separately for the Northern Hemisphere, AN, and the Southern Hemisphere, AS, during the period from 1983 to 2009, approximately spanning solar cycles 22 and 23. The solar north-south asymmetry is calculated by the definition as (ASAN) / (AS + AN). For the time-series of the SOI index, we have extracted daily data from the National Oceanic and Atmospheric Administration (NOAA) website (https://www.ncdc.noaa.gov/teleconnections), which provides various teleconnection indices representing large-scale anomalies of the pressure system in space and time that influence atmospheric circulation (Kim & Chang 2019). Having downloaded all the daily data, we averaged those data yearly so that we obtained in the normalized cross-correlations in yearly values. In Fig. 1, we show the yearly total sunspot area and the yearly solar north-south asymmetry in the top and middle panels, respectively. In the bottom panel, we show the yearly SOI index.

jass-36-4-225_F1
Fig. 1. Yearly total sunspot area, yearly solar north-south asymmetry, and yearly SOI index from top to bottom panels, respectively. SOI, Southern Oscillation Index.
Download Original Figure

The International Satellite Cloud Climatology Project (ISCCP) began in 1982 as a part of the World Climate Research Program (WCRP) to compile weather satellite radiance measurements, and it developed a climatology of cloud radiative properties (Schiffer & Rossow 1983). The legacy dataset known as ISCCP D series ended in 2009 and was upgraded to the H series. The resolution of the D2 product is 280 km in equal-area grid. Taking advantage of the classification of cloud type in the latest ISCCP D2 product, we downloaded cloud quantities in percentages from the ISCCP website (https://isccp.giss.nasa.gov/ products/browsed2.html). In Fig. 2, we show the yearly coverage of cloud in terms of total cloud, high cloud, middle cloud, and low cloud for the whole globe, and cloud in tropical regions (15°N–15°S) during the period from 1983 to 2009, from top to bottom panels, respectively. Cloud over fraction (%) represents the fractional area covered by clouds as observed by satellites. The altitude ranges covered by cloud are: greater than 680 mbar (0–3.2 km) for low cloud, 440–680 mbar (3.2–6.5 km) for mid-level cloud, 50–440 mbar (6.5–16 km) for high cloud. In Fig. 3, we show the total column water vapor in cm used to calculate the column mass density of the cloud, the top pressure in mbar representing the location of the radiating top of the cloud, the top temperature in K from which the top pressure can be calculated, the cloud optical depth representing the optical thickness of the cloud at visible wavelengths (approximately 0.6 μm), the global surface air temperature in K during the period from 1983 to 2009, from the top to the bottom panel, respectively.

jass-36-4-225_F2
Fig. 2. Yearly coverage of clouds in terms of total cloud, high cloud, middle cloud, low cloud, and cloud in tropical regions (15°N–15°S) from top to bottom panels, respectively.
Download Original Figure
jass-36-4-225_F3
Fig. 3. Total column water vapor of cloud, cloud top pressure, cloud top temperature, cloud optical depth, and global surface air temperature from top to bottom panels, respectively.
Download Original Figure

3. NORMALIZED CROSS-CORRELATION

As a measure of similarity between two physical quantities, one may calculate the normalized cross-correlation which is essentially similar to a convolution in nature. That is, the normalized cross-correlation is defined by

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

where

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

and

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

where τ is the displacement and u* (τ) is the complex conjugate of u (τ). For all t, |z(t)|1. By this normalization, the cross-correlation becomes a time-dependent Pearson correlation coefficient, that is, z(0) becomes the Pearson correlation coefficient. Note that the argument of the function at the extremum is called the lag if it is positive, or the lead if it is negative. In other words, 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.

In Fig. 4, we show the normalized cross-correlation resulting from the coverage of total cloud, cloud in the tropics, high-level cloud, medium-level cloud, and lowlevel cloud, in columns from left to right, respectively. In the top row, the normalized cross-correlations are calculated between the solar north-south asymmetry and those parameters mentioned just above, respectively. In the middle row, the normalized cross-correlations are calculated between the total sunspot area and those shown in Fig. 2, respectively. In the last row, for comparison, the normalized cross-correlations are calculated between the SOI index and those data shown in Fig. 2, respectively. For the solar north-south asymmetry, the coverage of highlevel cloud shown in panel (c) is positively correlated with approximately +7 years and the cloud coverage of low-level cloud in panel (e) is marginally correlated with approximately –2 years. Hence, the coverage of highlevel cloud is at a maximum when the solar north-south asymmetry is approximately at a minimum, or when more sunspots appear in the northern hemisphere of the Sun, in the sense that the periodicity of the solar north-south asymmetry is ~12 years. The coverage of low-level cloud is maximum around one or two years after the solar northsouth asymmetry is at a maximum, or when more sunspots appear in the southern hemisphere of the Sun. For solar activity based on the total sunspot area, only the coverage of low level-cloud shown in panel (j) is positively correlated with approximately +3–4 years. That is, the coverage of lowlevel cloud is at a maximum during the descending phase of solar activity. For the SOI index shown for comparison, surprisingly, any cloud coverage is unlikely to be correlated (e.g., see Wang & Su 2015).

jass-36-4-225_F4
Fig. 4. Normalized cross-correlation resulting from the coverage of total cloud, cloud in the tropics, high-level cloud, medium-level cloud, and low-level cloud in columns from left to right, respectively. In the top, middle, and bottom rows, the normalized cross-correlations are calculated using the solar north-south asymmetry, total sunspot are, and the SOI index, respectively. SOI, Southern Oscillation Index.
Download Original Figure

In Fig. 5, we show the normalized cross-correlation resulting from the total column water vapor, global surface air temperature, top temperature, top pressure, and cloud optical depth in columns from left to right, respectively. From top to bottom, as in Fig. 4, the normalized crosscorrelations are calculated using the solar north-south asymmetry, the total sunspot area, and the SOI index, respectively. For the solar north-south asymmetry, the global surface air temperature shown in panel (b) is positively correlated with approximately –5 years and the cloud optical depth in panel (e) is marginally anticorrelated with approximately –0.5 years. That is, the global surface air temperature reaches a maximum five years after the solar north-south asymmetry is at a maximum, as Cho et al. (2012) claimed earlier. The optical depth is at a minimum around when the solar north-south asymmetry is at a maximum. For solar activity based on the total sunspot area, the total column water vapor shown in panel (f ) is correlated with approximately 3 years, the global surface air temperature in panel (g) is anticorrelated with approximately 4 years and the cloud optical depth in panel (j) is marginally anticorrelated with approximately 4 years. Hence, the total column water vapor is at a maximum about 3 years after the solar maximum. This finding may correspond to what Svensmark et al. (2009) found, i.e., low clouds contain less liquid water following the Forbush decrease. Furthermore, the global surface air temperature and the cloud optical depth are at a minimum during the descending phase of solar activity, rather than exactly during the solar minimum. Any of the physical properties of cloud are unlikely to be correlated with the SOI index shown for comparison.

jass-36-4-225_F5
Fig. 5. Normalized cross-correlation resulting from the total column water vapor, global surface air temperature, top temperature, top pressure, and cloud optical depth in columns from left to right, respectively. In the top, middle, and bottom rows, the normalized crosscorrelations are calculated using the solar north-south asymmetry, total sunspot are, and the SOI index, respectively. SOI, Southern Oscillation Index.
Download Original Figure

4. CONCLUSIONS

GCRs modulated by solar activity are one of principal agents of solar variability in the changes to Earth’s climate. Motivated by the fact that GCR influx may modulate the coverage of low-level cloud by involving cloud processes such as forming CCN, and in turn may modulate the Earth’s radiation budget, we investigated associations between solar variabilities and the characteristics of clouds using normalized cross-correlations. We also explored the influence of the solar north-south asymmetry, comparing it with that of the SOI index. Unlike previous studies, we employed cloud data extracted from the ISCCP website as a whole time-series without detrending. We believe this to be crucial because the relative difference only yields a meaningful result when the relative change is similar to the overall trend, which is clearly not the case as shown in Figs. 2 and 3. The properties of cloud, including the cloud coverage, vary primarily by the distribution of cloud over the continents and oceans. The next most important agents of change are seasons and climatic oscillations, such as, ENSO, the North Atlantic Oscillation (NAO).

In this paper, we calculated the normalized crosscorrelations of the total sunspot area, the solar northsouth asymmetry, and the SOI index with cloud coverage at different altitude ranges, total column water vapor in cloud, global surface air temperature, cloud top temperature, cloud top pressure, and cloud optical depth. As a result, our main findings are:

  • (1) The coverage of high-level cloud is at a maximum when the solar north-south asymmetry is approximately at a minimum. The coverage of low-level cloud is at a maximum around one or two years after the solar northsouth asymmetry is at a maximum. The global surface air temperature reaches a maximum 5 years after the solar north-south asymmetry is at a maximum. The optical depth is at a minimum when the solar north-south asymmetry is at a maximum.

  • (2) The coverage of low-level cloud is at a maximum during the descending period of solar activity. The total column water vapor is at a maximum one or two years after the solar maximum. The global surface air temperature and the cloud optical depth are at a minimum during the descending phase of solar activity.

  • (3) Surprisingly, the SOI index, shown for comparison, is unlikely to be correlated with any cloud coverage, or any of the physical properties of cloud. Thus, our findings can be considered significant.

ACKNOWLEDGMENTS

The authors thank the anonymous referees for 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.

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).

2.

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).

3.

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).

4.

Carslaw KS, Harrison RG, Kirkby J, Cosmic rays, clouds, and climate, Science 298, 1732-1737 (2002).

5.

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

6.

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).

7.

Cho IH, Kwak YS, Chang HY, Cho KS, Park YD, et al., Dependence of GCRs influx on the solar North-South asymmetry, J. Atmos. Sol.-Terr. Phys. 73, 1723-1726 (2011).

8.

Cohen TJ, Sweetser EI, The spectra of the solar cycle and of data for Atlantic tropical cyclones, Nature 256, 295-296 (1975).

9.

Dickinson RE, Solar variability and the lower atmosphere, Bull. Am. Meteorol. Soc. 56, 1240-1248 (1975).

10.

Dunne EM, Gorden H, Kurten A, Almeida J, Duplissy J, et al., Global atmospheric particle formation from CERN CLOUD measurements, Science 354, 1119-1124 (2016).

11.

Dunne EM, Lee LA, Reddington CL, Carslaw KS, No statistically significant effect of a short-term decrease in the nucleation rate on atmospheric aerosols, Atmos. Chem. Phys. 12, 11573- 11587 (2012).

12.

Eddy JA, The maunder minimum, Science 192, 1189-1202 (1976).

13.

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

14.

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

15.

Forbush SE, World-wide cosmic ray variations, 1937–1952, J. Geophys. Res. 59, 525-542 (1954).

16.

Gordon H, Kirkby J, Baltensperger U, Bianchi F, Breitenlechner M, et al., Causes and importance of new particle formation in the present-day and preindustrial atmospheres, J. Geophys. Res. Atmos. 122, 8739-8760 (2017).

17.

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

18.

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

19.

Haigh JD, The impact of solar variability on climate, Science 272, 981-984 (1996).

20.

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

21.

Harrison RG, Stephenson DB, Empirical evidence for a nonlinear effect of galactic cosmic rays on clouds, Proc. R. Soc. A, 462, 1221-1223 (2006).

22.

Hodges RE, Elsner JB, Evidence linking solar variability with US hurricanes, Int. J. Climatol. 31, 1897-1907 (2011).

23.

Kane RP, Short-term periodicities in solar indices, Sol. Phys. 227, 155-175 (2005).

24.

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-6701 (2005).

25.

Kazil J, Zhang K, Stier P, Feichter J, Lohmann U, O’Brien K, The present-day decadal solar cycle modulation of Earth’s radiative forcing via charged H2SO4/H2O aerosol nucleation, Geophys. Res. Lett. 39, L02805 (2012).

26.

Kim JH, Chang HY, Association between solar variability and teleconnection index, J. Astron. Space Sci. 36, 149-157 (2019).

27.

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

28.

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).

29.

Kirkby J, Beam measurements of a CLOUD (Cosmics Leaving OUtdoor Droplets) chamber, European Organization for Nuclear Research (CERN), CERNOPEN-2001-028 (2001).

30.

Kolomeets EV, Mukanov JB, Shvartsman YE, 1973. Long-term variations in cosmic rays and solar activity, Proceedings of the 13th International Conference on Cosmic Rays, Denver, CO, 17-30 Aug 1973.

31.

Kopp G, Fehlmann A, Finsterle W, Harber D, Heuerman K, et al., Total solar irradiance data record accuracy and consistency improvements, Metrologia, 49, S29-S33 (2012).

32.

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).

33.

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

34.

Maunder EW, Note on the distribution of sun-spots in heliographic latitude, 1874 to 1902, Mon. Not. R. Astron. Soc. 64, 747-761 (1904).

35.

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

36.

Ney EP, Cosmic radiation and the weather, Nature 183, 451-452 (1959).

37.

Oey LY, Chou S, Evidence of rising and poleward shift of storm surge in western North Pacific in recent decades, J. Geophys. Res. Oceans, 121, 5181-5192 (2016).

38.

Pallé Bagó E, Butler CJ, The influence of cosmic rays on terrestrial clouds and global warming, Astron. Geophys. 41, 4.18-4.22 (2000).

39.

Pérez-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).

40.

Pierce JR, Adams PJ, Efficiency of cloud condensation nuclei formation from ultrafine particles, Atmos. Chem. Phys. 7, 1367- 1379 (2007).

41.

Pierce JR, Adams PJ, Can cosmic rays affect cloud condensation nuclei by altering new particle formation rates?, Geophys. Res. Lett. 36, L09820 (2009).

42.

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

43.

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).

44.

Rawal A, Tripathi SN, Michael M, Srivastava AK, Harrison RG, Quantifying the importance of galactic cosmic rays in cloud microphysical processes. J. Atmos. Sol.-Terr. Phys. 102, 243- 251 (2013).

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.

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

47.

Schiffer RA, Rossow WB, The International Satellite Cloud Climatology Project (ISCCP): the first project of the World Climate Research Programme, Bull. Am. Meteorol. Soc., 64, 779-784 (1983).

48.

Singh D, Singh RP, The role of cosmic rays in the Earth’s atmospheric processes, Pramana. J. Phys. 74, 153-168 (2010).

49.

Snow-Kropla EJ, Pierce JR, Westervelt DM, Trivitayanurak W, Cosmic rays, aerosol formation and cloud-condensation nuclei: sensitivities to model uncertainties, Atmos. Chem. Phys. 11, 4001-4013 (2011).

50.

Solanki SK, Krivova NA, Haigh JD, Solar irradiance variability and climate, Ann. Rev. Astron. Astrophys. 51, 311-351 (2013).

51.

Solanki SK, Schüssler M, Fligge M, Secular variation of the Sun’s magnetic flux, Astron. Astrophys. 383, 706-712 (2002).

52.

Solomon S, Qin D, Manning M, Chen Z, Marquis M, et al., 2007. Climate Change 2007: The Physical Science Basis. Contribution of Working Group i to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. (Cambridge University Press, Cambridge, UK).

53.

Spörer FWG, Vierteljahrsschr. Astron. Ges. Leipzig 22, 323 (1887).

54.

Svensmark H, Cosmoclimatology: A new theory emerges, Astron. Geophys. 48, 1.18-1.24 (2007).

55.

Svensmark H, Bondo T, Svensmark J, Cosmic ray decreases affect atmospheric aerosols and clouds, Geophys. Res. Lett. 36, L15101 (2009).

56.

Svensmark H, Enghoff MB, Shaviv NJ, Svensmark J, Increased ionization supports growth of aerosols into cloud condensation nuclei, Nat. Comm. 8, 2199 (2017).

57.

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).

58.

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).

59.

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).

60.

Todd MC, Kniveton DR, Short-term variability in satellite-derived cloud cover and galactic cosmic rays: an update, J. Atmos. Sol.-Terr. Phys. 66, 1205-1211 (2004).

61.

Usoskin IG, A history of solar activity over millennia, Living Rev. Sol. Phys. 10, 1 (2013).

62.

Wang H, Su W, The ENSO effects on tropical clouds and top-of-atmosphere cloud radiative effects in CMIP5 models, J. Geophys. Res. Atmos. 120, 4443-4465 (2015).

63.

Yu F, Luo G, Effect of solar variations on particle formation and cloud condensation nuclei, Environ. Res. Lett. 9, 045004 (2014).

64.

Yu F, Luo G, Liu X, Easter RC, Ma X, et al., Indirect radiative forcing by ion-mediated nucleation of aerosol, Atmos. Chem. Phys. 12, 11451-11463 (2012).