We have investigated drought periodicities occurred in Seoul to find out any indication of relationship between drought in Korea and solar activities. It is motivated, in view of solar-terrestrial connection, to search for an example of extreme weather condition controlled by solar activity. The periodicity of drought in Seoul has been re-examined using the wavelet transform technique as the consensus is not achieved yet. The reason we have chosen Seoul is because daily precipitation was recorded for longer than 200 years, which meets our requirement that analyses of drought frequency demand long-term historical data to ensure reliable estimates. We have examined three types of time series of the Effective Drought Index (EDI). We have directly analyzed EDI time series in the first place. And we have constructed and analyzed time series of histogram in which the number of days whose EDI is less than -1.5 for a given month of the year is given as a function of time, and one in which the number of occasions where EDI values of three consecutive days are all less than -1.5 is given as a function of time. All the time series data sets we analyzed are periodic. Apart from the annual cycle due to seasonal variations, periodicities shorter than the 11 year sunspot cycle, ~ 3, ~ 4, ~ 6 years, have been confirmed. Periodicities to which theses short periodicities (shorter than Hale period) may be corresponding are not yet known. Longer periodicities possibly related to Gleissberg cycles, ~ 55, ~ 120 years, can be also seen. However, periodicity comparable to the 11 year solar cycle seems absent in both EDI and the constructed data sets.
1. INTRODUCTION
Timely and accurate forecasts that can be used to minimize their damage are highly desirable since droughts are among the greatest natural disasters threatening the human population in the world. Unfortunately, however, the causes of drought are various, and even different for a separate region and each season, and, moreover, most elements that trigger droughts cannot be regulated. This is partly why many countries rely on the drought monitoring system as the most effective system for allaying drought damage by early drought detection and issuing warnings (e.g., Boken et al. 2005, Oh et al. 2010). Factors including relationships between solar activity and terrestrial weather that affect the amount of rainfall eventually controlling drought attributes have been under careful scrutiny to understand its characteristics, namely its possible duration, intensity/severity, and frequency (Ponce et al. 2000, Kim et al. 2009, Chang & Oh 2012).
In the meantime there is growing interest in roles that the Sun plays on terrestrial weather and climate, not to mention space weather (Svensmark & Friis-Christensen 1997, Tinsley 2000, Scafetta & West 2006, Cho & Chang 2008, Kniveton et al. 2008, Cho et al. 2012). Changes in the temperature and dynamics in the troposphere are suggested to correlate with the Earth's magnetic/electric changes corresponding to solar activities (Burns et al. 2007, 2008). The study of responses of atmospheric ionization and the global electric circuit to varying GCRs demonstrates that the atmospheric transparency is associated with GCR flux which itself is modulated by the solar magnetic field (Roldugin & Tinsley 2004). The response of daily sea-level pressure (SLP) over Korean Peninsula to the high-speed solar wind stream (HSS) is examined, resulting in 0.81 with day lag in the correlation coefficient between the average profiles of superposed SLP and HSS (Cho et al. 2012).
The solar activity cycles excite periodical and pseudo periodical variations of all Earth systems, such as, climate, water cycles, rainfalls, streamflows, groundwater, and soon. This rhythm has been associated with possible forcing by the Hale solar magnetic cycle and/or the 18.6-year lunar nodal cycle (Hodell et al. 2001, Hwang et al. 2005, Byun et al. 2008, Penalba & Vargas 2004, Leal-Silva & Velasco Herrera 2012, Yoo et al. 2012). For instance, Cook et al. (1997) made an assessment of these possible solar/lunar forcings on drought area index (DAI), finding a reasonably strong statistical association between the bidecadal drought area rhythm in the Western United States and years of Hale solar cycle minima and 18.6-year lunar tidal maxima. Analysis of the historical time series of rainfall data in Italy has revealed significant 11-year, 22-year and 44-year variabilities that are possibly linked to coherent fluctuations in the high-speed solar wind precipitation into the atmosphere (Mazzarella & Palumbo 1992). Periods corresponding to the 11 year sunspot cycle, the Hale cycle and others have been detected in various hydrologic data (Rao et al. 1992). The Southern Oscillation Index (SOI) is also shown to track the pairing of sunspot cycles in ~88 year periods (Landscheidt 2000, Baker 2008, Bettolli et al. 2010).
As mentioned above, drought occurs due to various causes. In some places, the solar forcing may thus play a role on one hand. In other places, on the other hand, it can be overshadowed by other reasons, such as meteorological causes, even if the Sun might have a playing role indeed, so that the association of the solar forcing with drought is negligible in a drought modeling. One may be interested in finding such a relation with data collected in Korea when especially tailored model for Korean drought is under development. Min et al. (2003) have found that the frequency of occurrence of droughts in Korea has significant time intervals of 2-3 and 5-8 years and has been increasing since the 1980s using the Standardized Precipitation Index (SPI) as a drought index. Wang et al. (2007) have also analyzed monthly precipitation data from the Seoul observatory and showed that spectral peaks exist in 4-year cycles. Recently, Byun et al. (2008) have claimed that drought periodicities correlated with the solar cycles can be found in records of precipitation of Seoul and in other available data including records of ritual praying for rain, drought damage, and years of deficient precipitation for the duration of 536 years. In their analysis, drought intensity is quantified by using the Effective Drought Index (EDI), which is an intensive measure that considers daily water accumulation with a weighting function for time passage (Byun & Whilhite 1999). EDI is one of widely used drought indices, including the Palmer drought severity index (PDSI; Palmer 1965), rainfall anomaly index (RAI; van Rooy 1965), deciles (Gibbs & Maher 1967), crop moisture index (CMI; Palmer 1968), standardized precipitation index (SPI; McKee et al. 1993, 1995), and cropspecific drought index (CSDI; Meyer & Hubbard 1995), and surface water supply index (SWSI; Shafer & Dezman 1982). Compared to other drought indices, the EDI has several advantages (see, e.g., Mishra & Singh 2010).
In this contribution the periodicity of drought in Seoul, where daily precipitation was recorded from 1778 to 2012, has been re-examined using the wavelet transform technique (e.g., Chang 2006). By doing so, the solar influences on local variations of drought can be investigated by means of available time series of observed meteorological data. This analyzing technique is used in many fields of physics and engineering, such as, acoustics, geophysics, helioseismology, image processing, and soon. It is a very powerful tool to obtain information on the temporal behavior of an oscillatory signal. This paper is organized as follows. We begin with brief descriptions of data in Section2. We present and discuss results in Section 3. Finally, we summarize and conclude in Section 4
2. DATA
We have used for the present analysis EDI of Seoul in Korea during the period from 1778 to 2012. The data set has been taken from the Meteorological Disaster Research Laboratory website , Department of Environmental Atmospheric Sciences, Pukyong National University. Time series of EDIs for other stations distributed in South Korea, North Korea, Japan, and China can be also obtained as ASCII text files. Computer codes for EDI are provided so that one may easily appreciate how EDI is defined and computed. Typically, EDI varies in the range from -2.0 to 2.0. The drought range of EDI indicates extremely dry conditions at EDI < -2.0, severe drought at -2.0 < EDI < -1.5 and moderate drought at -1.5 < EDI < -1.0. In
Fig. 1
, we show the daily EDI as a function of time for Seoul. One may identify the year with relative ease when the severe drought occurred.
3. RESULTS AND DISCUSSION
In
Fig. 2
, we show result from the wavelet transform
Daily EDI as a function of time for Seoul during the period from 1778 to 2012.
analysis of the data shown in
Fig. 1
. The horizontal axis represents time in day counting from the beginning of the data for 235 years, and the vertical axis the cyclic frequency in
The power is shown in gray scale. The power is scaled by the mean value of the noise. Several high power peaks appear at ~0.00274, ~0.001, ~0.00075, ~0.0005, ~0.0002
whose periods correspond to ~1, ~2.7, ~3.7, ~5.5, ~15 years, respectively. The most dominant peak is the one of frequency ~0.00274
It is due to the fact that drought has a seasonal variation. That is, in Korea basically we have rainy summer seasons and dry spring and winter seasons, which leads an annual modulation. We have found some marginal peaks corresponding to periods of ~70-100 years. But considering that the total data length is 235 years peaks should be regarded as statistically insignificant. We, therefore, conclude that periodicities associated with the solar cycles cannot be found in time series of EDI data. Another point to note is that most power is concentrated in 19th century and that in modern times there seems no severe drought occurred. Unlike the Fourier analysis, the wavelet transform analysis does also show where the modes of the signal exist in time. This is a strong point of the wavelet transform analysis over the traditional Fourierbased analysis. This finding can be understood in the sense the average value of EDI is more negative in 19th century implying severity of drought, as seen in
Fig. 1
.
In
Fig. 3
, we show result from the wavelet transform analysis of time series of histogram in which the number of days whose EDI is less than -1.5 for a given month of
Result of wavelet transform analysis from time series shown in Fig. 1. The horizontal axis represents time in day counting from the beginning of the data for 235 years, and the vertical axis the cyclic frequency in day/1. The power is shown in gray scale. The power is scaled by the mean value of the noise.
Result of the wavelet transform analysis of time series in which the number of days whose EDI is less than -1.5 for a given month of the year is given as a function of time. The horizontal axis represents time in month counting from the beginning of the data for 235 years, and the vertical axis the cyclic frequency in 1/month The power is shown in gray scale. The power is again scaled by the mean value of the noise.
the year is given as a function of time. The horizontal axis represents time in month counting from the beginning of the data for 235 years, and the vertical axis the cyclic frequency in
The power is shown in gray scale. The power is again scaled by the mean value of the noise. High peaks can be seen at ~0.083, ~0.028, ~0.022, ~0.014
whose periods correspond to ~1, ~2.9, ~3.8, ~6 years, respectively. The most dominant peak is obviously the one of frequency ~0.083
In
Fig. 4
, we show similar plot to
Fig. 3
except the vertical scale is different. Here, we see peaks at ~0.0015, ~0.0007
whose periods correspond to ~55, ~120 years, respectively. Even though a peak at ~0.007
corresponding to the 11 year solar cycle does not appear in the power spectrum, peaks possibly related to Gleissberg cycles can be found. We have repeated the
Similar as Fig. 3, except the frequency range is different.
same analysis with the critical EDI setting to -2.0 instead of -1.5. We find peaks at the same frequencies. Yet, any peaks relating 11 year solar cycle cannot be found. We point out that, despite apparent similarity, searching for a periodicity in a power spectrum resulting from the histogram as in
Figs. 3
and 4 is rather a proper way than from time series of EDI itself as in
Fig. 2
. Since the time series like one shown in
Fig. 1
commonly shows periodic behaviors of droughts (negative EDIs) and floods (positive EDIs) at the same time, the resulting power spectrum is a sum of signals due to pulses of positive and negative EDIs. Hence, power spectrum as shown in
Fig. 2
should be interpreted with due care.
In
Fig. 5
, we show result from the wavelet transform analysis of time series in which the number of occasions where EDI values of three consecutive days are all less than -1.5 is given as a function of time. By this, we meant to count severe or prolonged droughts. Here we show the power spectrum of occurrence of severe droughts in a sense. The horizontal axis represents time in month, and the vertical axis the cyclic frequency in
The power is shown in gray scale normalized by the mean value of the noise. High peaks can be seen at ~ 0.083, ~0.039, ~0.028, ~0.022, ~0.013, ~0.01
whose periods correspond to ~1, ~2.1, ~2.9, ~3.8, ~6.4, ~8.3 years, respectively. In
Fig. 6
, we show similar plot to
Fig. 5
except the vertical scale is different. Here, we see peaks at ~0.0015, ~0.0007
whose periods correspond to ~55, ~120 years, respectively. Even though a peak at ~0.007
corresponding to the 11 year solar cycle does not appear in the power spectrum, peaks possibly related to Gleissberg cycles can be found. We have repeated the same analysis with the critical EDI setting to -2.0 instead of -1.5. We also repeated a case of five consecutive day long drought with critical EDI of -1.5 and -2.0. Generally speaking, same conclusions are drawn so that our conclusion is robust.
Result of the wavelet transform analysis of time series in which the number of occasions where EDI values of three consecutive days are all less than -1.5 is given as a function of time. The horizontal axis represents time in month counting from the beginning of the data for 235 years, and the vertical axis the cyclic frequency in 1/month. The power is shown in gray scale. The power is again scaled by the mean value of the noise.
Similar as Fig. 5, except the frequency range is different.
4. SUMMARY AND CONCLUSION
The main purpose of study on periodicities of droughts is to provide a scientific tool to predict their occurrences in the future. Another scientific motivation in view of solarterrestrial connection is to search for an example of extreme weather condition controlled by solar activity. We have investigated drought periodicities occurred in Seoul to find out any indication of a relationship between drought in Korea and solar activities. The periodicity of drought in Seoul, where daily precipitation was recorded from 1778 to 2012, has been re-examined using the wavelet transform technique as the consensus is not achieved yet. The reason we have chosen Seoul is because daily precipitation was recorded for longer than 200 years, which meets our requirement that analyses of drought frequency demand long-term historical data to ensure reliable estimates.
Supposed the site is a random choice from homogeneous stations in Korea, we have investigated the solar influences on local variations of drought in Korea using time series data sets constructed with EDI data. We have examined three types of time series of EDI. We have directly analyzed EDI time series in the first place. And we have constructed and analyzed time series of histogram in which the number of days whose EDI is less than -1.5 for a given month of the year is given as a function of time, and one in which the number of occasions where EDI values of three consecutive days are all less than -1.5 is given as a function of time. We have repeatedly generated and tested similar time series data with different criteria too. All the time series data sets we analyzed are periodic. Apart from the annual cycle, periodicities shorter than the 11 year sunspot cycle reported have been confirmed. Longer periodicities possibly related to Gleissberg cycles can be also seen. However, periodicity comparable to the 11 year sunspot cycle seems absent in both EDI and the constructed data sets.
Acknowledgements
HYC was supported by the National Research Foundation of Korea Grant funded by the Korean government (NRF2011-0008123).
Baker RGV
2008
Exploratory Analysis of Similarities in Solar Cycle Magnetic Phases with Southern Oscillation Index Fluctuations in Eastern Australia.
Geographical Research
46
380 -
398
DOI : 10.1111/j.1745-5871.2008.00537.x
Bettolli ML
,
Penalba O
,
Rivera J
2010
Trends and periodicities in the annual amount of dry days over Argentina.
looking towards the climatic change
95
27 -
33
Boken VK
,
Cracknell AP
,
Heathcote RH
2005
Monitoring and Predicting Agricultural Drought: A Global Study
Oxford University Press
New York
472 -
Burns GB
,
Tinsley BA
,
Frank-Kamenetsky AV
,
Bering EA
2007
Interplanetary magnetic field and atmospheric electric circuit influences on ground-level pressure at Vostok.
JGR
112
D04103 -
DOI : 10.1029/2006JD007246
Burns GB
,
Tinsley BA
,
French WJR
,
Troshichev OA
,
Frank-Kamenetsky AV
2008
Atmospheric circuit influences on ground-level pressure in the Antarctic and Arctic.
JGR
113
D15112 -
DOI : 10.1029/2007JD009618
Byun HR
,
Lee SJ
,
Mori S
,
Choi K S
,
Lee S M
(2008)
Study on the Periodicities of Droughts in Korea
APJAS
44
417 -
441
Chang HY
2006
Time/Frequency Analysis Of Terrestrial Impact Crater Records.
JASS
23
199 -
208
Cho IH
,
Chang HY
2008
Long term variability of the sun and climate change.
JASS
25
395 -
404
Cho IH
,
Kwak YS
,
Marubashi K
,
Kim YH
,
Park YD
2012
Changes in sea-level pressure over South Korea associated with high-speed solar wind events.
Adv. Space Res.
50
777 -
782
DOI : 10.1016/j.asr.2011.06.025
Gibbs WJ
,
Maher JV
1967
Rainfall Deciles as Drought Indicators, In: Bureau of Meteorology Bull. 48, Commonwealth of Australia
Bureau of Meteorology
Melbourne
Hodell DA
,
Brenner M
,
Curtis JH
,
Guilderson T
2001
Solar forcing of drought frequency in the Maya lowlands.
Science
292
1367 -
1370
DOI : 10.1126/science.1057759
Hwang SH
,
Kim J
,
Won YI
,
Cho HK
,
Kim JS
2005
Statistical characteristics of secondary ozone density peak observed in Korea.
Adv. Space Res.
36
952 -
957
DOI : 10.1016/j.asr.2005.05.080
Kim DW
,
Byun HR
,
Choi KS
2009
Evaluation, modification, andapplication of the Effective Drought Index to 200-Yeardrought climatology of Seoul
Korea. J. Hydro.
378
1 -
12
DOI : 10.1016/j.jhydrol.2009.08.021
Kniveton DR
,
Tinsley BA
,
Burns GB
,
Bering EA
,
Troshichev OA
2008
Variations in global cloud cover and the fair-weather vertical electric field.
JASTP
70
1633 -
1642
Leal-Silva MC
,
Velasco Herrera VM
2012
Solar forcing on the ice winter severity index in the western Baltic region
JASTP
89
98 -
109
McKee TB
,
Doesken NJ
,
Kleist J
1993
The Relationship of Drought Frequency and Duration to Time Scales, In: 8th Conference on Applied Climatology
American Meteorological Society
Anaheim
179 -
184
McKee TB
,
Doesken NJ
,
Kleist J
1995
Drought Monitoring with Multiple Time Scales, In: 9th Conference on Applied Climatology
American Meteorological Society
Dallas
233 -
236
Mazzarella A
,
Palumbo F
1992
Rainfall Fluctuations over Italy and Their Association with Solar Activity.
Theor. Appl. Climatol.
45
201 -
207
DOI : 10.1007/BF00866193
Meyer SJ
,
Hubbard KG
1995
Extending the Crop-specific Drought Index to Soybean, In: Ninth Conf. on Applied Climatology
American Meteorological Society
Dallas
258 -
259
Min SK
,
Kwon WT
,
Park EH
,
Choi YG
2003
Spatial And Temporal Comparisons Of Droughts Over Korea With East Asia, Int
J. Climatol.
23
223 -
233
DOI : 10.1002/joc.872
Oh SB
,
Kim DW
,
Choi KS
,
Byun HR
2010
Introduction of East Asian Drought Monitoring System.
SOLA
6A
9 -
12
DOI : 10.2151/sola.6A-003
Palmer WC
(1965)
Meteorologic Drought, US Department of Commerce, Weather Bureau
Research Paper
45
1 -
58
Penalba OC
,
Vargas WM
(2004)
Interdecadal and interannual variations of annual and extreme precipitation over central-northeastern Argentina.
In: International Journal of Climatology
24
1565 -
1580
DOI : 10.1002/joc.1069
Rao AR
,
Jeong GD
,
Chang FJ
1992
Estimation of periodicities in hydrologic data.
Stochastic Hydrol. Hydraul.
6
270 -
288
DOI : 10.1007/BF01581621
Roldugin VC
,
Tinsley BA
2004
Atmospheric transparency changes associated with solar wind-induced atmospheric electricity variations.
JASTP
66
1143 -
1149
Scafetta N
,
West BJ
2006
Phenomenological solar contribution to the 1900-2000 global surface warming.
GRL
33
L05708 -
Shafer BA
,
Dezman L
(1982)
Development of a Surface Water Supply Index (SWSI) to Assess the Severity of Drought Conditions in Snowpack Runoff Areas.
In: Preprints, Western Snow Conf., Reno, NV, Colorado State University
164 -
175
Svensmark H
,
Friis-Christensen E
1997
Variation of cosmic ray flux and global cloud coverage-a missing link in solar-climate relationships.
JASTP
59
1225 -
1232
Tinsley BA
2000
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
DOI : 10.1023/A:1026775408875
van Rooy MP
1965
A rainfall anomaly index independent of time and space.
Notos
14
43 -
Wang B
,
Jhun JG
,
Moon BK
2007
Variability and singularity of Seoul, South Korea, Rainy Season
J. Climate
20
2572 -
2580
DOI : 10.1175/JCLI4123.1
Yoo J
,
Kwon HH
,
Kim WT
,
Ahn JH
2012
Drought frequency analysis using cluster analysis and bivariate probability distribution.
J. Hydrol.
420-421
102 -
111
DOI : 10.1016/j.jhydrol.2011.11.046