The Impacts of Water Cycle Components on Streamflow in a Changing Climate of Korea: Historical and Future Trends

This paper investigates the historical and future trends in water balance components and their impacts on streamflow. The trend analyses were applied to the daily climatic and hydrologic variables from 109 subbasins in Korea during the historical period and future period obtained by a multimodel ensemble of 13 global circulation models (GCMs) of the Coupled Model Intercomparison Project, Phase 5 (CMIP5). A calibrated hydrologic model, the precipitation-streamflow modeling system (PRMS) model, was applied to obtain hydrologic data. The results revealed apparent trends in streamflow, with increases in spring and decreases in the other seasons during the historical period. The reduction (or increase) in the amount of streamflow was counterbalanced by the reduction (or increase) in precipitation, groundwater, and soil moisture, which was mainly impacted by the increase (or reduction) in actual evapotranspiration. However, opposite trends are projected for the future period for streamflow and water cycle components, in which spring and winter are projected to have increasing trends mostly counterbalanced by the decreasing trends in precipitation and groundwater. The reasons for the reduction in streamflow include elevated evapotranspiration compared to precipitation, reduced soil moisture, and a significant decrease in groundwater recharge. In addition, the results of the seasonal variability among basins revealed higher variability in summer for the historical period and in winter for the future period, with maximum variability in the Sumjin River basin, indicating that streamflow fluctuated more strongly in the Sumjin River basin during the historical and future periods.


Introduction
Water demand is increasing worldwide; however, a reliable water supply is not always available. This concern becomes more challenging based on the increasing water demand due to population growth and climate change. Recent studies reported almost negative effects on water resources [1,2]. Rivers provide more than 50% of water demands; however, their streamflow volume is highly impacted by precipitation and temperature variations [3]. Therefore, it is necessary to understand how the water cycle components affect streamflow or are expected to impact streamflow in the future to overcome probable water scarcity by intelligent water resource management.
There is general agreement regarding the increasing trend of temperature by using different climate change scenarios anywhere during the 21st century [3,4]. However, higher temperatures cause an increase in the capacity of water holding and therefore evaporation capacity in the atmosphere, August, autumn from September to November, and winter from December to February. However, the climatic conditions vary significantly among the seasons and are dependent on both continental and oceanic effects. The temperature in summer ranges from 22-25 • C due to Pacific high pressure, while in winter, it ranges from −5 • C to −3 • C depending on the elevation. The mean annual precipitation generally increases from north to south with approximately 1500 mm in the middle of Korea [20]. Additionally, approximately 70% of the annual precipitation occurs during the summer, which is the rainy season in Korea; winter is the driest season and constitutes less than 10% of the total annual precipitation. Figure 2 shows the mean monthly values for the minimum temperature, maximum temperature, and precipitation in South Korea.

Climate and Hydrologic Data
Daily climate data were used to calibrate the model parameters as well as downscale the GCM scenarios. The observed climate data were extracted from the Ministry of Land, Infrastructure and Transport (MOLIT) and the Korea Meteorological Administration (KMA) from 1966 to 2016. Moreover, the observed daily dam inflow data were extracted from the Korean Water Management Information System (WAMIS) for the six gauged dam basins during the historical period from 1966 to 2016. In addition, the required GIS layers, such as the vegetation type, digital elevation models (DEMs), land use, and soil maps, were obtained from the WAMIS website.
Sustainability 2020, 12, x FOR PEER REVIEW 3 of 21 and oceanic effects. The temperature in summer ranges from 22-25 °C due to Pacific high pressure, while in winter, it ranges from −5 °C to −3 °C depending on the elevation. The mean annual precipitation generally increases from north to south with approximately 1500 mm in the middle of Korea [20]. Additionally, approximately 70% of the annual precipitation occurs during the summer, which is the rainy season in Korea; winter is the driest season and constitutes less than 10% of the total annual precipitation. Figure 2 shows the mean monthly values for the minimum temperature, maximum temperature, and precipitation in South Korea.

