Effect of Climate Change on Annual Precipitation in Korea Using Data Screening Techniques and Climate Change Scenarios

Precipitation is essential for understanding hydrological processes and identifying the characteristics that must be considered to protect human lives and property from natural disasters. Hydrological analyses assume that precipitation shows stationarity. However, because of the recent changes in climate, the stationarity of climate data has been widely debated, and a need has arisen to analyze its nonstationary nature. In this study, we reviewed a method to analyze the stationarity of annual precipitation data from 37 meteorological stations that have recorded data for more than 45 years. Six stations that showed abnormal precipitation during the previous year were selected to evaluate the normality of future precipitation. The results showed that a significant trend was present in four out of 37 stations with unstable precipitation in 22 stations and persistent precipitation in 4 stations. The stationarity analysis of future annual precipitation using climate change scenarios suggested that no trend would be present in 11 stations and that unstable precipitation would be present in six stations. Persistent precipitation was identified in four stations. A comparison between the historical and predicted precipitation data conducted with the climate change scenarios showed that an increasing number of stations presented nonstationarity. Therefore, both stationarity and nonstationarity should be considered when performing hydrological analyses using annual precipitation data in Korea. Accordingly, prior to conducting any such analyses, the effect of climate change on annual precipitation should also be considered.


Introduction
Climate change is altering precipitation patterns in South Korea during summer by increasing the intensity and frequency of precipitation events. In particular, an increase in precipitation has been observed in August, following the monsoon season. Ref. / Oh et al., suggested that the occurrence of precipitation events with ≥50 mm of rainfall per hour has increased 2.37-fold on average each year since 1998, [1]. Thus, it can be expected that the present and future climate would show nonstationarity, making it important to infer various statistics from current hydrological data. However, conventional hydrological data analyses assume stationarity and use data from extreme precipitation events to estimate a design scale of hydro-infrastructure for each of the frequencies. As climate change has altered the characteristics of extreme precipitation events, it is important to employ hydrological statistical analyses that do not assume stationarity to identify changes in extreme precipitation events. In particular, statistical probability estimations of extreme events from the analysis of the hydrological using the RCP 8.5 scenario, the temperature is expected to rise by about 4.6 • C and for precipitation, it is expected to increase by about 0.5 mm/day −1 , [12].
Moreover, a wide range of analytical methods and modeling techniques were applied, depending on the purpose of the study, to produce accurate and meaningful results.
Furthermore, a range of analytical methods and modeling techniques were applied depending on the objectives of the study for producing more accurate and meaningful results. However, the aforementioned studies either analyzed trends and stationarity separately or used a trend analysis to determine stationarity. This indicates that there were limitations with regard to the way in which these authors determined stationarity without considering persistence and stability, which also influences stationarity. Therefore, when analyzing the stationarity of hydrological time series data, persistence and stability must also be considered in the trend analysis.

Overview
In this study, we addressed the necessity of determining stationarity based on a time series analysis of the precipitation data of South Korea by considering persistence, trends, and stability. The analysis was conducted using annual precipitation data from 1971 to 2018 from 37 weather stations across the country. Stationarity was defined as the state in which persistence and trends were not observed yet stability existed. In addition, to evaluate climate variability with respect to annual precipitation, the coefficients of variation (CV) for all the stations were used. The results of a previous study [13] that used data from 1971 to 1995 were compared with those of the present study to evaluate change in the time series of each hydrological dataset.
The spatial disaggregation with quantile delta mapping (SDQDM) method interpolates with the inverse distance weighted using the values of the surrounding global climate model (GCM) grids before performing a bias correction, and then the rate of change for the future period is estimated from the original GCM for each quantile of the data. This statistical method maintains convenience. In addition, the study used climate change scenarios with SDQDM, a model for statistical downscaling, targeting the stations for which nonstationarity was observed. With this, the study aimed to identify changes in the characteristics of future annual precipitation caused by climate change by analyzing the stationarity of such data. Figure 1 shows the process followed in this study.
Atmosphere 2020, 11, x FOR PEER REVIEW 3 of 27 using the RCP 8.5 scenario, the temperature is expected to rise by about 4.6 °C and for precipitation, it is expected to increase by about 0.5 mm/day , [12]. Moreover, a wide range of analytical methods and modeling techniques were applied, depending on the purpose of the study, to produce accurate and meaningful results.
Furthermore, a range of analytical methods and modeling techniques were applied depending on the objectives of the study for producing more accurate and meaningful results. However, the aforementioned studies either analyzed trends and stationarity separately or used a trend analysis to determine stationarity. This indicates that there were limitations with regard to the way in which these authors determined stationarity without considering persistence and stability, which also influences stationarity. Therefore, when analyzing the stationarity of hydrological time series data, persistence and stability must also be considered in the trend analysis.

