Trend Analysis of Long-Term Reference Evapotranspiration and Its Components over the Korean Peninsula

: In this study, the spatial and temporal trends of reference evapotranspiration (ET o ) and its components consisting of the energy term (EN o ) and the aerodynamic term (AE o ) were considered over the Korean Peninsula. The T-test and Mann–Kendall (MK) test were used to detect parameter trends after removing the effect of serial correlation from annual and seasonal time series between 1980 and 2017. Due to the lack of solar-radiation data for North Korea (NK), a regionally calibrated model based on South Korea (SK) weather data was developed to estimate daily solar radiation in NK. The results showed that spatial distribution of the ET o increased southward in the range from 705 mm/year in the northeast to 1195 mm/year in the southeast of the Korean Peninsula. The spatial patterns of the EN o and AE o varied from the minimum in the north and increased southward, reaching their maximum values in the southern parts of the Korean Peninsula. The mean annual ET o values of SK and NK were also compared. Over the 37-year period, mean annual evapotranspiration in SK was approximately 18% higher than that in NK. Moreover, mean areal EN o and AE o in SK were higher than in NK. The trend of the EN o on annual and seasonal scales was also upward. In contrast, the trend of the AE o decreased over the Korean Peninsula through all seasons and annual scales. These opposite trends in the EN o and AE o parameters mitigated the signiﬁcant trends of the ET o . Finally, the stronger signiﬁcant upward trend of the energy term led to signiﬁcant increasing trends of ET o on the Korean Peninsula, with EN o being the dominant component in the increase of the ET o .


Introduction
Aside from precipitation, evapotranspiration (ET) is the most important component of the hydrological cycle. ET determines crop-water requirements and can play an important role in groundwater recharge, water distribution, seasonal water management, long-term estimates of water supply, and the design and management of hydraulic structures. A common procedure to estimate ET is to calculate reference evapotranspiration (ET o ) and multiply it by an appropriate crop coefficient [1]. Careful attention to ET o is necessary for the determination of water budget, efficient irrigation scheduling, and water-resource management [2]. Since the variability of ET o directly affects the availability of water resources, detecting the trends of ET o is necessary [3]. Additionally, this parameter can reflect the effects of climate change [4].
Many empirical methods have been suggested and utilized by researchers to estimate ET [5]. The Food and Agriculture Organization (FAO) of the United Nations recommends the FAO Penman-Monteith (FAO PM) method as a unique standard method for estimating ET o [6], which consists of energy and aerodynamic terms.
In recent decades, several studies have detected global or regional ET trends. On the global scale, both increasing and decreasing trends worldwide have been reported for ET [7][8][9]. Previous studies have demonstrated that the ET o trend varies regionally. Some studies reported increasing trends in ET o (e.g., Wang et al. (2014) in China [10], Tabari et al. (2012) in Iran [11], and Silva (2004) in Brazil [12]). Downward trends in ET have also been found in several countries worldwide, including China [13], Japan [14], New Zealand [15], Thailand [16], and the United States [17]. Moreover, several researchers have reported both positive and negative ET o trends in areas with different climate types. Gao et al. [18] found a decreasing trend in most areas of east China, while there was a positive trend in the west and the northern parts of northeast China. Moreover, Li et al. [19] indicated mixed trends of significantly increasing ET o in the lower Huai River Basin (HRB) and decreasing ET o in the upper-middle HRB and Yi-Shu-Si River Basin of China. In another study, Mosaedi et al. [20] reported mixed trends of significantly increasing and decreasing ET o in different parts of Iran (i.e., in mountainous, semidesert, coastal wetland, and semimountainous areas).
The above studies provide useful information on ET trends; however, these studies have obtained more or less similar results in different parts of the world. Interestingly, as of now, there is no study investigating which ET o term (energy term or aerodynamic term) has the main impact on the change in ET o trends and how these terms change over the seasons in an area characterized by heterogeneous weather and geographic conditions. The Korean Peninsula, which is exposed to subtropical monsoons, is covered by nearly 70% mountainous areas and, being a peninsula, is surrounded by water on three sides, resulting in complex spatial and temporal weather conditions. In this study, we focus on the analysis of ET o trends, as well the as aerodynamic and energy terms of ET o , to understand how it varies annually and seasonally over the Korean Peninsula. To this end, the spatial and temporal variation of ET o and the energy and aerodynamic terms of the FAO PM method were explored using the T-test and MK test over seasonal and annual time scales.