Climate and Hydrologic Data
Daily climate data were used to calibrate the model parameters as well as downscale the GCM scenarios. The observed climate data were extracted from the Ministry of Land, Infrastructure and Transport (MOLIT) and the Korea Meteorological Administration (KMA) from 1966 to 2016. Moreover, the observed daily dam inflow data were extracted from the Korean Water Management Information System (WAMIS) for the six gauged dam basins during the historical period from 1966 to 2016. In addition, the required GIS layers, such as the vegetation type, digital elevation models (DEMs), land use, and soil maps, were obtained from the WAMIS website.  and oceanic effects. The temperature in summer ranges from 22-25 °C due to Pacific high pressure, while in winter, it ranges from −5 °C to −3 °C depending on the elevation. The mean annual precipitation generally increases from north to south with approximately 1500 mm in the middle of Korea [20]. Additionally, approximately 70% of the annual precipitation occurs during the summer, which is the rainy season in Korea; winter is the driest season and constitutes less than 10% of the total annual precipitation. Figure 2 shows the mean monthly values for the minimum temperature, maximum temperature, and precipitation in South Korea.

Climate and Hydrologic Data
Daily climate data were used to calibrate the model parameters as well as downscale the GCM scenarios. The observed climate data were extracted from the Ministry of Land, Infrastructure and Transport (MOLIT) and the Korea Meteorological Administration (KMA) from 1966 to 2016. Moreover, the observed daily dam inflow data were extracted from the Korean Water Management Information System (WAMIS) for the six gauged dam basins during the historical period from 1966 to 2016. In addition, the required GIS layers, such as the vegetation type, digital elevation models (DEMs), land use, and soil maps, were obtained from the WAMIS website.

Hydrologic Model
In this study, a semidistributed model, namely precipitation-streamflow modeling system (PRMS), was used. The minimum data of the PRMS model are daily precipitation, the maximum temperatures, and the minimum temperatures [26,27]. The PRMS model simulates the water and the energy balance for each hydrologic response unit (HRU). In fact, the basins were divided into HRUs according to physiographic characteristics, such as slope, elevation, land use, and geology [10]. The sum of the water and energy balance of the HRUs weighted by area defines the daily streamflow [23,25].
The water balance components in the PRMS model consist of precipitation, evapotranspiration, soil moisture, groundwater, and streamflow. In the PRMS model, the amount of precipitation for each HRU is adjusted by a monthly correlation factor. The potential evapotranspiration (PET) is computed by the Hamon method [28,29]. There are different methods for estimating PET; however; Hamon method is simple and only require the temperature as input and compare to other methods provide reliable estimation of PET [30][31][32][33]. The Hamon method is the temperature-based formula that uses the daily mean temperature (T mean ) and total sunshine hours (L d ): e s = 6.108 × exp 17.26939 × T mean T mean + 273. 3 (1) where e s is the saturated vapor pressure (mb), ρ v is the saturated vapor density (g/m 3 ) at T mean ( • C), and C h is the monthly coefficient, which has a value of 0.0040 for July and August and a value of 0.0055 for other months, which was obtained using parameter estimation. Moreover, the actual evapotranspiration (AET) is calculated as a function of PET and soil moisture (SM) conditions, which depend on the soil type.

Climate Change Scenarios
In the present study, first we used 26 CMIP5 GCMs (Table 1) of RCP4.5 for the period between 1976 and 2099. The data of 26 GCMs were collected within 119 • E-135 • E and 29 • N-43 • N and clipped to the Korean peninsula to provide daily meteorological data under the same initial conditions of the observed period (from 1976 to 2005). Using the climate projections at the grid points of each GCM, three statistical downscaling methods were applied to downscale the climate projections into the 60 automated synoptic observing system (ASOS) stations as shown in Figure 1. These GCM scenarios were prepared for the Fifth Assessment Report (AR5), which was done by the APEC (Asia-Pacific Economic Cooperation) Climate Center and the Climate Change Adaptation for Water Resources (CCAW). The procedure for selecting 13 GCMs from 26 GCMs is shown in Figure 3.  Table 1. The description and characteristics of the mentioned GCMs are shown in Ghafouri-Azar and Bae [34]. More information about this procedure was explained in Cannon et al. [35] and Eum and Cannon [13].

Projecting the Hydrologic Data
The five main river basins were divided into 109 HRUs (subbasins) based on the Korean water resources unit map. In this study, we use the PRMS model, whose parameters for all 109 subbasins were calibrated and verified at six gauging stations referring to the statistics of long-term observed streamflow, and the results are shown in Bae et al. [23]. The hydrologic data for each subbasin were obtained by running the PRMS model using daily downscaled data for the period 1976 to 2099 of maximum temperature, minimum temperature, and precipitation. Then, the daily time series of precipitation (P), actual evapotranspiration (AET), soil moisture (SM), groundwater (GW), and streamflow (R) were calculated for the same period and for all 13 GCMs. Next, the multimodel ensemble (MME) of 13 GCMs was calculated by computing the average of the GCMs. Then, the MME was applied for further analysis. The analysis in the present study was performed for the historical period from 1976 to 2005 and the future period from 2071 to 2099. The absolute change in the mean temperature computed by subtracting the mean temperature of the future period from the historical period was compared to the reference period. The relative changes were computed as the percentage

Projecting the Hydrologic Data
The five main river basins were divided into 109 HRUs (subbasins) based on the Korean water resources unit map. In this study, we use the PRMS model, whose parameters for all 109 subbasins were calibrated and verified at six gauging stations referring to the statistics of long-term observed streamflow, and the results are shown in Bae et al. [23]. The hydrologic data for each subbasin were obtained by running the PRMS model using daily downscaled data for the period 1976 to 2099 of maximum temperature, minimum temperature, and precipitation. Then, the daily time series of precipitation (P), actual evapotranspiration (AET), soil moisture (SM), groundwater (GW), and streamflow (R) were calculated for the same period and for all 13 GCMs. Next, the multimodel ensemble (MME) of 13 GCMs was calculated by computing the average of the GCMs. Then, the MME was applied for further analysis. The analysis in the present study was performed for the historical period from 1976 to 2005 and the future period from 2071 to 2099. The absolute change in the mean temperature computed by subtracting the mean temperature of the future period from the historical period was compared to the reference period. The relative changes were computed as the percentage change (%) of the future period with respect to the historical period. The negative change revealed a decrease, and the positive change indicated an increase in the future period compared to the historical period.

Trend Analysis
The Mann-Kendall (MK) statistical test is a nonparametric test to analyze the significance of trends in time series and is mostly applied to identify trends in hydrological variables [36][37][38]. Kendall's statistic S is calculated using Equations (4) and (5): where n is the number of data points, X j and X i are sequential data (j > i), and Sgn (X j − X i ) is a sign function that extracts the sign of X j − X i . The statistic S is expected to be normally distributed with zero mean (E(S) = 0) and standard deviation, which is computed as (Hamed and Rao, 1998): where m is the number of ties, t is the size of the ith tie, and Σ represents the summation of all ties. Additionally, the Mann-Kendall statistic Z is estimated using Equation (7): where Z is the standard normal variable; a positive value of Z shows an increasing trend, whereas a negative value indicates a decreasing trend. Additionally, the null hypothesis H 0 that Z with no significant trend is , where α is the significance level for the test and ±Z 1− α 2 are the standard normal deviates. The trend analysis was performed for the time series of water cycle components of all 109 subbasins for two separate periods, historical  and future periods (2071-2099).

Changes in Precipitation and Mean Temperature
Variation in temperature patterns caused by climate change could impact the temporal and spatial distribution of water resources [39]. The absolute changes in seasonal mean temperature of the future period compared to that of the historical periods for the five major river basins are shown in Figure 4.
In fact, this figure shows the increase in mean seasonal temperature for all subbasins (26 subbasins for the Han, 33 subbasins for the Nakdong, 21 subbasins for the Geum, 15 subbasins for the Sumjin, and 14 subbasins for the Youngsan River basins) with respect to the reference period. The results showed that all basins during all seasons are expected to exhibit increased mean temperature by the impact of climate change. In addition, the higher increase in mean temperature is projected for autumn and winter compared to that of spring and summer. Moreover, the Han and Geum River basins revealed higher increases in mean seasonal temperature than that of the other basins.

