PDF
OA 학술지
Latitudinal Distribution of Sunspots Revisited
ABSTRACT
Latitudinal Distribution of Sunspots Revisited
KEYWORD
Sun , Sunspots , data analysis
본문
• ### 1. INTRODUCTION

Sunspots represent one of the most obvious manifestations of magnetic fields on the solar surface. The complexity of the physical models that explain how the solar magnetic flux is generated and works is mainly due to observational constraints set by temporal and spatial distributions of sunspots (Parker 1955, Babcock 1961, Leighton 1969). That is, any potential model should explain the latitude distribution and drifts of sunspots, the tilt angle and the polarity reversal of the bipolar groups, and so on. The cyclic behavior of the solar large-scale magnetic field has also provided invaluable information about the physical processes involved. The sunspot cycle itself has two main aspects: the periodic variation in the number(including polarity) of sunspots and the migration of the appearance position of sunspots.

Sunspot latitude is an indication of the phase of the sunspot cycle: within a cycle, the higher the latitude, the earlier the phase in the cycle. The progressive change in latitude of sunspot groups is best presented by plotting their latitudes against appearance time. In 1904 Maunder published the first version of such a diagram displaying the sunspot group distribution in latitude for each synodic rotation of the Sun (Maunder 1904). That diagram was soon named ‘butterfly diagram’ and became one of the most popular tools for a compact description of the spot zone evolution. The migration of the sunspot activity belt toward the equator is a direct consequence of the passive transport of the toroidal magnetic field by the meridional circulation (Choudhuri et al. 1995, Nandy &Choudhuri 2001).

One thing that must be kept in mind here is that the diagram takes no account of the sunspot lifetime, nor the spatial size. Since all groups are given equal weight,regardless of their temporal and spatial extention, the diagram is dominated by small sunspots, which scatter over wider latitude ranges than larger ones. That makes the butterfly diagram appear noisy and hard to be understood.This is the first reason why we revisit the butterfly diagram in the present paper. By replotting the butterfly diagram by weighting the sunspot with its area, we would like to check general concepts, such as, 'the distance of the sunspot to the equator monotonically decreases as the solar cycle proceeds'. One may also wish to check its validity of some research on the latitudinal migration of solar activity. For example, some researchers have reported that latitudinal migration of sunspots asynchronously occur in the northern and southern hemispheres (Waldmeier 1971, Swinson et al. 1986, Zolotova & Ponyavin 2006, 2007, Donner & Thiel 2007, Zolotova et al. 2009, Li et al. 2009a, b, 2010). They have attempted to explain the solar North-South asymmetry with such a phase difference between paired wings of the butterfly diagram instead of appreciating the real asymmetry in two hemispheres. They have drawn the conclusion by taking a mean latitude from the butterfly diagram without taking due care into account. Furthermore, if such is the case, it is very important since it gives a hint that solar dynamo in two hemispheres is only weakly connected. As a consequence, solar dynamo theory should explain the relative phase shift of the paired wings of the butterfly diagram in the future (Goel & Choudhuri 2009). Until now, magnetic flux transport dynamo models of the Sun’s global magnetic field have been shown to reproduce the amplitude and duration fairly well, but not the relative phase difference in the solar activity.

The North-South asymmetry of activity phenomena in the solar atmosphere also represents an important detail of the solar dynamo action. This is why various solar activity indices, such as, sunspot number, sunspot area, flare index, are studied. Asymmetries between the northern and southern hemispheres have been previously found in various solar indices (Waldmeier 1971, Roy 1977, White & Trotter 1977, Ichimoto et al. 1985, Swinson et al. 1986, Ozguc & Ucer 1987, Verma 1987, Tritakis et al. 1988, Vizoso & Ballester 1989, 1990, Antonucci et al. 1990, Mouradian & Soru-Escaut 1991, Schlamminger 1991, Yi 1992, Carbonell et al. 1993, Verma 1993, Oliver & Ballester 1994, Javaraiah & Gokhale 1997, Li et al. 1998, 2002, Atac & Ozguc 1996, 2001, Temmer et al. 2001, 2002, Krivova & Solanki 2002, Vernova et al. 2002, Berdyugina & Usoskin 2003, Knaack et al. 2004, 2005, Ballester et al. 2005, Gigolashvili et al. 2005a, b, Javaraiah & Ulrich 2006, Zaatri et al. 2006, Chang 2007a, b, 2008, 2009). This is the second motivation for the present study. We would like to see if one can find the North-South asymmetry of the time series of the area-weighted latitudinal variation of sunspots and related periodicities.