Overview
In this study, we addressed the necessity of determining stationarity based on a time series analysis of the precipitation data of South Korea by considering persistence, trends, and stability. The analysis was conducted using annual precipitation data from 1971 to 2018 from 37 weather stations across the country. Stationarity was defined as the state in which persistence and trends were not observed yet stability existed. In addition, to evaluate climate variability with respect to annual precipitation, the coefficients of variation (CV) for all the stations were used. The results of a previous study [13] that used data from 1971 to 1995 were compared with those of the present study to evaluate change in the time series of each hydrological dataset.
The spatial disaggregation with quantile delta mapping (SDQDM) method interpolates with the inverse distance weighted using the values of the surrounding global climate model (GCM) grids before performing a bias correction, and then the rate of change for the future period is estimated from the original GCM for each quantile of the data. This statistical method maintains convenience. In addition, the study used climate change scenarios with SDQDM, a model for statistical downscaling, targeting the stations for which nonstationarity was observed. With this, the study aimed to identify changes in the characteristics of future annual precipitation caused by climate change by analyzing the stationarity of such data. Figure 1 shows the process followed in this study.   The data were obtained from observations of the hydrological time series, and the trend and periodicity characteristics are components. In addition, the time series can be roughly classified into two groups: one for deterministic components and the other for stochastic components. The composition of the time series is shown in Equation (1): where X t denotes time series data, D t is the definite component, T t is the trend component, P t is the periodic component, and t is the stochastic component (randomness). In general, a hydrological time series shows either stationarity or nonstationarity. When its average and variance do not change with time, the time series is considered to show stationarity. The opposite holds true for a time series that shows nonstationarity, wherein a trend or periodicity is present or a lack of homogeneity is observed in the hydrological data [13].

Estimation of Persistence
(1) Cumulative Deviation Cumulative deviation is a parametric test to assess whether the average value is different before and after a random time during an observation period. As it assumes that the data are normally distributed, this test aims to evaluate whether the average changes after the selected observation, as shown in Equations (2) and (3): where µ is the average before the change, ∆ is the change in the average, and E(x i ) is the expected average value before and after the time of change. The Q statistic, which is based on the average, is shown in Equation (4): Given that the change at the time of S * * k has the maximum test value, a negative value of S * k indicates that the average of the former part is higher than that of the later part [1]. S * k and S * * k are defined in Equations (5) and (6): where n is the number of data points, x i is the time series data value, x is the average of the sample group, s is the standard deviation of the sample group, s 2 is the variance of the sample group, S * k is cumulative deviation, and S * * k is the standardized cumulative value of annual precipitation. (2) Serial Correlation Coefficient Hydrological persistence indicates a condition in which the hydrological quantities that constitute a hydrological time series are not in a random and independent pattern but continuously maintain similar size. This feature of the time series should be tested prior to a frequency analysis or stochastic simulation. If time series data show persistence, they cannot be used in a frequency analysis, and thus a serial correlation coefficient should be used to determine persistence. The serial correlation coefficient is a statistic that evaluates the persistence of a time series. When a serial correlation of variables measures zero, no relationship is present between the observed values, which are independent of one another. However, when the serial correlation is skewed towards one variable, the observed values are correlated within the series, and the current value has an influence on future values. In this case, the consecutive correlation coefficient should be used and is defined in Equation (7) [14]: where x i is an observe value, x i+1 is the next value, x is the average of a time series, and n is the number of data points. The equation for the 5% significance interval is shown in Equation (8):

Spearman Rank Correlation
Trends in hydrological data are influenced by a combination of natural and artificial factors. The data will either have an increasing or decreasing trend with regard to the average and variance of the time series. Therefore, when statistics, such as the average and variance, remain constant, a time series can be considered to present stationarity. A trend can be identified with listed time series data; however, a statistical analysis is necessary if a clear trend is not observed. When a trend is present in annual and seasonal data, strong auto-correlation is generally present among data points. The statistical methods used to analyze trends include the T-test, Hotelling-Pabst test, nonlinear trend test, Mann-Kendall test, and Sen test. This study used a Spearman rank correlation analysis as it is widely used to evaluate correlations between two variables using a rank-based nonparametric scale [15]. The Pearson correlation analysis is most effective when the data are normally distributed, whereas the Spearman rank correlation has a statistical efficiency of 70% or higher for all possible variables. The Spearman rank correlation equation is shown in Equation (9): where n is the number of data points, D i for Kx i − Ky for the time sequence, Kx i for the data sequence with the size, Ky is the time sequence of the raw data of Kx i , and x is the data time sequence [13].

Simple T-Test and Simple F-Test
The stability of hydrological time series data can be tested by assessing if statistics, such as the average and variance of a time series, change over time. If the average, variance, or other statistics vary with time, the data set is said to present nonstationarity. In general, the methods to analyze the stability of time series data include the Mann-Whitney test, sign test, Abbe test, simple T-test, and simple F-test. This study used a simple T-test and simple F-test to analyze stability, owing to their ease and simplicity of use. The T-test is a parametric statistical analysis that is used to evaluate the difference between the averages of two groups, whereas the F-test is used to determine whether a difference is present between variances of two groups and is shown in Equation (11) [13]: Atmosphere 2020, 11, 1027 6 of 27 where s 2 is the variance. A time series is considered to show stationarity once the following condition is satisfied: v 1 , v 2 , 2.5% < F t < v 1 , v 2 , 97.5%. In this condition, v 1 = (n 1 −1) and v 2 = (n 2 −1) are the degrees of freedom for the numerator and denominator, respectively.
The T-test for average stationarity is shown in Equation (12): where x 1 and x 2 are the average values, s 2 1 and s 2 2 are the variances, and n 1 and n 2 are the number of data points for each subset.

Station Used for Analysis
The hydrological time series data used for the analysis required testing. However, in the case of the data from South Korea, the test results showed a low confidence interval because of the short observation period. Therefore, in this study, data were collected for each of the 37 stations located across the country that began recording observations in 1971. These stations recorded data for at least 45 years, and the characteristics of annual precipitation time series data were assessed. Figure 2 shows the locations of the 37 stations selected for this study.
Atmosphere 2020, 11, x FOR PEER REVIEW 6 of 27 where and are the average values, and are the variances, and and are the number of data points for each subset.

Station Used for Analysis
The hydrological time series data used for the analysis required testing. However, in the case of the data from South Korea, the test results showed a low confidence interval because of the short observation period. Therefore, in this study, data were collected for each of the 37 stations located across the country that began recording observations in 1971. These stations recorded data for at least 45 years, and the characteristics of annual precipitation time series data were assessed. Figure 2 shows the locations of the 37 stations selected for this study.
To identify the trends with regard to changes in annual precipitation for South Korea, we analyzed the characteristics of the time series for the 37 weather stations taking into account the spatial distribution of the annual precipitation averages and standard deviations of the data from these stations. The annual precipitation average was highest for the Jeju area (1905 mm/year for Seogwipo and 1850 mm/year for Seongsan) and lowest for the Daegu station (1052 mm/year). In addition, the study calculated the CVs to evaluate the variability in annual precipitation. The CVs were obtained by dividing the standard deviation by the numerical mean. A low CV indicated that the annual data had low variability and high stability. The annual precipitation CV was found to be lowest for the Sokcho station and highest for the Seongsan station, indicating that the annual precipitation data of the Sokcho station was stable and least variable while that of the Seongsan station was unstable and highly variable. Figure 3 shows the spatial distribution of averages, standard deviations, and CVs of the annual precipitation data in this study.  To identify the trends with regard to changes in annual precipitation for South Korea, we analyzed the characteristics of the time series for the 37 weather stations taking into account the spatial distribution of the annual precipitation averages and standard deviations of the data from these stations. The annual precipitation average was highest for the Jeju area (1905 mm/year for Seogwipo and 1850 mm/year for Seongsan) and lowest for the Daegu station (1052 mm/year). In addition, Atmosphere 2020, 11, 1027 7 of 27 the study calculated the CVs to evaluate the variability in annual precipitation. The CVs were obtained by dividing the standard deviation by the numerical mean. A low CV indicated that the annual data had low variability and high stability. The annual precipitation CV was found to be lowest for the Sokcho station and highest for the Seongsan station, indicating that the annual precipitation data of the Sokcho station was stable and least variable while that of the Seongsan station was unstable and highly variable. Figure 3 shows the spatial distribution of averages, standard deviations, and CVs of the annual precipitation data in this study.

Climate Change Scenarios
The GCM, a numerical climate model, has several types and produces different simulation results depending on the type of input data used. Meteorological data predicted by the GCM are statistically downscaled for practical use. The KMA operates a certification system to encourage the Atmosphere 2020, 11, 1027 9 of 27 use of diverse scenarios. The system recommends 13 scenarios at a regional level that are produced with the combination of GCMs scenarios and statistical down-scaling (Table 1). For such statistical downscaling, SDQDM has been recommended as it was designed to preserve the long-term trends of climate models and conducts a bias correction [16]. In this study, 13 regional scale scenarios produced by SDQDM, a statistically detailed technique, were used with the GCM scenarios provided by the Standard Scenario Certification System of the Korea Meteorological Administration. Figure 4 shows the climate change scenarios. Figure 5 shows the spatial distribution of averages, standard deviations, and CVs of the future annual precipitation data in this study. The results showed that the annual precipitation average will be highest for the Seogwipo station (2445 mm/year) and lowest for the Daegu station (1023 mm/year). The standard deviation was the highest in Seogwipo (733), the lowest in Ulleung (201), the highest in CV (Gunsan (0.32)), and the lowest in Pohang (0.21). The GCM, a numerical climate model, has several types and produces different simulation results depending on the type of input data used. Meteorological data predicted by the GCM are statistically downscaled for practical use. The KMA operates a certification system to encourage the use of diverse scenarios. The system recommends 13 scenarios at a regional level that are produced with the combination of GCMs scenarios and statistical down-scaling (Table 1). For such statistical downscaling, SDQDM has been recommended as it was designed to preserve the long-term trends of climate models and conducts a bias correction [16]. In this study, 13 regional scale scenarios produced by SDQDM, a statistically detailed technique, were used with the GCM scenarios provided by the Standard Scenario Certification System of the Korea Meteorological Administration. Figure 4 shows the climate change scenarios. Figure 5 shows the spatial distribution of averages, standard deviations, and CVs of the future annual precipitation data in this study. The results showed that the annual precipitation average will be highest for the Seogwipo station (2445 mm/year) and lowest for the Daegu station (1023 mm/year). The standard deviation was the highest in Seogwipo (733), the lowest in Ulleung (201), the highest in CV (Gunsan (0.32)), and the lowest in Pohang (0.21).

Performance of Climate Change Scenarios
To verify the performance of MME from the 13 climate change scenarios of GCMs, the study draws a box plot, using spatially calculated averages for monthly precipitation and mean temperature data provided by 80 weather stations for the reference period for each of the models. In addition, the coefficient of variation (CV) between models was calculated by station point for annual precipitation and average temperatures of the late 21st century (2071-2100).
The study statistically down-scaled the 13 national standard scenarios, using the SDQDM of the AIMS provided by the Asia Pacific Economic Cooperation Climate Center (APCC). In this chapter, the study compares probability density functions (PDFs) and spatial distribution between MME, (precipitation and temperature), which is down-scaled for the period of reference climate, and observation data. This aims to verify the applicability of statistical down-scaling results, for North Korea's climate observation.
For annual precipitation, the PDFs were compared between the observation data and before application or after application of the SDQDM. This is well matched to the bell-shaped, normal distribution for the average value, while variables were gathered around the median. Comparing average values (mode and median) and variances between before down-scaling (blue line in Figure 6), and after down-scaling (red line in Figure 6), the study found that they were well matched to those of observation data after down-scaling. For quantitative results, the average value of OBS is 1116 mm, the average of the raw scenario is 1196, the average of the downscaled scenario is 1162 mm, and it is confirmed that the observation and the after down-scaling are similar. The result of comparing spatial distribution between before/after down-scaling suggested that the latter simulated the spatial distribution very close (Figure 7a,c) to that of observation data (Figure 7a). To verify the performance of MME from the 13 climate change scenarios of GCMs, the study draws a box plot, using spatially calculated averages for monthly precipitation and mean temperature data provided by 80 weather stations for the reference period for each of the models. In addition, the coefficient of variation (CV) between models was calculated by station point for annual precipitation and average temperatures of the late 21st century (2071-2100).
The study statistically down-scaled the 13 national standard scenarios, using the SDQDM of the AIMS provided by the Asia Pacific Economic Cooperation Climate Center (APCC). In this chapter, the study compares probability density functions (PDFs) and spatial distribution between MME, (precipitation and temperature), which is down-scaled for the period of reference climate, and observation data. This aims to verify the applicability of statistical down-scaling results, for North Korea's climate observation.
For annual precipitation, the PDFs were compared between the observation data and before application or after application of the SDQDM. This is well matched to the bell-shaped, normal distribution for the average value, while variables were gathered around the median. Comparing average values (mode and median) and variances between before down-scaling (blue line in Figure  6), and after down-scaling (red line in Figure 6), the study found that they were well matched to those of observation data after down-scaling. For quantitative results, the average value of OBS is 1116 mm, the average of the raw scenario is 1196, the average of the downscaled scenario is 1162 mm, and it is confirmed that the observation and the after down-scaling are similar. The result of comparing spatial distribution between before/after down-scaling suggested that the latter simulated the spatial distribution very close (Figure 7a,c) to that of observation data (Figure 7a).

Hydrological Persistence Test
A persistence test was conducted using the cumulative deviation to evaluate annual precipitation data. When persistence exists, data averages show inconsistency and randomness. Therefore, persistence in data is indicative of nonstationarity. After evaluating persistence using a cumulative deviation analysis, one station showed increasing annual precipitation since 1978. Out of the 37 selected stations, 19 and 17 stations showed increasing trends in annual precipitation from 1982 to 1989 and from 1994 to 1997, respectively. The study also used the serial correlation coefficient to evaluate persistence. In a previous study, persistence was only explained for the Daegwallyeong station [13]. However, the results of this study suggested that four stations (i.e., Daegwallyeong, Ulleung, Seogwipo, and Ganhwa) showed persistence, as summarized in Table 2. In addition, the graphs of Figures 8 and 9 show the results of the persistence analysis using the serial correlation coefficient. Each graph shows an inflection point, annual precipitation, the averages before and after the inflection point, and a standardized cumulative value of annual precipitation. The graphs in Figure 8 show the persistence analysis results for the Daegwallyeong, Ulleung, Seogwipo, and Ganhwa stations. In contrast, the graphs in Figure 9 show the results for the Busan, Gwanju, Haenam, and Jeju stations, which did not show persistence. The map shown in Figure 10 indicates the stations with and without persistence.

Hydrological Persistence Test
A persistence test was conducted using the cumulative deviation to evaluate annual precipitation data. When persistence exists, data averages show inconsistency and randomness. Therefore, persistence in data is indicative of nonstationarity. After evaluating persistence using a cumulative deviation analysis, one station showed increasing annual precipitation since 1978. Out of the 37 selected stations, 19 and 17 stations showed increasing trends in annual precipitation from 1982 to 1989 and from 1994 to 1997, respectively. The study also used the serial correlation coefficient to evaluate persistence. In a previous study, persistence was only explained for the Daegwallyeong station [13]. However, the results of this study suggested that four stations (i.e., Daegwallyeong, Ulleung, Seogwipo, and Ganhwa) showed persistence, as summarized in Table 2. In addition, the graphs of Figures 8  and 9 show the results of the persistence analysis using the serial correlation coefficient. Each graph shows an inflection point, annual precipitation, the averages before and after the inflection point, and a standardized cumulative value of annual precipitation. The graphs in Figure 8 show the persistence analysis results for the Daegwallyeong, Ulleung, Seogwipo, and Ganhwa stations. In contrast, the graphs in Figure 9 show the results for the Busan, Gwanju, Haenam, and Jeju stations, which did not show persistence. The map shown in Figure 10 indicates the stations with and without persistence.

Hydrological Trend Test
To evaluate the annual precipitation data, we used the Spearman rank correlation test to analyze possible trends. This study also conducted a two-tailed statistical test with significance at the 5% level for each of the dataset. An increasing or decreasing trend with regard to the average and variance may be present in a time series. Therefore, when the statistical characteristics such as average and variance are constant with time, the data cannot be considered stationary. Therefore, when statistics, such as the average and variance, are constant over time, the data are considered to show stationarity. Table 3 shows the results of the Spearman rank correlation analysis. Out of the 37 stations, the annual precipitation data from the Uleungdo, Seongsan, Seogwipo, and Inje stations showed increasing trends with variations between stations. Among the four stations, Seongsan station showed the greatest trend. In a past study [13], data from 1972 to 1995 were used to analyze the trends of Daegwallyeong station. However, in this study, no trend was observed for this station as the decrease in annual precipitation since 1994 was offset by the current increasing trend. The map in Figure 11 shows the results of the trend analysis.
Atmosphere 2020, 11, x FOR PEER REVIEW 15 of 27 Figure 11. Trend analysis result for each station.