Linear Seasonal Trends of Water Cycle Components
Tables 2-6 reveal the statistical analysis of the linear trends of the five major basins for the historical and future periods. Linear trends were performed for each water cycle component during the annual time series of each season for the period of 1976 to 2005 for the historical period and the annual time series of 2071 to 2099 for the future period. This analysis was performed separately for all five major river basins of South Korea, namely the Han, Nakdong, Geum, Sumjin, and Youngsan River basins. The positive value of slope (B) defines an increasing trend, and a negative B represents a decreasing trend in a time series. In addition, the standard deviation (SD) of the seasonal water cycle components are shown in Tables 1-5, which revealed the variability in the water cycle components over the seasonal time series.
A comparison of the seasonal precipitation trends during the historical and future periods is shown in Table 2. The linear trend of mean seasonal precipitation during the historical period exhibits increasing trends in spring and decreasing trends in other seasons for all basins. However, different trends were projected for the future period in which spring and winter are projected to have decreasing trends and summer and autumn are expected to have increasing trends. Summer and autumn have the strongest decreasing trends for precipitation, and winter has the weakest trends for all basins. Moreover, SD showed higher variability for precipitation for summer in the historical period and winter in the future period, with maximum variability in the Sumjin River basin. In the present study, there are increasing trends for spring and decreasing trends for the other months, and a significant trend is only detected for summer in the Youngsan Basin. In a previous study (Bae et al., 2008), a decreasing trend in spring and an increasing trend in summer were reported for precipitation, and a significant trend was reported for spring in the Geum, Sumjing, and Youngsan river basins. The difference between the results is due to the different time periods, which are 1976 to 2005 in the present study and 1968 to 2001 for the previous study.
The standard deviation of precipitation for the future period shows extremely higher variability than that for the historical period. For example, SD for the winter in the Youngsan River basin exhibits SD = 15.79 mm for the historical period, while it is projected to be 236.84 mm for the future period, which shows almost 15 times higher variability than that in the historical period. This indicates that the degree of dispersion of precipitation in the winter during the considered time period and the higher seasonality of precipitation in the winter are due to the effects of climate change.

Linear Seasonal Trends of Water Cycle Components
Tables 2-6 reveal the statistical analysis of the linear trends of the five major basins for the historical and future periods. Linear trends were performed for each water cycle component during the annual time series of each season for the period of 1976 to 2005 for the historical period and the annual time series of 2071 to 2099 for the future period. This analysis was performed separately for all five major river basins of South Korea, namely the Han, Nakdong, Geum, Sumjin, and Youngsan River basins. The positive value of slope (B) defines an increasing trend, and a negative B represents a decreasing trend in a time series. In addition, the standard deviation (SD) of the seasonal water cycle components are shown in Tables 1-5, which revealed the variability in the water cycle components over the seasonal time series.
A comparison of the seasonal precipitation trends during the historical and future periods is shown in Table 2. The linear trend of mean seasonal precipitation during the historical period exhibits increasing trends in spring and decreasing trends in other seasons for all basins. However, different trends were projected for the future period in which spring and winter are projected to have decreasing trends and summer and autumn are expected to have increasing trends. Summer and autumn have the strongest decreasing trends for precipitation, and winter has the weakest trends for all basins. Moreover, SD showed higher variability for precipitation for summer in the historical period and winter in the future period, with maximum variability in the Sumjin River basin. In the present study, there are increasing trends for spring and decreasing trends for the other months, and a significant trend is only detected for summer in the Youngsan Basin. In a previous study (Bae et al., 2008), a decreasing trend in spring and an increasing trend in summer were reported for precipitation, and a significant trend was reported for spring in the Geum, Sumjing, and Youngsan river basins. The difference between the results is due to the different time periods, which are 1976 to 2005 in the present study and 1968 to 2001 for the previous study.
The standard deviation of precipitation for the future period shows extremely higher variability than that for the historical period. For example, SD for the winter in the Youngsan River basin exhibits SD = 15.79 mm for the historical period, while it is projected to be 236.84 mm for the future period, which shows almost 15 times higher variability than that in the historical period. This indicates that the degree of dispersion of precipitation in the winter during the considered time period and the higher seasonality of precipitation in the winter are due to the effects of climate change. Table 3 shows the linear trend in mean seasonal actual evapotranspiration during the historical and future periods. The slope of the trend showed an increasing trend during the historical period for all seasons; however, spring and winter are projected to have decreasing trends during the future period for all basins. In addition, higher variability in AET was estimated during the summer for the historic period, while the variability in AET increased significantly for the future period, with a maximum in summer for all basins. Table 4 shows the statistics regarding the linear trend in soil moisture for the historic and future periods. The slope of the time series of soil moisture during the historical period showed a decreasing trend for all seasons; however, during the future period, the Geum, Sumjin, and Youngsan River basins are projected to have an increasing trend in soil moisture. Table 5 shows the mean seasonal linear trends in groundwater for the historical and future periods. The slope of the trend showed a decreasing trend for summer, autumn, and winter during the historical period, while opposite trends were projected for the future period, in which spring showed a decreasing trend and the other seasons revealed increasing trends.
A comparison of changes in mean seasonal soil moisture between two periods showed a slight decrease for the future period compared to the reference period. In fact, a higher increase in temperature caused higher evapotranspiration and therefore a decrease in soil moisture. However, the groundwater is projected to slightly increase from the historical period to the future period as a result of the increase in precipitation. The increase in groundwater was mainly caused by a precipitation rate higher than the evapotranspiration rate, which caused recharge of the water deficit rate and increased soil moisture and groundwater.      Table 6 shows the results of the seasonal linear trend in streamflow for the five basins during the historical and future periods. The linear trend in streamflow showed an increasing trend in spring and a decreasing trend in other seasons. However, spring and winter are projected to have downward trends, while summer and autumn are projected to have upward trends during the future period. Similar to precipitation and actual evapotranspiration, the SD showed higher variability in summer for the historical period and in winter for the future period, with maximum variability in the Sumjin River basin, indicating that streamflow fluctuated more strongly in the Sumjin River basin during the historical and future periods.