This paper begins with descriptions of data and how the center-of-latitude (COL) is defined in Section 2. We present results obtained with the Lomb-Scargle periodo2gram and the wavelet transform in Section 3. Finally, we discuss and conclude in Section 4.

### 2. CENTER-OF-LATITUDE OF SUNSPOT

We have used for the present analysis the Greenwich sunspot group data during the period from 1874 to 1976, and the sunspot group data from the Solar Optical Observing Network (SOON) of the US Air Force (USAF)/ US National Oceanic and Atmospheric Administration (NOAA) during the period of from 1977 to 2009. The daily sunspot area and its latitude have been downloaded from the NASA website1.

In Fig. 1 we show the averaged COL as a function of time during the period from 1874 to 2009. The COL is basically defined by

where Ai is the i-th sunspot’s area and Li is its latitude. Thick and thin curves represent the yearly averaged COL and the monthly averaged COL, respectively. The summation is carried out over a calendar month or a calendar year, instead of the Carrington month. In the first and the second panels, we separately plot the averaged

COL for sunspots appeared in the northern and southern hemispheres, respectively. Positive and negative COL represents the northern and southern latitudes, respectively. In the third panel, the averaged COL for sunspots appeared in both hemispheres is plotted. If the mean latitude where sunspots appear in both hemispheres is symmetrical, COL should be in a form of the random noise. However, we observe that the difference is large every time when the solar cycle begins. What is also barely seen is what is expected when there is a relative phase lag between the two wings. What we see is rather that the wings of the butterfly diagram is not symmetric, nor steadily drifts to the equator. In the last panel we plot the averaged COL of average of two hemispheres defined by

such that the absolute value of the latitude is used rather than the latitude itself. It is shown that the solar cycle begins when COL is high. And as the solar cycle proceeds COL decreases in general. We note, however, that COL is not monotonically decreasing as commonly assumed. Small humps (or short plateaus) between solar minima can be seen. Actually humps appear around every solar maxima.

1http://solarscience.msfc.nasa.gov/greenwch.shtml

### 3. NORTH-SOUTH ASYMMETRY AND PERIODICITIES

In Fig. 2, we show the difference between the north-

rern and southern hemispheres both in sunspot area and in COL. In the upper panel, we plot the difference in the sunspot area between the northern and southern hemispheres, which is averaged over the whole period for a given solar cycle. Therefore, the upper panel shows the asymmetrical behavior of the solar cycle and thereby which one is the dominant hemisphere. The difference is defined by

where Σ AN and Σ AS are thus the whole sunspot area appearing in the northern and southern hemispheres during a given solar cycle, respectively. In the lower panel, on the other hand, we plot the scaled difference of COL between the northern and southern hemispheres, which is also averaged over the period for a given solar cycle. The difference is defined by

where ΣCONN and ΣCOLS are the averaged COL in the northern and southern hemispheres during a given solar cycle as a whole, respectively. We note that COL is positive (negative), when the northern (southern) hemisphere is magnetically dominant. The exception occurs in the cycle 23.

Generally, there is no obvious reason why these two plots look similar, since the first panel tells us how large areas are covered by sunspots for a given hemisphere during the particular solar cycle while the second panel tells us the latitude where large sunspots are distributed for a given hemisphere during the particular solar cycle. For instance, supposed a bigger sunspot is at low latitude in the northern hemisphere and smaller sunspot is at high latitude in the southern hemisphere, it is not impossible that one can observe the positive DA (northern dominance) and negative COL (southern dominance). Therefore, the fact that this kind of ‘odd’ thing does not happen for 11 solar cycles may give a clue that these two phenomena are regulated by a consistent mechanism.

