Changes in Reference Evapotranspiration and Its Contributing Factors in Jiangsu , a Major Economic and Agricultural Province of Eastern China

Reference evapotranspiration (ETref) is a key parameter of hydro-meteorological studies as well as water resource planning. In this study, we adopted the Penman–Monteith FAO 56 model to estimate ETref and through the differential equation and detrending method to determine sensitivities and the contributions of four meteorological parameters to ETref based on daily weather data from 60 stations of Jiangsu province during 1961–2015. Results reveal that ETref and its trends in the three sub-regions of the Jiangsu province had a significant spatial heterogeneity. A significant decreasing tendency of ETref (p < 0.001) was observed in the Huaibei region, while a slightly increasing tendency was identified in the Jianghuai and Sunan regions. These changes of ETref were caused by a significant increasing trend in air temperature (TA) and significant decreasing trends in wind speed (WS), sunshine duration (SD) as well as a non-significant change trend in actual vapor pressure (VP). However, the VP was the meteorological parameter to which ETref was most sensitive, whereas ETref was more sensitive to TA and SD in the summer but less so in the winter; the least sensitive factor, WS, had the opposite trend. Across the whole region, WS contributed most to ETref, followed by SD, while the positive contribution of TA to ETref could not offset the negative contributions of WS and SD. Although the effect of VP on changes in ETref is small, it could not be ignored, especially in the winter. The reverse relationship between increasing TA and decreasing ETref, namely the “evaporation paradox,” occurred in Jiangsu province. Thus, the outcomes of this study will contribute to thorough insight into the response to changes in ETref to the provincial water planning and management in eastern China.


Introduction
Reference evapotranspiration (ET ref ) plays a vital role in climatological and hydrological researches, where spatial and temporal variations are a key indicator of changing climate and water resource allocation [1,2].The study of spatial and temporal changes of ET ref has drawn much attention Water 2017, 9, 486 2 of 20 in recent decades to regional-scale hydro-climatic change because of its nexus in the water resource management and terrestrial ecosystems [3][4][5][6].Additionally, ET ref not only plays a role in estimating agricultural water requirements, but also in irrigation scheduling and management [7].Therefore, knowledge into the changes of ET ref and its precise estimation responses to climate change conditions has become more important ever [8,9], particularly in the context of climate warming starting in the 1970s [10].
Meteorological factors and ET ref have changed drastically in the recent decade, due to the anthropogenic activities and land use changes [11][12][13].In the last century, mean air temperature has risen by 0.6 • C globally on average, and mean temperature has also increased by 0.5-0.8• C in China [10].Despite an increasing trend in the air temperature, a significant number of works have found decreasing pan evaporation and ET ref in many parts of the world, including China [14][15][16], Australia [17], New Zealand [18], Canada [3], India [19], Thailand [20], and Mexico [21].This has led to the coinage of the term "evaporation paradox" introduced by Brutsaert and Parlange [22].However, the reason for its occurrence is still debatable [23].Hobbins et al. [24] showed that the pan evaporation paradox is an indicator of the close linkage between actual evaporation and potential evaporation.Meanwhile, other researchers have also reported some strong proofs in accordance with the opinion of Hobbins et al. [24] that decreasing rates of ET ref are responsible not only for the rise in air temperature but also for other meteorological factors such as the reduction in wind speed, solar radiation, and relative humidity [25,26].
The Penman-Monteith FAO 56 model suggested by the Food and Agriculture Organization (FAO) has been regarded as a robust and universally accepted method for calculating ET ref [27,28].From the Penman-Monteith method, it is evident that most climatic parameters are affected by ET ref , including air temperature, wind speed, relative humidity, and sunshine duration, as well as net radiation.In addition, various attribution analysis techniques have been applied to determine the contribution of climatic factors on ET ref changes.The most widely used technique is a sensitivity coefficient analysis, which has been employed to measure the changes of ET ref in response to an observed change of climatic parameters [23,29].However, this method does not investigate the evaporative process, but only shows the relative contributions of climatic factors on ET ref changes.Xu et al. [8] have developed a technique for analyzing the relative change of evaporation demand through statistical and mathematical modeling.This mathematical approach has been commonly used to calculate the relative changes in ET ref analysis.For instance, Xu et al. [8], Li et al. [30], Liu et al. [31], and Huo et al. [32] used this technique for an attribution analysis of change in ET ref in different regions of China.Similarly, in the present study, we successfully apply this technique to Jiangsu province in eastern China.
However, Liu et al. [15] showed a reverse relationship between rising temperature and reducing evaporation in southwest China, which appeared to represent a paradoxical situation and decreases the contribution of wind speed to the ET ref changes to compensate for the increased contribution of mean air temperature.Wang et al. [23] found that the contributions of solar radiation, wind speed, and relative humidity are the major driving factors for the decrease in ET ref for the three-river source region of China.Hence, further evaluation is essential to determine the relationship between climatic factors that causes changes in ET ref in Jiangsu province.Attribution analysis of changes in ET ref will involve a detrending method that will be discussed in detail in this study.
Evapotranspiration (ET), an important climatic parameter that mainly controls the water balance process, has a great impact on hydro-meteorology processes under a changing climate [33,34].In water-limited regions, changes in ET ref will be the result of changes in ET (and so the water balance), whereas in energy-limited regions they drive changes in ET.The change of ET ref is a key component of hydro-meteorological studies in Jiangsu province.Several regional studies have documented the ET ref in the Yangtze River basin [8,35,36].Furthermore, a few studies have assessed the spatial distribution and temporal changes in the ET ref trends and its magnitude and spatial patterns of ET ref trends in eastern China.However, so far, changes in ET ref and its contributing factors in Jiangsu province have yet to be investigated thoroughly, and the province is located in the Yangtze River Water 2017, 9, 486 3 of 20 Delta, which contributes to the highest GDP in economy and is an emerging developed region.As an important economic province, Jiangsu also ranks at the top in terms of agricultural systems in China.This province has distinctive monsoon climatic characteristics, which makes it inherently sensitive to climate change [37].It is therefore important to investigate changes in ET ref trends and its contributing factors, which can provide a theoretical basis and scientific guidance for the management and allocation of water resources in Jiangsu province.The goals of this study are to explore the spatial-temporal trends in ET ref in the Jiangsu province based on weather data from 60 stations during 1961-2015, and to estimate sensitivities and contributing factors to changes in ET ref through a detrending method applied to the Penman-Monteith FAO 56 model.It is anticipated that the outcomes of this study will provide better guidance for agricultural production and economic development in this vital region.

