Next Article in Journal
The Effect of Short-Term Upwelling Events on Fish Assemblages at the South-Eastern Coast of the Baltic Sea
Previous Article in Journal
Discriminant Analysis of Water Inrush Sources in the Weibei Coalfield, Shaanxi Province, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Meteorological Influences on Reference Evapotranspiration in Different Geographical Regions

Department of Civil and Environmental Engineering, Dongguk University, Seoul 04620, Republic of Korea
*
Author to whom correspondence should be addressed.
Water 2023, 15(3), 454; https://doi.org/10.3390/w15030454
Submission received: 30 November 2022 / Revised: 13 January 2023 / Accepted: 19 January 2023 / Published: 23 January 2023
(This article belongs to the Section Hydrology)

Abstract

:
It is critical to understand how meteorological variables impact reference evapotranspiration ( ET o ) since meteorological variables have a different effect on ET o . This study examined the impact of meteorological variables on ET o on the Korean Peninsula under complex climatic and geographic conditions in order to understand how ET o and meteorological variables have changed over the past 42 years. Different geographical conditions were analyzed, including plains, mountains, and coastal areas on a seasonal and annual basis. ET o was estimated using the Penman-Monteith method by the Food and Agriculture Organization (FAO) using daily relative humidity ( RH ) , solar radiation ( R s ) , maximum temperature ( T max ) , minimum temperature ( T min ) , and wind speed ( W s ) . According to the results, the maximum mean seasonal and annual ET o occurred on the southern coast, while the minimum occurred in the mountainous area and along the east coast. Seasonal ET o is highest in summer, and is lowest in winter for all regions. The investigation of meteorological variables on ET o revealed that the response varied by area, and the magnitudes of sensitivity varied by location and season. RH is the most critical meteorological variable to affect ET o in all seasons, except summer, when T min is the most sensitive parameter. The results revealed that different regions showed different responses to the change in ET o by changing the meteorological variables. Meteorological variables affecting ET o differ with different geologic conditions and seasons. in mountainous areas revealed almost similar responses to the change in RH, R s , and T max (±10% change in ET o ) during the spring season. However, for other regions, RH and T max caused changes to ET o throughout, ranging from −15% to +20% in the plain area, −20% to +15% in the west and east coast, and −20% to +10% in the south coast. In addition, there were significant differences in parameter responses between regions and seasons, which was confirmed by the results.

1. Introduction

Evaporation (ET) is one of the main components of the water cycle, which has a significant impact on soil moisture [1]. ET is a crucial hydrological variable used for climate change impact studies, flood and drought analysis, irrigation scheduling, food and water security decisions by policy-makers, optimal water use, and developing hydrological water balance models [2,3,4,5].
Generally, ET exhibits an integrated response to key meteorological parameters, such as maximum temperature ( T max ) , minimum temperature ( T min ) , relative humidity ( RH ) , solar radiation ( R s ) , and wind speed ( W s ) at 2 m height [6,7]. ET is estimated as the crop coefficient (Kc) for a particular land use and land cover (LULC) multiplied by the reference evapotranspiration (ETo), which is computed for a reference crop when the evaporation behavior is similar to the surface of green grass at a height of 0.12 m, an albedo of 0.23, and a constant surface resistance of 70 s/m [8,9]. It is possible to directly measure the ET0 using the water balance approach and the lysimeter approach, while indirect estimation is possible with meteorological data [10,11]. The Food and Agriculture Organization Penman-Monteith (FAO PM) method has been suggested as the standard method for estimating ET o by the FAO, the American Society of Civil Engineering (ASCE) in Irrigation and Hydrology Committee, and the International Commission for Irrigation (ICID) [7,12,13]. However, FAO-56 PM is not widely applicable due to the lack of required meteorological data in many regions. Furthermore, directly measuring ET o using a lysimetric study is challenging. Thus, it is necessary to evaluate current ET o estimation methods on the basis of the FAO-56 PM model as a benchmark, which is useful for improving the efficiency of localized water management and crop water use [13].
A study of the trends of meteorological variables could provide insight into the impacts of key factors on ET o , as well as climate change [14,15]. Over the last few decades, numerous studies have identified global or regional trends in ET o , and have shown that the ET o trend varies by region. The upward trend of ET o was reported in China [16], Iran [17], Senegal [18], and in Brazil [19], as well as downward trends in India [20] and Africa [21].
It is necessary to improve our knowledge of the connection between meteorological variables and ET o , and to determine the role of each meteorological parameter in the change in ET o [22]. ET o is impacted by forcing variables, and these variables are sensitive to climate change [23]. Furthermore, an understanding of the sensitivity of ET o is needed to determine the required accuracy for measuring the meteorological variables used to estimate ET o [24]. Goyal [25] investigated the sensitivity of ET o in an arid zone in India. The results showed that ET o was most sensitive to air temperature. Liqiao et al. [26] analyzed the sensitivity of ET o to weather parameters in the Tao’er River Basin in a semi-arid region of China. According to their results, relative humidity was the most sensitive parameter in this area. A study conducted by Estevez et al. [27] in semi-arid regions of Spain showed that temperature and solar radiation were the most sensitive variables. Song et al. [28] analyzed distributions and trends in ET o in the relatively flat region of the North China Plain. The results indicated that decreasing solar radiation and wind speeds were the main impacts on ET o . Bakhtiari and Liaghat [9] considered a station in a semi-arid region of Iran and concluded that the vapor pressure deficit was the most sensitive parameter for ET o in this area. Yang et al. [29] examined the sensitivity of ET o to climatic parameters in arid and semi-arid regions of the Yellow River Basin, China. They found that relative humidity, followed by mean temperature, were the most sensitive parameters. Porter et al. [30] studied the sensitivity of ET o to weather parameters in a semi-arid zone of Bushland, Texas. Their results indicated wind speed and air temperature to be the most sensitive factors, and that special care is warranted in the siting, sensor placement, and maintenance of these parameters. Hou et al. [31] conducted relative change analysis at one station in the Ejina Oasis of the Heithe River in one of the most arid regions in China. Their results showed that shortwave radiation, followed by air temperature, were the most sensitive variables. Liu et al. [32] concluded that driving climatic factors varied for sub-regions of urban areas in China. The changes in temperature and relative humidity were caused to increase ET o in Taishan, Zhongshan, and Shenzhen, whereas the variations in sunshine hours and wind speed were responsible for decreasing ET o in Guangzhou and Zengcheng in China. Sharifi and Dinpashoh [22] conducted sensitivity analysis of eight stations in the arid and semi-arid regions of Iran. The study indicated that major meteorological variables changed based on the different climate zones or geographic conditions, while most previous studies were limited to stations in arid and semi-arid regions.
Therefore, it is of interest to us to investigate how the main meteorological variables impact ET o over an area with a temperate climate and complex geographic conditions. Moreover, the response of ET o to changing climatic parameters may be more crucial when climatic and geographic characteristics dramatically vary over a large-scale region. Therefore, the present study considered the Korean Peninsula as the domain, with complex climatic and geographic characteristics throughout the entire peninsula. Therefore, the aims of this paper were to analyze the impacts of meteorological variables on ET o in different geographic conditions using trend analysis, sensitivity analysis, and probability density function to determine the most dominate meteorological variables affecting ET o on seasonal and annual scales over different geographical conditions (plain, coastal, and mountainous areas).

2. Study Area and Data