Estimation of ET o
The FAO PM method is recommended as a standard method for estimating ET o . The reference surface is defined as a hypothetical crop that has evaporation similar to an extension surface of green grass with a height of 0.12 m, an albedo of 0.23, and a constant surface resistance of 70 s m −1 , and that is actively growing and well-watered. The FAO PM method can be expressed as follows [6]: where ET o is the reference evapotranspiration (mm day −1 ), ∆ is the slope vapor pressure curve (KPa • C −1 ), R n is the net radiation (MJ m −2 day −1 ), G is the soil heat flux density (MJ m −2 day −1 ), γ is the psychrometric constant (KPa • C −1 ), T is the average daily air temperature ( • C), u 2 is the wind speed at 2 m above ground level (m s −1 ), e s is the saturation vapor pressure (KPa), and e a is the actual vapor pressure (KPa). Other relevant equations used in computing ET o via the FAO PM method are explained in FAO 56 [6] and summarized in Table S1. Because the FAO PM method is a combination of energy and aerodynamic components, it can be written as follows: short and well-watered reference vegetation is exposed to bright sunshine and light wind speed. The aerodynamic term reveals the ET o when short and well-watered reference vegetation is exposed to high wind speed and a large vapor-pressure deficit. In this study, the daily ET o was calculated for each station. Then, the mean seasonal and annual ET o values for each year were derived by adding the daily values during each season and the annual ET o was obtained by summing the seasonal values.

Estimation of Solar Radiation
Solar radiation, or total sunshine hours, is one of the inputs used in the FAO PM method to estimate ET o . However, when solar-radiation or sunshine-duration data are unavailable, FAO 56 recommends using Hargreaves' radiation formula to estimate solar radiation. Hargreaves' radiation formula is calibrated and validated for various stations in different climate conditions. This formula, proposed by Reference [21], is expressed as follows: where R s is the solar radiation (MJ m −2 day −1 ); TD is the diurnal temperature, which is the difference between the maximum temperature (T max ) and the minimum temperature (T min ) ( • C); K Rs is the adjustment coefficient; and R a is the extraterrestrial radiation (MJ m −2 day −1 ), which varies throughout the day of the year and latitude. Hargreaves (1994) [22] suggested values of K Rs = 0.16 for inland areas and K Rs = 0.19 for coastal regions. However, other studies have recommended calibrating K Rs regionally to obtain an accurate estimation of solar radiation using long-term mean temperature (T mean ) and diurnal temperature (TD) [23][24][25]. In the regional calibration of K Rs , K Rs values are calculated using the ratio of the daily R s to R a and the square root of TD for each meteorological station. Then, the different relations for K Rs are determined using multiple regression with the long-term mean T mean and TD ( T mean TD ), station elevation (EL), and latitude (ϕ) as the independent variables. After statistical analysis, the selected relation for K Rs was applied to estimate R s using Equation (3).

Trend-Analysis Method
Parametric T-tests [26] and nonparametric Mann-Kendall trend tests [27][28][29] were used to determine the trends in the seasonal and annual ET o and components of the FAO PM method (i.e., EN o and AE o ).