Study Area and Data Source
Jiangsu province is located on the eastern coast of China and at the lower reaches of the Yangtze River and Huai River, which across its south and north, respectively, cover an area of 1.03 × 10 5 km 2 .The terrain primarily consists of plains, including the southern Jiangsu plain (south of Yangtze River), the middle Jianghuai plain (between Yangtze River and Huai River) and the northern Huang-Huai plain (north of Huai River), which comprises about 70% of the provincial area (Figure 1).The study area belongs to the East Asian monsoon zone and the transition between a subtropical zone and a warm temperate zone, which has distinct monsoon climate characteristics.Generally, the large area of the southern Huai River belongs to subtropical humid monsoon climate, the northern area belongs to temperate semi-humid monsoon climate [38].There is a clear discrepancy between four seasons, the sunshine duration is 2000~2600 h, mean temperature is 13~16 • C, the frost-free period is about 200~240 days, and mean precipitation is 800 mm~1200 mm on an annual scale, with more than 60% of rainfall occurring during the summer season [39].
Water 2017, 9, 486 3 of 21 the Yangtze River Delta, which contributes to the highest GDP in economy and is an emerging developed region.As an important economic province, Jiangsu also ranks at the top in terms of agricultural systems in China.This province has distinctive monsoon climatic characteristics, which makes it inherently sensitive to climate change [37].It is therefore important to investigate changes in ETref trends and its contributing factors, which can provide a theoretical basis and scientific guidance for the management and allocation of water resources in Jiangsu province.The goals of this study are to explore the spatial-temporal trends in ETref in the Jiangsu province based on weather data from 60 stations during 1961-2015, and to estimate sensitivities and contributing factors to changes in ETref through a detrending method applied to the Penman-Monteith FAO 56 model.It is anticipated that the outcomes of this study will provide better guidance for agricultural production and economic development in this vital region.

Study Area and Data Source
Jiangsu province is located on the eastern coast of China and at the lower reaches of the Yangtze River and Huai River, which across its south and north, respectively, cover an area of 1.03 × 10 5 km 2 .The terrain primarily consists of plains, including the southern Jiangsu plain (south of Yangtze River), the middle Jianghuai plain (between Yangtze River and Huai River) and the northern Huang-Huai plain (north of Huai River), which comprises about 70% of the provincial area (Figure 1).The study area belongs to the East Asian monsoon zone and the transition between a subtropical zone and a warm temperate zone, which has distinct monsoon climate characteristics.Generally, the large area of the southern Huai River belongs to subtropical humid monsoon climate, the northern area belongs to temperate semi-humid monsoon climate [38].There is a clear discrepancy between four seasons, the sunshine duration is 2000~2600 h, mean temperature is 13~16 °C, the frost-free period is about 200~240 days, and mean precipitation is 800 mm~1200 mm on an annual scale, with more than 60% of rainfall occurring during the summer season [39].The locations of meteorological stations and the elevation are shown in Figure 1.The Jiangsu Province is divided into Huaibei, Jianghuai, and Sunan regions on the basis of two major rivers (Huai River and Yangtze River) and their regional climate characteristics.The basic climatic characteristics of the three sub-regions and the whole Jiangsu province are shown in Table 1.As seen from Table 1, the annual mean temperature (15.09 °C), the mean relative humidity (76.15%), wind speed (2.79 m•s −1 ), sunshine duration (2116.07h), and annual reference evapotranspiration (892.24mm) have been observed for the last 55 years  in the whole of Jiangsu province.The locations of meteorological stations and the elevation are shown in Figure 1.The Jiangsu Province is divided into Huaibei, Jianghuai, and Sunan regions on the basis of two major rivers (Huai River and Yangtze River) and their regional climate characteristics.The basic climatic characteristics of the three sub-regions and the whole Jiangsu province are shown in Table 1.As seen from Table 1, the annual mean temperature (15.09 • C), the mean relative humidity (76.15%), wind speed (2.79 m•s −1 ), sunshine duration (2116.07h), and annual reference evapotranspiration (892.24mm) have been observed for the last 55 years  in the whole of Jiangsu province.Ground-based daily meteorological data from 60 meteorological stations provided by the Jiangsu Meteorological Information Centre (JMIC) of the China Meteorological Administration (CMA) including the daily mean temperature (T, • C), maximum temperature (T max , • C) and minimum temperature (T min , • C), precipitation (P, mm), sunshine duration (SD, h), wind speed (WS, m•s −1 ) at 10 m height, and relative humidity (RH, %) during 1961-2015.The station dataset is checked by the JMIC.In addition, this study has performed routine quality checks and error correction for data following the techniques defined by Peterson et al. [40].The four seasons are spring (March-May), summer (June-August), autumn (September-November), and winter (December-February of the next year).

Penman-Monteith FAO 56 Model for ET ref Estimation
The Penman-Monteith FAO 56 model is a broadly accepted method for estimating ET ref suggested by the Food and Agriculture Organization (FAO) of the United Nations [27].The method employed in this study because it is physically based and integrates physiological characteristics and aerodynamic variables.According to the Penman-Monteith FAO 56 model, the ET ref is estimated by the following: where ET re f is the reference crop evapotranspiration (mm•d −1 ), ∆ is the slope of temperature variation with saturated vapor pressure (kPa ), and G is the soil heat flux density (MJ•m −2 •d −1 ), which is negligible relative to R n , especially for vegetated coverage and at daily time steps, thus it may be ignored and G = 0. γ is a psychrometric constant (kPa• • C −1 ), T is the daily mean air temperature ( • C), u 2 is the wind speed at 2 m height (m•s −1 ), e s is the saturation vapor pressure (kPa), and e a is the actual vapor pressure (kPa).
For the calculation of ET ref , wind speed at 2 m height can be converted from the wind speed measured at 10 m based on the Penman-Monteith FAO 56 model as follows: where z is the measured height above the ground surface (m) and u 2 and u z are the wind speed at 2 and z meters height, respectively.In this study, the total net solar radiation in the background is estimated based on the measured data of sunshine duration, combined with the Angstrom formula [41], which is given below: where R s is the daily total solar radiation reaching the ground (MJ•m −2 •d −1 ), n is the daily measured sunshine duration hours (h), N is the greatest possible sunshine duration hours (h), R a is the extraterrestrial solar radiation (MJ•m −2 •d −1 ), and a s and b s are the regression coefficients-according Water 2017, 9, 486 5 of 20 to Chen et al. [42], the values were set to 0.19 and 0.53, respectively.Equations ( 4)-( 6) are used to estimate incoming shortwave and outgoing long-wave radiation: where R ns is the net incoming shortwave radiation (MJ•m −2 •d −1 ), α (=0.23) is the albedo of the reference grassland, R nl is the net outgoing long wave radiation ), and R s /R so is the relative shortwave radiation.Detailed calculation procedures can be found in Chapter 3 of the FAO Irrigation and Drainage paper No. 56 [27].

Statistical Test for Trend Analysis
The Mann-Kendall (MK) statistical test [43,44] is a rank-based nonparametric method that has been extensively applied to the analysis of trends in reference evapotranspiration and related meteorological factors in hydro-meteorological time series data.One advantage of this MK test is that the dataset does not require the underlying data to follow a specific distribution.It can handle non-normality, missing values, or seasonality and is reliable for biased variables [45,46].A nonparametric MK test was employed to show the trend variation in annual and seasonal ET ref for each region in this study.Positive (negative) normalized Z values from the MK test indicate rising (declining) trends, respectively.An absolute Z value of the MK test exceeding 1.96 indicates that ET ref has a significant change trend at the 95% confidence level, while an extremely noteworthy change trend was identified at the 99% significant level for the absolute Z value of the MK test more than 2.56.Additionally, to identify whether a trend exists, the Theil-Sen slope estimator developed by Sen [47] was used to determine the extent of a trend.This technique has been generally applied to detect the slope of a trend line in a hydro-meteorological time series dataset [48].

Spatial Interpolation Method for ET ref
The IDW (Inverse Distance Weighted), a simpler numerical method, was used to investigate the spatial distribution of ET ref in Jiangsu province.It has been widely applied to the spatial interpolation of climatic and hydrological point data like ET ref [49], as well as rainfall [50].Among various interpolation techniques, the IDW method was employed for the spatial distribution in the present study due to its easiness and estimated accuracy in comparison with other interpolation models like kriging.So we calculated the annual average estimated ET ref at each station and then spatially interpolated it.All spatial interpolations were accomplished using ArcGIS 9.3 software (Environmental Systems Research Institute (ESRI), Redlands, CA, USA).The justification for using the IDW method is that this technique can compute the spatially interpolated values very fast and accurately.

Calculation of Sensitivity Coefficient
To identify the sensitivity of ET ref related to climatic factors, a mathematically based sensitivity coefficient proposed by McCuen [29] in 1974 was applied in this study.Equation ( 7) is used to calculate the sensitivity of ET ref as a mathematical basis: where S(v i ) is the non-dimensional sensitivity coefficient of ET re f with respect to the climate variable (v i ).Positive (negative) sensitivity coefficients indicate that ET re f increases or decreases with increases Water 2017, 9, 486 6 of 20 in the climate variable v i (where v i is TA, VP, WS or SD).The absolute value measures the effect a given climate variable has on ET re f .∂ET re f /∂v i was calculated by partial differentiation of Equation ( 1); the detailed calculations can be found in the appendix of this paper.

Detrending Method
To evaluate the contributions of key meteorological factors to the change trend of ET ref quantitatively, the detrending method, a statistical and mathematical operation of removing trend from the series, was employed in this study [8,30,31].This is a simple but effective method.The main step of this method calculated by simple linear regression is demonstrated as follows: where b and â are the fitted linear regression coefficient and constant, respectively.ĝt is the fitted trend in time t.Then we removed the trend by subtracting the value of the trend line from the original data.Furthermore, to avoid negative values (e.g., VP, WS, and SD), the detrended dataset is currently defined as below: where x t is the original time series of meteorological factors, ĝ1 , corresponds to the first value of ĝt , and y t is the ultimate detrended time series.After removing the change trend in key meteorological factors, recalculating the ET ref by using the detrended data series of one stationary meteorological factor and using the original data of remaining variables, comparing the recalculated ET ref with the original value, the difference between them is identified as the contributions of that variable to the change trend of ET ref .In order to quantify the contribution, an assessment index R is used to calculate by the following: where ET 0 0 and ET R 0 represent the original and recalculated ET ref , respectively.n denotes the length of the time series.The larger the absolute value of R, the greater the contribution of that meteorological variable to the change trend of ET ref .R > 0, R < 0 and R = 0 indicate that the trend in that meteorological variable has a positive, negative, or no contribution to the trend in ET ref , respectively.

Spatial Distribution of Seasonal and Annual ET ref
Different spatial distributions of seasonal and annual ET ref are shown in Figure 2. The ET ref was higher in summer (335 mm~381 mm), followed by spring (221 mm~297 mm), autumn (176 mm~221 mm), and winter (85 mm~105 mm).The mean annual ET ref over the whole Jiangsu province was about 892 mm during 1961-2015 (Table 1).In spring, the highest value was found in the Huaibei region.However, in the summer season, the spatial distribution of ET ref showed a discrepancy, where the highest values were noticed in the Sunan region.Compared to the two seasons like autumn and winter, which had a similar spatial distribution of ET ref , the higher value was observed in the southeast part of Jiangsu province (Sunan and Jianghuai region).In general, the annual spatial distribution of ET ref had a relatively high value in northwest and south parts and low value in the central part of Jiangsu province (Figure 2a).

Temporal Trends in Seasonal ETref and Climatic Factors
Table 2 presents the results of temporal trend for ETref.On the annual scale, the ETref values showed a slightly insignificant decreasing trend (−0.245 mm•a −2 ) in the whole of Jiangsu province.A significant decreasing tendency (at 99.9% confidence level) was identified in the Huaibei region (−1.468mm•a −2 ), whereas a slightly increasing trend of ETref was observed in the Jianghuai (0.207 mm•a −2 ) and Sunan (0.592 mm•a −2 ) regions.ETref showed a significant increasing tendency in spring and decreasing tendency in the summer, with rates of 0.526 mm•a −2 and −0.701 mm•a −2 , respectively, at the seasonal scale.The trend of ETref in winter was lower than in autumn.ETref in spring showed an increasing trend, especially significant in the Jianghuai and Sunan regions.The most significant decreasing trend of ETref was presented in the summer, particularly in Huaibei.A slight increasing trend was shown in the whole region in autumn, mainly due to the slight significant declining trend in Huaibei region combined with the slight significant increasing trends in Jianghuai and Sunan regions.Despite a decreasing trend shown in winter for the Huaibei region, the whole region exhibited a slight increasing trend, mainly due to the contributions of the slight increasing trends in the Jianghuai and Sunan regions.
The slopes of linear trends in climatic factors are shown in Table 3. Mean temperature in different seasons and regions revealed a significant increasing tendency, except for the summer   ET ref showed a significant increasing tendency in spring and decreasing tendency in the summer, with rates of 0.526 mm•a −2 and −0.701 mm•a −2 , respectively, at the seasonal scale.The trend of ET ref in winter was lower than in autumn.ET ref in spring showed an increasing trend, especially significant in the Jianghuai and Sunan regions.The most significant decreasing trend of ET ref was presented in the summer, particularly in Huaibei.A slight increasing trend was shown in the whole region in autumn, mainly due to the slight significant declining trend in Huaibei region combined with the slight significant increasing trends in Jianghuai and Sunan regions.Despite a decreasing trend shown in winter for the Huaibei region, the whole region exhibited a slight increasing trend, mainly due to the contributions of the slight increasing trends in the Jianghuai and Sunan regions.
The slopes of linear trends in climatic factors are shown in Table 3. Mean temperature in different seasons and regions revealed a significant increasing tendency, except for the summer season.On the contrary, wind speed and sunshine duration were shown to have significant declining trends, in addition to the sunshine duration in spring, which demonstrated a slight decreasing trend in the Huaibei region and a slight increasing trend in the Jianghuai and Sunan regions.However, compared to the other three climatic factors, the change trend of VP is too small, only exhibiting a significant downward trend in summer, except for the Huaibei region.

Sensitivity Coefficients of ET ref to Climatic Factors
The sensitivity coefficients of seasonal and annual trends of ET ref to each climatic factor by partial derivative are displayed in Figure 3.There are clear differences in the sensitivity characteristics for the four climatic factors and a slight discrepancy exists in the three regions of Jiangsu province.Generally, TA (air temperature), WS (wind speed), and SD (sunshine duration) produced a positive effect, whereas the VP (actual vapor pressure) had an adverse effect on seasonal and annual ET ref .The results show that S(TA) and S(SD) were larger in the summer, but smaller in the winter, indicating that ET ref was most sensitive to TA and SD in the summer compared to other seasons.The sensitivity coefficient of ET ref to climatic factors ranked in a sequence as VP > TA > SD > WS.It can be said that the actual vapor pressure (VP) and air temperature (TA) were the greatest sensitive factors, while wind speed (WS) was shown to be the least sensitive factor of ET ref , contrary to TA and SD, which had a relatively higher value in winter but lower in summer (Figure 3).

Contributions of Climatic Factors to the Trends in ETref
The original and detrended climatic factors over the whole Jiangsu province are presented in Figure 4. Noticeable differences can be seen between the original and detrended datasets for each climatic factor, with a negative trend for TA and positive trends for VP, WS, and SD.The detrended TA was lower than the original data series, whereas the detrended VP, WS, and SD were higher than the original data series, especially the most obvious difference in WS.However, the difference between the original and detrended data series of VP was smaller than for the other three climatic factors at an annual scale in the whole region.

Contributions of Climatic Factors to the Trends in ET ref
The original and detrended climatic factors over the whole Jiangsu province are presented in Figure 4. Noticeable differences can be seen between the original and detrended datasets for each climatic factor, with a negative trend for TA and positive trends for VP, WS, and SD.The detrended TA was lower than the original data series, whereas the detrended VP, WS, and SD were higher than the original data series, especially the most obvious difference in WS.However, the difference between the original and detrended data series of VP was smaller than for the other three climatic factors at an annual scale in the whole region.

Contributions of Climatic Factors to the Trends in ETref
The original and detrended climatic factors over the whole Jiangsu province are presented in Figure 4. Noticeable differences can be seen between the original and detrended datasets for each climatic factor, with a negative trend for TA and positive trends for VP, WS, and SD.The detrended TA was lower than the original data series, whereas the detrended VP, WS, and SD were higher than the original data series, especially the most obvious difference in WS.However, the difference between the original and detrended data series of VP was smaller than for the other three climatic factors at an annual scale in the whole region.To understand a change event in seasonal and annual ETref due to the contributing factors, it is essential to analyze the contribution of each climatic factor to the change in seasonal and annual ETref.Results in Figure 5 demonstrate that the most important climatic factor dominating the trends in annual ETref on the whole Jiangsu province followed the order WS > SD > TA > VP.In general, the ETref recalculated with detrended TA and detrended VP was smaller than that calculated from the original climatic factors.On the contrary, the ETref recalculated with detrended WS and detrended SD was larger than that calculated from the original climatic factors.In other words, TA and VP had positive contributions to ETref, except VP in the Huaibei region, which exhibited a negative contribution, while WS and SD showed the opposite ones.In the end, the negative contributions would offset the positive one.To understand a change event in seasonal and annual ETref due to the contributing factors, it is essential to analyze the contribution of each climatic factor to the change in seasonal and annual ETref.Results in Figure 5 demonstrate that the most important climatic factor dominating the trends in annual ETref on the whole Jiangsu province followed the order WS > SD > TA > VP.In general, the ETref recalculated with detrended TA and detrended VP was smaller than that calculated from the original climatic factors.On the contrary, the ETref recalculated with detrended WS and detrended SD was larger than that calculated from the original climatic factors.In other words, TA and VP had positive contributions to ETref, except VP in the Huaibei region, which exhibited a negative contribution, while WS and SD showed the opposite ones.In the end, the negative contributions would offset the positive one.On the annual scale, the contributions of WS and SD were greater than the total contributions of TA and VP, indicating a downward trend of ET ref in the whole of Jiangsu province.Nevertheless, the variations of the climatic factors in three sub-regions caused different contributions to ET ref .Similar results of ET ref calculated with detrended WS, SD, and TA can be found in each sub-region of Jiangsu province; however, a small but interesting difference can be obtained between the original ET ref and the recalculated ET ref with detrended VP in the three sub-regions.The recalculated ET ref with detrended VP was higher than that calculated from the original climatic factors in the Huaibei region; on the contrary, the result was opposite in the Sunan region, which exhibited a lower value compared with the original ET ref .The reverse results demonstrated the varying trend of VP in the annual scale, which had been detected in Table 3; this also means it plays an important role in the contribution of an ET ref trend across the different regions to some extent, although the variation is small and non-significant.From the results presented above, the contribution of ET ref was a synergistic effect of various climatic factors; for example, significant decreasing WS and SD were the most crucial factors for the decreasing trend in ET ref .TA was the main contributing factor to the increasing trend in ET ref in each region of Jiangsu province, followed by WS and SD.The non-significant trend in VP made the smallest contribution to the trend of ET ref ; however, it could not be ignored due to the variability and the effect of the contribution to ET ref , as well as being the most sensitive factor to ET ref .In contrast, the contribution of WS was greater than other factors.
On the seasonal scale, due to the limited space, we established an assessment index R that was used to quantify the contributions of TA, VP, WS, and SD to the change trends of ET ref .Similar findings like Figure 5 can be obtained from Table 4, where TA had a positive contribution to the trend of ET ref ; however, WS and SD had an adverse contribution in most regions except Jianghuai and Sunan in spring.In spring, the contribution of TA was the largest, followed by WS and SD except for the Huaibei region, whereas WS made a much greater contribution, followed by TA, VP, and SD.In summer, the contribution of SD played a dominant role in each region.In autumn, the contribution of WS was larger than other climatic factors, followed by TA, SD, and VP.Meanwhile, in winter, the contribution of the WS to ET ref was much stronger than in the other three seasons.TA was the unique positive contributing factor, followed by WS, which led to an increasing trend in ET ref .However, contrary to our expectations, the contribution of VP was much greater than in other seasons, and it also exhibited a negative contribution to ET ref .This also confirmed the importance of VP, which has been discussed above, especially in the winter for the whole region.

Evaporation Paradox Analysis with Trends of TA and ET ref
The evaporation paradox has become a debatable issue against the backdrop of global climate warming and decreasing evaporation scenarios [23].It can be seen from Figure 6 that the mean temperature (1961-2015) had a significant increase in the whole of Jiangsu province when the rate of climatic tendency ranged between 0.013 • C•a −1 and 0.047 • C•a −1 , with an especially high rate in the Jianghuai and Sunan regions.The climatic trend rate of ET ref ranged between −2.46 mm•a −2 and 2.93 mm•a −2 , indicating a distinctly various distribution on the spatial scale.A positive trend of ET ref was found in the Jianghuai and Sunan regions, while a significant overall downward trend was observed in Huaibei; the reason is the decrease of WS and SD, despite being accompanied by a significant increase in TA.
Meanwhile, in the spring, the distribution of the TA was the same as the annual-both strong increasing tendencies.ET ref increased significantly in Jianghuai and Sunan and had increasing and decreasing trends, both insignificant, in Huaibei.In summer, a significant increasing trend of TA showed in southeast Jiangsu province; moreover, the increasing trend in other regions was insignificant.ET ref showed a downward trend, especially significant in Huaibei and along the lower Yangtze River Basin.In autumn, both the TA and ET ref were similar to the annual scale distribution.In the winter, a significant downward trend of ET ref existed only in Huaibei, although a strong upward trend of TA was exhibited in the whole of Jiangsu province.These results indicated that the evaporation paradox was apparent in Jiangsu province, especially in summer and in the Huaibei region in autumn, and also on an annual scale.observed in Huaibei; the reason is the decrease of WS and SD, despite being accompanied by a significant increase in TA.
Meanwhile, in the spring, the distribution of the TA was the same as the annual-both strong increasing tendencies.ETref increased significantly in Jianghuai and Sunan and had increasing and decreasing trends, both insignificant, in Huaibei.In summer, a significant increasing trend of TA showed in southeast Jiangsu province; moreover, the increasing trend in other regions was insignificant.ETref showed a downward trend, especially significant in Huaibei and along the lower Yangtze River Basin.In autumn, both the TA and ETref were similar to the annual scale distribution.In the winter, a significant downward trend of ETref existed only in Huaibei, although a strong upward trend of TA was exhibited in the whole of Jiangsu province.These results indicated that the evaporation paradox was apparent in Jiangsu province, especially in summer and in the Huaibei region in autumn, and also on an annual scale.

Change Trend of ETref and Its Contributing Factors
This study is driven by the need to know why ETref has changed in Jiangsu province and what factors are responsible.This study, therefore, applied physical modeling to determine the sensitivity of ETref to different climatic factors and a detrending method to identify the main contributions of climatic factors to ETref changes.In the current research, decreasing and increasing trends were found in Jiangsu province.A decreasing trend showed in annual and summer ETref, while an increasing trend was observed in spring, autumn, and winter ETref for the whole province.The results in this study highlighted the greater sensitivity of ET ref in Jiangsu province to TA, followed by SD during the summer, when ET ref reached its peak position.VP had a much higher impact on ET ref than the other three climatic factors, especially in the winter, when ET ref reached its lowest.WS had the smallest influence on the ET ref estimation, with a higher sensitivity in the winter but lower in the summer.The sensitivity of seasonal and annual ET ref to climatic factors was variable in different regions.Furthermore, the declining ET ref trend was more sensitive to VP in the present study.Our findings provide a new perspective for further study, due to the fact that previous research usually chose relative humidity as the influence factor, such as in the basins of the Yangtze River, Tao'er River, and Wei River in China [4,8,[51][52][53].However, the variation of ET ref was most sensitive to WS in the basins of the Shiyang River and Jinghe River [11,54].In addition, TA and R n were the utmost sensitive parameters in the Haihe River basin and Poyang Lake catchment in China [55,56].It can be said that ET ref has various responses to meteorological factors in various regions of China.
The contributions of climatic factors to the changes in ET ref were also diverse in various regions.In this study, WS was the primary variable, followed by SD to annual ET ref across the whole of Jiangsu province, which was similar to the outcomes reported by Wang et al. [36] in the Yellow River basin.Zhang et al. [57] indicated that the decreasing trend in ET ref was due to a decrease in WS and R n across the Tibetan Plateau.However, WS and TA were the predominant factors contributing to its inter-annual variation in Northwest China [58].Results of changes in ET ref and its contributing factors in Jiangsu province agreed with the conclusion of Gao et al. [59], who showed that WS was the chief factor, followed by R n in the west of the Liao River basin.In addition, we found that VP is also an important influence factor, especially in winter in the whole region, although the change trend of VP is very small.
As the decreasing WS had been widely reported in various regions, Zhang et al. [60] proposed that the weakening atmospheric circulation may be the explanation for the WS decreasing.However, the significant weakened East Asian Monsoon is also closely linked to decreasing WS [26,61].In addition, the increase of terrestrial surface roughness due to the increased vegetation coverage, such as the afforestation put forward by Vautard et al. [62], could also explain the decreasing WS.Meanwhile, the significantly decreasing SD also made a negative contribution to the decreasing ET ref , closely followed by WS, especially in the summer, while the functional mechanism of the decreasing SD is still a matter for debate.A significant number of studies have shown that decreasing solar radiation had occurred in different regions before 1980s [63,64]; however, this decreasing trend has been mitigated since the 1990s, via the phenomenon we call "From dimming to Brightening" [65][66][67].Combined with the sunshine duration in Figure 4, this phenomenon is seemingly not evident in Jiangsu province.Meanwhile, some recent studies have shown that the decreasing SD is mainly caused by the increased regional cloud cover [65,68] and increasing aerosol loading from anthropogenic emissions of pollutants [69,70], coupled with significant declining WS.These phenomena were confirmed by Zhang et al. [71] in east and south China.Furthermore, due to the absence of measured solar radiation, the parameters (a and b) in the Angstrom formula might vary across the past 55 years as the aerosol optical depths have changed in this region.This is an important potential constraint on the observed effects of changes in R n on ET ref , which will surely have been underestimated as a result of using this approximation of R s from the sunshine duration.

Evaporation Paradox
The "evaporation paradox" is a controversial topic that has earned attention among researchers worldwide [22,72].The results indicated that the evaporation paradox was apparent in Jiangsu province, mainly in the summer and in Huaibei region in the autumn.Furthermore, the significant increase in TA in Jiangsu province appeared to be in agreement with the global warming reported in recent decades [10].However, interestingly, TA showed a significant result with the decrease of ET ref (Tables 2 and 3).These findings disagree with the earlier results of Ye et al. [56], Liu and Zeng [73], and Li et al. [74].The "evaporation paradox" occurred due to the significant decrease of WS and SD in the whole of Jiangsu province.The strong influences on ET ref from these contributing factors compensate for the decrease of ET ref from the significant increase in TA, resulting in an overall decrease in ET ref .Such findings are echoed by the similar works of Roderick and Farquhar [17,18], Xu et al. [8], Cong et al. [14], and Zuo et al. [52].However, it is evident that ET ref has decreased in many parts of the world in the context of global warming.Xu et al. [75] summarized the major reason: higher air temperature leads to more powerful atmospheric water holding capacity, resulting in an increase of cloud cover and a decline of solar radiation reaching the Earth's surface.

Influence on Agricultural Production and Water Management
Jiangsu is a major agricultural province in China, and main crops such as winter wheat and single-season rice have been greatly affected by climate change, especially in regions with variable climate [38].Decreasing wind speed (WS), sunshine duration (SD), actual vapor pressure (VP), and increasing air temperature (TA) were the driving forces behind a slightly decreasing ET ref in the Jiangsu province, and climate change may complicate impacts on agricultural production and water resource management.The significantly decreasing ET ref in summer may be responsible for greater water reserves and an increase in regional vegetation cover and crop production, especially the cultivation of single-season rice across the whole of Jiangsu province, which further assists with irrigation water resource distribution and management.Meanwhile, excessive rainfall caused frequent floods in summer; this should be paid more attention, especially in the Huaibei region, which belongs to the Huai River Basin [76].However, the increasing ET ref in autumn would not be suitable for the seeding of winter wheat in Jianghuai and Sunan, and could affect crop growth and development.Moreover, the slight increase in ET ref in the winter would reduce the soil moisture and be harmful for the overwintering of winter wheat.More seriously, the significant increase in ET ref in spring may inhibit the reviving and jointing of winter wheat and would likely increase the frequency of spring droughts, despite the abundant rainfall in this province of eastern China.However, the significant decrease i ET ref in Huaibei region in summer, autumn, and on an annual scale would promote single-season rice and winter wheat growth, whereas the slight decreasing trend of ET ref in the winter and the increasing trend in the spring could not be ignored.
For the past few years, with the rapid increase of population and intensified land utilization, the rising problems of regional water resources and ecological environment have been increasingly highlighted.These issues deserve more attention in future studies on how the changes of annual and seasonal ET ref are associated with potential and actual evapotranspiration.Integration with remote sensing techniques will be important to determine the subsequent effect on the regional water balance of Jiangsu province in the light of population growth, urbanization, and climate change.

Conclusions
This research examined the spatiotemporal changes in reference evapotranspiration (ET ref ) based on the Penman-Monteith FAO 56 model in Jiangsu province, a major economic and agricultural province of eastern China, by studying the daily weather data of 60 meteorological stations from 1961 to 2015.In addition, physical and mathematical modeling and the detrending method were applied to determine the sensitivity of ET ref and the main contributing climatic factors to ET ref .The following findings were drawn from this study: (1) Decreasing and increasing trends were identified in Jiangsu province at seasonal and annual scales.A significant declining trend of ET ref was observed in the Huaibei region, while other regions demonstrated a slight rising trend of ET ref .
ET ref also showed a significant increasing tendency in spring and a decreasing tendency in summer.(2) VP was the most sensitive climatic factor, followed by TA, SD, and WS in each region.ET ref was more sensitive to TA and SD in the summer but less so in the winter; however, the sensitivities of WS and VP to ET ref were the opposite in Jiangsu province.
(3) The contribution of ET ref was a synergistic effect of various climatic factors; significantly decreasing WS and SD were the most crucial factors contributing to the decreasing trend in ET ref , while TA was the main contributing factor for the increasing trend in ET ref , followed by WS and SD.The non-significant trend in VP made the smallest contribution to the trend of ET ref ; however, it could not be ignored due to the variability and the effect of the contribution to ET ref , as well as being the most sensitive factor to ET ref , especially in winter.
(4) It was evident that the phenomenon of evaporation paradox was apparent in Jiangsu province, especially in the summer and in the Huaibei region in the autumn.(5) In summary, these results will provide a theoretical background and practical guidance for future water resource distribution and management and lay the basis for the future study of agricultural water requirements and economic development in the whole of Jiangsu province in the scenario of climate change. ET

Figure 2 .
Figure 2. Spatial distribution of seasonal and annual ET ref in Jiangsu Province during 1961-2015.Note: (a)-(e) represent ET ref in annual, spring, summer, autumn, and winter seasons respectively.

3. 2 .
Temporal Trends in Seasonal ET ref and Climatic Factors Table 2 presents the results of temporal trend for ET ref .On the annual scale, the ET ref values showed a slightly insignificant decreasing trend (−0.245 mm•a −2 ) in the whole of Jiangsu province.A significant decreasing tendency (at 99.9% confidence level) was identified in the Huaibei region (−1.468mm•a −2 ), whereas a slightly increasing trend of ET ref was observed in the Jianghuai (0.207 mm•a −2 ) and Sunan (0.592 mm•a −2 ) regions.

Figure 3 .
Figure 3. Sensitivities of seasonal and annual ETref in relation to climatic factors in (a) Whole, (b) Huaibei, (c) Jianghuai, and (d) Sunan regions of the Jiangsu province from 1961 to 2015.

Figure 3 .
Figure 3. Sensitivities of seasonal and annual ET ref in relation to climatic factors in (a) Whole, (b) Huaibei, (c) Jianghuai, and (d) Sunan regions of the Jiangsu province from 1961 to 2015.

Figure 3 .
Figure 3. Sensitivities of seasonal and annual ETref in relation to climatic factors in (a) Whole, (b) Huaibei, (c) Jianghuai, and (d) Sunan regions of the Jiangsu province from 1961 to 2015.

Figure 4 .
Figure 4.The original and detrended climatic factors (a) TA, (b) VP, (c) WS, and (d) SD for the whole of Jiangsu province during 1961−2015.

Figure 4 .
Figure 4.The original and detrended climatic factors (a) TA, (b) VP, (c) WS, and (d) SD for the whole of Jiangsu province during 1961−2015.

Figure 5 .
Figure 5. Contributions of climatic factors to the trends in ET ref in (a) Whole, (b) Huaibei, (c) Jianghuai, and (d) Sunan regions of the Jiangsu province during 1961-2015.

Figure 6 .
Figure 6.Spatial distributions of climate tendency rates of the TA and ETref from 1961 to 2015 in Jiangsu province.Note: S and NS represent significant and non-significant trends at the 0.05 level of the TA and ETref, respectively.(a), (c), (e), (g), and (i) represent the trends of TA in annual, spring, summer, autumn, and winter seasons, respectively, (b), (d), (f), (h), and (j) represent the trends of ETref in annual, spring, summer, autumn and winter seasons, respectively.

Figure 6 .
Figure 6.Spatial distributions of climate tendency rates of the TA and ET ref from 1961 to 2015 in Jiangsu province.Note: S and NS represent significant and non-significant trends at the 0.05 level of the TA and ET ref , respectively.(a), (c), (e), (g), and represent the trends of TA in annual, spring, summer, autumn, and winter seasons, respectively, (b), (d), (f), (h), and (j) represent the trends of ET ref in annual, spring, summer, autumn and winter seasons, respectively.

1 .
Change Trend of ET ref and Its Contributing FactorsThis study is driven by the need to know why ET ref has changed in Jiangsu province and what factors are responsible.This study, therefore, applied physical modeling to determine the sensitivity of ET ref to different climatic factors and a detrending method to identify the main contributions of climatic factors to ET ref changes.In the current research, decreasing and increasing trends were found in Jiangsu province.A decreasing trend showed in annual and summer ET ref , while an increasing trend was observed in spring, autumn, and winter ET ref for the whole province.

Table 1 .
Basic climatic characteristics of three sub-regions and the whole of Jiangsu province during 1961-2015.

Table 1 .
Basic climatic characteristics of three sub-regions and the whole of Jiangsu province during 1961-2015.

Table 2 .
Temporal trends for seasonal and annual ETref (mm•a −2 ) in each region of Jiangsu province from 1961 to 2015.

Table 2 .
Temporal trends for seasonal and annual ET ref (mm•a −2 ) in each region of Jiangsu province from 1961 to 2015.Note: + , *, ** and *** indicate the significance level of 0.1, 0.05, 0.01, and 0.001, respectively.β is the estimated slope of ET ref trends, β > 0 and β < 0 indicate an increasing and a decreasing trend, respectively.

Table 3 .
Temporal trends for climatic factors in Jiangsu province from 1961 to 2015.

Table 4 .
Values of assessment index R for quantifying the contributions of climatic factors to the change trends of seasonal ET ref in each region of Jiangsu province.
Note: Values in bold face indicate the dominant contributing factors.