Weather and Climate in Monsoon RegionsView this Special Issue
Spatiotemporal Characteristics of Evapotranspiration Paradox and Impact Factors in China in the Period of 1960–2013
Downward trend of potential evaporation accompanied with upward of air temperature which is denoted as evaporation paradox has been reported in many regions over the past several decades in the world. In this paper, evaporation paradox and key factors attributed to ET0 changes are systematically analyzed based on data from 599 meteorological stations during 1960–2013. Results show that (1) Evaporation paradox exists in all regions in1960–2013 and 1960–1999 except SWRB in 1960–2013 but no evaporation paradox in 2000–2013. (2) Evaporation paradox exists in large areas in spring and summer, the extent and range fall in autumn, and there is no evaporation paradox in winter. (3) The evaporation paradox area accounts for 73.7% of China in 1960–2013 and 91.2% in 1969–1999. (4) Sunshine hours, humidity, wind speed, and maximum temperature appear to be the most important variables which contributed to ET0 change in China.
Climate change characterized by global warming has been the focus of diversified research fields such as water resource, agriculture, ecosystem, and human health. It is widely accepted that global air temperature had been increasing in recent decades, it has risen by about 0.85 (0.65–1.06)°C from 1951 to 2012, and the average rising rate was 0.12 (0.08–0.14)°C (IPCC ). In China, temperature has increased by about 0.5–0.8°C and precipitation has large regional fluctuation but no significant trend in the recent 100 years (Wang et al. ).
Potential evapotranspiration (ET0) is one of the most important components of the hydrological system which refers to “the quantity of water evaporated per unit area, per unit time from an idealized, extensive free water surface under existing atmospheric conditions.” It is an important indicator of atmospheric evaporative demand for estimating terrestrial evaporation and crop water requirements. There have been many discussions on methods of calculating ET0 (Penman , Hargreaves and Samani , Pereira and Pruitt ), the spatial-temporal variations (Irmak et al. , Dinpashoh et al. , Liang et al. , Croitoru et al. ), and its influencing factors (Feng et al. , Liu and Yang , Harmsen et al. , Tang et al. ). Declining trends in both pan evaporation (McVicar et al. ) and potential evaporation (ET0) have been reported to be occurring simultaneously in many regions with increasing trends of air temperature, which has been denoted as the evaporation paradox (Roderick and Farquhar ) and it has been one of the hot issues of hydrological system. Over the past several decades evaporation paradox had been verified in many regions of the world such as the former Soviet Union (Peterson et al. ), the United States (Golubev et al. ), China (Thomas , Ma et al. ), India (Chattopadhyay et al. ), Thailand (Tebakari et al. ), Italy (Moonen et al. ), Romania (Croitoru et al. ), Australia, New Zealand (Roderick and Farquhar ), Canada (Burn and Hesch ). But there existed exception to this rule (Cohen et al. ). In China both at national scale (Yin et al. , Han et al. ) and at regional scale such regions as the Northwest China (Liang et al. ), the YeRB (Wang et al. ), the HaRB (Xing et al. ), the YaRB (Xu et al. ), the Northwest China (Yang et al. ), the Loess Plateau (Li et al. ), and the Tibet Plateau (Liu et al. ) evaporation paradox had been found.
In fact pan evaporation observations mostly ended in 2001 in China; evaporation paradox was concluded based on annual pan evaporation of 1960–2000 (H. Yang and D. Yang ) or potential evaporation from 1960–2010 without no data segment. However, change of temperature and potential evaporation transformed around 2000. In this paper observed meteorological variables are divided into two parts taking year 2000 as the boundary and the objectives of this study are to investigate changes in ET0 and temperature in China since 1960s; to examine the existence of evaporation paradox in different periods and regions; to determine potential key factors attributed to ET0 changes in the whole country as well as different river basins.
2. Data and Methodology
Daily meteorological data were obtained from 754 stations from the China Meteorological Administration (CMA) and National Meteorological Information Center of China (NMIC); 599 stations (Figure 1) of these had complete records of all climatic factors calculating ET0 in time series of 1960–2013. The daily meteorological data included precipitation, relative humidity, sunshine hours, vapor pressure, wind speed, maximum, minimum, and mean air temperature. A few missing data (mainly in 1967, 1968, 1969) were estimated by averaging the value of the other years observed at the same station.
In the data set, the 10 river basins are the first-order basin in China (Figure 1). 56 stations are in the Songhua River basin (SRB), 36 are in the Liao River basin (LRB), 33 are in the Hai River basin (HaRB), 67 are in the Yellow River basin (YeRB), 38 are in the Huai River basin (HuRB), 143 are in the Yangtze River basin (YaRB), 28 are in the southeast rivers basin (SERB), 67 are in the Pearl River basin (PRB), 35 are in the southwest rivers basin (SWRB), and 97 are in the Northwest Rivers Basin (NWRB). In the 599 stations, the Taiwan Island is the one that we could not collect observation data from; therefore, it is excluded from the study region.
2.2.1. Penman-Monteith Method
In this paper, potential evapotranspiration (ET0) was estimated using the Penman-Monteith (PM) method (Allen et al. ); the formula is given as where ET0 is the daily potential evapotranspiration (mm d−1), and the yearly and monthly value of ET0 will be used in this paper; is the net radiation at the top surface (MJ·m−2·d−1); G is the soil heat flux density (MJ·m−2·d−1); is the mean daily air temperature at 2 m height (°C); is daily average wind speed at 2 m height (m·s−1); is the saturation vapor pressure (kPa); is the actual vapor pressure (kPa); is the slope of the vapor pressure curve (kPa°C−1); is the psychrometric constant (kPa°C−1). In the model the radiation term was calculated by experience formula and its accuracy depends on the experience coefficients which were often only effective in particular regions. In this paper, ET0 was calculated by corrected radiation. The correcting net radiation is as follows (Yin et al. ): where is constant of Stefan-Boltzmann (4.903 × 10−9MJ·K−4·m−2·d−1), , is the absolute maximum and minimum temperature (K), is the actual sunshine hours (h), is the duration of possible sunshine (h), and is the Sunny radiation (MJ·m−2). Soil heat flux is small compared with the relative net radiation and in the day time scale.
2.2.2. Trend Analysis Method
The simple linear regression method was used to estimate the trend magnitudes (slope) in ET0 and other climatic variables. The linear equation is where is the simulated value of climatic variables; is the trend which denoted the change trend of climatic variables per decade; and is the time series (Yang et al. ). Meanwhile, the nonparametric Mann-Kendall (M-K) method (Mann , Kendall et al. ) is highly recommended by the World Meteorological Organization for analyzing hydrological series as it did not need any distributional assumption for the data and it was used to detect the significance of the trend.
2.2.3. Stepwise Regression
The basic idea is to introduce the influencing factors into regression equation one by one. Significant test is carried out when introducing one variable into the model, retaining the significant factors and rejecting the insignificant ones until there are no variables introduced into the model and no one rejected. This method can eliminate the variables which contribute little to principal component or those existing linear relations and can overcome the multicollinearity based on guaranteeing the regression effects.
2.2.4. Region Average of Variables
In previous researches, regional value was obtained by using an arithmetic mean method from meteorological station. However, meteorological stations are not distributed evenly but dense in the east and sparse in the west in China. Therefore, it is necessary to assign different weights for different stations when evaluating climate change accurately for different regions. When calculating the average value of an area, the weight of a station is determined by the percentage of the Thiessen polygon in the whole area. Thiessen polygon method was more accurate than simple mean method and less workload grid data set method.
3.1. Observed Changes of Temperature and ET0
In 1960–2013, 98.2% of the 599 stations show upward trend (91.2% of all stations are at 95% significance level). The average daily temperature in China as a whole (Figure 2) rises at the rate of 0.24°C per decade (95% significance level). Corresponding with significant warming trend, the mean national ET0 declines at the rate of −3.9 mm per decade (95% significance level), so there exists evaporation paradox in China as a whole.
3.2. Temporal Trends of Evaporation Paradox
According to Figure 2, the mean annual temperature climaxed in around year 2000 and then decreased slowly, and ET0 reached the lowest value around 1993 and then rose slowly. Taking the change into account comprehensively, this paper took the year 2000 as the dividing line. At the same time in order to compare with the proceeding results of other researchers, we analyzed the characteristics of evaporation paradox in the period of 1960–2013, 1960–1999, and 2000–2013 (Table 1). In 1960–2013 and 1960–1999, the annual temperature increased significantly, while ET0 decreased significantly at the rate −3.9 mm per decade (58.4% of all stations) and −14.78 mm per decade (75.1% of all stations); in 2000–2013, temperature decreased (58.6% of all stations) with the rate being −0.08°C per decade while ET0 increased (56.3%) with the rate being 10.08 mm per decade. In 1960–1999, ET0 of 75.1% stations dropped and temperature of 90% stations rose and the opposite change of temperature and ET0 between 1960–1999 and 2000–2013 made the range of temperature rise and ET0 dropping moderate and evaporation paradox weaken in 1960–2013.
The mean annual temperature in the 10 river basins all rose at 95% significance level and ET0 all showed downward trend except in SERB, YeRB, and SWRB in 1960–2013; all river basins indicated upward in temperature and downward in ET0 in 1960–1999. The maximum downward in ET0 and upward in temperature appeared in NWRB and SRB with values being −27.65 mm per decade and 0.44°C per decade, respectively, in 1960–1999. In 2000–2013, ET0 in YaRB, SERB, SWRB, NWRB, and PRB increased while temperature in SWRB and PRB decreased. In other basins the trend in temperature and ET0 was the same. In this period temperature only in YaRB and SWRB increased, whether the increase was a fluctuation in the whole upward process or the beginning of decrease needs further investigation.
3.3. Seasonal Change of Evaporation Paradox
Figure 3 showed trends of evaporation paradox in 4 seasons.
Spring (March to May): in 1960–2013, evaporation paradox existed in HaRB, LRB, SRB, NWRB, SWRB, PRB, and China as a whole. In 1960–1999, evaporation paradox existed in all regions except SRB; in 2000–2013, there was no evaporation paradox; change of temperature and ET0 was the same in eight river basins and temperature in 5 river basins dropped. The opposite changes of temperature and ET0 in 1960–1999 and 2000–2013 weakened the evaporation paradox of 1960–2013.
Summer (June to August): in 1960–2013, ET0 and temperature were the highest values in a whole year; the slope of ET0 and percent of downward stations were the highest values too. ET0 descended in YaRB, HaRB, HuRB, YeRB, LRB, NWRB, and China as a whole; temperature rose in all regions, except HuRB which was at 99% confidence level. In 1960–1999, evaporation paradox existed in all river basins except SWRB and HuRB; in China as a whole the percent of ET0 downward climaxed 77% and 8 river basins were more than 70%, and so evaporation paradox was the most prominent in all statistical periods. In 2000–2013, the variation of ET0 and temperature was the same except the HuRB.
Autumn (September to November): the slope and range of ET0 decline reduced in autumn comparing with that of spring and summer. In 1960–2013, ET0 in HaRB, LRB, SERB, SRB, and NWRB decreased and temperature increased significantly. The phenomenon existed in the 5 river basins in 1960–1999 too. In 2000–2013, HuRB, LRB, SRB, and PRB showed evaporation paradox.
Winter (December to February next year): change in temperature was the most severe compared with the other seasons. Except in PRB in 1960–1999, temperature in 1960–2013 and 1960–1999 all rose significantly. Opposite to the severe increase in temperature, decrease of ET0 in winter was moderate. In 1960–2013, ET0 only in HuRB, HaRB, LRB, and NWRB declined slightly and other regions showed upward trend. In 1960–1999, ET0 in HaRB, HuRB, PRB, YaRB, NWRB, and China as a whole insignificantly decreased. In 2000–2013, temperature showed biggest fall and only SWRB showed upward trend. In winter ET0 changed the smallest in the four seasons and evaporation paradox was moderate.
3.4. Spatial Distribution of Evaporation Paradox
In 1960–2013 and 1960–1999, the percent of rising stations in temperature exceeded 90%, so stations in which ET0 decreased can be judged as where evaporation paradox existed (Figure 4). The evaporation paradox distribution can be obtained from the interpolation of statistic of ET0. In 1960–2013, 57.6% of the site of annual ET0 decreased in China, the regions where ET0 increased were mainly located in the northeast of the NWRB, northwest of SRB, three rivers sources regions, middle reach in YeRB, northeast and southeast of HuRB, SWRB, coastal area of PRB, middle of YaRB, and so on. Overall coastal areas in the south of 37°N, most of the regions between 30°N and 40°N, 90°E–110°E, northwestern of the northeast China, and southeastern of SWRB were the areas where evaporation paradox does not exist. The evaporation paradox area accounted for 73.7% of the 10 river basins. In 1960–1999, ET0 of 75.1% stations showed downward trend; northwestern of SRB, Ningxia and middle Shaanxi section of the YeRB, three rivers sources regions and northeastern of HuRB, there is no evaporation paradox in such areas. The evaporation paradox area accounted for 91.2%. In 2000–2013, ET0 and temperature of 223 stations change oppositely and 322 stations were the same; 54 stations of statistics of ET0 or temperature were 0.
4.1. Relationship between ET0 and Precipitation
Precipitation and ET0 were two important segments of the hydrologic cycle; ET0 will decline with the increase of precipitation according to the Bouchet assumption. The annual average precipitation is 808 mm in China during 1960–2013; it rose insignificantly and the rising percent was 48.4%. The annual average precipitation was 811 mm in 1960–1999 and the insignificant rising percent was 55.9%; in 2000–2013 the average value was 799.5 mm; change of ET0 and precipitation in China were consistent with the Bouchet hypothesis and Figure 5 showed the relationship of them. In 1960–2013, reverse trend stations were 301 and were located in NWRB, SERB, three rivers sources regions, lower reaches of YaRB, northwestern of SRB, HuRB, and PRB where precipitation increased. Precipitation and ET0 drop sites were mainly located in LRB, HaRB, YeRB, upper reaches of PRB, southeastern of SRB, and the middle reaches of the YaRB. Decreased stations of precipitation in 1960–1999 were more than that of 1960–2013, but its distribution was basically the same. The two factors both indicated upward trend in 2000–2013, whose precipitation rose in 330 stations, unchanged in 45 sites. 352 stations showed opposite change of them which accounted for 58.8% of the whole stations. Regions where precipitation decreased were mainly located in NWRB, SRB, LRB, HaiRB, and the middle reaches of the YeRB.
4.2. Impacts of Meteorological Factors on ET0
Meteorological factors change had profound impacts on ET0; in this paper stepwise regression was used to extract the influencing factors of ET0. Yearly ET0 and 7 meteorological factors such as mean temperature (), maximum temperature (), minimum temperature (), relative humidity (RH), sunshine hours (), average wind speed (), and average water pressure () were firstly normalized in order to remove the impacts of inconsistent units. The entering order of climate variability was showed in Table 2. was the primary contributor which caused ET0 change in China as a whole, SERB, YeRB, and SWRB and the standardized coefficients were 0.75, 0.54, 0.87, and 0.54, respectively. contributed most to ET0 change in YaRB, HuRB, LRB, and PRB with the standardized coefficients of 0.56, 0.66, 0.82, and 0.70. had maximum impact on ET0 in HaRB with standardized coefficients being 0.80. The largest contribution in SRB was RH which was negative with ET0.
Table 2 indicated that , , , and RH were the most important factors influencing ET0. Table 3 showed slope of the listed meteorological elements and ET0 in different statistical period. In 1960–2013, ET0 decreased in all regions expect SERB, YeRB, and SWRB; in NWRB, LRB and HaRB ET0 decreased at 99% level of confidence. In 1960–1999, ET0 decreased in all regions and the decline rate was much more than that of 1960–2013; in 2000–2013, ET0 in 5 river basins was downward trend. which was in positive relationship with ET0 decreased at the rate of −0.11 m s−1 per decade significantly in China as a whole and it was found to be the primary contributor which caused ET0 to decrease in the past 54 years; except PRB in 1960–2013, decreased in all regions significantly. which was in positive relationship with ET0 too had decreased with a significant trend of −0.11 h per decade in 1960–2013 and −0.13 h per decade in 1960–1999 at 99% confidence level and the decline led ET0 to decrease in China as a whole; in the 10 river basins, decreased significantly except SWRB. The range and scope of , decline both diminished in 2000–2013. contributing positive impact on ET0 growth increased in all regions in 1960–1999 and the increase was significant except HuRB in 1960–2013; the most dramatic change occurred in in 2000–2013 among the 4 variables; it changed into decrease from increase in SERB, HaRB, HuRB, YeRB, LRB, SRB, PRB, and China as a whole. RH only increased in YaRB, NWRB, and SWRB in 1960–1999 and LRB and SRB in 2000–2013. The change of ET0 was influenced comprehensively by all these factors; ET0 was in a positive relationship with , , and negative relationship with RH. In China, the decline of , made ET0 reduce and decline RH and ascension of made ET0 ascend. The comprehensive effect of the four elements was the decline of ET0. In 1960–1999, the decline rate of , strengthened corresponding to the weakness in rising and RH decreasing strengthened the decline rate of ET0 to −14.78 mm per decade. In 2000–2013, the trend of meteorological factors changed compared with that of 1960–1999, decline rate of , reduced greatly, switched from increase into decrease, and decline rate of RH increased substantially. The combination caused ET0 switch from decrease to increase.
(1)In 1960–2013, temperature in 98.2% stations of 599 stations increased in China. The decline rate of annual national ET0 was −3.9 mm per decade so evaporation paradox existed. In 1960–1999, ET0 of 75.1% stations was downward and temperature of 90% stations was upward which indicated the most prominent evaporation paradox. In 2000–2013 there was no evaporation paradox. The opposite change of temperature and ET0 in 1960–1999 and 2000–2013 weaken the evaporation paradox in 1960–2013 compared with that of 1960–1999.(2)In 1960–2013, evaporation paradox existed in spring, summer, and autumn in China as a whole; it existed in 6, 6, 5, and 4 river basins in spring, summer, autumn, and winter; the decline rate of ET0 and percent of temperature downward climaxed in summer. In 1960–1999, except SRB in spring, SWRB in summer, HaRB, HuRB, and YeRB in autumn, HuRB, YeRB, LRB, and SRB in winter evaporation paradox exists in other times. In 2000–2013, there was no evaporation paradox.(3)There was no evaporation paradox in the southeastern coastal areas south of 37°N, most of areas in 30°N–40°N, 90°E–110°E, northwestern in SRB, and southeastern of SWRB in China in 1960–2013; the area accounted for 26.3% of the 10 river basins. No evaporation paradox area was only in NWRB, northeastern of SRB, middle reach of YeRB, three river source regions, and northeastern of HuRB which accounted for 8.8% merely.(4)Precipitation in NWRB, SERB, three river source regions, lower reaches of YaRB, northwestern of SRB, northwestern of HuRB, and lower reaches of PRB increased and in such regions ET0 decreased in 1960–2013. Most of stations in which ET0 and precipitation change inversely were located south of 27°N and north of 32°N; the number of the stations was 346 in 1960–1999. In 2000–2013, the stations which precipitation increased were located in north of 32°N and the number of stations in which ET0 and precipitation change inversely was 352.(5), , , RH were the most important variations affecting ET0 change; , , were positive and RH was negative relationship with ET0; , , RH mainly decreased and mainly increased in China and the comprehensive function of them made ET0 decrease in 1960–2013 and 1960–1999 and increase in 2000–2013.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This study was supported by the National Natural Science Foundation of China (Grant no. 51279063); Program for Innovative Research Team (in Science and Technology) in University of Henan Province (15IRTSTHN030); and Program for New Century Excellent Talents in University (Grant no. NCET-13-0794).
IPCC, Climate Change 2013: The Physical Science Basis, Working Group I Contribution to Fifth Assessment Report of Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, UK, 2013.
S. Wang, J. Ye, D. Gong et al., “Construction of mean annual temperature series for the last one hundred years in China,” Quarterly Journal of Applied Meteorology, vol. 9, no. 4, pp. 392–401, 1998.View at: Google Scholar
H. L. Penman, “Natural evaporation from open water, bare soil and grass,” Proceedings of the Royal Society of London. Series A, vol. 193, pp. 120–145, 1948.View at: Publisher Site | Google Scholar
G. H. Hargreaves and Z. A. Samani, “Reference crop evapotranspiration from temperature,” Applied Engineering in Agriculture, vol. 1, no. 2, pp. 96–99, 1985.View at: Google Scholar
A. R. Pereira and W. O. Pruitt, “Adaptation of the Thornthwaite scheme for estimating daily reference evapotranspiration,” Agricultural Water Management, vol. 66, no. 3, pp. 251–257, 2004.View at: Publisher Site | Google Scholar
S. Irmak, I. Kabenge, K. E. Skaggs, and D. Mutiibwa, “Trend and magnitude of changes in climate variables and reference evapotranspiration over 116-yr period in the Platte River Basin, central Nebraska–USA,” Journal of Hydrology, vol. 420-421, pp. 228–244, 2012.View at: Publisher Site | Google Scholar
Y. Dinpashoh, D. Jhajharia, A. Fakheri-Fard, V. P. Singh, and E. Kahya, “Trends in reference crop evapotranspiration over Iran,” Journal of Hydrology, vol. 399, no. 3-4, pp. 422–433, 2011.View at: Publisher Site | Google Scholar
L. Q. Liang, L. J. Li, and Q. L. Liu, “Temporal variation of reference evapotranspiration during 1961–2005 in the Taoer River basin of Northeast China,” Agricultural and Forest Meteorology, vol. 150, no. 2, pp. 298–306, 2010.View at: Publisher Site | Google Scholar
A.-E. Croitoru, A. Piticar, C. S. Dragotǎ, and D. C. Burada, “Recent changes in reference evapotranspiration in Romania,” Global and Planetary Change, vol. 111, pp. 127–132, 2013.View at: Publisher Site | Google Scholar
J. Feng, D. Yan, C. Li, F. Yu, and C. Zhang, “Assessing the impact of climatic factors on potential evapotranspiration in droughts in North China,” Quaternary International, vol. 336, pp. 6–12, 2014.View at: Publisher Site | Google Scholar
Q. Liu and Z. Yang, “Quantitative estimation of the impact of climate change on actual evapotranspiration in the Yellow River Basin, China,” Journal of Hydrology, vol. 395, no. 3-4, pp. 226–234, 2010.View at: Publisher Site | Google Scholar
E. W. Harmsen, N. L. Miller, N. J. Schlegel, and J. E. Gonzalez, “Seasonal climate change impacts on evapotranspiration, precipitation deficit and crop yield in Puerto Rico,” Agricultural Water Management, vol. 96, no. 7, pp. 1085–1095, 2009.View at: Publisher Site | Google Scholar
B. Tang, L. Tong, S. Kang, and L. Zhang, “Impacts of climate variability on reference evapotranspiration over 58 years in the Haihe river basin of north China,” Agricultural Water Management, vol. 98, no. 10, pp. 1660–1670, 2011.View at: Publisher Site | Google Scholar
T. R. McVicar, M. L. Roderick, R. J. Donohue et al., “Global review and synthesis of trends in observed terrestrial near-surface wind speeds: implications for evaporation,” Journal of Hydrology, vol. 416-417, pp. 182–205, 2012.View at: Publisher Site | Google Scholar
M. L. Roderick and G. D. Farquhar, “The cause of decreased pan evaporation over the past 50 years,” Science, vol. 298, no. 5597, pp. 1410–1411, 2002.View at: Publisher Site | Google Scholar
T. C. Peterson, V. S. Golubev, and P. Y. Groisman, “Evaporation losing its strength,” Nature, vol. 377, no. 6551, pp. 687–688, 1995.View at: Google Scholar
V. S. Golubev, J. H. Lawrimore, P. Y. Groisman et al., “Evaporation changes over the contiguous United States and the former USSR: a reassessment,” Geophysical Research Letters, vol. 28, no. 13, pp. 2665–2668, 2001.View at: Publisher Site | Google Scholar
A. Thomas, “Spatial and temporal characteristics of potential evapotranspiration trends over China,” International Journal of Climatology, vol. 20, no. 4, pp. 381–396, 2000.View at: Publisher Site | Google Scholar
X. Ma, M. Zhang, S. Wang, Q. Ma, and S. Pan, “Evaporation paradox in the Yellow River Basin,” Acta Geographica Sinica, vol. 67, no. 5, pp. 645–656, 2012.View at: Google Scholar
T. Tebakari, J. Yoshitani, and C. Suvanpimol, “Time-space trend analysis in pan evaporation over Kingdom of Thailand,” Journal of Hydrologic Engineering, vol. 10, no. 3, pp. 205–215, 2005.View at: Publisher Site | Google Scholar
A. C. Moonen, L. Ercoli, M. Mariotti, and A. Masoni, “Climate change in Italy indicated by agrometeorological indices over 122 years,” Agricultural and Forest Meteorology, vol. 111, no. 1, pp. 13–27, 2002.View at: Publisher Site | Google Scholar
M. L. Roderick and G. D. Farquhar, “Changes in New Zealand pan evaporation since the 1970s,” International Journal of Climatology, vol. 25, no. 15, pp. 2031–2039, 2005.View at: Publisher Site | Google Scholar
D. H. Burn and N. M. Hesch, “Trends in evaporation for the Canadian Prairies,” Journal of Hydrology, vol. 336, no. 1-2, pp. 61–73, 2007.View at: Publisher Site | Google Scholar
S. Cohen, A. Ianetz, and G. Stanhill, “Evaporative climate changes at Bet Dagan, Israel, 1964–1998,” Agricultural and Forest Meteorology, vol. 111, no. 2, pp. 83–91, 2002.View at: Publisher Site | Google Scholar
Y. Yin, S. Wu, G. Chen, and E. Dai, “Attribution analyses of potential evapotranspiration changes in China since the 1960s,” Theoretical and Applied Climatology, vol. 101, no. 1, pp. 19–28, 2010.View at: Publisher Site | Google Scholar
S. Han, D. Xu, and S. Wang, “Decreasing potential evaporation trends in China from 1956 to 2005: accelerated in regions with significant agricultural influence?” Agricultural and Forest Meteorology, vol. 154-155, pp. 44–56, 2012.View at: Publisher Site | Google Scholar
W. Wang, Q. Shao, S. Peng et al., “Reference evapotranspiration change and the causes across the Yellow River Basin during 1957–2008 and their spatial and seasonal differences,” Water Resources Research, vol. 48, no. 5, pp. 27–35, 2012.View at: Publisher Site | Google Scholar
W. Xing, W. Wang, Q. Shao et al., “Changes of reference evapotranspiration in the Haihe River Basin: present observations and future projection from climatic variables through multi-model ensemble,” Global and Planetary Change, vol. 115, pp. 1–15, 2014.View at: Publisher Site | Google Scholar
C.-Y. Xu, L. Gong, T. Jiang, D. Chen, and V. P. Singh, “Analysis of spatial distribution and temporal trend of reference evapotranspiration and pan evaporation in Changjiang (Yangtze River) catchment,” Journal of Hydrology, vol. 327, no. 1-2, pp. 81–93, 2006.View at: Publisher Site | Google Scholar
Y. Yang, Z. Feng, H. Q. Huang, and Y. Lin, “Climate-induced changes in crop water balance during 1960–2001 in Northwest China,” Agriculture, Ecosystems & Environment, vol. 127, no. 1-2, pp. 107–118, 2008.View at: Publisher Site | Google Scholar
Z. Li, F.-L. Zheng, and W.-Z. Liu, “Spatiotemporal characteristics of reference evapotranspiration during 1961–2009 and its projected changes during 2011–2099 on the Loess Plateau of China,” Agricultural and Forest Meteorology, vol. 154-155, pp. 147–155, 2012.View at: Publisher Site | Google Scholar
X. Liu, H. Zheng, M. Zhang, and C. Liu, “Identification of dominant climate factor for pan evaporation trend in the Tibetan Plateau,” Journal of Geographical Sciences, vol. 21, no. 4, pp. 594–608, 2011.View at: Publisher Site | Google Scholar
H. Yang and D. Yang, “Climatic factors influencing changing pan evaporation across China from 1961 to 2001,” Journal of Hydrology, vol. 414-415, pp. 184–193, 2012.View at: Publisher Site | Google Scholar
R. G. Allen, L. S. Pereira, D. Raes, and M. Smith, “Crop evapotranspiration guidelines for computing crop water requirements,” FAO Irrigation and Drainage Paper 56, FAO, Rome, Italy, 1998.View at: Google Scholar
Y. Yin, S. Wu, D. Zheng, and Q. Yang, “Radiation calibration of FAO56 Penman-Monteith model to estimate reference crop evapotranspiration in China,” Agricultural Water Management, vol. 95, no. 1, pp. 77–84, 2008.View at: Publisher Site | Google Scholar
Y. Yang, X. Zhou, and B. Wu, “Urban expansion prediction for Zhangzhou city based on GIS and spatiotemporal logistic regression model,” Journal of Geo-Information Science, vol. 13, no. 3, pp. 374–382, 2011.View at: Publisher Site | Google Scholar
H. B. Mann, “Nonparametric tests against trend,” Econometrica, vol. 13, pp. 245–259, 1945.View at: Publisher Site | Google Scholar | MathSciNet
M. G. Kendall, Rank Correlation Methods, Griffin, London, UK, 1975.