Indication of the Significant Trend in Seasonal Water Cycle Components
The MK test was performed to detect the significant seasonal trends in water cycle components to understand how the trends in water cycle components changed during the historical period  and how they are projected to change during the future period (2071 to 2099) by the impact of climate change. The analysis was performed for all 109 subbasins of the Korean river basins. The results are shown in Figures 3-7.

Precipitation
The spatial distribution of the MK test for the seasonal precipitation trends in the historical and future periods are shown in Figure 5. The results of the trend test for the historical period revealed an increasing trend in precipitation during the spring with a significant increasing trend on the eastern coast. Other seasons showed a decreasing trend in precipitation with a significant decreasing trend for the two subbasins of the Han River during the autumn season. However, a different pattern is projected for the trend in precipitation for the future period. Opposite trends in precipitation are projected for spring and all basins affected by the overall decreasing trend in precipitation. During summer, both increasing and decreasing trends are projected, while autumn showed an increasing trend with a significant increasing trend in the eastern subbasins of Nakdong and the central basin of the Han River. Winter showed the same decreasing trend in precipitation with significant decreasing trends in the southern parts.  show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively.

Actual Evapotranspiration
The trends in actual evapotranspiration are shown in Figure 6 for the historical and future periods. The results of the MK test revealed an increasing trend in actual evapotranspiration during all seasons except for some southern subbasins in summer and autumn. In particular, spring and winter in most of the subbasins revealed significantly increasing trends. Conflicting behavior was The basins bordered by red show an increasing trend, while the basins outlined by blue show a decreasing trend. The white color indicates no significant trends. Light and dark blue areas show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively.

Actual Evapotranspiration
The trends in actual evapotranspiration are shown in Figure 6 for the historical and future periods. The results of the MK test revealed an increasing trend in actual evapotranspiration during all seasons except for some southern subbasins in summer and autumn. In particular, spring and winter in most of the subbasins revealed significantly increasing trends. Conflicting behavior was projected for spring and winter during the future period. Spring is projected to have increasing trends in the northern and southern subbasins and decreasing trends in the eastern and central subbasins. Winter is projected to have a decreasing trend in actual evapotranspiration in most of the basins, with significant decreasing trends in the southern part. An increasing trend in actual evapotranspiration indicates a reduction in the available water; however, the decreasing trend in actual evapotranspiration caused an increasing trend in the water availability.

Soil Moisture
The trend in soil moisture during the historical and future periods is shown in Figure 7. The decreasing trend in soil moisture, similar to that of evapotranspiration, revealed a reduction in the water availability in the system. The results of the trend test for the historical period revealed a decreasing trend for all seasons, with striking decreasing trends in summer for all subbasins except for the northern subbasin of the Han and Geum River basins, whose decreasing trends are not significant. Additionally, autumn revealed a significant decreasing trend for the southern subbasins (Nakdong, Sumjin, and Youngsan). The decreasing trend in soil moisture is projected for the future periods, which revealed that the amount of soil moisture still has a decreasing trend for these seasons; however, the amount of soil moisture in autumn is projected to increase for the future period. In winter, the central part of Korea is projected to have an increasing trend for soil moisture during 2071 to 2099; however, the western and eastern parts are expected to have reduced soil moisture.
The results of the study conducted by Jung et al. [40] showed that ET is affected by solar radiation or vapor pressure deficit (VPD) when there is enough soil moisture. However; when the soils are too dry, ET is mostly impacted by soil moisture and the lower amount of ET can cause to increase the atmospheric dryness. Moreover, reduction in soil water content leading to high evaporative demand and reduction of plant production, and possibly limiting the amount of water cycle [41].   and future period (2071-2099) were identified using the Mann-Kendall test. The analysis was performed for 109 subbasins. The basins bordered by red show an increasing trend, while the basins outlined by blue show a decreasing trend. The white color indicates no significant trends. Light and dark blue areas show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively.

Soil Moisture
The trend in soil moisture during the historical and future periods is shown in Figure 7. The decreasing trend in soil moisture, similar to that of evapotranspiration, revealed a reduction in the water availability in the system. The results of the trend test for the historical period revealed a decreasing trend for all seasons, with striking decreasing trends in summer for all subbasins except for the northern subbasin of the Han and Geum River basins, whose decreasing trends are not significant. Additionally, autumn revealed a significant decreasing trend for the southern subbasins (Nakdong, Sumjin, and Youngsan). The decreasing trend in soil moisture is projected for the future periods, which revealed that the amount of soil moisture still has a decreasing trend for these seasons; however, the amount of soil moisture in autumn is projected to increase for the future period. In winter, the central part of Korea is projected to have an increasing trend for soil moisture during 2071 to 2099; however, the western and eastern parts are expected to have reduced soil moisture.
The results of the study conducted by Jung et al. [40] showed that ET is affected by solar radiation or vapor pressure deficit (VPD) when there is enough soil moisture. However; when the soils are too dry, ET is mostly impacted by soil moisture and the lower amount of ET can cause to increase the atmospheric dryness. Moreover, reduction in soil water content leading to high evaporative demand and reduction of plant production, and possibly limiting the amount of water cycle [41]. show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively. Figure 8 shows the results of the MK test for the groundwater during the historical and future periods. The increasing trend in the groundwater during the historical period was only estimated in spring for the coastal areas, while other subbasins in spring and other seasons revealed decreasing trends. However, a significant decrease in soil moisture was detected for most of the subbasins in the eastern and northern parts and in winter for the southern subbasins. However, for the future period, opposite trends were projected for the groundwater, and the subbasins that showed a significant decreasing trend exhibited increasing trends in autumn and winter. The basins bordered by red show an increasing trend, while the basins outlined by blue show a decreasing trend. The white color indicates no significant trends. Light and dark blue areas show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively. Figure 8 shows the results of the MK test for the groundwater during the historical and future periods. The increasing trend in the groundwater during the historical period was only estimated in spring for the coastal areas, while other subbasins in spring and other seasons revealed decreasing trends. However, a significant decrease in soil moisture was detected for most of the subbasins in the eastern and northern parts and in winter for the southern subbasins. However, for the future period, opposite trends were projected for the groundwater, and the subbasins that showed a significant decreasing trend exhibited increasing trends in autumn and winter.

Spring
Summer Autumn Winter Historical period Future period show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively.

Streamflow
The historical and future trends in streamflow that revealed the combined impacts of all water cycle components are shown in Figure 9. The results of the MK test for the historical period exhibited increasing trends in 80 subbasins out of 109 subbasins for the spring season. However, other seasons revealed a decreasing trend in streamflow with a significant decrease in some catchments of Nakdong in summer, Han in autumn, and Geum, Sumjin, and Youngsan in winter.
However, the future period is expected to have different trends in streamflow. The increasing trend in spring for the historical period is expected to change to a decreasing trend for the future period, with a significant decreasing trend in some basins of the Nakdong and Geum River basins. Summer showed a decreasing trend in streamflow except for the northern subbasins. Autumn in the future period is expected to have increased streamflow except for the Nakdong River basin, which is projected to have a decrease in the amount of streamflow. Streamflow in all subbasins in Korea is expected to increase for the period 2071 to 2099, with a significant increase in the Nakdong River subbasins. The basins bordered by red show an increasing trend, while the basins outlined by blue show a decreasing trend. The white color indicates no significant trends. Light and dark blue areas show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively.

Streamflow
The historical and future trends in streamflow that revealed the combined impacts of all water cycle components are shown in Figure 9. The results of the MK test for the historical period exhibited increasing trends in 80 subbasins out of 109 subbasins for the spring season. However, other seasons revealed a decreasing trend in streamflow with a significant decrease in some catchments of Nakdong in summer, Han in autumn, and Geum, Sumjin, and Youngsan in winter.
However, the future period is expected to have different trends in streamflow. The increasing trend in spring for the historical period is expected to change to a decreasing trend for the future period, with a significant decreasing trend in some basins of the Nakdong and Geum River basins. Summer showed a decreasing trend in streamflow except for the northern subbasins. Autumn in the future period is expected to have increased streamflow except for the Nakdong River basin, which is projected to have a decrease in the amount of streamflow. Streamflow in all subbasins in Korea is expected to increase for the period 2071 to 2099, with a significant increase in the Nakdong River subbasins.

Spring
Summer Autumn Winter Historical period Future period show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively. Table 7 shows the relative change in MME projected by 13 GCMs for the future period relative to the historical period for water cycle components. The positive value of relative change revealed an increasing change rate, while the negative values showed a decreasing change rate compared to that of the historical period.

Relative Changes in Seasonal Water Cycle Components
The comparison of the mean seasonal change rate of precipitation of the future period with respect to the historical period showed the increasing change rate for all basins, with higher increases in spring and lower increases in winter. Additionally, the comparison of the mean seasonal change rate of the actual evapotranspiration of basins showed a higher change rate in winter for the Han, in autumn for the Nakdong, and in spring for the Geum, Sumjin, and Youngsan river basins. However, the lowest mean seasonal change rate of precipitation occurred in winter for the Nakdong, Sumjin, and Youngsan river basins. The comparison of the mean seasonal precipitation and actual evapotranspiration exhibited a higher change rate of precipitation than actual evapotranspiration during spring and summer seasons. However, the amount of actual evapotranspiration is higher in winter for all basins. In addition, the mean seasonal change rate of soil moisture during the historical and future period revealed that the amounts of soil moisture in all basins are expected to decrease in relation to the historic period, with a greater decrease in the autumn and winter except for the Han River basin, where the higher decrease in soil moisture is projected to be in summer and autumn.
The mean seasonal change rate of groundwater showed a decreasing change rate for spring and an increasing change rate for summer and autumn in all basins. Moreover, winter showed a decreasing change rate for groundwater in all basins except for the Han River basin. Additionally, Table 7 shows the change rate of precipitation relative to the reference period. The results of the mean seasonal change rate of the Han River basin exhibited a different tendency compared to that of the other basins. The Han River basin is expected to have the minimum change rate in spring and the maximum change rate in winter, while in other basins, the maximum change rate is projected to be in summer and autumn, with the lowest change rate in winter. The basins bordered by red showed an increasing trend, while the basins outlined by blue show a decreasing trend. The white color indicates no significant trends. Light and dark blue areas show decreasing trends at 90% and 95% confidence levels, respectively. The light and dark red areas show increasing trends at 90% and 95% confidence levels, respectively. Table 7 shows the relative change in MME projected by 13 GCMs for the future period relative to the historical period for water cycle components. The positive value of relative change revealed an increasing change rate, while the negative values showed a decreasing change rate compared to that of the historical period.

Relative Changes in Seasonal Water Cycle Components
The comparison of the mean seasonal change rate of precipitation of the future period with respect to the historical period showed the increasing change rate for all basins, with higher increases in spring and lower increases in winter. Additionally, the comparison of the mean seasonal change rate of the actual evapotranspiration of basins showed a higher change rate in winter for the Han, in autumn for the Nakdong, and in spring for the Geum, Sumjin, and Youngsan river basins. However, the lowest mean seasonal change rate of precipitation occurred in winter for the Nakdong, Sumjin, and Youngsan river basins. The comparison of the mean seasonal precipitation and actual evapotranspiration exhibited a higher change rate of precipitation than actual evapotranspiration during spring and summer seasons. However, the amount of actual evapotranspiration is higher in winter for all basins. In addition, the mean seasonal change rate of soil moisture during the historical and future period revealed that the amounts of soil moisture in all basins are expected to decrease in relation to the historic period, with a greater decrease in the autumn and winter except for the Han River basin, where the higher decrease in soil moisture is projected to be in summer and autumn.
The mean seasonal change rate of groundwater showed a decreasing change rate for spring and an increasing change rate for summer and autumn in all basins. Moreover, winter showed a decreasing change rate for groundwater in all basins except for the Han River basin. Additionally, Table 7 shows the change rate of precipitation relative to the reference period. The results of the mean seasonal change rate of the Han River basin exhibited a different tendency compared to that of the other basins. The Han River basin is expected to have the minimum change rate in spring and the maximum change rate in winter, while in other basins, the maximum change rate is projected to be in summer and autumn, with the lowest change rate in winter. In general, the reason for the decreasing change rate for streamflow in the Han River basin during the spring occurs due to the almost similar amount of increasing change rate in precipitation and actual evapotranspiration, which cause a decrease in the amount of soil moisture and therefore streamflow. However, the change rate of streamflow in the Han River basin is projected to increase by 23.64%. This increasing change rate occurred due to the lower decrease in soil moisture and higher change rate of groundwater in this season. In contrast, the decreasing change rate of streamflow for the Youngsan River basin in winter occurred due to the elevated amount of actual evapotranspiration compared to precipitation, which caused reduced soil moisture and therefore decreased groundwater recharge. Additionally, the increase in the change rate of streamflow in summer for the Nakdong River basin (+15.90%) was mainly caused by a precipitation (+12.96%) change rate higher than the evapotranspiration change rate (+10.46%), which caused recharge of the water deficit rate and increased soil moisture and groundwater. Soil moisture has a strong influence on evapotranspiration and the amount of water that percolates into the underlying aquifer or contributes directly to surface streamflow, which causes an increase in the amount of streamflow.

Conclusions
In this study, the trends of water cycle components in the time series of precipitation, actual evapotranspiration, soil moisture, groundwater, and streamflow were analyzed using the multimodel ensemble of 13 GCMs of RCP 4.5 for two separate periods: the historical period from 1976 to 2005 and the future period from 2071 to 2099. The parameters of the hydrologic model were calibrated and verified to produce the hydrological variables. To assess the trends in the water cycle components, two statistical tests were applied: the linear trend and the Mann-Kendall test. Linear trend analysis was performed for five major river basins, and the Mann-Kendall test was performed for all 109 subbasins of the Korean river basins.
The results of the linear trend test during the historical period revealed both increasing and decreasing trends for water cycle components during seasons with similar tendencies among basins. During the historical period, the spring season showed increasing trends for streamflow, while other seasons showed decreasing trends mainly affected by similar trends in precipitation and groundwater.
Overall, the increasing or decreasing trend in streamflow occurs due to the combined effects of all water cycle components. For example, during the future period, a decreasing trend in streamflow during the spring season occurred due to a decreasing trend in precipitation and an increasing trend in actual evapotranspiration in the northern and southern parts, and a decreasing trend in soil moisture and groundwater caused a decreasing trend in streamflow. However, although the trend of precipitation is decreasing in winter for the future period, the decreasing trend in actual evapotranspiration, increasing trend for soil moisture (especially in the central subbasins), and increasing trend for groundwater caused an increasing trend for streamflow in this period. In general, the results of the trend analysis for the streamflow revealed opposite trends for streamflow in spring (increasing for the historical period and increasing for the future period), autumn (significant decreasing in the northern and southern subbasins for the historical period while increasing trend in the northern and southern part for the future period), and winter (decreasing for the historical and increasing for the future period).