Parametric T-Test
The T-test is a parametric test considering the assumption that data are normally distributed and uses a statistical technique for detecting trends in a time series by fitting a linear equation to observations. The estimated equation uses a linear relationship between X as an independent variable and Y as a dependent variable, and is written as follows: where A and B denote the intercept and the slope of the regression line, respectively. Positive slope A represents an upward trend and a negative slope indicates a downward trend in a time series.
The T-tests apply the hypothesis test on the standard slope coefficient, and the result defines as T o and is calculated with the following equation: where B 1 is the standard slope of the regression line and SE is the standard error of the slope. SE is calculated as follows: where, x i is an independent variable and y i is an independent variable in the time series.
x is the average of x i and y is the average of y i in the time series, and n is the number of observed data. T o has a t distribution with n-2 degrees of freedom, where n is the total number of observations. Null hypothesis has no significant trend for where t α 2 ,n−2 and −t α 2 ,n−2 represent the standard normal deviates and α determines the level of significance. In this test, 5% and 10% significance levels are considered.

Nonparametric Mann-Kendall test
The Mann-Kendall (MK) test is a nonparametric test used to analyze the null hypothesis of no change in the mean of a series (no trend) against the alternative hypothesis of an increase or decrease over time (existence of a trend) [30]. This method is frequently used to determine trends in hydrological variables. Kendall's statistic S is obtained using the following equations: where n is the number of observations, 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 . Statistic S has normal distribution with zero mean and variance and is computed as (10) where m is the number of ties, t is the size of the ith tie, and Σ is the summation of all ties. MK statistic Z is estimated as follows: where Z is the standard normal variable, with positive values of Z indicating increasing trends, whereas negative values indicate decreasing trends. The null hypothesis H o that Z has no significant trend is , where α is the significance level for the test (5% and 10% significance levels) and +Z 1− α 2 and −Z 1− α 2 are the critical values for the two-sided hypothesis. The slope of the trends in the nonparametric MK test can be calculated as follows: where X j and X i are sequential data (j > i) and β is Theil-Sen's estimator, which is a robust estimate of the trend magnitude in the MK test. A positive Theil-Sen slope indicates an increasing trend and a negative slope represents a decreasing trend in a time series.

Impact of Serial Correlation in a Time Series
The parametric T-test and MK test are designed for serially independent series. Therefore, the presence of incorrect rejection of the null hypothesis when there is no trend (Type I error) is correct only when the data are independent. As serial correlation increases, the variance of the test statistic and the rejection rate of the null hypothesis increase, specifically when there is a positive serial correlation [30,31]. In this study, the prewhitening approaches proposed by von Storch and Navarra [32] were used to eliminate the effect of serial correlation by removing the lag-one autoregressive component (r 1 ) of the series with serial correlation before applying trend tests. To this end, r 1 was computed and compared to the 95% confidence limits. If r 1 ranged between −1−1.96 , r 1 was not significant; otherwise, r 1 is significant (n is the number of data). The trend tests are applied to the original time series (x 1 , x 2 , . . . , x n ) if r 1 is not significant. However, if r 1 is significant, then the time series should be prewhitened as (x 2 − r 1 x 1 , x 3 − r 1 x 2 , . . . , x n − r 1 x n−1 ) to remove the effect of serial correlation prior to application of the trend tests [31,33].

