Short Periodicities in Latitudinal Variation of Sunspots

Journal of Astronomy and Space Sciences.
2011.
Jun,
28(2):
103-108

This is an Open Access article distributed under the terms of theCreative Commons Attribution Non-Commercial License(http://creativecommons.org/licenses/by-nc/3.0/)which permits unrestrictednon-commercial use, distribution, and reproduction in any medium,provided the original work is properly cited.

- Received : April 08, 2011
- Accepted : April 22, 2011
- Published : June 15, 2011

Download

PDF

e-PUB

PubReader

PPT

Export by style

Share

Article

Metrics

Cited by

TagCloud

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}
. 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
A_{i}
is the
i
-th sunspot's area and
L_{i}
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.
^{1}http://solarscience.msfc.nasa.gov/greenwch.shtml
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.

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
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

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
Lager Image

Lager Image

Lager Image

Lager Image

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**

Basu S
,
Antia HM
2001
A study of possible temporal and latitudinal variations in the properties of the solar tachocline
MNRAS
324
498 -
508
** DOI : 10.1046/j.1365-8711.2001.04364.x**

Chang H-Y
2007
A new method for north-south asymmetry of sun spot area
JASS
24
261 -
268

Chang H-Y
2008
Stochastic properties in North-South asymmetryof sunspot area
NewA
13
195 -
201
** DOI : 10.1016/j.newast.2007.08.007**

Chang H-Y
2009
Periodicity of North-South asymmetry of sunspot area revisited: cepstrum analysis
NewA
14
133 -
138
** DOI : 10.1016/j.newast.2008.07.001**

Cho I-H
,
Chang H-Y
2011
Latitudinal distribution of sunspots revisited
JASS
28
1 -
7
** DOI : 10.5140/JASS.2011.28.1.00**

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**

Javaraiah J
2003
Long-term variations in the solar differential rotation
SoPh
212
23 -
49
** DOI : 10.1023/A:1022912430585**

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
2010
Latitude migration of solar filaments
MNRAS
405
1040 -
1046
** DOI : 10.1111/j.1365-2966.2010.16508.x**

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**

Li KJ
,
Gao PX
,
Zhan LS
2009b
Synchronization of sunspot numbersand sunspot areas
SoPh
255
289 -
300
** DOI : 10.1007/s11207-009-9328-7**

Li KJ
,
Gao PX
,
Zhan LS
,
Shi XJ
,
Zhu WW
2008
On the asynchronization of hemispheric high-latitude solar activity
MNRAS
391
L34 -
L38
** DOI : 10.1111/j.1745-3933.2008.00551.x**

Li KJ
,
Liang H-F
,
Feng W
2010a
Phase shifts of the paired wings of butterfly diagrams
RAA
10
1177 -
1185
** DOI : 10.1088/1674-4527/10/11/008**

Li KJ
,
Liang HF
,
Yun HS
,
Gu XM
2002
Statistical behavior of sunspot groups on the solar disk
SoPh
205
361 -
370
** DOI : 10.1023/A:1014288424727**

Li KJ
,
Liu XH
,
Gao PX
,
Zhan LS
2010b
The North-South asymmetry of filaments in solar cycles 16-21
NewA
15
346 -
350
** DOI : 10.1016/j.newast.2009.09.009**

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

Schwabe HS
1843
Die Sonne
AN
20
283 -
286
** DOI : 10.1002/asna.18430201706**

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
2007a
The butterfly diagram fine structure
SoPh
240
153 -
164
** DOI : 10.1007/s11207-006-0261-8**

Ternullo M
,
Looking inside the butterfly diagram
2007b
AN
Looking inside the butterfly diagram
328
1023 -
1032
** DOI : 10.1002/asna.200710868**

Ternullo M
2010a
The butterfly diagram internal structure
Ap&SS
328
301 -
305
** DOI : 10.1007/s10509-010-0270-9**

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

Waldmeier M
1971
The asymmetry of solar activity in the years 1959-1969
SoPh
20
332 -
344
** DOI : 10.1007/BF00159763**

Wolf R
1862
Schreiben des Herrn Prof. R. Wolf Directors der Sternwarte in Zurich an den Herausgeber
AN
56
355 -
358
** DOI : 10.1002/asna.18620562303**

Zolotova NV
,
Ponyavin DI
2006
Phase asynchrony of the North-South sunspot activity
A&A
449
L1 -
L4
** DOI : 10.1051/0004-6361:200600013**

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**

Zolotova NV
,
Ponyavin DI
,
Marwan N
,
Kurths J
2009
Long-term asymmetry in the wings of the butterfly diagram
A&A
503
197 -
201
** DOI : 10.1051/0004-6361/200811430**

Citing 'Short Periodicities in Latitudinal Variation of Sunspots
'

@article{ OJOOBS_2011_v28n2_103}
,title={Short Periodicities in Latitudinal Variation of Sunspots}
,volume={2}
, url={http://dx.doi.org/10.5140/JASS.2011.28.2.103}, DOI={10.5140/JASS.2011.28.2.103}
, number= {2}
, journal={Journal of Astronomy and Space Sciences}
, publisher={The Korean Space Science Society}
, author={Kim, Bang-Yeop
and
Chang, Heon-Young}
, year={2011}
, month={Jun}