Hydrological Stable Test
A simple T-test and F-test were used to test the stability of annual precipitation. Two-tailed tests with significance at the 5% level were conducted for each dataset. The data in which statistics, such as the average and variance, were not consistent were considered nonstationary. Table 4 summarizes

Hydrological Stable Test
A simple T-test and F-test were used to test the stability of annual precipitation. Two-tailed tests with significance at the 5% level were conducted for each dataset. The data in which statistics, such as the average and variance, were not consistent were considered nonstationary. Table 4 summarizes the simple T-test and F-test results. For the simple T-test, 23 of the 37 stations showed unstable annual precipitation, including the Uljin, Cheongju, Daejeon, Chupungnyeong, Pohang, Gunsan, Daegu, Jeonju, Busan, and Ulsan stations. The simple F-test results show that the annual precipitation data from two stations, Seongsan and Haenam, were unstable. According to [13], out of the stations that showed unstable data, the Uljin, Cheongju, Daejeon, and Chupungnyeong stations showed the greatest changes in trends, and data from the Daegwallyeong, Ganghwa, and Imsil stations were found to be stable. However, in this study, data from all stations, except Imsil, were unstable, whereas the data from 23 stations, including Uljin, Cheongju, and Imsil, were analyzed for instability. Figure 12 shows the simple T-test and F-test results.

Estimation of Hydrological Stationarity
As shown in Table 5, data from the 14 stations that showed stability without trends or persistence were deemed to show stationarity. The data from 23 stations that showed abnormality in two out of the three features (i.e., persistence, trends, and stability) were deemed nonstationary. According to [13], data from the Daegwallyeong, Uleungdo, Yangpyeong, and Jecheon stations showed nonstationarity. However, in this study, we found that the data from 23 stations, including Uleungdo station, showed nonstationarity. Moreover, the data from Daegwallyeong station were considered to show stationarity. Figure 13 shows the results of the stationarity analysis for the annual precipitation data of the 37 stations. Table 5. Estimation of the stationarity of annual precipitation data.

No.
Station Stationarity   90  Sokcho  NP  S  NT  O  156  Gwangju  NP  US NT  X  100  Daegwanryeong P  S  NT  X  159  Busan  NP  US NT  X  101  Chuncheon  NP  S  NT  O  162  Tongyeong NP  US NT  X  105  Gangneung  NP  S  NT  O  165  Mokpo  NP  US NT  X  108  Seoul  NP  S  NT  O  168  Yeosu  NP  US NT  X  112  Incheon  NP  S  NT  O  170  Wando  NP  US NT  X  114  Wonju  NP  S  NT  O  184  Jeju  NP  S  NT  O  115  Ulleung  P  S  T  X  188  Sungsan  NP  US  T  X  119 Suwon

Selected Stations
We selected 25 stations that showed nonstationarity based on the results of the stationarity analysis using current annual precipitation data. With the exception of Chupungnyeong and

Selected Stations
We selected 25 stations that showed nonstationarity based on the results of the stationarity analysis using current annual precipitation data. With the exception of Chupungnyeong and Seongsan stations, for which climate change scenarios were not available, all stations were used for the stationary analysis using future annual precipitation scenarios for South Korea with SDQDM, a statistical downscaling model. Figure 14 shows the spatial distribution of the observations for this study. The blue dot is the entire station, and the red dot is the station analyzed using climate change scenario data.

Hydrological Persistence Test Using Climate Change Scenarios
To determine the hydrological persistence of future annual precipitation data, we used a persistence analysis with the cumulative deviations and serial-correlation coefficients of future annual precipitation data of the 23 aforementioned stations. [13] suggested that persistence was only present for Daegwallyeong station. However, in this study, the analysis of the current annual precipitation indicated the presence of persistence in the Jeongeup, Uleungdo, Seogwipo, and Seongsan stations, whereas the analysis using climate change scenarios suggested that this feature did not appear in any station. Figure 15 shows the cumulative graphs of the stations for which the data did not show persistence, namely the Jeongeup, Daejeong, Seogwipo, and Jeonju stations. Table  6 summarizes the results of the serial correlation coefficient analysis, and Figure 16 shows the results of the persistence analysis.

Hydrological Persistence Test Using Climate Change Scenarios
To determine the hydrological persistence of future annual precipitation data, we used a persistence analysis with the cumulative deviations and serial-correlation coefficients of future annual precipitation data of the 23 aforementioned stations. [13] suggested that persistence was only present for Daegwallyeong station. However, in this study, the analysis of the current annual precipitation indicated the presence of persistence in the Jeongeup, Uleungdo, Seogwipo, and Seongsan stations, whereas the analysis using climate change scenarios suggested that this feature did not appear in any station. Figure 15 shows the cumulative graphs of the stations for which the data did not show persistence, namely the Jeongeup, Daejeong, Seogwipo, and Jeonju stations. Table 6 summarizes the results of the serial correlation coefficient analysis, and Figure 16 shows the results of the persistence analysis.

Hydrological Trend Test Using Climate Change Scenarios
This study used a Spearman rank correlation analysis to analyze trends in annual precipitation using climate change scenarios. A two-tailed statistical test with significance at the 5% level was used for each dataset. Table 6 shows the trends in annual precipitation data for 11 out of the 23 stations, including the Busan, Seogwipo, Yangpyeong, and Jeongeup stations. In a previous study [13], data from 1972 to 1995 was used to analyze the trends of Daegwallyeong station. However, in this study, the analysis of this data indicated no significant trend for this station as the decrease in the annual precipitation since 1995 was offset by the current increasing trend. A comparison between the current data and that of the Future 3 scenario (2071 to 2100), the latter of the climate change scenarios, indicated that Seogwipo station showed the highest rate of increase (28.7%). The lowest rate of decrease was present at Daegu station (10.1%). In addition, Busan and Yangpyeong stations showed increasing trends of 8.6% and 3.8%, respectively. In addition, Uljin and Inje stations showed decreasing trends of 1.3% and 2.3%, respectively. These results suggest that variations in the trends of annual precipitation of the target stations. The results of the trend analysis are summarized in Table 7, while Figure 17 provides the graphs of six stations and depicts changes in current and future annual precipitation. Figure 18 shows the results of the trend analysis.

Hydrological Stability Analysis Using the Climate Change Scenarios
For the annual precipitation stability analysis using climate change scenarios (2011-2099), we used simple T-test and F-test analyses. With the results, we used a two-tailed test with significance at the 5% level for each data. Table 8 shows the Simple T-test and F-test results. The Uljin, Busan, Seogwipo, Yangpyeong, Inje, and Jeongeup stations showed instability. [13] indicated that the data showed instability for the Daegwallyeong, Ganghwa, and Imsil stations. In this study, the results suggested that the data were unstable for all 21 stations (Table 7), which was also supported by the simple T-test and F-test results. Figure 19 shows the results of the stability test.

Hydrological Stability Analysis Using the Climate Change Scenarios
For the annual precipitation stability analysis using climate change scenarios (2011-2099), we used simple T-test and F-test analyses. With the results, we used a two-tailed test with significance at the 5% level for each data. Table 8 shows the Simple T-test and F-test results. The Uljin, Busan, Seogwipo, Yangpyeong, Inje, and Jeongeup stations showed instability. [13] indicated that the data showed instability for the Daegwallyeong, Ganghwa, and Imsil stations. In this study, the results suggested that the data were unstable for all 21 stations (Table 7), which was also supported by the simple T-test and F-test results. Figure 19 shows the results of the stability test.

Hydrological Stability Analysis Using the Climate Change Scenarios
For the annual precipitation stability analysis using climate change scenarios (2011-2099), we used simple T-test and F-test analyses. With the results, we used a two-tailed test with significance at the 5% level for each data. Table 8 shows the Simple T-test and F-test results. The Uljin, Busan, Seogwipo, Yangpyeong, Inje, and Jeongeup stations showed instability. [13] indicated that the data showed instability for the Daegwallyeong, Ganghwa, and Imsil stations. In this study, the results suggested that the data were unstable for all 21 stations (Table 7), which was also supported by the simple T-test and F-test results. Figure 19 shows the results of the stability test.

Estimation of Hydrological Stationarity Using Climate Change Scenarios
We analyzed persistence, trends, and stability to determine the stationarity of annual precipitation data based on climate change scenarios. The results are summarized in Table 9 for the 23 stations that showed nonstationary with current data, including the Daegwallyeong, Uleungdo, Yangpyeon, and Jecheon stations. However, the analysis of future data found that nonstationary was Figure 19. Stability test results using the climate change scenarios.

Estimation of Hydrological Stationarity Using Climate Change Scenarios
We analyzed persistence, trends, and stability to determine the stationarity of annual precipitation data based on climate change scenarios. The results are summarized in Table 9 for the 23 stations that showed nonstationary with current data, including the Daegwallyeong, Uleungdo, Yangpyeon, and Jecheon stations. However, the analysis of future data found that nonstationary was present in 23 stations, including those of the current study. Figure 20 shows the results of the stationarity analysis.

Conclusions
The study collected annual precipitation data from 37 weather stations across Korea that had collected data for 45 years or longer with the aim to analyze the stationarity of the data. We also conducted persistence, trend, and stability analyses and then screened out stations that showed nonstationarity of annual precipitation data. For the remaining stations, we used climate change scenarios and analyzed the stationarity of future annual precipitation data. The results can be summarized as follows.
(1) After analyzing the persistence of data from 37 stations, annual precipitation at Haenam station was found to have increased since 1978. In addition, 19 stations, including Gangneung, Ulsan, Jeonju, Ganghwa, and Seongsan, showed increasing precipitation trends from 1982 to 1989, whereas

Conclusions
The study collected annual precipitation data from 37 weather stations across Korea that had collected data for 45 years or longer with the aim to analyze the stationarity of the data. We also conducted persistence, trend, and stability analyses and then screened out stations that showed nonstationarity of annual precipitation data. For the remaining stations, we used climate change scenarios and analyzed the stationarity of future annual precipitation data. The results can be summarized as follows.
(2) The trend analysis for the annual precipitation of the target stations suggested that no trends were present for most regions. However, the Uleungdo, Seongsan, Seogwipo, and Inje stations showed some evidence of trends. Among these, Seongsan station showed the distinct trend of increasing annual precipitation.
(3) The stability analysis of 23 stations, including that of the Uljin, Cheongju, Daejeon, and Chupungnyeong stations, showed that the annual precipitation data was unstable. In particular, the Seogwipo station registered the highest annual precipitation (more than 2000 mm) from 2009 to 2012 and 2014 to 2016. Therefore, it can be concluded that the influence of such periods caused data instability.
(4) Based on the results of the persistence, trend, and stability analyses, stationarity was detected. In particular, 14 stations showed stationarity, whereas 25 stations showed nonstationarity. The stations with nonstationarity accounted for 67% of all stations. Most of the stations in the central region, except the Inje, Hongcheon, and Yangpyeong stations, and most of the stations in the southern region, except the Jinju and Jeju stations, showed stationarity. The principle reason for the nonstationarity of these stations in the southern region was the influence of typhoons that led to persistence and instability in the annual precipitation data.
(5) The stationarity analysis for future annual precipitation using climate change scenarios showed that out of 23 stations, 12 stations showed trends, 23 stations had unstable data, and no station showed persistence as climate change caused changes in annual precipitation characteristics. Therefore, future hydrological analysis should consider the influence of climate change on the characteristics of hydrological time series.
In addition, the results obtained were different from those of a previous study [13], wherein nonstationarity was observed for multiple stations, including Daegwallyeong, Uleungdo, Yangpyeong, and Jecheon.
However, in this study, 25 stations, including Uleungdo, showed nonstationarity. Moreover, the results of the stationarity analysis of future annual precipitation using climate change scenarios suggests that climate change will influence the characteristics of hydrological time series data. For example, the stations that showed a trend in the present also showed persistence without the trend. This underlines the necessity of reconsidering present hydrological analyses with the variable characteristics of annual precipitation. Therefore, in hydrological analyses that use annual precipitation in South Korea, nonstationarity needs to be considered beyond the current assumption of stationarity. In addition, current analyses should be reviewed by considering the current characteristics of annual precipitation and the influence of climate change on future annual precipitation. Moreover, it is also necessary to statistically identify climate characteristics in diverse conditions by analyzing the correlations among a range of hydro-meteorological factors with regard to annual precipitation in the future as time series characteristics with respect to trends, stability, and persistence are likely to be affected by climate change. Furthermore, the results of this study are expected to contribute to long-term plans regarding the preventive management measures that must be undertaken with regard to natural phenomena, including those for flooding, as well as additional studies that are conducted on time series that consider the influence of climate change on diverse hydro-meteorological factors.