In Fig. 3, we show Lomb-Scargle periodograms (Press & Rybicki 1989) of COL in units of month-1. The power shown here is in an arbitrary unit. In the upper left panel, we show the result of COL in the northern hemisphere. In the upper right one, we show the result of COL variation in the southern hemisphere. Not surprisingly, main peak corresponds to ~11 years. The secondary, yet significant, peak we also see here corresponds to ~5.5 years (confidence level is greater than 90%). We consider this period-

icity is due to humps between the minima we have seen in Fig. 1. In the lower left and right ones, we show results from COL of sunspots appeared in both hemispheres obtained by taking the unsigned COL, and the signed latitude into account, respectively. The power spectrum in the lower left panel is somewhat similar to power spectra in upper panels. In the lower right panel, one might expect a peak corresponding to ~ 9 year periodicity if s/he expects the solar North-South asymmetry (Chang 2007a, b). Instead, the periodicity of ~12 years (frequency ？0.007 month-1) can be found, which appears in the periodogram of the reciprocal of the total sunspot number,

(Ballester et al. 2005). This can be understood since we take the sunspot area as the weight function in COL calculations. It should be noted that the the highest peak corresponds to ~49 year periodicity (frequency ~0.0017 month-1). It can also be understood in terms of a long term variation seen in Fig. 2. A similar periodicity of of ~44 years has been reported in the solar North-South asymmetry by Javaraiah (2003) and Zolotova & Ponyavin (2007), who suggested the existence of a '44-yr' cycle or 'double Hale cycle' and further ~80 year periodicity (Gleissberg 1971).

To see how its phase keeps consistently we analyze the same time series data using the wavelet transform technique (Chang 2006), which shows the power distribution in time. The power is normalized in terms of the mean level of noise so that the plot contains the signal-to-noise ratio by itself. As shown in Fig. 4, the power is apparently maintained for the almost entire period of data set. We have to conclude therefore that the oscillating component of ~49 year periodicity is significant and consistent. Moreover, one may also see another component with ~167 year periodicity (frequency ~0.0005 month-1), which is also found in Javaraiah (2003) where the Sun’s motion about the center of mass of the solar system is related to the long term modulation in the differential solar rotation velocity derived by the sunspot observation.

### 4. DISCUSSION AND CONCLUSIONS

Since the butterfly diagram is introduced to the field of solar study, the accepted paradigm has pretended that spots are scattered around a mean latitude, which monotonically drifts equatorward. For the sunspot activity a second order polynomial curve is widely used to give a fit to the monthly mean latitudes of sunspot activity. Accordingly, most of theoretical calculations have been aimed at developing models predicting the location of such a monotonic latitudinal migration of sunspots. On the other hands, some nonstandard solar magnetic models are introduced to describe non-steady equatorward progression (Bell 1960, Antalova & Gnevyshev 1983, Solanki et al. 2008).

We have examined the latitudinal variation of COL of sunspots in each hemisphere. As a fine structure of the butterfly diagram interests some researchers (Solanki et al. 2008, Ternullo 2010) an issue of how one may take care of different size sunspots will become more important. In this work we have suggested one way of doing this. By taking the area of sunspots into account, one may draw a more reliable conclusion on the evolution and distribution of sunspots, such as, 'active latitude' recently suggested. Another example of use of COL is a new diagnostic of the solar cycle and of models of the generation of the Sun’s magnetic field. As Solanki et al. (2008) demonstrated the latitudinal variation of position of sunspots can be used in distinguishing some dynamo models. On doing so, we may take care of the area of sunspots by using COL.

The summary of what we have found is as follows:(1) It is shown that the solar cycle begins when COL is high and that as the solar cycle proceeds COL decreases in general. COL in both hemispheres is not, however, monotonically decreasing as commonly pictured in theoretical calculations. Small humps (or short plateaus) between every solar minima can be seen around solar maxima.