The Korean Peninsula was chosen as the representative region for this study. The Korean Peninsula is a part of East Asia and consists of South Korea (45% of land area) and North Korea (55% of land area). Subtropical monsoons affect the Korean Peninsula, which is nearly 70% mountainous and bounded by three seas: the East Sea, the East China Sea, and the Yellow Sea, resulting in complex atmospheric conditions. The Korean Peninsula can be divided into five distinct regions: a mostly mountainous area in the northern part, the west coast in the western part, the east coast in the eastern part, the south coast in the south, and the plains in the rest of the peninsula. There are four distinct seasons on the Korean Peninsula: spring from March to May, summer from June to August, autumn from September to November, and winter from December to February. The climate significantly varies between seasons due to both continental and oceanic influences, with extremes in summer, impacted by pacific high pressure from the south, and winter, influenced by the cold Siberian air mass [33]. The mean annual precipitation in the Korean Peninsula increases southward, ranging from 1000 mm to 1800 mm. The mean annual temperature ranges from 22.5 °C to 25 °C in summer and from −5 °C to −2.5 °C in winter [34].
This study utilized the Korea Meteorological Administration website to obtain daily climatic data from 21 stations in South Korea and 27 stations in North Korea from 1980 to 2021. Table 1 shows the information on the downloaded meteorological data from the KMA website for both South and North Korea. Only 21 stations in South Korea contained data for observed solar radiation and were selected, and the stations in North Korea had no data for solar radiation. The location and description of the meteorological station is shown in Figure 1 and Table 2.

3. Methods

Figure 2 shows the methodology used in this study. The required meteorological data for South Korea and North Korea were extracted for the period 1980 to 2021. In order to estimate Rs over the Korean Peninsula, a regression equation was obtained for K Rs , and then Evapotranspiration was calculated over the stations of the Korean Peninsula. The temporal variations of ET o and meteorological variables were analyzed in different geographic conditions to determine how they had changed during the past 42 years in different seasons. Moreover, the Mann-Kendall test was conducted to find significant trends in ET o and meteorological variables. The sensitivity analysis and probability density function of the relative change in ET o were used to determine the most dominate meteorological variables affecting ET o on seasonal and annual scales. In the following sections, the method is explained in detail.

3.1. Estimation of R S

Whenever data on solar radiation or sunshine hours are unavailable, Allen et al. [12] recommend estimating solar radiation data using the formula developed by Hargreaves and Samani [35]:
R s = K Rs   T max   T min     R a
where T max (°C) is the maximum temperature, T min (°C) is the minimum temperature, R a represents the extraterrestrial solar radiation (MJ m−2 d−1), and K Rs is a coefficient of adjustment, which is calculated using a regression equation. In order to estimate an accurate value of solar radiation, K Rs needs to be calibrated for each station in the Statistical Package for the Social Sciences (SPSS) software as an independent variable, and the long-term observations of T mean TD ( T mean is the mean temperature, and TD is the difference between maximum and minimum temperatures) and EL (elevation) as dependent variables. The results of this analysis were reported in the Coefficient Table of SPSS software to determine the most suitable equation for K Rs . In SPSS, the significant value should be compared to 0.05 when testing the hypothesis. If the significant value is greater than 0.05, it means the null hypothesis is accepted. If the significant value is less than 0.05, it shows a rejection of the null hypothesis.
Since there are no R s data for North Korea, a regional calibrated model for K Rs was constructed for South Korea based on the observed R s data in South Korea. This formula was then used to estimate R s in both South Korea and North Korea. Using multiple regression analysis, Equation (2) was derived using the mean values over the long term of T mean TD and EL as the independent variables, and the values of K Rs as the dependent variables.
K Rs = 0.032 T mean TD + ( 2.29 × 10 5 EL ) + 0.105
By having the K Rs , Equation (1) was used to estimate daily R s over the Korean Peninsula. It should be noted that K Rs was regionally calibrated using observed solar radiation data in South Korea, and therefore, it may be subject to uncertainty in the estimated R s for stations in North Korea. Table 3 shows the comparison between calculated daily ET o by using observed and estimated daily R s during 1980–2021 for the stations in South Korea. Based on the results, the root mean square error (RMSE) for all stations was low and the determination coefficient (R2) was high, which determined the good correlations for ET o computed by observed and estimated R s .

3.2. Estimation of E T o

There are different empirical formulas available for estimations of ET o depending on the availability of meteorological data in the past [36]. When all required data were available, the FAO PM method described by Allen et al. [12] was directly used to estimate daily ET o . The input data required for the FAO PM method were determined using daily maximum air temperature ( T max ) , minimum air temperature ( T min ) , solar radiation ( R s ) , relative humidity ( RH ) , and wind speed ( W s ) . The FAO PM method is expressed in Equation (3):
ET o = 0.408 Δ ( R n G ) + γ 900 T + 273 W s ( e s e a ) Δ + γ ( 1 + 0.34 W s )
where ET o is the reference evapotranspiration [mm day−1], R n is the net radiation at the crop surface [MJ m−2 day−1], G is the soil heat flux density [MJ m−2 day−1], T is the mean daily air temperature at 2 m height [°C], W s is the wind speed at 2 m height [m s−1], e s is the saturation vapor pressure [kPa], e a is the actual vapor pressure [kPa], ( e s e a ) is the saturation vapor pressure deficit [kPa], Δ is the slope of the vapor pressure curve [kPa °C−1], and γ is a psychrometric constant [kPa °C−1]. Whenever there are missing data or they cannot be calculated, it is recommended that the user estimates the missing climatic data using one of the procedures described in Chapter 3, by Allen et al. [4]. While there are some alternative methods for calculating ET o that require fewer meteorological parameters, they are less recommended. The data on solar radiation was only available for South Korean stations, so we developed a regression equation to estimate the data for North Korea as explained in the next section.

3.3. The Impact of Meteorological Variables on E T o

In order to identify significant trends and to understand how ET o and meteorological variables have changed over the past 42 years on the Korean Peninsula, seasonal and annual trends of ET o and meteorological variables were analyzed. Sensitivity analysis was used to detect the possible change in ET o caused by variation in meteorological variables. In addition, the probability density function was used to provide the probability density function of ET o change due to the change in meteorological variables. Analysis was completed for both annual and seasonal scales. The standard slope in the linear regression (as defined by B) reflects the percentage of increase or decrease in ET o and meteorological variables in a time series.

3.3.1. Mann-Kendall Test

