Spatiotemporal Rainfall Variability and Drought Assessment during Past Five Decades in South Korea Using SPI and SPEI

About 41% of the earth is drought-affected, which has impacted nearly 2 billion people, and it is expected that more than 90% of terrestrial areas will be degraded by 2050. To evade and mitigate the harmful impacts of drought, it is necessary to study the rainfall variability and assess the drought trend at a global and regional level. This study utilized 70 meteorological stations in South Korea to evaluate the rainfall variability, drought, and its trend during the past five decades using the standardized precipitation evapotranspiration index (SPEI) and the standardized precipitation index (SPI). Rainfall data normality was assessed with mean, standard deviation, skewness, and kurtosis. The highest amount of rainfall was observed in the months of June, July, and August. The SPI and SPEI 12-month results revealed that 1982, 1988, 2008, 2015, and 2017 were dry years throughout the country, while from 2013 to 2017 mixed drought events were observed for the 6-month time series. The Mann-Kendall trend test was applied to the 1- and 12-month time series, and the results revealed that the months of January, March, April, May, June, and August had a significant negative trend, which means drought is increasing in these months, while the months of September, October, and December had a significant positive trend, which means wetter conditions prevailed in these months during the study period. It was observed in the 12-month time series that only two met stations had a significant negative trend, while only one had a significant positive trend. It was found that January and March were the driest months, and October was the wettest month. The detected drought events in this research are consistent with ENSO events. We have observed differences in drought characteristics (duration and frequency) for both indices. Climatic data revealed that South Korea has faced drought conditions (rainfall deficit) due to a shortened monsoon season. This study can provide guidance on water management strategies under the changing pattern of drought in South Korea.


Introduction
Drought is one of the most frequently occurring and severe hydrometeorological phenomena, resulting from a decrease in rainfall from normal conditions, causing water scarcity for various water-dependent activities [1,2]. Drought is a costly and destructive hazard with an extensive impact on the environment, socioeconomics, and populations around the globe [3,4]. Due to the shift in climatological patterns (rainfall and temperature), droughts are occurring more frequently [5,6]. Dai, (2013) [7] has projected that the intensity and frequency of drought will increase in the 21st century. Drought is a slow onset event with effects can be seen even after the end of the event, which makes it difficult to identify the start and end of it [8][9][10]. Drought is characterized into four types, deficit of rainfall from its average, initiating meteorological drought [3,11], which leads to a dearth of streamflow developed various drought indices to measure drought severity, i.e., the standardized precipitation index (SPI) [37], the standardized precipitation evapotranspiration index (SPEI) [38], the Palmer drought severity index (PDSI) [39], the reconnaissance drought index (RDI) [40], and the crop moisture index (CMI) [41]. SPI requires only precipitation data, while RDI, PDSI, and SPEI require additional potential evapotranspiration (PET) data to compute drought. Thus, the key aim of this research is to study meteorological drought and its characteristics using the SPI and SPEI drought indexes and to assess the trend of drought using the Mann-Kendall test (MK test). The efforts made through this study will significantly help the authorities to mitigate the impacts of meteorological drought at an early stage before the situation worsens and leads to hydrological and agricultural drought conditions.

Study Area and Data
South Korea is a part of East Asia located between Japan and China at 33-43 • N latitude and 124-131 • E longitude, and the total area of the country is 100,210 km 2 . Figure 1 shows the administrative boundaries, elevation, and meteorological stations (70) used in this study. South Korea experiences four different seasons: winter (Dec-Feb), spring (March-May), summer (June-August), and autumn (September-November) [42,43]. South Korea lies in the Asian monsoon region, which experiences wet summers with more than 50% of annual rainfall with heavy rain and typhoons [27]. Many parts of South Korea receive more than 1200 mm of annual rainfall, which is greater than the global average rainfall [28,44], but some parts receive less than 1000 mm rainfall and have a shortage of rainfall in the crop growing season (April-October); thus, the country experiences droughts [45]. Due to topographic variability and mountainous terrain the climate of South Korea has complex spatial variations. The average temperature in South Korea ranges between 10-15 • C [46,47].
Atmosphere 2021, 12, x FOR PEER REVIEW 3 of 25 on a regional and global scale. To study drought, meteorological scientists have developed various drought indices to measure drought severity, i.e., the standardized precipitation index (SPI) [37], the standardized precipitation evapotranspiration index (SPEI) [38], the Palmer drought severity index (PDSI) [39], the reconnaissance drought index (RDI) [40], and the crop moisture index (CMI) [41]. SPI requires only precipitation data, while RDI, PDSI, and SPEI require additional potential evapotranspiration (PET) data to compute drought. Thus, the key aim of this research is to study meteorological drought and its characteristics using the SPI and SPEI drought indexes and to assess the trend of drought using the Mann-Kendall test (MK test). The efforts made through this study will significantly help the authorities to mitigate the impacts of meteorological drought at an early stage before the situation worsens and leads to hydrological and agricultural drought conditions.

Study Area and Data
South Korea is a part of East Asia located between Japan and China at 33-43° N latitude and 124-131° E longitude, and the total area of the country is 100,210 km 2 . Figure 1 shows the administrative boundaries, elevation, and meteorological stations (70) used in this study. South Korea experiences four different seasons: winter (Dec-Feb), spring (March-May), summer (June-August), and autumn (September-November) [42,43]. South Korea lies in the Asian monsoon region, which experiences wet summers with more than 50% of annual rainfall with heavy rain and typhoons [27]. Many parts of South Korea receive more than 1200 mm of annual rainfall, which is greater than the global average rainfall [28,44], but some parts receive less than 1000 mm rainfall and have a shortage of rainfall in the crop growing season (April-October); thus, the country experiences droughts [45]. Due to topographic variability and mountainous terrain the climate of South Korea has complex spatial variations. The average temperature in South Korea ranges between 10-15 °C [46,47].

Data Collection
The Korean Meteorological Administration (KMA) has installed more than 70 meteorological stations throughout the country for detailed information on precipitation and temperature. The total annual rainfall and mean monthly, minimum, and maximum tem-  (Table 1). The data of meteorological stations were not consistent, as some of the stations were established earlier than others. Therefore, to avoid inconsistency in the analysis, we have chosen 1979 for the beginning of this study. Table 2 shows the geographic location and altitude of all the meteorological stations. To control data quality, the missing data were replaced with average values of the data of nearby meteorological stations [48]. For drought calculations, the SPI and SPEI indices were applied using monthly rainfall data in the SPI calculation and rainfall and temperature data for SPEI calculation in this study. The Mann-Kendall trend test was applied for the calculation of the drought trend in the study region.

Methods
The detailed methodological framework of this study is shown below in Figure 2. Initially, the daily precipitation and temperature data of 70 meteorological stations were acquired from the Korean Meteorological Administration (KMA) for the period of 1979-2020 (42 years). The detailed summary of the meteorological stations is presented in Table 2. The daily precipitation and temperature datasets were summed to create the monthly and annual datasets. The standardized normalized homogeneity test (SNHT) was used to find the peaks and outliers in the data [49]. During data quality evaluation, it was found that only 2-3% of the data were missing which would not have a substantial impact on the results [49]; however, we still replaced the missing data with the average values of the nearby meteorological stations to maintain the quality [48]. Next, we calculated the descriptive statistics in XLSAT Adin of Excel software, the details of which are presented in Table 3. Then, we created the monthly and annual rainfall maps (Figures 3 and 4). We used SPI and SPEI to calculate the drought events on 1-, 6-, and 12-month time scales, which are presented in Figures 5 and 6. The trend of drought was assessed using the Mann-Kendall trend test, and maps were prepared (Figures 7 and 8). Finally, the drought frequency and duration was calculated (Figures 9 and 10).

Methods
The detailed methodological framework of this study is shown below in Figure 2. Initially, the daily precipitation and temperature data of 70 meteorological stations were acquired from the Korean Meteorological Administration (KMA) for the period of 1979-2020 (42 years). The detailed summary of the meteorological stations is presented in Table  2. The daily precipitation and temperature datasets were summed to create the monthly and annual datasets. The standardized normalized homogeneity test (SNHT) was used to find the peaks and outliers in the data [49]. During data quality evaluation, it was found that only 2-3% of the data were missing which would not have a substantial impact on the results [49]; however, we still replaced the missing data with the average values of the nearby meteorological stations to maintain the quality [48]. Next, we calculated the descriptive statistics in XLSAT Adin of Excel software, the details of which are presented in Table 3. Then, we created the monthly and annual rainfall maps (Figures 3 and 4). We used SPI and SPEI to calculate the drought events on 1-, 6-, and 12-month time scales, which are presented in Figures 5 and 6. The trend of drought was assessed using the Mann-Kendall trend test, and maps were prepared (Figures 7 and 8). Finally, the drought frequency and duration was calculated (Figures 9 and 10).

Standardized Precipitation Index (SPI)
The standardized precipitation index (SPI) was developed by Mc Kee in 1993 [37]. It is one of the simplest and most commonly used drought indexes and compares the normalized rainfall with average rainfall to express the deficit and surplus of rainfall for a particular time period, location, and climate [50]. SPI is globally recommended by the World Meteorological Organization (WMO) for drought assessment [51]. SPI can be calculated based on long-term rainfall data for each meteorological station with gamma distribution using the maximum likelihood estimation method for different time scales, i.e., 1, 3, 6, 9, 12, 24, and 48 months. SPI can be calculated on both short-term and long-term droughts, where its short term application is to assess the drought impact on moisture and precipitation, while the long term SPI is calculated to assess the drought impact on agriculture and water resources, i.e., surface water flow, reservoir, and groundwater supplies [11]. In this study we have discussed both short term SPI (1 and 6 months) and long-term SPI (12 months) to assess and characterize drought in South Korea.
SPI was calculated using the SPEI package in R programming. The specific formula for SPI is given below.
To calculate SPI, subtract monthly rainfall (y) from the mean (z) and divide by the standard deviation (σ) of rainfall calculated from the monthly time series. Data were normalized using the gamma distribution method, which is suggested by many researchers to calculate the SPI [52,53].

Standardized Precipitation Evapotranspiration Index (SPEI)
The standardized precipitation evapotranspiration index (SPEI) is similar to SPI except for the addition of the mean minimum and mean maximum temperature to rainfall. In the SPEI calculation, we need the monthly water balance, which is based on the difference between rainfall and potential evapotranspiration (PET). PET can be calculated using various parameters, i.e., temperature, relative humidity, solar radiation, air water vapor, and sensible and latent heat flux. However, most of the meteorological stations in the world do not record this kind of data [3]. So as a proxy, we calculated the PET using the Hargreaves method [54], which requires temperature data and the latitude of the concerned meteorological station. The PET was calculated using R programming, and the monthly PET was subtracted from the monthly precipitation data to obtain the water balance.
Here, B denotes the water balance deficit or surplus at the ith month, R indicates monthly rainfall at month i, and PET is the potential evapotranspiration at the ith month. Now, B i are accumulated for different time scales The SPEI was calculated for monthly time scales i, for n number of calculations using the SPEI package in R programming.

Mann-Kendall (MK)Trend Test
In the literature, we found many statistical methods [48] for trend detection of time series data. Every method has its own advantages and disadvantages. We used the Mann-Kendall (MK) test to identify the trend in SPI and SPEI for dry and wet conditions. It is a nonparametric statistical method [48,55] recommended by the WMO for trend assessment of hydrological and meteorological variables [55,56].

Spatiotemporal Rainfall Variability
The rainfall variability was analyzed based on 70 meteorological stations in South Korea. Most of the meteorological stations in South Korea are located at low elevation with the lowest Pohang (4 m) and the highest Daegwallyeong (772 m). We performed statistical analysis on all the meteorological stations to examine the central tendency and dispersion. The highest standard deviation was found for Geoje (525.68 mm), followed by Daegwallyeong (496.69), Seogwipo (454 mm), Namhae (444.52 mm), Changwon (419.42 mm), and Sancheong (419.75), which revealed high amount of annual rainfall variability at each met station. The lowest amount of annual rainfall variability was found for Yeongcheon (238.6 mm) followed by Daegu (245.2 mm), Uiseong (251.5 mm), Chupungryung (253.42 mm), Andong (260.54 mm), and Yeongdeok (263.7 mm). The data normality was checked using skewness and kurtosis. Skewness shows the symmetry of the data; Daegu, Gunsan, Gwangju, Inje, Cheonan, Jeongeup, Haenam, and Namhae showed a normal distribution, while Cheongju, Longevity, Andong, Jeonju, Mokpo, Pearl, Boryeong, Geumsan, Imsil, Namwon, Jangheung, Geochang, Hapcheon, Miryang, and Sancheong met stations showed negative skewness, which described negative symmetry from the normal. All other met stations showed a positive skewness, which indicated an increase in the amount of rainfall from the long-term normal. Kurtosis was used to measure the flatness and peaks of the data compared to their normal distribution, in which high (positive) values represents peaks of the data, while low (negative) values indicate the flatness of the data near the mean. Table 3 shows that 43 met stations were flat, while only 27 met stations had a high peak from the normal data.
The inverse distance weightage (IDW) interpolation technique was applied on monthly and annual rainfall data to spatially analyze the rainfall in the study area. The results revealed that in South Korea, 50-60% of the annual rainfall occurs in summer during June-August, whereas the spring season (March-May) receives 25-30% of the total annual rainfall ( Figure 3). On average, only 20% of rainfall occurs from September to February. Areas receiving less than 250 mm of rainfall are arid regions, while rainfall occurring between 250 and 750 mm and greater than 750 mm come under the umbrella of semi-arid and humid regions, respectively [57]. Thus, according to this classification, most of South Korea is a humid region; only North Gyeongsang province is arid and semi-arid ( Figure 4). From the mean monthly and mean annual rainfall, we can observe that the northern part of the country receives less rainfall compared to the south and coastal regions of South Korea [26] (Figures 3 and 4). Coastal regions of South Korea receive intense amounts of rainfall in the summer season due to typhoon-induced changes [21]. It was also observed that the western part of the country receives more rainfall compared to east, which may be due to the foehn warming phenomena [58]. The eastern and central part of the country have the Taebaek and Sobaek mountain ranges, respectively (encircled part in Figure 4). In mountain ranges, the windward side receives more rainfall compared to the leeward side (mostly on the east side) [58][59][60]. NGII, (2020) has also discussed the effects of the windward and leeward sides in Korea.

12-Month SPI and SPEI
The SPI and SPEI were applied to analyze the dry and wet situations in the study area during 1979-2020 using 70 meteorological stations' data. Long term SPI and SPEI are the accumulated effect of short-term droughts, thus, we incorporated both short-and long-term SPI and SPEI in this study for the time scale of 1, 6, and 12 months. The 12-month SPI and SPEI revealed drought conditions in South Korea in different years. It was found that in 1982 all provinces except Jeju and North Gyeongsang province experienced moderate drought in both indexes [23] (Figure 5a-i). In the same year, only North Gyeongsang experienced a severe drought condition ( Figure 4). Moreover, Jeju province experienced a severe rainfall deficit in 1984 (Figure 5i). Extreme drought conditions were observed for South Jeolla in 1988 in the SPI; however, the SPEI a showed moderate drought in South Chungcheong and severe drought conditions in South Jeolla (Figure 5b,h). In the same year, South Gyeongsang, North Jeolla, and Jeju experienced severe drought conditions, according to the SPI and SPEI (Figure 5d,g,i) while North Chungcheong Gangwon, and North Gyeongsang underwent moderate drought, according to the SPI and SPEI (Figure 5a,c,d). Moderate drought conditions were observed in 1994 for Gangwon, Gyeonggi, North and South Chungcheong, and North Jeolla province, with severe drought for North and South Gyeongsang (Figure 5d,e) and South Jeolla (Figure 5h) [23]. In 1995, moderate drought was observed for North Jeolla (Figure 5g) and North and South Gyeongsang (Figure 5d,e), while drought was severe for South Jeolla (Figure 5h

6-Month SPI and SPEI
In the present study, we calculated the 6-month SPI and SPEI for the months of June and December for South Korea. Figure 6 shows the 6-month drought for June and December for the SPI and SPEI. Gangwon province experienced a moderate drought in the month of June, according to the SPI in 1985, but it was not detected by the SPEI. In 1988, both the SPEI and SPI detected moderate drought, but the severity of the SPI was higher than the SPEI. In 1995, moderate drought was observed, according to the SPI, but it was not detected by the SPEI. However, in 2014, the SPEI spotted moderate drought, but the SPI did not observe it. The SPEI observed severe and extreme drought conditions for the month of June in Gangwon in 2015 and 2017, but the SPI identified moderate and severe drought conditions, respectively (Figure 6a). In the month of December, extreme and severe drought conditions were observed by the SPEI and SPI in Gangwon province for the years of 1979 and 2015, respectively (Figure 6b).
Gyeonggi province experienced three moderate droughts in the month of June for 1988, 1995, and 2000, according to both indices. In 2017, the SPEI spotted an extreme drought condition, but the SPI identified a severe drought condition (Figure 6c). For the month of December, Gyeonggi province went through moderate drought (1988, 1996, and 2016) and severe drought (1979, 2014, and 2015), according to both indices (Figure 6d). Jeju province experienced five years of moderate drought (1982, 1984, 2002, 2005, and 2013) and two years of severe drought (2000 and 2017), while it was found that the province did not undergo extreme drought for the month of June, according to both indices (Figure 6e). In the month of December, the SPEI showed a severe drought in 1988, but the SPI detected an extreme drought condition. In 2013, the SPEI and SPI revealed extreme drought conditions. Moreover, the province went through five moderate drought years in the month of December, according to both indices (1984,1994,1996,2005, and 2008) (Figure 6f).
North Chungcheong province underwent moderate drought in the month of June 1981, according to the SPI, but it was not detected by the SPEI. The SPI found severe drought in 1982, 1988, and 1995, while the SPEI observed moderate drought. Two more drought years (2014 and 2015) were revealed by both indices (Figure 6g). The month of December suffered moderate drought in 1982, 1983, and 2008, according to both indices, while three severe drought years were also detected by the SPEI and SPI (1994, 1996, and 2015). An extreme drought condition was observed by the SPI in 2001, but it was not seen in the SPEI (Figure 6h). South Gyeongsang province went through two severe droughts (1982 and 2000) and one extreme drought condition (2017) for the month of June, according to both indices. While in 1981 and 2000, the province experienced moderate drought (Figure 6i). Both indices revealed five years of moderate drought (1986, 1990, 2001, 2013, and 2017) in the month of December, while one severe drought (1996) was also detected. The SPEI discovered three severe drought years (1988, 1994, and 2008), but the SPI showed extreme drought in the same years (Figure 6j). For the month of June, North Jeolla province showed two moderate drought years (1992 and 2019) and two severe drought years (2012 and 2014), according to both indices (Figure 6k). However, December revealed five moderate drought years (1982, 1988, 1994, 1996, and 2017) and one severe drought year (2015). While in 2008, the SPI detected an extreme drought condition, but the SPEI revealed a severe drought condition (Figure 6l). South Chungcheong did not experience any extreme drought conditions in the months of June and December (Figure 6m,n). However, three moderate drought years (1982, 1988, and 1992) in the month of June (Figure 6m) and one moderate drought year (2008) in the month of December (Figure 6n) were found by both indices. Furthermore, one severe drought year was (2017) detected for the month of June ( Figure 6m) and four severe drought years (1994, 1996, 2001, and 2015) were found for the month of December (Figure 6n). Not a single year was detected as an extreme drought condition in both indices for the month of June and December (Figure 6m,n). South Gyeongsang was affected by two moderate droughts (1981 and 2007) and two severe droughts (1982 and 2000) for the month of June, according to both indices (Figure 6o). For the month of December, the province revealed four moderate drought years (1986, 1990, 2013, and 2017). In 1996, the province was affected by a severe drought condition. It was also shown that in 1988, 1994, and 2008, the SPEI detected a severe drought, but the SPI showed an extreme drought condition (Figure 6p). South Jeolla was affected by two severe drought years (1982 and 1994) and one extreme drought year (2017) for the month of June, according to both indices. In 2000, the SPEI detected severe drought conditions in South Jeolla province, but the SPI failed to detect this (Figure 6q). For the month of December, the SPEI detected three severe drought years (1988, 1995, and 2008), but the SPI found extreme drought in the same years. Two moderate drought years (1994 and 1996) were also observed by both indices (Figure 6r). In this study, a different pattern of SPI and SPEI was found for the 6-month time scale, which is in contrast to another study [25]. The over and underestimation of drought by SPI may be due to a single variable (rainfall) [27]; however, SPEI performed well in the detection of drought events, due to the additional PET parameter, which is also discussed by Uddin et al., (2020) [25].

Mann-Kendall Trend Test
In the present research, we utilized the MK test to analyze the trend in the 1-and 12-month SPEI and SPI series. However, we noticed a similar trend, which is why we did not incorporate the SPI trend results. The trend results were interpolated using IDW technique. We calculated the trend using a 95% significance level with a p value of ≤0.05 for both the 1-and 12-month SPEI series. In the trend test, positive values indicate wet conditions, while negative values indicate dry condition.
In the one-month SPEI series, a positive trend was found in the months of February (Northern part) and March (Southern part); whereas, in September, October, and December there was a positive trend in whole country. A significant positive trend (wet conditions) was revealed in the months of September, October, and December in 2, 20, and 11 meteorological stations, respectively. In the month of February, there was no significant trend at any meteorological station, and most parts of the country had wet conditions, except some parts in the southeast and west. In the month of July, there was also no significant trend, but there were wet conditions except in Jeju province. In the month of November, no significant trend was detected, but it also showed nominal wet conditions in the country. It was observed that the country has faced a negative trend in the month of January, (North), March (except the Southern part), April (North-East), June (whole country, except Jeju province), and August (North-East and South). A significant negative trend was found for the month of January (8 met stations), March (22), April (1), May (2), June (12), and August (2) [61]. Azam et al., (2018a) [26] also found the highest negative autocorrelation in the summer season and the least in the winter months.
More than 50% of the country's rainfall occurs in the months of June, July, and August [21,26]; however, in the present study, we found that the months of June and August had significant droughts, which is an alarming situation due to increasing temperature and the shortened monsoon period [21,26,62]. Many parts of South Korea receive precipitation in the form of snowfall from December-February, but in the month of January a significant negative trend was observed (dry condition), which is also a frightening situation because melt snow contributes to river flow in the summer season; thus, dry conditions in winter can ultimately cause water shortages in the summer season for the agriculture sector.
The 12-month SPEI results revealed drying conditions in the South-West, South-East, and East-Central part of the country [26], while wet conditions were detected in North, West, and in Jeju province. Jeonju and Daegwallyeong met stations encountered a significant decreasing trend (dry condition), while Seongsan experienced a significant increasing trend (wet condition), and all other met stations did not show any significant trend.