(2) When the northern (southern) hemisphere is magnetically dominant, COL is positive (negative), except the solar cycle 23. The fact that this kind of ‘odd’ thing does not happen for 11 solar cycles may give a clue that these two phenomena are consistently regulated by a mechanism.(3) The periodicities of ~5.5 and ~11 years are found in each hemisphere. The periodicity of ~49 years are found in the power spectrum of the averaged COL. A possible periodicity of ~167 years is also noticed.

참고문헌
• 1. Antalova A, Gnevyshev MN 1983 Latitudinal distribution of sunspot areas during the period 1874-1976 [CoSka] Vol.11 P.63-93
• 2. Antonucci E, Hoeksema JT, Scherrer PH 1990 Rotation of the photospheric magnetic fields: a North-South asymmetry [ApJ] Vol.360 P.296-304
• 3. Atac T, Ozguc A 1996 North-South asymmetry in the solar flare index [SoPh] Vol.166 P.201-208
• 4. Atac T, Ozguc A 2001 Flare index during the rising phase of solar cycle 23 [SoPh] Vol.198 P.399-407
• 5. Babcock HW 1961 The topology of the Sun’s magnetic field and the 22-year cycle [ApJ] Vol.133 P.572
• 6. Ballester JL, Oliver R, Carbonell M 2005 The periodic behaviour of the North-South asymmetry of sunspot areas revisited [A&A] Vol.431 P.L5-L8
• 7. Bell Bell 1960 On the structure of the Sunspot zone [SCoA] Vol.5 P.17-28
• 8. Berdyugina SV, Usoskin IG 2003 Active longitudes in sunspot activity: century scale persistence [A&A] Vol.405 P.1121-1128
• 9. Carbonell M, Oliver R, Ballester JL 1993 On the asymmetry of solar activity [A&A] Vol.274 P.497
• 10. Chang H-Y 2006 Time/frequency analysis of terrestrial impack crater records [JASS] Vol.23 P.199-208
• 11. Chang H-Y 2007a Variation in North-South asymmetry of sun spot area [JASS] Vol.24 P.91-98
• 12. Chang H-Y 2007b A new method for North-South asymmetry of sun spot area [JASS] Vol.24 P.261-268
• 13. Chang H-Y 2008 Stochastic properties in North South asymmetry of sunspot area [NewA] Vol.13 P.195-201
• 14. Chang H-Y 2009 Periodicity of North South asymmetry of sunspot area revisited: cepstrum analysis [NewA] Vol.14 P.133-138
• 15. Choudhuri AR, Schussler M, Dikpati M 1995 The solar dynamo with meridional circulation [A&A] Vol.303 P.L29
• 16. Donner R, Thiel M 2007 Scale-resolved phase coherence analysis of hemispheric sunspot activity: a new look at the north-south asymmetry [A&A] Vol.475 P.L33-L36
• 17. Gigolashvili MSh, Japaridze DR, Kukhianidze VJ 2005a Variations of the solar differential rotation associated with polarity reversal [SoPh] Vol.231 P.23-28
• 18. Gigolashvili MSh, Japaridze DR, Mdzinarishvili TG, Chargeishvili BB 2005b N-S asymmetry in the solar differential rotation during 1957-1993 [SoPh] Vol.227 P.27-38
• 19. Gleissberg W 1971 The probable behaviour of sunspot cycle 21 [SoPh] Vol.21 P.240-245
• 20. Goel A, Choudhuri AR 2009 The hemispheric asymmetry of solar activity during the last century and the solar dynamo RAA Vol.9 P.115-126
• 21. Ichimoto K, Kubota J, Suzuki M, Tohmura I, Kurokawa H 1985 Periodic behaviour of solar flare activity [Nature] Vol.316 P.422-424
• 22. Javaraiah J 2003 Predictions of strengths of long-term variations in sunspot activity [A&A] Vol.401 P.L9-L12
• 23. Javaraiah J, Gokhale MH 1997 Periodicities in the North-South asymmetry of the solar differential rotation and surface magnetic field [SoPh] Vol.170 P.389-410
• 24. Javaraiah J, Ulrich RK 2006 Solar-cycle-related variations in the solar differential rotation and meridional flow: a comparison [SoPh] Vol.237 P.245-265
• 25. Knaack R, Stenflo JO, Berdyugina SV 2004 Periodic oscillations in the North-South asymmetry of the solar magnetic field [A&A] Vol.418 P.L17-L20
• 26. Knaack R, Stenflo JO, Berdyugina SV 2005 Evolution and rotation of large-scale photospheric magnetic fields of the Sun during cycles 21-23. Periodicities north-south asymmetries and r-mode signatures [A&A] Vol.438 P.1067-1082
• 27. Krivova NA, Solanki SK 2002 The 1.3-year and 156-day periodicities in sunspot data: wavelet analysis suggests a common origin [A&A] Vol.394 P.701-706
• 28. Leighton RB 1969 A Magneto-Kinematic model of the solar cycle [ApJ] Vol.156 P.1
• 29. Li K-J, Gao PX, Zhan LS 2009a Synchronization of hemispheric sunspot activity revisited: wavelet transform analyses [ApJ] Vol.691 P.537-546
• 30. Li K-J, Gao PX, Zhan LS 2009b Synchronization of sunspot numbers and sunspot areas [SoPh] Vol.255 P.289-300
• 31. Li K-J, Liang HF, Feng W 2010 Phase shifts of the paired wings of butterfly diagrams [RAA] Vol.10 P.1177-1185
• 32. Li K-J, Schmieder B, Li Q-Sh 1998 Statistical analysis of the X-ray flares (M >= 1) during the maximum period of solar cycle 22 [A&AS] Vol.131 P.99-104
• 33. Li K-J, Wang JX, Xiong SY, Liang HF, Yun HS 2002 Regularityof the North-South asymmetry of solar activity [A&A] Vol.383 P.648-652
• 34. Maunder EW 1904 Note on the distribution of sun-spots in heliographiclatitude 1874-1902 [MNRAS] Vol.64 P.747-761
• 35. Mouradian Z, Soru-Escaut I 1991 On the dynamics of the large-scale magnetic fields of the sun and the sunspot cycle [A&A] Vol.251 P.649-654
• 36. Nandy D, Choudhuri AR 2001 Toward a mean field formulation of the babcock-leighton type solar dynamo. i. α-coefficient versus Durney’s double-ring approach [ApJ] Vol.551 P.576-585
• 37. Oliver R, Ballester JL 1994 The North-South asymmetry of sunspot areas during solar cycle 22 [SoPh] Vol.152 P.481-485
• 38. Ozguc A, Ucer C 1987 North-South asymmetries in the green corona brightness between 1947 and 1976 [SoPh] Vol.114 P.141-146
• 39. Parker EN 1955 Hydromagnetic dynamo models [ApJ] Vol.122 P.293
• 40. Press WH, Rybicki GB 1989 Fast algorithm for spectral analysis of unevenly sampled data [ApJ] Vol.338 P.277-280
• 41. Roy JR 1977 The North-South distribution of major solar flare events sunspot magnetic classes and sunspot areas 1955-1974 [SoPh] Vol.52 P.53-61
• 42. Schlamminger L 1991 Hemispherical asymmetries in sunspot areas and auroral frequencies [SoPh] Vol.135 P.407-413
• 43. Solanki SK, Wenzler T, Schmitt D 2008 Moments of the latitudinal dependence of the sunspot cycle: a new diagnostic of dynamo models [A&A] Vol.483 P.623-632
• 44. Swinson DB, Koyama H, Saito T 1986 Long-term variations in North-South asymmetry of solar activity [SoPh] Vol.106 P.35-42
• 45. Temmer M, Veronig A, Hanslmeier A 2002 Hemispheric sunspot numbers Rn and Rs: catalogue and N-S asymmetry analysis [A&A] Vol.390 P.707-715
• 46. Temmer M, Veronig A, Hanslmeier A, Otruba W, Messerotti M 2001 Statistical analysis of solar Hα flares [A&A] Vol.375 P.1049-1061
• 47. Ternullo M 2010 The butterfly diagram internal structure [MSAIS] Vol.14 P.202
• 48. Tritakis VP, Petropoulos B, Mavromichalaki H 1988 Asymmetric variations of the coronal green line intensity [SoPh] Vol.115 P.367-384
• 49. Verma VK 1987 On the increase of solar activity in the Southern hemisphere during solar cycle 21 [SoPh] Vol.114 P.185-188
• 50. Verma VK 1993 On the North-South asymmetry of solar activity cycles [ApJ] Vol.403 P.797-800
• 51. Vernova ES, Mursula K, Tyasto MI, Baranov DG 2002 A new pattern for the North-South asymmetry of sunspots [SoPh] Vol.205 P.371-382
• 52. Vizoso G, Ballester JL 1989 Periodicities in the North-South asymmetry of solar activity [SoPh] Vol.119 P.411-414
• 53. Vizoso G, Ballester JL 1990 The North-South asymmetry of sunspots [A&A] Vol.229 P.540-546
• 54. Waldmeier M 1971 The asymmetry of solar activity in the years 1959-1969 [SoPh] Vol.20 P.332-344
• 55. White OR, Trotter DE 1977 Note on the distribution of sunspots between the North and South solar hemispheres and its variation with the solar cycle [ApJS] Vol.33 P.391
• 56. Yi W 1992 The North-South asymmetry of sunspot distribution [JRASC] Vol.86 P.89-98
• 57. Zaatri A, Komm R, Gonza lez-Herna ndez I, Howe R, Corbard T 2006 North South asymmetry of zonal and meridional flows determined from ring diagram analysis of gong ++ data [SoPh] Vol.236 P.227-244
• 58. Zolotova NV, Ponyavin DI 2006 Phase asynchrony of the North-South sunspot activity [A&A] Vol.449 P.L1-L4
• 59. Zolotova NV, Ponyavin DI 2007 Synchronization in sunspot indices in the two hemispheres [SoPh] Vol.243 P.193-203
• 60. Zolotova NV, Ponyavin DI, Marwan N, Kurths J 2009 Long-term asymmetry in the wings of the butterfly diagram [A&A] Vol.503 P.197-201
OAK XML 통계
이미지 / 테이블
• [ Fig.1. ]  The averaged center-of-latitude (COL) as a function of time during the period from 1874 to 2009. Thick and thin curves represent the yearly averaged COL and the monthly averaged COL respectively. Positive and negative COL represents the northern and southern latitudes respectively. In the first and the second panels we separately plot the averaged COL for sunspots appeared in the northern and southern hemispheres respectively. In the third and fourth panels we show fitting results from the distributions of COL of sunspots appeared in both hemispheres with signed and unsigned latitudes respectively.
• [ Fig.2. ]  The difference between the northern and southern hemispheres both in sunspot area and in center-of-latitude (COL). In the upper panel we plot the difference in the sunspot area between the northern and southern hemispheres. In the lower panel we plot the scaled difference of COL between the northern and southern hemispheres.
• [ Fig. 3. ]  Lomb-Scargle periodograms of center-of-latitude (COL) in units of month-1. The power shown here is in an arbitrary unit. In the upper left panel we show the result of COL in the northern hemisphere. In the upper right one we show the result of COL variation in the southern hemisphere. In the lower left and right ones we show results from COL of sunspots appeared in both hemispheres obtained by taking the unsigned COL and the signed latitude into account respectively.
• [ Fig. 4. ]  The power distribution in time obtained by the wavelet transform technique. The power is apparently maintained for the almost entire period of data set. One may find peaks corresponding to periodicities of ~49 and ~167 years.
(우)06579 서울시 서초구 반포대로 201(반포동)
Tel. 02-537-6389 ｜ Fax. 02-590-0571 ｜ 문의 : oak2014@korea.kr