A Mann-Kendal (MK) test is a nonparametric statistical test which is recommended to determine if the trend in the dataset is significant over time (increasing or decreasing). The MK test was used in this study to assess the significance of seasonal and annual trends of ET o and meteorological variables in different geographic conditions. This method is commonly used to identify trends in meteorological variables, and can be found in [37] and [1]. The statistics of the Mann-Kendall test are determined as follows:
S = k = 1 n 1 j = k + 1 n sign ( x j x k )
sign ( x j x k ) = { + 1       if ( x j x k ) > 0 0       if ( x j x k ) = 0 1       if ( x j x k ) < 0  
Var   ( S ) = [ n ( n 1 ) ( 2 n + 5 ) i = 1 m t i ( t i 1 ) ( 2 t i + 5 ) ] 18
Z = { S 1 Var   ( S )       if   S > 0 0                   if   S = 0 S + 1 Var   ( S )         if   S < 0  
where S denotes the test statistics, n indicates the length of the data set, x j and x k indicate the sequential values, t i is the number of ties of extent i, m is the number of tied groups, and Z is the standardized Mann-Kendall statistic. A pre-whitening method was also used to eliminate serial correlations in the time series of ET o and its climate variables [38,39].

3.3.2. Sensitivity Analysis

Sensitivity analysis is widely used to identify the changes in the dependent variable ( ET o ) caused by the change in an independent meteorological variable (e.g., [40,41,42,43]). For compound models (such as the FAO PM model), it is difficult to compare sensitivity based on partial derivatives, since meteorological variables have different dimensions [23]. As a result, a dimensionless index is derived from the partial derivative SC i [44]:
SC i = lim Δ X i 0 ( Δ ET o / ET o Δ X i / X i ) = ET o X i · X i ET o
where SC i represents the sensitivity coefficient, Δ X i is the actual change of meteorological variable X i , and Δ ET o is the actual change in ET o induced by Δ X i . The sign of SC i determines how ET o responds to codirectional changes in input parameters. The sign of SC i determines how ET o would codirectionally react to the input parameter change. For example, a positive SC i of a meteorological parameter indicates that ET o will increase as the variable increases. In addition, the absolute value of SC i indicates the magnitude of ET o response to that meteorological variable [42,43,44,45].
In this study, sensitivity analysis of ET o was calculated for meteorological variables (i.e.,   RH ,   R s , T max , T min , and W s ) for ±5%, ±10%, ±15%, and ±20% for each of the variables, while keeping other meteorological variables constant. In order to calculate the mean seasonal and annual SC i , daily values were averaged for different geographical conditions on seasonal and annual scales.
The average values were calculated using geospatial interpolation instead of simple averaging. Spatial averaging involves different methods, including inverse distance weighting (IDW) and kriging. It is used for spatial averaging along surfaces. The IDW combines a set of sample points linearly weighted. Weight is determined by inverse distance. A kriging method, however, generates an estimated surface from a scattering of points by computing kriging weights. A kriging weight value was obtained for each station in this study using the geographical information system (GIS) Spatial Analyst. In this study, the kriging method [46] was used to interpolate the mean seasonal and annual sensitivity coefficients of all stations. This method is widely used by researchers for interpolating data to have the mean value or spatial distribution of a station [47,48,49]. A higher mean value of the sensitivity coefficient indicates that ET o is more sensitive to changes in meteorological variables.

3.3.3. Probability Density Function

The kernel distribution produces the probability density function (PDF) by summing the smooth curves for each data value, which creates a smooth, continuous probability density function for the dataset. In fact, in a non-parametric distribution, the density is entirely determined by the data without any strict distributional assumptions [50]. In other words, nonparametric statistical procedures do not estimate parameters regarding the shape or form of the PDF. Kernel distributions provide nonparametric probability density estimates rather than selecting a density with a specific parametric form and estimating its parameters. The density curve is generated by a kernel density estimator, a smoothing function that determines the shape of the PDF, along with a bandwidth value that controls its smoothness. In the present study, Kernel Density Estimation, which is a nonparametric technique for density estimation, was used to provide the probability density function of percent ET o change due to the percent change in meteorological variables. Kernel density estimation for the relative change of ET o from −20% to +20% change in each meteorological variables was constructed using following equation. Given a random X 1 ,   ,   X n with a continuous univariate density f, the kernel density estimator is [51]:
f ^ ( x , h ) = 1 nh i = 1 n K ( x X i h )
where K is kernel, f ^ ( x , h ) is the height of the curve at x (percent change of ET o ), x is relative change of meteorological variable from −20% to 20%, K (.) is the standard normal density, and h is the bandwidth of the density curve.

4. Results

4.1. Mean Seasonal and Annual E T o and Meteorological Variables

Temporal variation in ET o and meteorological variables over the different regions (i.e., mountainous area, plain area, west coast, east coast, and south coast) are plotted in Figure 3, and the statistical results are reported in Table 4. The trend in the slope of the mean annual ET o is highest in the east coast compared to the other regions. There is an increasing trend in the linear trend of mean annual ET o over the Korean Peninsula from 1980 to 2021 at a rate of B = 1.43   mm / yr 2 in the east coast, 1.09   mm / yr 2 in the plain area, and 0.78   mm / yr 2 in the mountainous area; however, the south coast revealed a slightly decreasing trend for annual ET o . By neglecting the slightly increasing trend for ET o in the winter season in the east coast, it could be concluded that the linear trend on the seasonal scale increases in all seasons except winter, which exhibits a decreasing trend range from 0.24   mm / yr 2 in the south coast to 0.19   mm / yr 2 in the west coast. Moreover, summer has the greatest seasonal rate of change in ET o for all regions. Increasing trends in ET o in South Korea are consistent with the results of a study conducted by Aydin et al. [52], which defined the increasing trend for pan evaporation and potential evapotranspiration in South Korea from 1980 to 2009.
According to Table 4, the highest mean annual ET o occurred in the south coast (1000.39 mm/yr), west coast (934.47 mm/yr), and then the plain area (912.28 mm/yr); however, the mean annual ET was lowest along the east coast (813.50 mm/yr) and in mountainous areas (912.28 mm/yr). In all regions, summer had the highest seasonal ET o value, and winter had the lowest.
There was a higher mean value of RH in coastal regions, particularly on the west and east coasts, while the lowest relative humidity was found in mountainous areas. In the mountainous area and on the east coast, the linear trend showed a decreasing trend for seasonal and annual relative humidity. Water availability decreased in the northern and eastern parts of this area, while it increased in other areas of the Korean Peninsula. According to these results, the RH trend over the Korean Peninsula displayed conflicting behavior.
The results of the linear trend of R s showed increasing trends on seasonal and annual scales over the Korean Peninsula (0.01–0.02 MJm−2d−1yr−1). The maximum increasing trends of R s occurred in summer and then in the spring over the Korean Peninsula, with the higher value in the mountainous area. The increasing trend in R s was reported by Russak [53], which concluded that the increase in R s trend could be a result of decreasing trends in cloudiness and aerosol thickness.
Moreover, the time series of T max and T min exhibited an increasing trend, with a higher slope for T min than that for T max in the plain area (autumn and winter) and west coast (during all seasons except spring). Based on this pattern, TD (amount of difference between a maximum and a minimum temperature) declined, suggesting that global warming has become more influential in this region, due to increased urbanization and higher aerosol concentrations in this area. However, in other areas, a higher slope for T max than T min caused the increasing trend for TD, and therefore, the decreasing trend for RH . Because TD is related to RH [54] and cloud effects, accordingly, TD decreases as cloudiness increases, and is correlated with a deficit in vapor pressure [55]. Therefore, decreasing trends in TD caused an increasing trend in RH in the west coast and plain area, while the increasing trend in TD caused a decreasing trend in RH , especially in the mountainous area.
According to the results of the trend analysis of wind speed over the Korean Peninsula, W s showed almost constant downward trends on both an annual and seasonal basis. W s showed the same slope of decreasing trend (−0.01 m s−1), with a higher decreasing rate in the south coast (−0.02 m s−1) during the spring and winter seasons.
Overall, the trend in ET o reflects the combined impact of all meteorological variables. Generally, the highest ET o has been recorded in the southern part of the Korean Peninsula. The rising ET o trend on Korea’s peninsula over the past 42 years may have been primarily caused by significant increases in T max and decreases in RH. Furthermore, the outweighing increases in climate parameters in the southern part of the peninsula result in a higher increase in ET o in comparison to the northern part.

4.2. Significant Trends of E T o and Meteorological Variables on Seasonal and Annual Scale

Mann-Kendall tests were conducted on pre-whitened time series in order to assess how ET o trends and meteorological variables have changed over the past 42 years. Figure 4 depicts the spatial distribution of trends at stations on an annual scale. The results for the seasonal and annual scales are also shown in Table 5, which shows how many stations in each region show significant increases, significant decreases, and non-significant trends.
ET o showed significant increases mainly along the east coast, plain area, and south coast; however, the mountainous region did not show any significant changes. The significant increasing ET o trends on the east coast were mainly driven by a significant decrease in RH . Furthermore, the significant increasing trends in the plain area and south coast occurred due to the most striking trends in R s , T max , and T min . Moreover, most stations in the south coast of the peninsula indicated a significant decreasing trend, which mainly occurred because of significant decreasing trends in R s , T max , and W s .

4.3. Sensitivity Analysis of Meteorological Variables to Change E T o over the Korean Peninsula

The daily sensitivity coefficients of ET o with respect to the key meteorological variables was calculated for all stations, and then the station value sensitivity coefficient was obtained for all meteorological variables on both a seasonal and annual basis. The sensitivity coefficient can be used to quantitatively analyze changes in ET o by changing one variable, while other variables remain constant. Percentage changes in ET o relative to changes in meteorological variables for each season and annual were obtained.
The bar charts in Figure 5 show the seasonal and annual sensitivity coefficients of meteorological variables over the Korean Peninsula, and the statistics are reported in Table 6. The spatial distribution of sensitivity coefficients indicated that the sensitive parameters are different in each season. According to the spatial distribution of the sensitivity coefficient, each season has different sensitive parameters. In spring and autumn, ET o responds almost in the same way to changes in meteorological variables, but, in summer and winter, this response is significantly different.
The most sensitive parameters to change in spring and autumn are T max and RH ; however, in summer, T max , T min , and R s have the greatest impact on ET o . In winter, RH is the key parameters to affect ET o and the northern part of the Korean Peninsula is more sensitive to change in RH than the southern part. However, some stations in the mountainous area are more sensitive to T min because the northern part of the peninsula is more impacted by Siberia’s cold weather in winter.
In general, the following results were obtained by sensitivity analysis: (i) T max , followed by RH , are the most sensitive parameters in spring for all regions, with very high sensitivity in the south coast and west coast. (ii) In summer, the coastal areas have high sensitivity to the change in T min (especially the south coast, with SC= 1.12); however, the mountainous and plain areas are the most sensitive to R S and T max . (iii) In autumn, T max is the main parameter that influences ET o , followed by RH , with very high sensitivity in the south coast and the west coast. However, in other regions, the increase in T max   followed by a decrease in RH is the key change that affects ET o in autumn. (iv) Winter shows the highest sensitivity among the seasons for all regions. In fact, the winter season is the driest season in the Korean Peninsula and ET o is very sensitive to the change in RH , with high sensitivity in the coastal regions. (v) On an annual scale, the average value of sensitivity analysis coefficients can be ranked as S RH S T max > S Rs > S T min S Ws in the mountainous and plain areas, and S RH > S T max > S Rs > S T min > S Ws in the coastal areas. Overall, analysis showed that the two major components to change ET o are RH and T max over the Korean Peninsula.
Figure 6 illustrates the amount of change in ET o by changing the meteorological variables from −20% to +20 with a 5% interval. Based on the relative change in ET o , different geographical conditions respond differently to changes in meteorological variables, with coastal areas exhibiting higher variability than mountainous and plain areas. Additionally, meteorological variables vary from season to season. In general, RH responds most strongly to changes in ET o except in the summer season. A high-pressure system over the western North Pacific and the Indian summer monsoon contribute to high relative humidity during summer monsoons. As a result, ET o is not sensitive to changes in RH in summer; however, other seasons, especially winter, which is the driest season on the Korean Peninsula, are highly sensitive to changes in RH.

4.4. Probability Density Function of Meteorological Variables to Change E T o over the Korean Peninsula

The probability density function (PDF) of percent change of ET o due to changes in meteorological parameters on seasonal and annual scales in different regions was obtained using a Kernel Density Estimator and is illustrated in Figure 7. The wider shape of the distribution from the normal distribution indicated more variability of meteorological parameters to change ET o . The relative change in ET o caused by percent change of each meteorological variables from −20% to +20% with intervals of 5% was shown as the probability density function during the seasonal and annual scales (Figure 7), which demonstrated the impact of meteorological variables varied during the seasons and throughout regions. In fact, changes in meteorological variables caused changes in the mean, variance, and shape of the distribution, which were different in each season and each geographic condition. According to the results, ET o responded differently to changes in meteorological variables based on geographic location. The shape of PDF near the normal distribution shows lower sensitivity, while dispersion from the normal distribution (mean) indicates the higher sensitivity of ET o by changing meteorological variables. In spring, the shape of the distribution was shifted from a normal distribution for T max and RH in spring, especially for the coastal area, indicating higher sensitivity of ET o to change those variables in this season. In summer, T max , T min , and R S , and in autumn, T max and RH , are more dispersed from the mean, which shows the higher sensitivity of ET o to change them. However, in winter, RH is the most sensitive parameter for all geographic conditions (with a higher sensitivity in coastal areas).

5. Discussion

The Korean Peninsula is affected by surrounded seas. The peninsula is also impacted by the Pacific Ocean over the Japanese archipelago. Therefore, the climatic conditions of the peninsula are impacted by both the ocean and the continent, leading to significant seasonal variability, such as monsoons [56]. In fact, the summer season in the Korean Peninsula is warm and wet, which is driven by the ocean, while the winter season is cold and dry, driven by the Siberian cold air mass (continental effect). Temporal variations of the ET o showed that the southern coast recorded the highest mean seasonal and annual values. Similar to the mountainous area, the east coast had the lowest ET o on a seasonal and annual basis. The results showed that meteorological variables have different tendencies to change ET o over the Korean Peninsula.
The results of the study by Nie et al. [57] showed that the sensitivity coefficient and contribution rate of ET o to meteorological factors vary by climatic region in China. Sensitivity coefficients and contribution rates to ET o spatially vary even within the same climatic region. In this study, the slope of the temporal trends indicates a decreasing trend for the south coast and an increasing trend for the east coast. The reason for the decreasing trend in ET o on the south coast is the increasing trend in RH, which is affected by TD. In fact, on the south coast, the rising trend of T min is higher than T max , leading to a decreasing trend in TD, and therefore, an upward trend in RH. The increasing trend in ET o due to a deficiency in RH was reported in Iran [9] and China [26,29]. In mountainous areas, affected by continental influence in the northern part of the Korean Peninsula, the mean RH is the lowest. While in coastal areas, particularly on the east coast and south coast, RH is the highest. R s in spring and summer has the highest value in the mountainous area in the northern Korean Peninsula, while autumn has the maximum value in the south coast and the minimum value in the east coast. The mean annual and seasonal temperature increases from north to south. R s and W s have the main influence on changing ET o in the North China plain [28]. Furthermore, the mean T min and T max are highest in the south coast and decline toward the west coast and east coast. The higher temperature in the southern part of the Korean Peninsula is related to the Tsushima Warm Current. However, the mountainous area has the minimum value of maximum and minimum temperature, since this area is affected by the North Korea Cold Current. The sensitivity of ET o to change T mean was also reported in a semi-arid region in India [25]. Moreover, the W s showed decreasing trends among all seasons and regions. Significant decreases in wind speed have recently been reported in different countries. Land use changes and increased surface roughness near meteorological stations have been attributed to the decrease in W s trend. Similar results were reported in [28,30]. The weakening trend in W s may occur due to temperature increases due to global change. According to the findings of Xiaomei et al. [58], the increasing trends in T max and T min cause decreases in W s , and it may be concluded that changes in temperature indirectly influence other meteorological variables.
This study was based on a 42-year historical dataset, and the impact of climate change was not taken into consideration. Therefore, we suggest following the same methodology by using global climate models to examine the impact of the main meteorological variables on ET o under different geographical conditions. Furthermore, we recommend assessing the risk and vulnerability by taking into account the change in ET o , and implementing the necessary strategy and deployment plan to ensure a sustainable water supply. Furthermore, we recommend conducting a similar study to compare the impact of meteorological variables to change water budgets, including evapotranspiration, precipitation, and discharge changes in different geographical conditions.

6. Conclusions

In this study, we examined the impacts of meteorological variables on ET o across different geographic conditions on the Korean Peninsula, using a 42-year daily dataset. It was found that seasonal and regional changes were influenced differently by different meteorological variables. In the summer, coastal areas showed higher sensitivity to the increase in T min . In contrast, autumn and spring were the most sensitive to the increase in T max . This indicated that, for analyzing the impact of the meteorological variables, T mean is not the appropriate variable and the influence of both T max and T min should be considered. In all regions, winter was the most sensitive to changes in RH , particularly on the east and west coasts. Among all Korean Peninsula regions, winter and summer were the most sensitive seasons. The magnitudes of sensitive parameters during each season varied depending on location, according to the spatial distribution of those parameters. As a result of the study, it can be determined that different seasons have different responses to changes in maximum and minimum temperatures, and, for a better understanding of the change of ET o to changes in temperature, both the maximum and minimum temperatures should be considered instead of T mean . Moreover, comparison between the sensitive coefficients of the seasonal and annual scales showed that it is not enough to only consider the impact of sensitive parameters for the annual scale. Although it revealed the key sensitive parameters in each region, different parameters have different responses in each season and in each region.

Author Contributions

Conceptualization, M.G.-A.; Data curation M.G.-A.; Formal analysis, M.G.-A.; Funding acquisition, S.-I.L.; Methodology, M.G.-A.; Project administration, S.-I.L.; Supervision, S.-I.L.; Visualization, M.G.-A.; Writing—original draft, M.G.-A.; Writing—review and editing, M.G.-A. and S.-I.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant by the Korea government (2021R1A2C2011193).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gocic, M.; Trajkovic, S. Analysis of trends in reference evapotranspiration data in a humid agricultural water management climate. Hydrol. Sci. J. 2014, 59, 65–180. [Google Scholar] [CrossRef]
  2. Valipour, M.; Bateni, S.M.; Gholami Sefidkouhi, M.A.; Raeini-Sarjaz, M.; Singh, V.P. Complexity of Forces Driving Trend of Reference Evapotranspiration and Signals of Climate Change. Atmosphere 2020, 11, 1081. [Google Scholar] [CrossRef]
  3. Luo, Y.; Gao, P.; Mu, X. Influence of Meteorological Factors on the Potential Evapotranspiration in Yanhe River Basin, China. Water 2021, 13, 1222. [Google Scholar] [CrossRef]
  4. McNally, A.; Arsenault, K.; Kumar, S.; Shukla, S.; Peterson, P.; Wang, S.; Funk, C.; Peters-Lidard, C.D.; Verdin, J.P. A land data assimilation system for sub-Saharan Africa food and water security applications. Sci. Data 2017, 4, 170012. [Google Scholar] [CrossRef] [Green Version]
  5. Muhammad Adnan, R.; Chen, Z.; Yuan, X.; Kisi, O.; El-Shafie, A.; Kuriqi, A.; Ikram, M. Reference Evapotranspiration Modeling Using New Heuristic Methods. Entropy 2020, 22, 547. [Google Scholar] [CrossRef]
  6. Zhang, X.; Kang, S.; Zhang, L.; Liu, J. Spatial variation of climatology monthly crop reference evapotranspiration and sensitivity coefficients in Shiyang river basin of northwest China. Agric. Water Manag. 2010, 97, 1506–1516. [Google Scholar] [CrossRef]
  7. Sahoo, B.; Walling, I.; Deka, B.C.; Bhatt, B.P. Standardization of Reference Evapotranspiration Models for a Subhumid Valley Rangeland in the Eastern Himalayas. J. Irrig. Drain. Eng. 2012, 138, 880–896. [Google Scholar] [CrossRef]
  8. Srivastava, A.; Sahoo, B.; Raghuwanshi, N.S.; Singh, R. Evaluation of Variable-Infiltration Capacity Model and MODIS-Terra Satellite-Derived Grid-Scale Evapotranspiration Estimates in a River Basin with Tropical Monsoon-Type Climatology. J. Irrig. Drain. Eng. 2017, 143, 04017028. [Google Scholar] [CrossRef] [Green Version]
  9. Bakhtiari, B.; Liaghat, A.M. Seasonal Sensitivity Analysis for Climatic Variables of ASCE Penman-Monteith Model in a Semi-Arid Climate. J. Agric. Sci. Technol. 2011, 13, 1135–1145. [Google Scholar]
  10. Pandey, P.K.; Dabral, P.P.; Pandey, V. Evaluation of reference evapotranspiration methods for the northeastern region of India. Int. Soil Water Conserv. Res. 2016, 4, 52–63. [Google Scholar] [CrossRef] [Green Version]
  11. Pandey, P.K.; Nyori, T.; Pandey, V. Estimation of reference evapotranspiration using data driven techniques under limited data conditions. Model. Earth Syst. Environ. 2017, 3, 1449–1461. [Google Scholar] [CrossRef]
  12. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; Food and Agriculture Organization of the United Nations: Rome, Italy, 1998. [Google Scholar]
  13. Vicente-Serranoa, S.M.; Azorin-Molinaa, C.; Sanchez-Lorenzoa, A.; Revuelto, J.; López-Morenoa, J.I.; González-Hidalgoc, J.C.; Moran-Tejedaa, E.; Espejod, F. Reference evapotranspiration variability and trends in Spain, 1961–2011. Glob. Planet. Chang. 2014, 121, 26–40. [Google Scholar] [CrossRef] [Green Version]
  14. Nouri, M.; Homaee, M.; Bannayan, M. Spatiotemporal reference evapotranspiration changes in humid and semi-arid regions of Iran: Past trends and future projections. Theor. Appl. Climatol. 2018, 133, 361–375. [Google Scholar] [CrossRef]
  15. Nouri, M.; Homaee, M.; Bannayan, M. Quantitative trend, sensitivity and contribution analyses of reference evapotranspiration in some arid environments under climate change. Water Resour. Manag. 2017, 31, 2207–2224. [Google Scholar] [CrossRef]
  16. Yan, X.; Mohammadian, A.; Ao, R.; Liu, J.; Chen, X. Spatiotemporal Variations in Reference Evapotranspiration and Its Contributing Climatic Variables at Various Spatial Scales across China for 1984–2019. Water 2022, 14, 2502. [Google Scholar] [CrossRef]
  17. Dinpashoh, Y.; Babamiri, O. Trends in reference crop evapotranspiration in Urmia Lake basin. Arab. J. Geosci. 2020, 13, 372. [Google Scholar] [CrossRef]
  18. Ndiaye, P.M.; Bodian, A.; Diop, L.; Deme, A.; Dezetter, A.; Djaman, K.; Ogilvie, A. Trend and Sensitivity Analysis of Reference Evapotranspiration in the Senegal River Basin Using NASA Meteorological Data. Water 2020, 12, 1957. [Google Scholar] [CrossRef]
  19. Da Silva, V.D.P.R. On climate variability in northeast of Brazil. J. Arid Environ. 2004, 58, 575–596. [Google Scholar] [CrossRef]
  20. Darshana; Pandey, A.; Pandey, R.P. Analysing trends in reference evapotranspiration and weather variables in the Tons River Basin in Central India. Stoch. Environ. Res. Risk Assess. 2013, 27, 1407–1421. [Google Scholar] [CrossRef]
  21. Obada, E.; Alamou, E.A.; Chabi, A.; Zandagba, J.; Afouda, A. Trends and Changes in Recent and Future Penman-Monteith Potential Evapotranspiration in Benin (West Africa). Hydrology 2017, 4, 38. [Google Scholar] [CrossRef] [Green Version]
  22. Sharifi, A.; Dinpashoh, Y. Sensitivity Analysis of the Penman-Monteith reference Crop Evapotranspiration to Climatic Variables in Iran. Water Res. Manag. 2014, 28, 5465–5476. [Google Scholar] [CrossRef]
  23. Seong, C.; Sridhar, V.; Billah, M.M. Implications of potential evapotranspiration methods for streamflow estimations under changing climatic conditions. Int. J. Climatol. 2018, 38, 896–914. [Google Scholar] [CrossRef]
  24. Irmak, S.; Payero, J.O.; Martin, D.L.; Irmak, A.; Howell, T.A. Sensitivity analyses and sensitivity coefficients of standardized daily ASCE-Penman-Monteith equation. J. Irrig. Drain. Eng. 2006, 132, 564–578. [Google Scholar] [CrossRef]
  25. Goyal, R.K. Sensitivity of evapotranspiration to global warming: A case study of arid zone of Rajasthan (India). Agric. Water Manag. 2004, 69, 1–11. [Google Scholar] [CrossRef]
  26. Liqiao, L.; Li, L.; Zhang, L.; Li, J.; Li, B. Sensitivity of Penman-Monteith reference crop evapotranspiration in Tao’er River Basin of northeastern China. Chin. Geogr. Sci. 2008, 18, 340–347. [Google Scholar] [CrossRef]
  27. Estévez, J.; Gavilán, P.; Berengena, J. Sensitivity analysis of a Penman–Monteith type equation to estimate reference evapotranspiration in southern Spain. Hydrol. Process. 2009, 23, 3342–3353. [Google Scholar] [CrossRef]
  28. Song, Z.; Zhang, H.; Snyder, R.; Anderson, F.; Chen, F. Distribution and trends in reference evapotranspiration in the North China Plain. J. Irrig. Drain. Eng. 2010, 136, 240–247. [Google Scholar] [CrossRef]
  29. Yang, Z.; Liu, Q.; Cui, B. Spatial distribution and temporal variation of reference evapotranspiration during 1961–2006 in the Yellow River Basin, China. Hydrol. Sci. J. 2011, 56, 1015–1026. [Google Scholar] [CrossRef] [Green Version]
  30. Porter, D.; Gowda, P.; Marek, T.; Howell, T.; Moorhead, T.J.; Irmak, S. Sensitivity of grass-and alfalfa-reference evapotranspiration to weather station sensor accuracy. Appl. Eng. Agric. 2012, 28, 543–549. [Google Scholar] [CrossRef]
  31. Hou, L.G.; Zou, S.B.; Xiao, H.L.; Yang, Y.G. Sensitivity of the reference evapotranspiration to key climatic variables during the growing season in the Ejina oasis northwest China. SpringerPlus 2013, 2 (Suppl. 1), S4. [Google Scholar] [CrossRef] [Green Version]
  32. Liu, Y.; Tang, G.; Wu, L.; Wu, Y.; Yang, M. Variations in reference evapotranspiration and associated driving forces in the Pearl River Delta of China during 1960–2016. J. Meteorol. Soc. Jpn. 2019, 97, 467–479. [Google Scholar] [CrossRef]
  33. Hwang, J.H.; Van, S.P.; Choi, B.-J.; Chang, Y.S.; Kim, Y.H. The physical processes in the Yellow Sea. Ocean Coast. Manag. 2014, 102, 449–457. [Google Scholar] [CrossRef]
  34. Teng, S.K.; Yu, H.; Tang, Y.; Tong, L.; Choi, C.I.; Kang, D.; Liu, H.; Chun, Y.; Juliano, R.O.; Rautalahti-Miettinen, E.; et al. Yellow Sea; GIWA Regional Assessment 34; University of Kalmar: Kalmar, Sweden, 2005. [Google Scholar]
  35. Hargreaves, G.H.; Samani, Z. Estimating Potential Evapotranspiration. J. Irrig. Drain. Div. 1982, 108, 225–230. [Google Scholar] [CrossRef]
  36. Xu, J.; Haginoya, S.; Saito, K.; Motoya, K. Surface heat balance and pan evaporation trends in eastern Asia in the period 1971–2000. Hydrol. Process. 2005, 19, 2161–2186. [Google Scholar] [CrossRef]
  37. Hirsch, R.M.; Slack, J.R.; Smith, R.A. Techniques of trend analysis for monthly water quality. Water Resour. Res. 1982, 18, 107–121. [Google Scholar] [CrossRef] [Green Version]
  38. Von Storch, H.; Navarra, A. Analysis of Climate Variability-Applications of Statistical Techniques; Springer: Berlin/Heidelberg, Germany, 1999. [Google Scholar]
  39. Jhajharia, D.; Dinpashoh, Y.; Kahya, E.; Singh, V.P.; Fakheri-Fard, A. Trends in reference evapotranspiration in the humid region of northeast India. Hydrol. Process. 2012, 26, 421–435. [Google Scholar] [CrossRef]
  40. Coleman, G.; DeCoursey, D.G. Sensitivity and model variance analysis applied to some evaporation and evapotranspiration models. Water Resour. Res. 1976, 12, 873–879. [Google Scholar] [CrossRef]
  41. Hupet, F.; Vanclooster, M. Effect of the sampling frequency of meteorological variables on the estimation of the reference evapotranspiration. J. Hydrol. 2001, 243, 192–204. [Google Scholar] [CrossRef]
  42. Yin, Y.; Wu, S.; Chen, G.; Dai, E. Contribution analyses of potential evapotranspiration changes in China since the 1960s. Theor. Appl. Climatol. 2010, 101, 19–28. [Google Scholar] [CrossRef]
  43. Patle, G.T.; Singh, D.K. Sensitivity of annual and seasonal reference crop evapotranspiration to principal climatic variables. J. Earth Syst. Sci. 2015, 124, 819–828. [Google Scholar] [CrossRef] [Green Version]
  44. McCuen, R.H. A sensitivity and error analysis CF procedures used for estimating evaporation. J. Am. Water Resour. Assoc. 1974, 10, 486–497. [Google Scholar] [CrossRef]
  45. Gong, L.; Xu, C.Y.; Chen, D.; Halldin, S.; Chen, Y.D. Sensitivity of the Penman–Monteith reference evapotranspiration to key climatic variables in the Changjiang (Yangtze River) basin. J. Hydrol. 2006, 329, 620–629. [Google Scholar] [CrossRef]
  46. Journel, A.G.; Huijbregts, C.J. Mining Geostatistics; Academic Press: New York, NY, USA, 1978. [Google Scholar]
  47. Oliver, M.A. Kriging: A Method of Interpolation for Geographical Information Systems. Int. J. Geogr. Inf. 1990, 4, 313–332. [Google Scholar] [CrossRef]
  48. Xu, C.Y.; Gong, L.; Jiang, T.; Chend, D.; Singh, V.P. Analysis of spatial distribution and temporal trend of reference evapotranspiration and pan evaporation in Changjiang (Yangtze River) catchment. J. Hydrol. 2006, 327, 81–93. [Google Scholar] [CrossRef]
  49. Liu, Q.; Yang, Z.; Cui, B.; Sun, T. The temporal trends of reference evapotranspiration and its sensitivity to key meteorological variables in the Yellow River Basin, China. Hydrol. Process. 2010, 24, 2171–2181. [Google Scholar] [CrossRef]
  50. Murphy, K. Machine Learning: A Probabilistic Perspective; MIT Press: Cambridge, MA, USA, 2012; p. 16. [Google Scholar]
  51. Wand, M.P.; Jones, M.C. Kernel Smoothing, Monographs on Statistics and Applied Probability; Chapman & Hall: Boca Raton, FL, USA, 1995. [Google Scholar]
  52. Aydin, M.; Jung, Y.S.; Yang, J.E.; Kim, S.J.; Kim, K.D. Sensitivity of soil evaporation and reference evapotranspiration to climatic variables in South Korea. Turk. J. Agric. For. 2015, 39, 652–662. [Google Scholar] [CrossRef]
  53. Russak, V. Changes in solar radiation and their influence on temperature trend in Estonia (1955–2007). J. Geophys. Res. 2009, 114, D00D01. [Google Scholar] [CrossRef]
  54. Verhoef, A.; Feddes, R.A. Preliminary Review of Revised FAO Radiation and Temperature Methods; FAO Report 16; Landbouwuniversiteit Wageningen: Wageningen, The Netherlands, 1991; pp. 1–116. [Google Scholar]
  55. Hargreaves, G.H.; Allen, R.G. History and evaluation of Hargreaves evapotranspiration equation. J. Irrig. Drain. Eng. 2003, 129, 53–63. [Google Scholar] [CrossRef]
  56. Youn, Y.H. The climate variabilities of air temperature around the Korean Peninsula. Adv. Atmos. Sci. 2005, 22, 575–584. [Google Scholar] [CrossRef]
  57. Nie, T.; Yuan, R.; Liao, S.; Zhang, Z.; Gong, Z.; Zhao, X.; Chen, P.; Li, T.; Lin, Y.; Du, C.; et al. Characteristics of Potential Evapotranspiration Changes and Its Climatic Causes in Heilongjiang Province from 1960 to 2019. Agriculture 2022, 12, 2017. [Google Scholar] [CrossRef]
  58. Xiaomei, Y.; Zongxing, L.; Qi, F.; Yuanqing, H.; Wenlin, A.; Wei, Z.; Weihong, C.; Tengfei, Y.; Yamin, W.; Theakstone, W.H. The decreasing wind speed in southwestern China during 1969–2009, and possible causes. Quat. Int. 2012, 263, 71–84. [Google Scholar] [CrossRef]
Figure 1. Study area.
Figure 1. Study area.
Water 15 00454 g001
Figure 2. Workflow diagram.
Figure 2. Workflow diagram.
Water 15 00454 g002
Figure 3. Temporal variation of ET o and meteorological variables on seasonal and annual scales in different geographic conditions (1980–2021).
Figure 3. Temporal variation of ET o and meteorological variables on seasonal and annual scales in different geographic conditions (1980–2021).
Water 15 00454 g003
Figure 4. Spatial distribution of the annual trend of stations for ET o and meteorological variables over the Korean Peninsula (1980–2021). Upward and downward triangles define increasing and decreasing trends, respectively. The pink color indicates 95% confidence level, the green color indicates 90% confidence level, and white color indicates the trend is not significant.
Figure 4. Spatial distribution of the annual trend of stations for ET o and meteorological variables over the Korean Peninsula (1980–2021). Upward and downward triangles define increasing and decreasing trends, respectively. The pink color indicates 95% confidence level, the green color indicates 90% confidence level, and white color indicates the trend is not significant.
Water 15 00454 g004
Figure 5. Spatial distribution of the sensitivity coefficients of ET o to meteorological variables over the Korean Peninsula.
Figure 5. Spatial distribution of the sensitivity coefficients of ET o to meteorological variables over the Korean Peninsula.
Water 15 00454 g005
Figure 6. Relative change in ET o with respect to the percent change of meteorological variables on seasonal and annual scales in different geographic conditions.
Figure 6. Relative change in ET o with respect to the percent change of meteorological variables on seasonal and annual scales in different geographic conditions.
Water 15 00454 g006
Figure 7. Probability density of changes ET o with respect to changes in meteorological variables on seasonal and annual scales in different geographic conditions.
Figure 7. Probability density of changes ET o with respect to changes in meteorological variables on seasonal and annual scales in different geographic conditions.
Water 15 00454 g007
Table 1. The list of the extracted data for this study.
Table 1. The list of the extracted data for this study.
Weather Observation Data by CountryVariablesPeriodLink
South Korea RH 1980–2021https://data.kma.go.kr/data/grnd/selectAsosRltmList.do?pgmNo=36
(accessed on 20 May 2022).
R s 1980–2021
T max 1980–2021
T min 1980–2021
W s 1980–2021
North Korea RH 1980–2021https://data.kma.go.kr/data/grnd/selectNkRltmList.do?pgmNo=58
(accessed on 20 May 2022).
T max 1980–2021
T min 1980–2021
W s 1980–2021
Table 2. Characteristics of the meteorological stations used in this study.
Table 2. Characteristics of the meteorological stations used in this study.
Station CodeStation NameLat. (°)Lon. (°)GeographyEle. (m)Station CodeStation
Name
Lat. (°)Lon. (°)GeographyEle. (m)
3Seonbong42.3130.4EC369Haeju38.0125.7WC81
5Samjiyeon41.8128.3M138670Gaeseong38.0126.6P70
8Cheongjin41.8129.8EC4375Pyeonggang38.4127.3P371
14Junggang41.8126.9P332100Daegwalryeong37.7128.7M842.5
16Hyesan41.4128.2M714101Chuncheon37.9127.7P75.6
20Ganggye41.0126.6P306105Gangnung37.8128.9EC26
22Pungsan40.8128.2M1206108Seoul37.6127.0P85.5
25Gim-haeg40.7129.2EC23112Incheon37.5126.6WC68.2
28Supung40.5124.9P83114Wonju37.3127.9P148.6
31Cheongjin40.4127.3M1081129Seosan36.8126.5P28.9
35Sinuiju40.1124.4P7131Cheongju36.6127.4P57.2
37Guseong40.0125.3P99133Daejeon36.4127.4P68.9
39Huicheon40.2126.3P155135Chupungyon36.2128.0P244.7
41Hamheung39.9127.6P38136Andong36.6128.7P139.4
46Sinpo40.0128.2EC19138Pohang36.0129.4EC2.3
50Anju39.6125.7P27143Daegu35.9128.6P53.4
52Yangdeog39.2126.7P279146Jeonju35.8127.2P62.4
55Wonsan39.2127.4EC36156Kwangju35.2126.9P72.4
58Pyeongyang39.0125.8P38159Busan35.1129.0SC69.6
60Nampo38.7125.4P47165Mokpo34.8126.4WC38
61Jangjeon38.7128.2EC35169Heuksan34.7125.5WC76.5
65Saliwon38.5125.8P52184Jeju33.5126.5SC20.45
67Singye38.5126.5P100185Gosan33.3126.2SC71.5
68Yongyeon38.2124.9WC5192Jinju35.2128.0SC21.3
Ele. = elevation, Lat. = latitude, Lon. = longitude, M = mountain, P = plain, WC = west coast, EC = east coast, and SC = south coast.
Table 3. Comparison of ET o computed by observed and estimated Rs.
Table 3. Comparison of ET o computed by observed and estimated Rs.
Station CodeStation NameR2RMSE
(MJ/m2)
100Daegwalryeong0.970.28
101Chuncheon0.950.35
105Gangnung0.960.33
108Seoul0.950.39
112Incheon0.950.38
114Wonju0.870.55
129Seosan0.950.34
131Cheongju0.950.40
133Daejeon0.950.33
135Chupungyong0.940.40
136Andong0.950.34
138Pohang0.960.34
143Daegu0.970.29
146Jeonju0.960.33
156Gwangju0.960.30
159Busan0.950.33
165Mokpo0.970.28
169Heuksando0.840.34
184Jeju0.960.29
185Gosan0.950.27
192Jinju0.960.30
Table 4. The seasonal and annual variation of ET o and meteorological variables from 1980 to 2021 in different geographic conditions.
Table 4. The seasonal and annual variation of ET o and meteorological variables from 1980 to 2021 in different geographic conditions.
PeriodMountainPlainWest CoastEast CoastSouth Coast
MEANBMEANBMEANBMEANBMEANB
ETo (unit mm yr−1 for the mean and mm yr−2 for the rate)
Spring227.140.31253.630.32243.250.18214.590.53269.32−0.09
Summer386.740.41419.950.65413.990.60370.700.65425.410.35
Autumn149.300.06175.120.09199.08−0.24165.520.17209.52−0.22
Winter50.340.0063.580.0278.15−0.1964.050.0796.16−0.24
Annual813.500.78912.281.09934.470.36814.851.431000.39−0.21
RH (unit % for the mean and % yr−1 for the rate)
Spring57.59−0.0662.58−0.0270.570.0070.46−0.0666.19−0.02
Summer68.43−0.0573.940.0079.370.0280.35−0.0478.000.00
Autumn59.66−0.0266.060.0472.450.0771.23−0.0169.000.05
Winter52.74−0.0563.98−0.0172.200.0568.95−0.0564.590.00
Annual59.68−0.0466.720.0073.720.0472.81−0.0469.540.01
RS (unit MJm−2d−1for the mean and MJm−2d−1yr−1 for the rate)
Spring18.060.0116.730.0215.750.0114.640.0217.150.01
Summer24.110.0123.420.0223.280.0222.020.0123.560.01
Autumn12.310.0111.970.0112.190.0110.940.0113.270.01
Winter9.070.006.820.006.26−0.015.810.007.510.01
Annual15.930.0114.770.0114.410.0113.400.0115.420.01
Tmax (unit °C for the mean and °C yr−1for the rate)
Spring11.150.0317.210.0516.200.0314.280.0618.720.03
Summer23.470.0428.380.0527.460.0424.860.0528.390.03
Autumn12.290.0418.560.0419.570.0217.520.0421.480.02
Winter−3.860.022.960.034.610.013.340.058.820.04
Annual10.820.0316.840.0417.040.0215.060.0519.420.03
Tmin (unit °C for the mean and °C yr−1 for the rate)
Spring−2.180.015.050.046.840.025.040.047.590.02
Summer12.820.0219.380.0420.400.0418.250.0420.740.03
Autumn−0.440.037.580.0510.650.048.400.0411.190.05
Winter−17.230.01−7.440.04−3.470.03−5.520.03−2.030.04
Annual−1.710.026.210.048.660.036.600.049.430.04
WS (unit ms−1for the mean and ms−1 yr−1 for the rate)
Spring1.640.001.44−0.012.34−0.011.550.001.90−0.02
Summer1.09−0.011.120.001.85−0.011.140.001.69−0.01
Autumn1.26−0.011.05−0.011.91−0.011.360.001.55−0.01
Winter1.56−0.011.26−0.012.35−0.011.530.001.88−0.02
Annual1.39−0.011.21−0.012.11−0.011.390.001.75−0.02
B = standard slope.
Table 5. Identification of the significant increasing, decreasing, and non-significant trends of ET o and meteorological variables on seasonal and annual scales.
Table 5. Identification of the significant increasing, decreasing, and non-significant trends of ET o and meteorological variables on seasonal and annual scales.
PeriodMountainPlainWest CoastEast CoastSouth Coast
Inc.Dec.NoInc.Dec.NoInc.Dec.NoInc.Dec.NoInc.Dec.No
ETo
Spring2038018212404112
Summer3021808302701103
Autumn0055417005125022
Winter0054220005017112
Annual20311114212503112
RH
Spring0143320104134103
Summer0234220104161202
Autumn01410115203134202
Winter1134418104035103
Annual0057217104134202
RS
Spring10414012203701301
Summer3022204410503400
Autumn20312014221503301
Winter0053221005305202
Annual2031907212701400
T max
Spring5002105203800301
Summer5002204203701202
Autumn4012006212404400
Winter1040026005305202
Annual1041808212503400
T min
Spring2031808203602202
Summer4012402401800400
Autumn5002105302602301
Winter0053023005206004
Annual4011808302503400
WS
Spring1131916023026022
Summer0140818032026022
Autumn0141916032026022
Winter1132915005017112
Annual0231817023017022
Inc. = increasing trend, Dec. = decreasing trend, No = no trend.
Table 6. The results of the sensitivity analysis on seasonal and annual scales in different geographic conditions.
Table 6. The results of the sensitivity analysis on seasonal and annual scales in different geographic conditions.
PeriodMountainPlainWest CoastEast CoastSouth Coast
Spring
RH−0.41−0.50−0.98−0.68−1.02
RS0.530.510.470.500.51
T max 0.550.861.040.841.10
T min 0.030.040.050.070.13
WS0.100.120.100.100.09
Summer
RH−0.19−0.24−0.57−0.34−0.66
RS0.740.760.730.760.74
T max 0.720.540.460.170.28
T min 0.220.580.870.911.12
WS0.010.030.010.000.01
Autumn
RH−0.49−0.55−1.15−0.85−1.27
RS0.300.360.380.340.43
T max 0.510.791.030.871.06
T min 0.020.110.240.170.41
WS0.180.190.190.210.15
Winter
RH−0.73−0.99−1.90−1.32−2.03
RS0.000.050.100.020.14
T max −0.310.180.410.290.95
T min −0.17−0.08−0.03−0.06−0.02
WS0.250.250.240.290.24
Annual
RH−0.46−0.57−1.15−0.80−1.24
RS0.390.420.420.410.45
T max 0.370.590.740.540.85
T min 0.020.160.280.280.41
WS0.140.160.140.160.13
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ghafouri-Azar, M.; Lee, S.-I. Meteorological Influences on Reference Evapotranspiration in Different Geographical Regions. Water 2023, 15, 454. https://doi.org/10.3390/w15030454

AMA Style

Ghafouri-Azar M, Lee S-I. Meteorological Influences on Reference Evapotranspiration in Different Geographical Regions. Water. 2023; 15(3):454. https://doi.org/10.3390/w15030454

Chicago/Turabian Style

Ghafouri-Azar, Mona, and Sang-Il Lee. 2023. "Meteorological Influences on Reference Evapotranspiration in Different Geographical Regions" Water 15, no. 3: 454. https://doi.org/10.3390/w15030454

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop