Estimation of Actual Evapotranspiration Distribution in the Huaihe River Upstream Basin Based on the Generalized Complementary Principle

College of Atmospheric Science, Nanjing University of Information Science and Technology, Nanjing 210044, China Inner Mongolia Meteorological Information Center, Huhehot 010000, China College of Applied Meteorology, Nanjing University of Information Science and Technology, Nanjing 210044, China Climate Centre of Jiangsu Province, Jiangsu Provincial Meteorological Bureau, Nanjing 210009, China Nanjing Broadcasting System, Nanjing 210001, China Jinjiang Meteorological Bureau, Quanzhou 362200, China College of International Students, Nanjing University of Information Science and Technology, Nanjing 210044, China


Introduction
Evapotranspiration from the land surface plays a key role in maintaining the balance of land surface water-lakes-reservoirs and the energy balance of the earth's surface [1][2][3][4]. Evapotranspiration is a very important step in the water cycle. During the evaporation process, water is transported from the surface to the atmosphere. e water circulation, energy balance, and carbon cycle are closely linked. Many studies have shown that at least 70% of surface precipitation will return to the atmosphere through evapotranspiration, and the value that rises to over 90% in drought areas. e accurate estimation of basin evaporation has very important applications for water resources assessment, vegetation drought monitoring, and ecological water use. e complementary principle was proposed by Bouchet [1]. is model, based on the principle of complementary correlation, is widely promoted and used because it only needs conventional meteorological data as inputs to calculate evapotranspiration rather than hydrological and soil data which are hard to obtain.
Actual evapotranspiration is defined as a combined process of both evaporation from soil and plant surfaces and transpiration through plant canopies [5]. e model for estimating actual evapotranspiration based on the complementary principle mainly includes the following: the advection-aridity (AA) [6], the complementary relationship areal evapotranspiration (CRAE) [7], and the Granger and Gray (GG) [8] model. Jian et al. [9] applied the AA model to the Tarim River to estimate actual evapotranspiration and verified the complementary correlation theory, which proved that the calculation results of the AA model in this area are reasonable and credible. Han et al. [10] compared the applicability of the AA, GG, and PM-Katerji models, and the results showed that the AA model had good estimation ability in generally wet conditions, whereas the estimated actual evapotranspiration in dry regions was relatively low and the estimation results under very wet conditions were relatively high. e Granger model was approximately equivalent to the AA model under neither dry nor wet conditions, and the limits of variation of the evapotranspiration ratio were practical over a wide range. e PM-Katerji model had poor estimation results when the variation of the evapotranspiration ratio was larger. Han et al. [11] first proposed boundary conditions in the formulation of a complementary relationship model. Based on the model by Han et al. [11], Brutsaert [2] set physical constraints for the ratio of actual evapotranspiration to possible evapotranspiration, as well as the ratio of wet evapotranspiration to evapotranspiration, and proposed a generalized complementary correlation principle. Crago et al. [12] revised the generalized complementary principle of Brutsaert [2] and concluded that when the actual evapotranspiration was close to 0, the ratio of wet evapotranspiration to possible evapotranspiration is not 0, but between 0 and 1. e revised model was tested using depth and precipitation data from the United States. Szilagyi et al. [13] took into account the limitations of water balance data and then derived a CR model without calibration that was a unique estimation Priestley-Taylor α parameter model with minimum data requirements and no dependencies for the calculation of the parameter α. e works of Crago et al. [12] and Szilagyi et al. [13] have been modified and updated based on Brutsaert [2], and new uncertainties have been discussed. Studies show that models based on the principle of complementary correlation are less effective in estimating actual evapotranspiration in humid and arid environments, and actual evapotranspiration cannot be calculated in areas where drought is severe or in very humid areas. However, the principle of generalized complementary correlation does not have this limitation.
Water resources per unit area and per capita in China are below the world average, and freshwater resources per capita in China are only a quarter of the world's per capita freshwater resources. e distribution of water resources in China is extremely uneven, and the evapotranspiration process is an important process of the water cycle. e accurate estimation of evapotranspiration has a very important value for water resources assessment and vegetation drought monitoring. In China, scarcely any research has been done based on the principle of generalized complementary correlation. In this study, we selected the basin above Wangjiaba in the upper Huaihe River as the study area and calculated the annual evapotranspiration over 40 years according to the generalized complementary principle of Brutsaert [2] using spatial simulation results of meteorological elements; we then verified the simulation accuracy of the actual evapotranspiration in the basin. Regarding the determination of model parameters and basic data inputs, the simulation results for evapotranspiration in different regions were different. Research on the coupling of spatial simulation methods of meteorological elements with an actual evapotranspiration model was also an important purpose of this study.

Data.
e Huaihe River Basin is in eastern China between the Yangtze River and the Yellow River. e Huaihe River is divided into three parts: the upper, middle, and lower reaches. e upper reach is Honghekou and above and is 360 km long between longitudes 113°27′ and 115°62′E and latitudes 31°47′ and 33°53′N. e watershed control area is approximately 31,000 km 2 .
e Huaihe River Basin is in the transitional climate zone between North China and South China and is characterized by a warm temperate and semihumid monsoon climate. Additionally, the winter and spring are characterized by droughts with rare rainfall, and the summer and autumn are usually hot with frequent rain. e rainfall in the basin is sufficient, and the annual average rainfall is approximately 920 mm. e spatial distribution of rainfall roughly decreases from the south to the north, with the mountainous areas receiving more rain than the plains areas. e annual average temperature is generally 11-16°C, and the monthly average temperature of July is the highest, usually approximately 25°C. e monthly average temperature of January is the lowest, usually approximately 0°C [14]. In this study, we chose the upper reach of the Huaihe River Basin as a research area not only because we have sufficient historical meteorological data for this area but also because the climate characteristics of this area represent the transitional climate between China's north and south. Figure 1 shows the distribution of the water system, meteorological stations, and hydrological stations in the upper reach of the Huaihe River Basin and its surrounding areas. e meteorological data between 1956 and 2000, which were obtained from the National Meteorological Information Center, include the daily average temperature, daily minimum temperature, daily maximum temperature, daily range of temperature, percentage sunshine, daily average wind speed, and daily precipitation data. e level-2 water resources data were obtained from the Water Resources, Hydropower Planning and Design Institute of the Ministry of Water Resources. e level-2 river network data and the digital elevation model (DEM) data source were derived from the National Geomatics Center of China. We estimated the distributed results of the monthly albedo from 1961 to 2000 using the data from NOAA/AVHRR satellite channels 1 and 2 from 1961 to 2000 and the estimation formula for the surface albedo.

Estimation of Actual Evapotranspiration.
In this study, the actual evapotranspiration is calculated based on the generalized complementary principle proposed by Brutsaert. He introduced nondimensional E TW and E TA via dividing them by E TP , x E TW /E TP , and y E TA /E TP , where E TA represents the actual evapotranspiration, E TP represents the potential evapotranspiration, and E TW represents the wet potential evapotranspiration under the "potential" conditions that the evapotranspiration water is not a limiting factor with the same solar radiation input and the same meso-alpha scale pressure gradient pattern [2]. e polynomial implementation can be written by expressing the relationship of x and y as follows: where a i is a coe cient and n is the degree of the polynomial.
To develop a signi cant equation, Brutsaert [2] imposed four physically determined boundary conditions that he used with (1) to develop a fourth-order polynomial solution. Speci cally, he made the following assumptions: Under very moist conditions, (i) shows that y 1 as x 1 and a 0 0, a 1 1, a 2 −2, and a 3 1; this yields the main result: or in the equation with the original variables: where E TW is given by the Priestley and Taylor [15] where α is the Priestley-Taylor parameter, related to the underlying surface conditions (here, we assume 0.72 ≤ α ≤ 1.57); R n is the net radiation (mm/d); G is the heat ux into the ground (mm/d); and Δ is the slope of the saturation vapor pressure curve (hPa). Additionally, where c is the psychrometric constant and E a is the drying power of the air in mm/d. e actual evapotranspiration can be expressed by (4) and (5) as follows: Since the formula [17] is based on the assumption that there is no advection, the inhomogeneity of the underlying surface of the basin will lead to the occurrence of advection. ere can be no full advection in the actual environment, and α thus has great variability [18] and needs to be calibrated. Based on previous studies [19], the actual evapotranspiration calculated by the water balance equation (7) is calibrated with steps of 0.01 in the model's α, which increases from 0.72 to 1.57 [20], depending on di erent surface conditions. e actual evapotranspiration is calculated with (6). Figure 2 shows a ow chart of the spatial distribution result of actual evapotranspiration based on the generalized complementary principle. e actual evapotranspiration is calculated from the water balance method as follows [17]: Advances in Meteorology 3 where P is the average precipitation in the basin; R is the runoff depth in the basin; and ΔW is the changing value of soil water storage in the basin, averaged for many years; the changing value of soil water storage in the basin is approximately zero.

Calculation of Model Components
(1) Slope of the Saturation Vapor Pressure Curve. e equation for calculating the slope of the saturation vapor pressure curve is given as follows [21]: where T is the temperature in°C, e * is the saturated vapor pressure (hPa), and a and b are coefficients. a � 17.23 and b � 237.3 at T ≥ 0°C, and a � 21.88 and b � 265.5°C at T < 0°C.
(2) Psychrometric Constant. e expression of psychrometric constants is defined as follows [15]: where T is the temperature in°C and p is the vapor pressure (hPa).
(3) Net Radiation. e expression of net radiation under the complex terrain is defined [4] as follows: where Q αβ is the total radiation under the complex terrain (MJ/(m 2 ·d)), α s is the surface albedo, and F αβ is the net surface longwave radiation (NSLR) (MJ/(m 2 ·d)). e unit needs to be unified in the calculation, with 1 mm/d � 0.408 MJ/(m 2 ·d) at T ≥ 0°C and 1 mm/d � 0.355 MJ/(m 2 ·d) at T < 0°C.
(1) Calculation of total radiation under the complex terrain e method of calculating the total solar radiation received by the surface under the complex terrain is as follows: where Q bαβ is the direct solar radiation energy under the complex terrain (MJ/(m 2 ·d)), Q dαβ is the solar scattered radiation energy under the complex terrain (MJ/(m 2 ·d)), and Q rαβ is the solar shortwave radiation reflected from the topography around the complex terrain (MJ/(m 2 ·d)). For the specific calculation formula of Q bαβ , Q dαβ , and Q rαβ , see the study by Zeng [22].
(2) Estimation of the surface albedo by remote sensing e surface albedo can be calculated by climatological calculation or satellite remote sensing inversion. e surface albedo of this study is calculated by the method given by

Advances in Meteorology
Valiente [23], by using the NOAA/AVHRR channel 1 and 2 data. e surface albedo inversion is calculated using monthly measurements over the 1961-2000 period [23]. e surface albedo is mainly used to calculate the simulation parameters of the net radiation in the model. Previous studies [24,25] carried out the calculations using some fixed values added to the net radiation model. For example, the albedo of the water surface is 7, the albedo of the grassland is 15, and the desert albedo is 28. e addition of the surface albedo using remote sensing makes the spatial distribution of net radiation finer. e spatial resolution is 8 km: where α s is the average monthly surface albedo and P CH1 and P CH2 are the results of the AVHRR channels 1 and 2, respectively.
(3) Calculation of NSLR under the complex terrain e NSLR under the complex terrain is calculated as follows [26]: where ε s is the surface emissivity, generally taken as 0.95; σ is the Stefan-Boltzmann constant (4.9 × 10 −9 MJ/(m 2 ·d·K)); T is the temperature (K); a 0 , a 1 , and b 0 are the coefficients; e is the actual vapor pressure (hPa); and S is the percentage sunshine. e simulation parameters of NSLR are shown in Table 1.
(4) Soil Heat Flux. e equation for estimating the soil heat flux G is defined as follows [15]: where G is the soil heat flux (MJ/(m 2 ·d)), T i+1 is the monthly average temperature of the next month (°C), and T i−1 is the monthly average temperature of the previous month (°C). e units need to be unified when calculating G, with 1 mm/d � 0.408 MJ/(m 2 ·d) at T ≥ 0°C and 1 mm/d � 0.355MJ/(m 2 ·d) at T < 0°C.  [15]. e expression for estimating the drying power E a is as follows [4,15,16]: where u 2 is 2-meter average wind speed in m/s. In this study, meteorological data generally only provide the wind speed at 10 meters from the surface, and the wind speed thus needs to be changed. e conversion formula is defined as follows [21]:

Verification.
e first evaluation of water resources in China occurred from 1956 to 1979, and a large number of reservoirs were built after 1980. ese greatly affected the water balance and are not conducive to estimating the actual evapotranspiration in the river basin. erefore, the data from 1956 to 1976 were chosen for estimation. e α in the upper Huaihe River Basin is 1.25. When using multiyear average data, another method of calculating the actual evapotranspiration from (7) can be used to determine that the changing value of soil water storage in the basin is approximately zero. e methods of this study can determine the actual evapotranspiration for each year. To verify the correctness of the actual evapotranspiration using these methods, we used the water balance equation (7) in the verification component of this study to calculate the actual evapotranspiration method for comparison. Due to the limited data and the necessity for verification, we chose the 10-year moving average method, which calculated the actual evapotranspiration using the annual average data from 1967 to 1977, 1968 to 1978, and 1969 to 1979. We chose the 10-year average to calculate the actual evapotranspiration by referring to the Priestley and Taylor [17] method, which mainly considered the stability of the calculation results. Table 2 shows the verification results. e average absolute error is 16.64 mm, the relative error is within ±5%, the average bias is 9.57 mm, the root-meansquared error (RMSE) is 20.13 mm, and the average relative error is 2.25%. e R-square of validation is 0.986 (Figure 3). Figure 3 is a square graph of water balance and generalized complementary simulation and shows their bias by the square axis and 45-degree line.  Advances in Meteorology

Spatial Estimation Results of Actual Evapotranspiration.
Based on the simulation method in this study, we obtained the spatial distribution of annual average actual evapotranspiration in the upper reach of the Huaihe River Basin from 1961 to 2000 (Figure 4). e annual average actual evapotranspiration in the Huaihe River Basin is high in the south and low in the north. e high evapotranspiration values are concentrated in the southwest of the basin, with the highest values within the vicinity of the Dabie Mountains and the low values in the northwest watershed. e average annual evapotranspiration of the whole basin over 40 years is 665.11 mm, and that of the water balance data is 677.38 mm, which is consistent by an order of magnitude. e absolute error is 12.27 mm, and the relative deviation is 1.81%. ese veri cation results illustrate the applicability of the generalized complementary correlation principle in the region. Figure 4 shows the annual average change curve of actual evapotranspiration in the upper reach of the Huaihe River Basin from 1961 to 2000. Figure 5 shows that the actual evapotranspiration is a signi cant downward trend over time, with a rate of decline of 24.1 mm/10 years.

Analysis of In uencing Factors.
Pearson correlations between actual evapotranspiration and the percentage sunshine, maximum temperature, minimum temperature, average temperature, daily temperature, 2-meter average wind speed, relative humidity, daily range of temperature, actual vapor pressure, and other variables are shown in Table 3. As seen from the table, the percentage sunshine, maximum temperature, daily range of temperature, 2-meter average wind speed, and actual evapotranspiration were positively correlated, and they all passed the 0.05 signi cance test. erefore, we mainly analyzed the in uence of the percentage sunshine, maximum temperature, daily range of temperature, and 2-meter average wind speed on the actual evapotranspiration. e t-test was carried out on the trend of the percentage sunshine, maximum temperature, daily range of temperature, and 2-meter average wind speed, and the test passed the signi cance test at the 99% con dence level. As seen from Table 3 and Figure 6, the decrease in the percentage sunshine (a), the decrease in the daily range of temperature (c), and the decrease in the 2-meter average wind speed (d) are the main reasons for the decrease of actual evapotranspiration. e increase in the maximum temperature (b) also slows down the downward trend of the actual evapotranspiration.

Discussion
In this study, we estimated actual evapotranspiration based on the generalized complementary principle. e model components (drying power, soil heat ux, net radiation, constant of wet and dry tables, etc.) and surface albedo all  affect the accuracy of the actual evapotranspiration. Using the Priestley-Taylor formula as a possible evapotranspiration formula, the parameter α is rated based on the influence of underlying surface factors. However, the α value changes as the range of the underlying surface changes. e smaller the calibration area, the more accurate the α parameter and the more reliable the possible evapotranspiration. e most accurate calibration method for the α parameter is a comparison of the estimation results with the measured data. e method of obtaining the actual evapotranspiration through the lysimeter is the only method that is closest to the true value; however, the lysimeter data for China are still scarce at large scales [27] and further development and study are needed. To ascertain whether the generalized complementary theory is applicable to other basins in China, further research and verification are needed.
Since the formula is based on the assumption that there is no advection, the inhomogeneity of the underlying surface of the basin will lead to the occurrence of advection. ere can be no full advection in the actual environment; thus, the α parameter has great variability [18], and estimating it accurately requires calibration. Eagleson [18] provided values from 0.72 to 1.57 that depended on different surface conditions. Relying on instruments for observation is the most reliable way to obtain continuous evapotranspiration data. As the Chinese economy develops, more evapotranspiration instruments will be deployed. e verification part of this study could only be carried out using the water  1980 Years (d) Figure 6: Variation trend graphs of maximum temperature, daily range of temperature, 2-meter average wind speed, and percentage sunshine with actual evapotranspiration.
balance method, and the data were limited. In the near future, evapotranspiration measured by instruments will be used to verify the results and will make the results more reliable. e sensitivity coefficient is a dimensionless index S calculated by the radio of the change in output to input on the condition that the other variables remain constant. e sensitivity of dependent E TA to a particular independent X (slope of the saturation vapor pressure curve, psychrometric constant, net radiation, soil heat flux, or drying power) can be calculated from the derivative of E TA with X, zE TA /zX [28,29]. e first-order Taylor series approximation was applied to calculate S [30][31][32]. To more clearly assess the sensitivity, S was ranked into four classes [31] (Table 4). e higher the sensitivity coefficient value, the more the impact of this variable on evapotranspiration. e sign of S determines if the model would react codirectionally to the input parameter change, that is, an increase (a decrease) of the parameter would lead to an increase (a decrease) of the output variable. Table 5 shows that the sensitivity coefficient of the variables is different and R n > Δ > c > G > E a . e sensitivity levels of the net radiation, slope of the saturation vapor pressure curve, and psychrometric constant are high, which indicates high sensitivity. e sensitivity impact of the slope of the saturation vapor pressure curve and net radiation is codirectional.
Many researchers [12,13] have attempted to correct Brutsaert et al.'s [2] generalized complementary correlation principle. Crago et al. [12] rescaled the x-axis by calculation and proposed a new version of the general CR. ey also believed that the new version of the general CR better predicted the overall measured actual evapotranspiration than did Brutsaert's [2] equation and maintained considerable continuity with the reasoning behind the general CR derived by Brutsaert et al. [2]. Szilagyi et al. [13] introduced a novel approach to calculate the value of the Priestley-Taylor parameter and made a new version of the CR equation with free calibration. e objective of this study was to apply the generalized complementary principle to evapotranspiration estimation at the watershed level in China and to verify its viability. New methods [12,13] have been introduced and have brought new uncertainties, which require further discussion. In our future research, we will focus on different updated CR models based on Brutsaert et al.'s [2] method and choose the best performance for one or some Chinese basins.

Conclusions
e theoretical model of the generalized complementary correlation principle was used to calculate the actual evapotranspiration using the spatial simulation results of meteorological elements in the upper reach of the Huaihe River. Estimation was successful, and the spatial distribution of the average actual evapotranspiration was well characterized by the distributed simulation. e estimated average annual evapotranspiration absolute error was 16.64 mm, and the average relative error was 2.25%. Pearson's correlation coefficient analysis showed that the reduction in the percentage sunshine and the decreases in the daily range of temperature and 2-meter average wind speed were the main reasons for the decrease in actual evapotranspiration. erefore, for vegetation and drought monitoring and ecological modeling, estimating evapotranspiration and evaluating water resources in the basin are very important.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this article.