Study Area and Data
The Korean Peninsula is in northeastern Asia and is divided into two parts: South Korea (SK) and North Korea (NK). Figure 1 shows the spatial distribution of stations (name and code) of topographic and land-cover characteristics of the study area. The Korean Peninsula is connected to the northeastern part of mainland China and is surrounded by oceans along the eastern, western, and southern coasts. The peninsula is mountainous in the northern and eastern parts and has well-developed plains in the western and southern parts. The climate of the peninsula is affected by the East Asian monsoon in summer and typhoons in autumn. Nearly two-thirds of the annual rainfall occurs during the flooding season in the summer months [34]. Winters are cold, with typical January temperatures below freezing and minimal precipitation and little snow accumulation outside of the mountainous areas. However, climate differs dramatically from North to South Korea. The northern regions have a colder climate in the more mountainous area (in the north and east), while the southern regions have a warm and wet climate in the flat plains of the southern and western peninsula. Monthly average rainfall is approximately 112.5 mm in SK and 74.5 mm in NK. Mean annual temperature is 13.1 • C in SK, from 1.2 • C in winter to 24.1 • C in summer, and 9.0 • C in NK, from -4.7 • C in winter to 21.5 • C in summer.
The detection of ET o trends requires the use of long-term meteorological data, including the daily observed R s ,T max , T min , dew-point temperature (T dew ), and wind speed (W s ). Daily meteorological data from 21 meteorological stations in SK and 27 stations in NK were collected from the Korea Meteorological Administration (https://data.kma.go.kr) within a sufficient meteorological time period . Statistics for the observed daily R s are available for the 21 weather stations in SK but not for those in NK. The spatial distribution of the stations in SK and NK is shown in Figure 1.
The data quality and integrity assessments were performed using the procedures explained by Reference [6] prior to the analysis. The missing data were filled with averaged values from the neighboring sites using an inverse-distance weighting method [35]. The numbers of missing data for each parameter of stations in SK and NK are provided in Table S2. For the quality check for temperature, the observed T max and T min data were compared against the long-term temperature extremes for each year [24]. To evaluate the integrity of solar-radiation data, the observed solar-radiation data of each station were compared with the clear-sky solar-radiation envelope (R so ) for each year. R so is a function of a day of the year, latitude, and elevation. The calculated R so can be plotted as an upper envelope of observed Rs. This method is especially useful for areas with high air pollution or airborne dust without clear skies. In these regions, there were differences between the measured R s and calculated R so , and R s was significantly below or above the computed R so . This issue was further investigated by using multiplicative factor ((R s ) cor = a.R s ), additive factor ((R s ) cor = a + R s ), or a combination of these variables to force the measured R S to reach the computed R so (where a is the corrective factor). Additional explanation about the integrity of solar-radiation data using clear-sky comparisons can be found in FAO 56 [6,24,36]. There were differences between R s and R so for some years at stations in SK. These mismatches were found and modified for all stations in SK. Figure 2 shows the time series of solar radiation before and after this prescreening for the Incheon station in 1980 and Seaosan station in 1982.
and North Korea (NK). Figure 1 shows the spatial distribution of stations (name and code) of topographic and land-cover characteristics of the study area. The Korean Peninsula is connected to the northeastern part of mainland China and is surrounded by oceans along the eastern, western, and southern coasts. The peninsula is mountainous in the northern and eastern parts and has welldeveloped plains in the western and southern parts. The climate of the peninsula is affected by the East Asian monsoon in summer and typhoons in autumn. Nearly two-thirds of the annual rainfall occurs during the flooding season in the summer months [34]. Winters are cold, with typical January temperatures below freezing and minimal precipitation and little snow accumulation outside of the mountainous areas. However, climate differs dramatically from North to South Korea. The northern regions have a colder climate in the more mountainous area (in the north and east), while the southern regions have a warm and wet climate in the flat plains of the southern and western peninsula. Monthly average rainfall is approximately 112.5 mm in SK and 74.5 mm in NK. Mean annual temperature is 13.1 °C in SK, from 1.2 °C in winter to 24.1 °C in summer, and 9.0 °C in NK, from -4.7 °C in winter to 21.5 °C in summer.

Regional Calibration of K Rs and Estimation of Solar Radiation
Due to the lack of solar-radiation data in NK, R s was estimated based on Hargreaves' formula and the regionally calibrated formula for K Rs . The different relationships were determined by multiregression of long-term mean values of T mean TD , EL, and ]ϕ as the independent variables and station values of K Rs as the dependent variable. The following equation was selected as the best model after statistical analysis: Figure 3 depicts the actual K Rs against the computed K Rs using Equation (13) in SK. The estimated K Rs values are well-correlated with the calculated values, with a high correlation coefficient (R 2 = 0.89) and small root-mean-square error (RMSE = 0.007). Using the K Rs values for all stations in SK and NK, R s was estimated using Hargreaves' formula (Equation (3)) for all stations in SK and NK. Table 1 shows the statistical analysis for the comparison of observed and calculated daily solar-radiation data for 21 stations in SK. The results show good agreement between observed and estimated solar-radiation data in SK. However, because K Rs was determined using regional calibration data in SK, there is uncertainty in the estimated R s in NK. estimated solar-radiation data in SK. However, because was determined using regional calibration data in SK, there is uncertainty in the estimated in NK.

Estimation of ET o Using the Penman-Monteith Method
The time series of the annual ET o for 27 stations in NK and 21 stations in SK are illustrated as a box-and-whisker plot in Figure 4a and b, respectively. In NK, the maximum and minimum estimated ET o values were at Huicheon station (1005.5 mm) and Seonbong station (660.3 mm), respectively. In SK, the maximum estimated ET o was 1192.8 mm for Daegu station, whereas the minimum was 826 mm for Daegwalryeong station, indicating that the mean annual ET o in SK is higher than that in NK. The time series of the annual ET for 27 stations in NK and 21 stations in SK are illustrated as a box-and-whisker plot in Figure 4a and b, respectively. In NK, the maximum and minimum estimated ET values were at Huicheon station (1005.5 mm) and Seonbong station (660.3 mm), respectively. In SK, the maximum estimated ET was 1192.8 mm for Daegu station, whereas the minimum was 826 mm for Daegwalryeong station, indicating that the mean annual ET in SK is higher than that in NK. The Kriging method was used to interpolate the mean annual and seasonal ET , EN , and AE over the Korean Peninsula by including the data of all stations in SK and NK. The Kriging method is widely used for spatial analysis by using geographic information systems (GIS). This method applies the distance-weighting approach to interpolating the data by including information from specific

Spatial Distribution of Mean Annual and Seasonal ET o , EN o , and AE o
The Kriging method was used to interpolate the mean annual and seasonal ET o , EN o , and AE o over the Korean Peninsula by including the data of all stations in SK and NK. The Kriging method is widely used for spatial analysis by using geographic information systems (GIS). This method applies the distance-weighting approach to interpolating the data by including information from specific locations [37]. The spatial distribution pattern of the mean annual ET o averaged from 1980 to 2017 is shown in Figure 5, which reflects the combined effect of meteorological variables and topographic characteristics across the Korean Peninsula. The annual ET o variation exhibits a clear spatial pattern (southward increase), ranging from 705 mm/year to 1195 mm/year over the entire basin. A lower ET o is present in the northeast, while higher values are found in the southeast of the Korean Peninsula.
due to the differences in the position of the sun [6], and it is higher near the equatorial region than the polar region. Thus, higher were observed in the southern part (higher latitude) of the Korean Peninsula than the northern part (lower latitude). Additionally, T in North Korea mainly affects the propagation of Siberian cold weather, and the movement of this cold stream from Siberia to North Korea produces temperature drops in the region. However, T in South Korea is higher, which is mostly affected by the subtropical warm current. Therefore, the higher magnitude of and T in South Korea causes a higher value of EN in SK than NK. AE primarily affects the relative humidity and wind speed. The intense value of relative humidity and wind speed in the southern part causes higher AE in SK relative to NK. In general, the higher value of EN and AE in SK than NK results in higher ET in SK.
Moreover, the seasonal spatial distribution of ET shows that the highest amount of ET occurs in the summer and then in spring, while winter has the lowest value of ET . This pattern exists due to the high amount of EN in the summer and the seasonal spatial distribution of EN . Furthermore, the comparison of the spatial distribution of ET with EN and AE in different seasons shows that ET in spring and summer is mostly affected by EN , whereas in autumn and winter it is mostly caused by AE . This finding indicates that ET in spring and summer has a stronger effect resulting from the superimposed effect of and T in autumn and winter by the combined influence of RH and WS.

Differences in the Mean Annual and Seasonal , , and between SK and NK
The time series of seasonal and annual ET , EN , and AE from 1980 to 2017 averaged over the entire catchment in SK and NK are graphically presented in Figure 6, and the comparison of the results is reported in Table 2. The mean annual ET in SK is 971.2 mm/year, whereas that for NK is 823.5 mm/year, which indicates that the mean annual ET in SK is 18% higher than that in NK over the 37-year period. In addition, the annual EN and AE terms are 708.4 and 262.7 mm/year in SK and 648.0 and 175.5 mm/year in NK, respectively. It follows that mean areal EN and AE are higher, at 9.3% and 49.7%, respectively, in SK than in NK.

The spatial distribution of mean annual EN o and AE o show that there is a distinct pattern of minimum values in the northern area and maximum values in the southern area of the Korean
Peninsula with different directions of variation. The spatial distribution of EN o decreases from the southwest to the northeast direction of the Korean Peninsula, while regarding the AE o , the maximum was established in the southeastward direction and the minimum in the northward. This pattern occurs because the energy term mainly affects the R s and T. R s is different based on the latitude due to the differences in the position of the sun [6], and it is higher near the equatorial region than the polar region. Thus, higher R s were observed in the southern part (higher latitude) of the Korean Peninsula than the northern part (lower latitude). Additionally, T in North Korea mainly affects the propagation of Siberian cold weather, and the movement of this cold stream from Siberia to North Korea produces temperature drops in the region. However, T in South Korea is higher, which is mostly affected by the subtropical warm current. Therefore, the higher magnitude of R s and T in South Korea causes a higher value of EN o in SK than NK. AE o primarily affects the relative humidity and wind speed. The intense value of relative humidity and wind speed in the southern part causes higher AE o in SK relative to NK. In general, the higher value of EN o and AE o in SK than NK results in higher ET o in SK.
Moreover, the seasonal spatial distribution of ET o shows that the highest amount of ET o occurs in the summer and then in spring, while winter has the lowest value of ET o . This pattern exists due to the high amount of EN o in the summer and the seasonal spatial distribution of EN o . Furthermore, the comparison of the spatial distribution of ET o with EN o and AE o in different seasons shows that ET o in spring and summer is mostly affected by EN o , whereas in autumn and winter it is mostly caused by AE o . This finding indicates that ET o in spring and summer has a stronger effect resulting from the superimposed effect of R s and T in autumn and winter by the combined influence of RH and WS.

Differences in the Mean Annual and Seasonal ET o , EN o , and AE o between SK and NK
The time series of seasonal and annual ET o , EN o , and AE o from 1980 to 2017 averaged over the entire catchment in SK and NK are graphically presented in Figure 6, and the comparison of the results is reported in   The most seasonal ET occurs in summer, which constitutes 44% and 47.5% of the total ET in SK and NK, respectively ( Table 2). The lowest ET occurs in winter in both SK and NK (8.3% in SK and 6.3% in NK). However, the annual ET exhibits similar ratios in spring and autumn.
To understand which of the components of ET (i.e., EN and AE ) result in higher ET in SK than NK, the percentages of mean EN and AE were compared to the mean ET for each season (see Table 2). The results reveal that spring and autumn have similar percentages of EN and AE in SK and NK, while the energy term (EN ) constitutes almost 70% and the aerodynamic term covers (AE ) 30% of the total ET on the Korean Peninsula. The EN in summer has the greatest impact  The most seasonal ET o occurs in summer, which constitutes 44% and 47.5% of the total ET o in SK and NK, respectively ( Table 2). The lowest ET o occurs in winter in both SK and NK (8.3% in SK and 6.3% in NK). However, the annual ET o exhibits similar ratios in spring and autumn.
To understand which of the components of ET o (i.e., EN o and AE o ) result in higher ET o in SK than NK, the percentages of mean EN o and AE o were compared to the mean ET o for each season (see Table 2). The results reveal that spring and autumn have similar percentages of EN o and AE o in SK and NK, while the energy term (EN o ) constitutes almost 70% and the aerodynamic term covers (AE o ) 30% of the total ET o on the Korean Peninsula. The EN o in summer has the greatest impact on the ET o when compared to the other seasons and contains approximately 85% of the total ET o in SK and NK. Furthermore, in winter, EN o and AE o have almost the same portion of ET o (approximately 50%). These results indicate that the EN o term plays an important role in ET o during the spring, summer, and autumn, while in the winter, both terms contribute equally to the amount of ET o .

Spatial Variation of Annual and Seasonal ET o , EN o , and AE o Trends
The results of prewhitening for ET o , EN o , and AE o on an annual scale are reported in Tables S3-S5 as supplementary documents. The spatial distributions of stations in SK and NK from 1980 to 2017 with significant increasing or decreasing trends according to the MK test are presented in Figure 7. As shown, 43% of stations in SK reveal significant increasing trends, and 14% show significant decreasing trends in annual ET o . The significant decreasing trend in SK is located in the southeastern area of the Korean Peninsula. However, 22% of stations in NK detected significant increasing trends, mainly located in the coastal areas, while 4% of the stations revealed significant decreasing trends. The annual spatial distribution of EN o showed that 66% of stations in SK, as well as NK, detected significant increasing trends, while no station showed significant decreasing trends for annual energy trends. In contrast, for AE o , significant decreasing trends were detected at 19% of stations in SK and 33% stations in NK; however, there were no significant increasing trends for this variable. These results indicated that significant increasing trends on the Korean Peninsula are caused by the significant increasing trends in the energy term. As described above, the EN o is primarily related to solar radiation and temperature. It follows that the trends of R s and T were increased on the Korean Peninsula. Similar results were also found using the T-test, and the results are shown in Figure 8.
The results showed that ET o increased more significantly on the plains than in the mountainous areas. EN o increased in both areas, while AE o decreased more significantly in the plains than in the mountainous areas.
In addition, the spatial distributions of the seasonal trends in ET o and its components were investigated. The spatial trends of ET o in spring revealed significant increasing trends for 29% and 44% of stations in SK and NK, respectively. Additionally, in summer, 33% of stations in SK and 59% of stations in NK detected significant increasing trends. This finding indicates that there is a more intense significant trend during spring and summer in NK. Moreover, the spatial distribution of ET o during autumn and winter revealed both significant upward and downward trends on the Korean Peninsula. Regarding the EN o , only significant increasing trends were detected for stations, and no station detected a significant downward trend. In contrast, almost all significant detected trends for AE o are decreasing for all seasons over the Korean Peninsula. These opposite trends in EN o and AE o caused the mitigation of the significant trends of ET o . Furthermore, the stronger significant upward trends of energy terms led to significant increasing trends of ET o on the Korean Peninsula. results indicated that significant increasing trends on the Korean Peninsula are caused by the significant increasing trends in the energy term. As described above, the EN is primarily related to solar radiation and temperature. It follows that the trends of and T were increased on the Korean Peninsula. Similar results were also found using the T-test, and the results are shown in Figure 8.
The results showed that ET increased more significantly on the plains than in the mountainous areas. EN increased in both areas, while AE decreased more significantly in the plains than in the mountainous areas.
ET  In addition, the spatial distributions of the seasonal trends in ET and its components were investigated. The spatial trends of ET in spring revealed significant increasing trends for 29% and 44% of stations in SK and NK, respectively. Additionally, in summer, 33% of stations in SK and 59% of stations in NK detected significant increasing trends. This finding indicates that there is a more intense significant trend during spring and summer in NK. Moreover, the spatial distribution of ET during autumn and winter revealed both significant upward and downward trends on the Korean Peninsula. Regarding the EN , only significant increasing trends were detected for stations, and no station detected a significant downward trend. In contrast, almost all significant detected trends for

Difference in Trends of Mean Annual and Seasonal ET o , EN o , and AE o in SK and NK
The time series of ET o , EN o , and AE o from 1980 to 2017 were averaged over the entire SK and NK using the kriging method. Then, the parametric T-test and nonparametric MK test were performed to detect the trends of ET o and its components in the seasonal and annual scales, and the results are reported in Table 3. The positive values of T and Z indicate increasing trends, whereas negative values indicate decreasing trends. B and β show the standard slope of the linear regression and Theil-Sen slope of the MK test, respectively, which reflect the rate of increase or decrease on seasonal and annual timescales. The results of the trends test for ET o showed that spring and summer trends in SK and NK exhibited significant increasing patterns, with stronger trends in summer. Furthermore, NK had more significant increasing trends for both seasons than those in SK. However, autumn showed increasing trends and winter decreasing trends, which were not significant on the Korean Peninsula. In addition, the annual trends of ET o indicated that although there are increasing trends in both SK and NK, significant trends were detected only in SK. A study conducted by Kim et al. [38] focused on the precipitation trends over the Korean Peninsula, and their results indicated that there was an increasing trend in precipitation in SK and a decreasing trend in NK during the summer months. The results of the present study showed that a significant increasing trend for ET o was observed in summers in NK. Lower precipitation and higher ET o in the summer months in NK could cause a reduction in soil moisture and vegetation cover and induce severe drought in NK. The results of the trend tests in EN o SK indicate that all trends are significantly positive in SK and NK on the seasonal and annual scales, and the trends are stronger in SK except for autumn. In contrast, AE o exhibited decreasing trends during all seasons and the annual time series in SK and NK. However, significant decreasing trends were shown in spring, summer, autumn, and the annual time scale in SK, while the other trends were not significant. The results of the trends test of EN o and AE o indicate that there are opposite trends for the energy term and aerodynamic term on the Korean Peninsula. Although the energy term exhibited significant increasing trends, the aerodynamic terms showed decreasing trends on seasonal and annual time scales. These trends are primarily affected by meteorological effects on the energy and aerodynamic terms.
Additionally, the comparison between the results of the two tests indicates that similar results were obtained by the parametric T-test and nonparametric MK test on the Korean Peninsula on both seasonal and annual scales, which is related to the degree of normality of the distribution [39,40].

Conclusions
This study investigated ET o , EN o , and AE o trends in seasonal and annual time series for the Korean Peninsula from 1980 to 2017 by considering the lack of solar-radiation data in NK. First, a regionally calibrated model was constructed based on SK weather data to estimate the daily solar radiation in NK. Meteorological data were collected from 21 stations in SK and 27 stations in NK. Then, ET o was estimated using the FAO PM method. Next, the effects of serial correlation in the time series with significant r 1 were removed using prewhitening, and trend analysis was performed to understand the temporal variations in ET o over the 37-year period in SK and NK using parametric linear regression and nonparametric MK tests on both seasonal and annual time series. The following results may be drawn from the present study:

•
The spatial distribution of the mean annual ET o exhibited increasing spatial variation in annual ET o from NK to SK. A lower ET o was evident for the northeastern Korean Peninsula, and higher values were found over the southeastern Korean Peninsula.

•
The spatial distribution of the EN o and AE o revealed a clear distribution from the minimum in NK and reached the peak in SK for both energy and aerodynamic terms. The energy term is influenced by R s and T, and the aerodynamic term is affected by RH and WS. The lower latitude in NK causes a lower amount of R s than that in SK, and mean annual T in NK is lower than in SK because T in the northern area of the Korean Peninsula is mostly affected by the influx of cold air from the Siberian high. These superimposed effects cause the higher energy term in SK. In general, the findings of the present study showed an increasing trend in ET o over the Korean Peninsula. However, this increasing trend is an undesirable outcome from water-resource perspective because it represents the increasing demand for water in the Korean Peninsula. Therefore, the results obtained from this research, which is the first study to determine ET o trends for the Korean Peninsula, could be helpful for the future planning of agriculture and water demands. Moreover, understanding how the trends of these terms change during seasons could provide insight into the possible impact of climate change on future water resources of the Korean Peninsula. Although this study is limited to the Korean Peninsula, these results might be applicable to East Asia as well. Moreover, further research on the impacts of meteorological parameters on ET o over the Korean Peninsula could be of interest.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4441/10/10/1373/ s1: Table S1: required equations for computing ET o using the FAO PM method; Table S2: numbers of missing data for each parameter for each station;