Spatial Distribution of Drought Characteristics
In this study, we calculated the characteristics of drought (frequency and duration) for both indices at different time scales; then, we interpolated it with an IDW technique (Figures 9 and 10). It was noticed that the drought duration for the one-month SPEI and SPI was different in spatial distribution. In the SPEI, the highest drought duration was detected for Jeju province (Figure 9a), while Jeju province experienced moderate drought duration, according to the SPI (Figure 8). It was also spotted that North Gyeongsang fell under the highest drought duration in the one-month SPI index (Figure 9d), but it fell under the lowest drought duration in the 1-month SPEI (Figure 9a). Similarly for the 6-month data, the drought duration was high for Jeju province, according to the SPEI (Figure 9b), but it fell under moderate drought duration, according to the SPI (Figure 9e). The 12-month SPI and SPEI drought duration were slightly different (Figure 9c,f).
For a drought frequency of 1-month, according to the SPEI, the highest frequency was distributed in the central part (North Chungcheong, North Jeolla) and in the south (Jeju province) ( Figure 10a); while, according to the SPI, the highest drought frequency was observed for North Gyeongsang province (Figure 10d), which is in contrast to a study conducted by Azam et al., (2018a) [26]. At the 6-month time scale in the SPEI, a moderate to highest drought frequency was centered on North and South Gyeongsang, North and South Chungcheong, and North Jeolla province (Figure 10b). However, in the SPI, the moderate to highest frequency was distributed in Gangwon, Gyeonggi, and the North Jeolla province (Figure 10e). At the 12-month time scale in the SPEI, the highest drought frequency was distributed in all provinces except Gyeonggi, Gangwon, and Jeju (Figure 10c), while in the SPI, it was only found for Gangwon province (Figure 10f).

Discussion
Shortening of the summer monsoon period or less rainfall for a longer period are major causes of drought in South Korea, because more than 50% of the total rainfall occurs in the summer season, which is shrinking [21,26]. In the present study, drought conditions were witnessed in the central, eastern, and western coastal regions of the country, according to the 12-month SPEI. Azam et al., (2018a) [26] also found increasing drought conditions in the mid-latitudes, southeast, and northeast coastal region, and west side of the country. The results revealed that the summer months had significant drought conditions as compared to the winter months, which is consistent with a study conducted by Azam et al., (2018a) [26]. This is due to the decreasing amount of rainfall in the shortened monsoon period.
In this study, we examined that in last decade the intensity of drought has increased in South Korea and; this was also discussed and projected by Kim et al., (2012) [63]. Lee et al., (2021) [64] also observed and discussed that in last ten years drought events have increased in South Korea and found severe drought conditions in 2013 and 2017; the current study also detected severe drought in those years, as well as increasing drought frequency. Hong et al. (2016) [23] discussed that from early 2014, the precipitation in South Korea has decreased extensively at record levels, which has triggered drought conditions in Korea, which lasted for more than 20 months. Bae et al., (2019) [65] found that the pattern of rainfall is changing every 4-5 years, which endorses the drought events occurring in South Korea every 5-6 years since 1960 [23]. Zhang and Zhou [66] revealed that drought conditions in Japan, China, and the Korean peninsula are usually dominated by East Asian monsoon variability. Choi et al., (2013) [67] explained that in the spring season (April-May), the amount of rainfall increases temporarily due to the northern hemisphere atmospheric circulation, but the rainfall in the spring season decreases due to the weakening of the subtropical western North Pacific high (SWNPH), which can cause droughts in the spring season [68]. Byun HR (1996) [69] found that deepening of cyclones and abnormally low and high temperatures can cause summer droughts in Korea. Choi et al., (2011) [70] observed that intensified drought in the summer season occurs in Northeast Asia due to the northeasterlies cyclonic and anticyclonic circulation. Park and Schubert (1997) [71] revealed that the 1994 drought in East Asia was forced as a result of the deficit of the summer/monsoon rainfall with an early development of upper level anticyclonic flow to east of the Tibetan plateau. Bae et al., (2018) [72] detected intensive drought in the 1-month SPEI for almost all met stations in Korea, which was also consistent with the 1-month SPEI and SPI results calculated in this study. According to the Korean climate change and drought survey report, Korea faced drought events in 1981-1982, 1994-1995, 2001-2007, and 2014-2015, which was also reflected in this study using SPI and SPEI for 1982, 1988, 1994, 2001, 2008. Almost the same years of drought have been detected in other studies conducted in the Asian region [4,11]. Hong et al., (2016) [23] discussed the 2014-2015 drought as the worst in the history of South Korea in which annual rainfall decreased 35-50% from its normal, and the same years were detected as severe to extreme drought in this study for six out of a total of nine provinces in South Korea in the results of the 12-month SPI and SPEI. Moreover, the 1994-1995 drought event was also found in the present study for all regions, except Jeju province. However, the 1981-1982 drought event was not clearly observed by both the indices. Comparing the drought events with the El Niño-Southern Oscillation (ENSO), it was found that most of the drought events were consistent with moderate, strong, and very strong ENSO events [73]. It is expected that in future, the severity of droughts will increase in South Korea due to climate change [53]. It was shown in the present study that droughts occur frequently throughout the year in different parts of the country with various intensities, which was also revealed by Azam et al., (2018a) [26]. It was noted that the frequency of drought was high in the 1-month SPEI for the southern and central regions of the country, which is consistent with a study conducted by Bae et al., (2018) [72]; however, the 12-month SPI drought frequency was high in the northeast region, which was in contrast with Azam et al., (2018a) [26]. Jang (2018) [27] projected that the 12-month drought intensity would increase in Daegu, which is in contrast with this study; however, we noticed that in the 1-month drought data, the duration and frequency were high for Daegu [27].
We compared the SPI and SPEI, and it was found that the SPEI performed better than the SPI, which is consistent with a study performed by Sung et al., (2020) [29], because SPEI incorporates an additional evapotranspiration parameter for drought identification. Jang (2018) [27] also discussed that the high temperature influences the evapotranspiration, leading to the possibility of drought conditions, which cannot be detected only with precipitation, as in case of SPI.

Conclusions
In this study, we evaluated the spatial and temporal changes of drought in South Korea using SPEI and SPI for the period of 1979-2020 using 70 meteorological stations. We analyzed the rainfall variability and found the highest variability for the Geoje met station and the lowest for Yeongcheon. It was observed that the northern part of the country received less rainfall than the southern part. The western part of the country also received more rainfall due to the mountain ranges in the eastern part of the country, which cause foehn warming phenomena. We calculated SPI and SPEI for 1, 6, and 12-month time series for dry and wet conditions and found that 1982, 1988, 1994, 2001 2008, 2015, and 2017 were the dry years in the 12-month time series, according to both indices. In the 6-month time series the highest number of mixed dry events (moderate, severe, and extreme drought) were found in the period of 2013-2017. In the one-month time series, a significant positive trend was observed for the months of September, October, and December, while a significant negative trend was detected in the months of January, March, April, May, June, and August. It was noticed that October was the wettest month, while March and June were the driest months in the one-month time series. In the 12-month time series, significant dry conditions were observed for Jeonju and Daegwallyeong met stations, and significant wet conditions were observed for Seongsan met station. The climatic variability needs longer time data records, but, unfortunately, consistent data for South Korea was only available from 1979 onwards. Trend analysis based on longer temporal data records will provide more authentic results. However, this study provided a reference to evaluate the different types of droughts for the long term. There is a dire need to formulate proper strategies and make policies to cope with the changing climate and its related disasters. We plan to study different types of droughts in more detail using future projection data (GCM) of precipitation and temperature. We also plan to study the spatiotemporal changes in aridity and its influence of land use using past and future projection data.  Data Availability Statement: The data acquired from Korean Meteorological Administration which can not be provided openly. The data can be requested from the first author.