The latitudinal variation of sunspots appearing during the period from 1874 to 2009 has been studied in terms of centerof-latitude (COL). The butterfly diagram has been used to study the evolution of the magnetic field and the dynamics at the bottom of the solar convection zone. Short-term periodicities have been of particular interest, in that they are somehow related to the structure and dynamics of the solar interior. We thus have focused our investigation on shortterm periodicities. We first calculated COL by averaging the latitude of sunspots with the weight function in area. Then,we analyzed the time series of COL using the wavelet transform technique. We found that a periodicity of ~5 years is the most dominant feature in the time series of COL, with the exception of the ~11 year solar cycle itself. This periodicity can be easily understood by considering small humps between the minima in the area-weighted butterfly diagram.However, we find that periodicities of ~1.3 (0.064), ~1.5 (0.056), or ~1.8 (0.046) years ( 1 / month ), which have been previously suggested as evidence of links between the changing structure of the sunspot zone and the tachocline rotation rate oscillations,are insignificant and inconsistent. We therefore conclude that the only existing short-term periodicity is of ~5 years, and that periodicities of ~1.3, ~1.5, or ~1.8 years are likely to be artifacts due to random noise of small sunspots.
1. INTRODUCTION
The solar magnetic cycle can be studied in terms of the latitudinal variation of the location of sunspots with time(Antalova & Gnevyshev 1965, 1983, Pelt et al. 2000, Temmer et al. 2002, 2006, Zolotova & Ponyavin 2006, 2007,Donner & Thiel 2007, Chang 2007, 2008, 2009, Ternullo 2007a, 2010a, Solanki et al. 2008, Goel & Choudhuri 2009,Li et al. 2008, 2009a, b, 2010a, b, Zolotova et al. 2009, Li 2010). For instance, Tuominen (1952) noticed the existence of a difference in the meridional motions of spot groups during the odd and the even cycles, and suggested the existence of a 22-year periodicity in the solar meridional flow. Pulkkinen et al. (1999) also studied the long-term variation of sunspot latitudes, and reported a period of about 90 years. Moreover, they demonstrated an idea that the solar North-South asymmetry phenomenon could be explained by the oscillation of a quadrupolar component of a solar-mixed parity mode magnetic field. In this way, the butterfly diagram may be used to study the solar activity, leading even to constraining a solar dynamo model (Waldmeier 1957, 1971, Li et al. 2002, Solanki et al. 2008).
Since sunspots are formed soon after the emergence of an active region and generally die long before the active region itself, they can be considered as tracers of freshly emerged magnetic flux. Thus, under the assumption that the sunspot butterfly diagram reflects the evolution of the magnetic field and the dynamics at the bottom of the solar convection zone, which has been shown as reasonable (Hathaway et al. 2003), the search for such a shortterm periodicity is intriguing, and calls for more careful analyses.
Since sunspots are formed soon after the emergence of an active region and generally die long before the active region itself, they can be considered as tracers of freshly emerged magnetic flux. Thus, under the assumption that the sunspot butterfly diagram reflects the evolution of the magnetic field and the dynamics at the bottom of the solar convection zone, which has been shown as reasonable(Hathaway et al. 2003), the search for such a shortterm periodicity is intriguing, and calls for more careful analyses.
Recently, there have been some efforts to search for short-term periodicities in the latitudinal distribution of sunspots (Krivova & Solanki 2002, Knaack et al. 2005) in addition to long-term periodicities, such as ~11-year (Schwabe 1843), ~44-year (Javaraiah 2003), ~53-year (Le & Wang 2003), ~78-year (Wolf 1862), ~80-90-year (Gleissberg 1971), and ~65-130 year quasi-periodic cycles (Nagovitsyn 1997). Interestingly, these short-term periodicities have been related to the structure and dynamics of the solar interior, such as oscillations at ~1.3 years observed in the rotation rate at the base of the convection zone (Howe et al. 2000), or a torsional wave pattern with alternating latitude zones of slow and fast rotation (Howard & LaBonte 1980). For instance, Krivova & Solanki (2002) attempt to demonstrate that the periodicity of the tachocline rotation rate oscillations can actually be traced in the sunspot activity.
In this paper, we take the area of sunspots into account in the latitudinal solar-cyclic behavior of sunspots (Cho & Chang 2011). This consideration is required, as it reduces the noise level of a power spectrum in the sense that small sunspots prevail in the butterfly diagram, which scatter over wider latitude ranges than larger ones. The butterfly diagram is quite noisy because all sunspots are given equal weight in the butterfly diagram, regardless of their temporal and spatial extension. For example, 65% of groups (the smallest ones) contribute to only ≈10% of the total spotted area, as was pointed out by Ternullo (2007b). Hence, periodical tests are apt to be biased by randomly scattered small sunspots unless the area of sunspots is properly managed. So far, this simple fact that the butterfly diagram appears to be noisy due to small sunspots has been ignored. For instance, Ternullo (2010a) has claimed that there is a fine structure in the butterfly diagram that evolves over time, and that the averaged latitude of sunspots apparently drifts alternatingly equatorward and poleward for short time intervals, i.e. ~2 years. He came to this conclusion through a visual inspection, without any periodogram analysis. More seriously,he averaged the latitudes of sunspots with a moving window, whose width is a decisive factor to the result, without taking care of the area of sunspots. We review the recently reported short-term periodicities in the butterfly diagram by replotting it by weighting the sunspot with its area (Krivova & Solanki 2002, Ternullo 2007a, 2010a,b). By doing so, we may or may not confirm the short-term periodicities recently reported.
This paper begins with descriptions of how the center-of-latitude (COL) is defined, and of the wavelet analysis in Section 2. We present and discuss the results obtained with the technique of wavelet transform in Section 3. Finally, we summarize and conclude in Section 4.
2. CENTER-OF-LATITUDE OF SUNSPOT AND WAVELET ANALYSIS
The data used for analysis were the daily sunspot area and its latitude obtained from the National Aeronautics and Space Administration (NASA) website
1
. The data set we analyze in the current analysis spans from 1874 to 2009. Here, we consider the distribution of the COL. We collectively compute COL for sunspots appearing in both hemispheres by taking the absolute value of the latitude, which is given by
where
Ai
is the
i
-th sunspot's area and
Li
is its latitude. We have carried out the summation over a period of one calendar month. It should be noted that the collective COL is not identical to the average of COLs, due to the fact that some sunspots separately appear in only one of the northern or southern hemispheres.
In
Fig. 1
, we show the averaged COL as a function of time. The yearly averaged COL is also shown with a thick curve for comparison. Note, however, that the thick curve is given simply for clarity, and the thin curve alone is used in the analysis. The solar cycle begins when COL is high, and COL gradually decreases in general as the solar cycle proceeds. We note, however, that small humps (or short plateaus) appear around almost all solar maxima. That
The averaged center-of-latitude (COL) as a function of time from 1874 to 2009. Thick and thin curves are yearly and monthly averaged COL respectively. Note however that the thick curve is given for a simple clarity and the thin curve alone is used in the analysis.
is, COL is not monotonically decreasing. A similar conclusion was drawn by Antalova & Gnevyshev (1983), who stated that the 11-year cycle consists of the superposition of two maxima which peak at different times and different latitudes, and that the second peak occurs at a lower latitude than the first. Recently, Ternullo (2010b) has also shown that the sunspot has the tendency to preferentially appear in some photospheric regions tightly limited in latitude,givinga 'leopard skin' appearance to the butterfly diagram. This is in contrast to the traditional view that the average latitude of sunspots in the butterfly diagram steadily drifts equatorward. One may be interested in further pursuing this non-monotonic behavior to relate the latitudinal preference to postulating 'active latitude' suggested by Solanki et al. (2008).
In the present analysis, we employ the wavelet analysis technique associated with the Gabor transform. This method is sometimes referred to as time/frequency analysis. Compared to the Fourier transform
FT
, which is defined by
the Gabor transform
GT
is effectively defined as
where
S(t)
is the signal as a function of time
t, ν
is the cyclic frequency, and
G(t)
is the window function. In the above form, the wavelet transform WT can be considered as a transform of a signal derived from an original wavelet by dilation in time, i.e., as
where
d
is a dilation coefficient. Put simply, the wavelet transform is the correlation of the signal
S(t)
and a running wavelet
W(t)
.
Several wavelet-forms can be selected, but in practice the most popular choice is the Morlet wavelet (Grossmann & Morlet 1984). It is a sinusoidal wave with a Gaussian window:
where
d
is a dilation coefficient. Thus, it will provide information around time
t
0
and around frequency
ν
=1/
d
. Its squared modulus is proportional to the power of the signal.This information is displayed in a two-dimensional plot, showing variations of the power with time and frequency. The shape of the window is an important parameter, asit will allow optimization in time and in frequency of the information extracted. In this paper, we adopt the Morlet wavelet transform throughout the analysis.
1http://solarscience.msfc.nasa.gov/greenwch.shtml
3. RESULTS AND DISCUSSION
In
Fig. 2
, as an example, we show a 2-dimensional power spectrum obtained by the wavelet analysis to demonstrate its robustness. As mentioned in Section 2, the power spectrum is shown in a 2-dimensional plot as a function of time in month and frequency in 1/month. The data set covers 1628 months, from May in 1874 to December in 2009. The power is normalized with noise, and shown in an arbitrary unit. Not surprisingly, the main peak in the whole power spectrum occurs where the frequency is ~0.00781/month, whose period corresponds to ~10.68 years. In
Fig. 3
, we show a part of the power spectrum where the second most conspicuous peak in the whole power spectrum can be found. The power is consistently concentrated at frequency ~0.016 1/month, which corresponds to ~5.21 years. Even though the peak can be understood in terms of a subharmonic of the main peak, we consider this periodicity to be due to the humps between the minima we have seen in
Fig. 1
. It can be considered as evidence that there are two populations of sunspots,as will be shown in a subsequent paper.
As stated above, we are only interested in short-term periodicities, and focus on these in the present analysis.In
Fig. 4
, we show the power spectrum obtained by
The 2-dimensional power spectrum obtained by the wavelet analysis. The power spectrum is shown as a function of time in month and frequency in 1/month. The data set covers 1628 months from May in 1874 to December in 2009. The main peak occurs where the frequency is ~0.0078 1/month whose period corresponds to ~10.68 years.
the wavelet analysis in frequency ranges from 0.03 to 0.09 1/month to check whether the tendency of sunspots
A part of the whole power spectrum where the second most conspicuous peak can be found. The power is consistently concentrated at frequency ≈0.0161/month which corresponds to ~5.21 years.
A part of the whole power spectrum in frequency ranges from 0.03 to 0.09 1/month. Any periodicities of ~1.3 (0.064) ~1.5 (0.056) or ~1.8 (0.046) years (1/month) can be hardly seen.
Zoom-in the whole power spectrum in duration ranges from September in 1986 to December in 2009 (solar cycles 22 and 23) which overlap the data period when Howe et al. (2000 2007) analyzed.
to appear in limited latitudes indeed has biennial periodicities,which has been suggested as evidence of links between the changing structure of the sunspot zone and the tachocline rotation rate oscillations. This frequency range corresponds to one from 11 (0.91 years) to 33 months (2.75 years) in period. As shown in
Fig. 4
, however,the power appears as randomly distributed in almost the entire region of the power spectrum. In other words, we cannot find any periodicities of ~1.3 (0.064), ~1.5 (0.056), or ~1.8 (0.046) years (1/month), as claimed previously(Krivova&Solanki 2002, Knaack et al. 2005, Ternullo 2010a). We, therefore, have to conclude that the oscillating component of ~2 year periodicity is insignificant and inconsistent. According to what we see in Figs. 3 and 4, the only reasonable short-term periodicity one can find in the latitudinal variation of sunspot positions is apparently~5.21 years, but not ~2 years. Previously reported periodicities of ~2 years seem to be biased results due to small sunspots.
As a matter of fact, the 1.3 year periodicity detected by Howe et al. (2000) in the rotation rate at the base of the convection zone is still an open issue, in that Antia & Basu (2000) and Basu & Antia (2001), with a slightly different analysis of the same MDI and GONG data, reported no significant variations. Basu & Antia (2001) saw only an insignificant signal. Moreover, the periodic signal disappears in the post-2001 data, even in Howe et al. (2007). In the later analysis by the original authors, the ~1 year periodicity was considered as variable. That is, even according to the original authors, the signal of ~1.3 year periodicity can be seen in the data set from 1995 to 2000 only, but cannot be seen in the data set from 2001 up to the near present. Even though we do not see any significant peaks around a frequency corresponding to ~1.3 years in the power spectrum in
Fig. 4
, we decide to zoom-in the power spectrum for the duration from September 1986 to December 2009 (solar cycles 22 and 23), whose duration overlaps the data period analyzed by Howe et al. (2000, 2007). If and only if we do see some similar temporal behaviors in the power distribution at frequency of ~0.064 1/month, which corresponds to the 16 month (~1.3 years) period, we can still struggle to clean the power spectrum around that frequency to obtain a signal. As shown in
Fig. 5
, unfortunately, we indeed find a consistent peak at frequency~0.063 1/month from ~1990 to 2009, which is one of the other peaks of noise. In other words, there is no behavior seen as expected (i.e., only strong peaks from ~1995 to ~2000 and no peaks otherwise). Hence, even if what Howe et al. (2000, 2007) have found is a real feature, and even if the component at frequency ~0.063 is significant, we do not see any evidence whatsoever supporting its association in the latitudinal variation of the position of sunspots in view of COL.
4. SUMMARY AND CONCLUSION
The variation of area-weighted latitude of sunspots appearing during the period from 1874 to 2009 has been studied. We focus our investigation on short-term periodicities. We first determine the COL by averaging the latitude with the weight function in area. Then, we analyzed the time series of COL using the wavelet transform technique. We find that a periodicity of ~5 years is the second most dominant feature in the time series of COL. The power is consistently concentrated at frequency~0.0161/month. We consider this periodicity to be due to humps between the minima. We show the power spectrum in frequency ranges from 0.03 to 0.09 1/month to check whether the sunspot tendency to appear in limited latitudes indeed has biennial periodicities, which have been suggested as evidence of links between the changing structure of the sunspot zone and the tachocline rotation rate oscillations. As shown, however, we cannot find any periodicities of ~1.3 (0.064), ~1.5 (0.056), or ~1.8 (0.046) years (1month), as claimed previously. We therefore have to conclude that the oscillating component of ~2 year periodicity is insignificant and inconsistent. Previously reported periodicities of ~2 years seem to be biased results due to small sunspots. Hence, even if what Howe et al. (2000, 2007) have found is a real feature, and even if the component at frequency ~0.063 is significant, we do not see any evidence whatsoever supporting its association in the latitudinal variation of the position of sunspots in view of COL.
Acknowledgements
HYC is grateful to Sasha Kosovichev and Philip Scherrer for constructive discussions and hospitality while visiting Hansen Experimental Physics Laboratory, Stanford University where most of work has been done. This work was supported by a Research Grant funded by the Korea Aerospace Research Institute.
Antalova A
,
Gnevyshev MN
1965
Principal characteristics of the 11-year solar activity cycle
SvA
9
198 -
201
Antalova A
,
Gnevyshev MN
1983
Latitudinal distribution of sunspot areas during the period 1874-1976
CoSka
11
63 -
93
Antia HM
,
Basu S
2000
Temporal variations of the rotation rate in the solar interior
ApJ
541
442 -
448
DOI : 10.1086/309421
Chang H-Y
2007
A new method for north-south asymmetry of sun spot area
JASS
24
261 -
268
Donner R
,
Thiel M
2007
Scale-resolved phase coherence analysis of hemispheric sunspot activity: a new look at the North-South asymmetry
A&A
475
L33 -
L36
DOI : 10.1051/0004-6361:20078672
Gleissberg W
1971
Bemerkungen uber einen vermuteten extrasolaren Einfluss auf die Sonnentatigkeit
MitAG
30
151 -
153
Goel A
,
Choudhuri AR
2009
The hemispheric asymmetry of solar activity during the last century and the solar dynamo
RAA
9
115 -
126
DOI : 10.1088/1674-4527/9/1/010
Grossmann A
,
Morlet J
1984
Decomposition of Hardy functions into square integrable wavelets of constant shape
SIAM J Math Anal
15
723 -
736
DOI : 10.1137/0515056
Hathaway DH
,
Nandy D
,
Wilson RM
,
Reichmann EJ
2003
Evidence that a deep meridional flow sets the sunspot cycle period
ApJ
589
665 -
670
DOI : 10.1086/374393
Howard R
,
LaBonte BJ
1980
The sun is observed to be a torsional oscillator with a period of 11 years
ApJ
239
L33 -
L36
DOI : 10.1086/183286
Howe R
,
Christensen-Dalsgaard J
,
Hill F
,
Komm RW
,
Larsen RM
2000
Dynamic variations at the base of the solar convection zone
Sci
287
2456 -
2460
DOI : 10.1126/science.287.5462.2456
Howe R
,
Christensen-Dalsgaard J
,
Hill F
,
Komm R
,
Schou J
2007
Temporal variations in solar rotation at the bottom of the convection zone: The current status
AdSpR
40
915 -
918
DOI : 10.1016/j.asr.2006.12.047
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
438
1067 -
1082
DOI : 10.1051/0004-6361:20042091
Krivova NA
,
Solanki SK
2002
The 1.3-year and 156-day periodicities in sunspot data: wavelet analysis suggests a common origin
A&A
394
701 -
706
DOI : 10.1051/0004-6361:20021063
Le G-M
,
Wang J-L
2003
Wavelet analysis of several important periodic properties in the relative sunspot numbers
ChJAA
3
L391 -
L394
DOI : 10.1088/1009-9271/3/5/391
Li KJ
,
Gao PX
,
Zhan LS
2009a
Synchronization of hemispheric sunspot activity revisited: wavelet transform analyses
ApJ
691
L537 -
L540
DOI : 10.1088/0004-637X/691/1/537
Nagovitsyn YA
1997
A nonlinear mathematical model for the solar cyclicity and prospects for reconstructing the solar activity in the past
AstL
23
742 -
748
Pelt J
,
Brooke J
,
Pulkkinen PJ
,
Tuominen I
2000
A new interpretation of the solar magnetic cycle
A&A
362
1143 -
1150
Pulkkinen PJ
,
Brooke J
,
Pelt J
,
Tuominen I
1999
Long-term variation of sunspot latitudes
A&A
341
L43 -
L46
Solanki SK
,
Wenzler T
,
Schmitt D
2008
Moments of the latitudinal dependence of the sunspot cycle: a new diagnostic of dynamo models
A&A
483
623 -
632
DOI : 10.1051/0004-6361:20054282
Temmer M
,
Rybák J
,
Bendik P
,
Veronig A
,
Vogler F
2006
Hemispheric sunspot numbers Rn and Rs from 1945-2004: catalogue and N-S asymmetry analysis for solar cycles 18-23
A&A
447
735 -
743
DOI : 10.1051/0004-6361:20054060
Temmer M
,
Veronig A
,
Hanslmeier A
2002
Hemispheric sunspot numbers Rn and Rs: catalogue and N-S asymmetry analysis
A&A
390
707 -
712
DOI : 10.1051/0004-6361:20020758
Ternullo M
,
Looking inside the butterfly diagram
2007b
AN
Looking inside the butterfly diagram
328
1023 -
1032
DOI : 10.1002/asna.200710868
Ternullo M
2010b
The butterfly diagram internal structure
MSAIS
14
202 -
205
Tuominen J
1952
On the dependence of the systematic drift of sunspots in heliographic latitude on phase in the 22-year cycle of the Sun. With 3 figures
ZA
30
261 -
274
Waldmeier M
1957
Der langeSonnenzyklus. Mit 3 Textabbildungen
ZA
43
149 -
160
Zolotova NV
,
Ponyavin DI
2007
Was the unusual solar cycle at the end of the XVIII century a result of phase asynchronization?
A&A
470
L17 -
L20
DOI : 10.1051/0004-6361:20077681