Intra-Annual Variability of Evapotranspiration in Response to Climate and Vegetation Change across the Poyang Lake Basin, China

: Improving understanding of changes in intra-annual variability (IAV) of evapotranspiration (ET) and the underlying drivers is an essential step for modeling hydrological processes in response to global change. Previous studies paid special attention to climatic regulations of IAV of ET. However, ignoring the role of landscape characteristics (e.g., vegetation coverage) can introduce great uncertainty in the explanation of ET variance. In this work, the Poyang Lake Basin, which is a typical humid basin in China, was taken as the study area. It has experienced an obvious climate change and revegetation since the 1980s. Here, trends of IAV of ET and their responses to four climatic variables (i.e., air temperature, precipitation, downward shortwave radiation and wind speed) and vegetation coverage were explored from 1983 to 2014. The results show that IAV of ET exhibited contrary trends during the past decades. It signiﬁcantly ( p < 0.05) declined with a signiﬁcant linear slope of − 0.52 mm/year before 2000, and then slightly increased (slope = 0.06 mm/year, p > 0.05) over the basin, which was generally consistent with the IAV of temperature and radiation. The proposed variables could well capture the change in IAV of ET, while their dominators were different during the two contrasting phases mentioned above. The IAV of radiation and temperature dominated the change of the IAV of ET over 77.82% and 35.14% of the basin, respectively, before and after the turning point. Meanwhile, the rapid increase in vegetation coverage, which was associated with afforestation, signiﬁcantly ( p < 0.05) reduced IAV of ET over about 35% of the study area. The achievements of this study should be beneﬁcial for a sophisticated awareness of responses of intra-annual variability of ET to climate and land cover changes at the basin scale.


Introduction
Evapotranspiration (ET) is the only way in which terrestrial water returns to the atmosphere in the earth system [1]. It is a crucial component to maintain land surface water and energy balance [2,3]. Terrestrial ET is mainly composed of two parts [4]. The first is plant transpiration, which is mainly controlled by leaf stomata, and the other is water evaporation from moist surfaces of vegetation and/or soil. Transpiration and evaporation are regulated by both abiotic (e.g., climate and soil) and biotic (e.g., leaf area and coverage) variables [5]. Hence, it is meaningful and challenging to understand comprehensive and complex responses of terrestrial ET to changes in climate and land cover.
Recently, variabilities of ET at the intra-annual time scale have been increasingly paid attention to [6][7][8][9]. It is a rhythmic characteristic of water loss, which is usually indicated by The watershed size is about 1.6225 × 10 5 km 2 , i.e., 9% of the total area of the Yangtze River Basin. The Poyang Lake Basin is a typical energy-limited basin with a sub-tropical monsoon climate [19]. That is, actual evaporation is limited by the rate of energy supply rather than water supply. The annual temperature and rainfall are in the range 16.3-17.5 °C and 1341-1943 mm, respectively. The vegetation coverage across the basin experienced an obvious change during the past decades [18]. The Poyang Lake Basin suffered an intensive deforestation after 1950s, where forest coverage decreased to 33.1% in 1978. Then, with the implementation of afforestation, greening was widely observed across the basin. The forest coverage continuously increased since 1990s and raised to 63.1% until 2011. That is, the coverage showed a fast increase mainly resulting from anthropogenic activity. This provides a natural laboratory for learning impacts of vegetation restoration on IAV of ET. In our present work, the open water bodies of the Poyang Lake Basin were excluded from the subsequent analysis to minimize the effects of water exchange to evaporation.

Evapotranspiration
Two long-term satellite-based ET datasets were adopted to calculate the IAV of ET in this study. The first was the Advanced Very High-Resolution Radiometer (AVHRR) monthly ET data (unit: mm/month) with a spatial resolution of 8 km from 1983 to 2006 [20,21]. The other was the Moderate-Resolution Imaging Spectroradiometer (MODIS) monthly ET product (unit: mm/month) with a resolution of 1 km covering the period of 2001-2014 [22,23]. For the relative high accuracy, these two datasets have been widely The vegetation coverage across the basin experienced an obvious change during the past decades [18]. The Poyang Lake Basin suffered an intensive deforestation after 1950s, where forest coverage decreased to 33.1% in 1978. Then, with the implementation of afforestation, greening was widely observed across the basin. The forest coverage continuously increased since 1990s and raised to 63.1% until 2011. That is, the coverage showed a fast increase mainly resulting from anthropogenic activity. This provides a natural laboratory for learning impacts of vegetation restoration on IAV of ET. In our present work, the open water bodies of the Poyang Lake Basin were excluded from the subsequent analysis to minimize the effects of water exchange to evaporation.

Evapotranspiration
Two long-term satellite-based ET datasets were adopted to calculate the IAV of ET in this study. The first was the Advanced Very High-Resolution Radiometer (AVHRR) monthly ET data (unit: mm/month) with a spatial resolution of 8 km from 1983 to 2006 [20,21]. The other was the Moderate-Resolution Imaging Spectroradiometer (MODIS) monthly ET product (unit: mm/month) with a resolution of 1 km covering the period of 2001-2014 [22,23]. For the relative high accuracy, these two datasets have been widely applied in studies on water cycle in response to climate change and human activities at multiple spatial scales [24,25]. In this study, both of these datasets were aggregated to grids of 1/12 • × 1/12 • to match the spatial resolution of the vegetation coverage data derived from remote sensing vegetation index mentioned below. Afterwards, the two ET datasets were fused according to the overlapping period to obtain a long time series from 1983 to 2014 [18]. Additionally, field observations of hydrological variables, including precipitation, soil moisture and runoff, were collected to validate our satellite-based ET data [26].

Climatic Variables
Four climatic variables, i.e., surface air temperature (AT, unit: • C), precipitation (PR, unit: mm h −1 ), downward shortwave radiation (SR, unit: W m −2 ) and wind speed (WD, unit: m s −1 ) were used in this study, which are recognized the principal regulators of ET [27]. They were obtained from a monthly dataset of the high-resolution Chinese Meteorological Forcing Dataset (CMFD) [28][29][30]. CMFD was made through fusion of satellite-, reanalysisand in-situ station-based data with a spatial resolution of 0.1 • . It is the first high spatialtemporal resolution gridded near-surface meteorological dataset and has been widely used in studies of land surface process in China [31,32]. In the present work, the monthly climatic data during 1983-2014 were resampled to the resolution of 1/12 • by the nearest neighbor method to match the resolution of the ET data.

Vegetation Coverage
The Normalized Difference Vegetation Index (NDVI) product of NOAA Global Inventory Monitoring and Modeling System (GIMMS, version number 3g.v1) was adopted to calculate vegetation coverage across the study area [33]. This dataset covers the period of 1983-2014 with a temporal resolution of half-month and a spatial resolution of 1/12 • (about 0.083 • ) [34]. After a series of corrections (e.g., orbital drift effects, calibration, viewing geometry, stratospheric volcanic aerosols and other errors unrelated to vegetation change) [35,36], GIMMS NDVI 3g.v1 is well consistent with other high-precision NDVI products (e.g., MODIS), and has been widely applied in detecting responses of vegetation to climate change over the world [37,38]. In this study, the maximum value composite was carried out for the NDVI dataset to further get rid of white noise points [39].
A long-term time series of vegetation coverage was derived following an improved approach described in Wang et al. [18]. This approach employed an optimized dimidiate pixel model, in which dynamic background values and Moderate-Resolution Imaging Spectroradiometer (MODIS) Vegetation Continuous Fields (VCF) (MOD44B) were used. By validation, the NDVI-based vegetation coverage of this work performed well against the field investigation data of vegetation cover which were collected from the Statistics Yearbook of Jiangxi Province and Statistics Yearbook of China [18,26]. Overall, the correlation coefficient (r) between them was 0.94, indicating that the satellite-based vegetation cover could appropriately capture the revegetation process in the Poyang Lake Basin.

Trend Analysis
To detect changes of IAV of ET over the Poyang Lake Basin, two regression models, i.e., linear regression model (Equation (1)) and piecewise regression model (Equation (2)) [40], were adopted and compared. The former assumed there was only one trend over the entire study period, while the latter one assumed that the trend of SD ET significantly changed during the past decades [41][42][43].
In the models, y indicates the standard deviation of monthly ET in a year (i.e., SD ET ), t indicates the year, and α indicates the turning point (TP) of the SD ET time-series data. k, k 1 and k 2 indicate the magnitudes of global SD ET trend, SD ET trend before the TP and SD ET trend after the TP, respectively. b, b 1 , b 2 and ε indicate the intercept and the residual random error of the regression model, respectively. The models were fitted to the SD ET timeseries by the least squares approach. The performances of the models were evaluated by coefficient of determination (R 2 ), root mean square error (RMSE) and Akaike information criterion (AIC). Generally, a model with a relatively lower AIC value performs better than another. In this study, the value of AIC was quantified following Peng et al. [43]. The difference between the AIC value of linear regression model and that of piecewise regression model (termed δAIC) was calculated. Here, if the δAIC was lower than −2, the piecewise regression model was significantly preferred, and vice versa. Additionally, changes in IAV of climatic variables and vegetation coverage were also detected over the study area.

Path Analysis
Relationships of IAV of ET to the climatic variables and vegetation coverage were explored by partial correlation analysis. Moreover, path analysis was applied to investigate direct effects of the environmental variables on IAV of ET pixel by pixel over the study area. In the analysis, path coefficients were standardized weights which can be used in examining the possible causal linkage between statistical variables [44][45][46]. In the present work, SD ET was employed as the dependent variable, indicating IAV of ET. IAV of climatic variables (also indicated by the intra-annual standard deviation of the monthly variable), as well as yearly vegetation coverage (VC), were employed as independent variables. Specifically, the climatic variables included temperature (SD AT ), solar radiation (SD SR ), precipitation (SD PR ) and wind speed (SD WD ).

Analysis of Dominant Variable
According to the significance of the path coefficients, response patterns of IAV of ET to the variables could be recognized. Furthermore, the environmental variable with the highest absolute value of path coefficient (direct effect) was identified as the dominant variable of IAV of ET. To explore differences in the effects of climatic variables and VC throughout the continuous vegetation restoration, all of the analysis mentioned above were applied during the entire study period and the phases before and after the turning point.

Changes of IAV of ET
Changes in IAV of ET, measured by SD ET , were detected by both linear regression and piecewise regression in the Poyang Lake Basin ( Figure 2). In general, SD ET significantly (p < 0.05) declined from 26.2 mm in 1983 to 20.8 mm in 2014 over the entire basin, with a linear slope of −0.14 mm/year (Figure 2a). Similarly, significant (p < 0.05) decreases in SD ET were observed in all of the sub-basins, despite variations in magnitude (from −1.0 to −1.7 mm/year) (Figure 2b-f). However, the piecewise regression model performed better, with a higher R 2 (0.56) and lower RMSE (1.73 mm) and AIC (44.56), than the linear regression model (R 2 = 0.25, RMSE = 2.27 mm, and AIC = 56.78) in capturing the change of SD ET the across the Poyang Lake Basin ( Table 1). The turning point of the SD ET timeseries was found in the year of 1999 ( Figure 1). That is, the trend of IAV of ET remarkably changed over the Poyang Lake Basin during the study period. This result was consistent with that in each of the sub-basins. The differences of AIC between the piecewise and linear regression models were negative (with an average value of −10.91) across all of the sub-basins. The turning points of the SD ET trend were found in 1999 across the sub-basins except the Ganjiang Basin (in 2002) ( Figure 1). To minimize spatial heterogeneity, the year of 1999 was selected as the turning point for the subsequent analysis. Before 1999, the SD ET significantly (p < 0.05) decreased by −0.52 mm/year, while it slightly (p > 0.05) increased by 0.06 mm/year after 2000 over the whole study area (Figure 1a). For the sub-basins, likely, the trends of the SD ET ranged from −0.42 to −0.56 mm/year and from 0. significantly (p < 0.05) decreased by −0.52 mm/year, while it slightly (p > 0.05) increased by 0.06 mm/year after 2000 over the whole study area (Figure 1a). For the sub-basins, likely, the trends of the SDET ranged from −0.42 to −0.56 mm/year and from 0.3 to 1.3 mm/year before and after 2000, respectively (Figure 1b-f). Moreover, the spatial distribution of IAV of ET across the study area was investigated during different phases (    Table 1.   Figure 4 shows the changes in IAV of the climatic variables and VC over the Poyang Lake Basin from 1983 to 2014. In general, IAV of radiation and temperature first decreased and then increased around 2000, which was very consistent with IAV of ET. In contrast, IAV of wind speed exhibited gradual decline, while that of precipitation varied without a significant trend. For VC, a continuous increase was observed over the study area. It rapidly increased to 60% from 1990 to 2000 due to revegetation and maintained a slow increase since 2000. Relationships of IAV of ET to IAV of the climatic variables and VC were investigated by partial correlation analysis over the Poyang Lake Basin and its sub-basins. Table 2a provides the average values of the coefficients of partial correlation (r) during 1983-2014. Significantly (p < 0.01) positive correlations between SDET and SDSR dominated the entire study area. That is, IAV of ET increased with IAV of radiation. On the contrary, an increase in VC reduced IAV of ET in most of the basin, showing significant (p < 0.05) negative correlation coefficients between SDET and VC in the study area except the Raohe Basin. The two phases before and after the turning point were taken into account in the analysis. During the first phase (1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999), both SDSR and VC were significantly (p < 0.05) correlated with SDET in most of the sub-basins (Table 2b). In other words, an increasing IAV of radiation and VC would substantially enhance and weaken IAV of ET, respectively, in the Poyang Lake Basin. Notably, the absolute values of correlation coefficient for   Figure 4 shows the changes in IAV of the climatic variables and VC over the Poyang Lake Basin from 1983 to 2014. In general, IAV of radiation and temperature first decreased and then increased around 2000, which was very consistent with IAV of ET. In contrast, IAV of wind speed exhibited gradual decline, while that of precipitation varied without a significant trend. For VC, a continuous increase was observed over the study area. It rapidly increased to 60% from 1990 to 2000 due to revegetation and maintained a slow increase since 2000. Relationships of IAV of ET to IAV of the climatic variables and VC were investigated by partial correlation analysis over the Poyang Lake Basin and its sub-basins. Table 2a provides the average values of the coefficients of partial correlation (r) during 1983-2014. Significantly (p < 0.01) positive correlations between SD ET and SD SR dominated the entire study area. That is, IAV of ET increased with IAV of radiation. On the contrary, an increase in VC reduced IAV of ET in most of the basin, showing significant (p < 0.05) negative correlation coefficients between SD ET and VC in the study area except the Raohe Basin. The two phases before and after the turning point were taken into account in the analysis. During the first phase (1983-1999), both SD SR and VC were significantly (p < 0.05) correlated with SD ET in most of the sub-basins (Table 2b). In other words, an increasing IAV of radiation and VC would substantially enhance and weaken IAV of ET, respectively, in the Poyang Lake Basin. Notably, the absolute values of correlation coefficient for VC during 1983-1999 were larger than those during the whole study period, indicating a greater effect of VC during a rapid revegetation. Generally, stronger effects of IAV of temperature, indicated by the significant correlation coefficients between SD AT and SD ET , were observed in all of the sub-basins from 2000 to 2014 (Table 2c). In contrast, SD SR was only significantly (p < 0.05) associated with SD ET in the Raohe Basin, where the effect of IAV of temperature was not significant anymore. Moreover, the relationship between VC and IAV of ET remarkably decoupled across the basin, where the correlation between VC and SD ET was not statistically significant. After that, the direct effects of the climate variables and VC on IAV of ET were quantitatively explored across the Poyang Lake Basin during the different periods. The spatial distribution of the significant path coefficients (p < 0.05) demonstrated IAV of solar radiation dominated IAV of ET from 1983 to 2014, followed by VC and IAV air temperature ( Figure 5). IAV of radiation enhanced IAV of ET at most (about 85.3%) of study area (Figure 5c). Likewise, IAV of temperature also promoted IAV of ET in the central of the Poyang Lake Basin (24.2% of the area) (Figure 5b). A significant suppression of VC on IAV of ET was widely observed in the basin (47.7% of the area) (Figure 5a). IAV of precipitation significantly reduced IAV of ET as well, mainly in the Xiushui and Fuhe Basin, accounting  After that, the direct effects of the climate variables and VC on IAV of ET were quantitatively explored across the Poyang Lake Basin during the different periods. The spatial distribution of the significant path coefficients (p < 0.05) demonstrated IAV of solar radiation dominated IAV of ET from 1983 to 2014, followed by VC and IAV air temperature ( Figure 5). IAV of radiation enhanced IAV of ET at most (about 85.3%) of study area (Figure 5c). Likewise, IAV of temperature also promoted IAV of ET in the central of the Poyang Lake Basin (24.2% of the area) (Figure 5b). A significant suppression of VC on IAV of ET was widely observed in the basin (47.7% of the area) (Figure 5a). IAV of precipitation significantly reduced IAV of ET as well, mainly in the Xiushui and Fuhe Basin, accounting for 24.2% of the area (Figure 5e). The significant effects of IAV of wind speed on that of ET were only observed in parts of the Xinjiang and Ganjiang Basin (about 5.6% of the area) (Figure 5d). 4.3% of the area, respectively) were mainly located in the Xiushui and Xinjiang Basin be-fore 2000. After 2000, the areas controlled by IAV of both radiation and VC remarkably reduced to 34.1% and 4.1%, respectively (Figure 7). The significant effects of IAV of temperature on the ET variance, in constant, widely expanded, which was observed in 40.7% of the study area. The response of IAV of ET to IAV of precipitation also enhanced in western of Ganjiang and Fuhe Basin. The effect of IAV of wind speed was still limited during 2000-2014.  Figure 6 shows spatial distribution of the significant (p < 0.05) path coefficients of the climatic variables and VC before the turning point of 1999. Similarly to the entire study period, IAV of solar radiation mainly controlled IAV of ET across the Poyang Lake Basin (88.3% of the area), indicating that an increasing fluctuation of energy would greatly stimulate the seasonal deviation of ET. Significantly negative effects of VC were widespread in the basin as well, accounting for 35% of the area. However, IAV of temperature hardly impacted IAV of ET during 1983-1999. Only 7.3% of the study area occurred significant positive effects of temperature, mainly in central of the Xinjiang Basin and southern of the Ganjiang Basin. The areas regulated by IAV of precipitation and wind speed (7.4% and 4.3% of the area, respectively) were mainly located in the Xiushui and Xinjiang Basin before 2000. After 2000, the areas controlled by IAV of both radiation and VC remarkably reduced to 34.1% and 4.1%, respectively (Figure 7). The significant effects of IAV of temperature on the ET variance, in constant, widely expanded, which was observed in 40.7% of the study area. , measured by SD ET , over the Poyang Lake Basin from 1983 to 1999. SD ET , SD AT , SD SR , SD WD and SD PR denote the intra-annual standard deviation of monthly ET, temperature, radiation, wind speed and precipitation, respectively. The colored areas indicate where the path coefficient of the variable was statistically significant (p < 0.05) in the path analysis.

Dominators of IAV of ET
According to the significant path coefficients of the climatic variables and VC, the response patterns of IAV of ET were categorized during the different study periods. Spatial distributions of the first six dominant patterns are shown in Figure 8. The result showed that Pattern I (SD SR ) and Pattern V (SD SR + VC) were widespread, accounting for about 31.99% and 29.78% of the basin, respectively, from 1983 to 2014 (Figure 8a). The former one was mainly located in the Raohe and Ganjiang Basin, meanwhile the latter one widely occurred in the Ganjiang and Fuhe Basin. Similarly to the entire study period, Pattern I dominated the most area (51.99%) of the Poyang Lake basin, followed by Pattern V which was observed in 11.36% of the area during 1983-1999 (Figure 8b). After the turning point (i.e., the year of 1999), the patterns of IAV of ET in response to IAV of the variables gave a quite different picture (Figure 8c). Only 12.58% and 8.86% of the area showed Pattern I and Pattern V, respectively, mainly in the Raohe and Fuhe Basin as well. Pattern II (SD AT ) were widely observed in the north of the Poyang Lake Basin (about 30.21% of the study area).

Dominators of IAV of ET
According to the significant path coefficients of the climatic variables and VC, the response patterns of IAV of ET were categorized during the different study periods. Spatial distributions of the first six dominant patterns are shown in Figure 8. The result showed that Pattern Ⅰ (SDSR) and Pattern Ⅴ (SDSR + VC) were widespread, accounting for about 31.99% and 29.78% of the basin, respectively, from 1983 to 2014 (Figure 8a). The former one was mainly located in the Raohe and Ganjiang Basin, meanwhile the latter one widely occurred in the Ganjiang and Fuhe Basin. Similarly to the entire study period, Pattern Ⅰ dominated the most area (51.99%) of the Poyang Lake basin, followed by Pattern Ⅴ which was observed in 11.36% of the area during 1983-1999 (Figure 8b). After the turning point (i.e., the year of 1999), the patterns of IAV of ET in response to IAV of the variables gave a quite different picture (Figure 8c). Only 12.58% and 8.86% of the area showed Pattern Ⅰ and Pattern Ⅴ, respectively, mainly in the Raohe and Fuhe Basin as well. Pattern Ⅱ (SDAT) were widely observed in the north of the Poyang Lake Basin (about 30.21% of the study area). Furthermore, the dominator of IAV of ET was further identified for each of the pixels in the basin in light of the maximum absolute value of the significant path coefficient (Figure 9). During the entire study period, IAV of ET was dominated by IAV of radiation (i.e., SD SR ) in the most (about 72.0%) of the Poyang Lake Basin (Figure 9a). VC and IAV of temperature (SD AT ) followed that, accounting for 12.92% and 12.15% of the basin, respectively. Notably, the roles of both SD SR and VC were more important before the year of 1999 (Figure 9b). The area dominated by SD SR and VC increased to 77.82% and 18.0%. In contrast, the control of IAV of temperature to IAV of ET dramatically enhanced in north and central of the study area (35.14%) from 2000 to 2014 (Figure 9c), while the area dominated by IAV of radiation just decreased to 16.27% of the basin. However, there was a 39.46% of the area where IAV of ET was significantly regulated by none of the selected variables during this period.
tively. Notably, the roles of both SDSR and VC were more important before the year of 1999 (Figure 9b). The area dominated by SDSR and VC increased to 77.82% and 18.0%. In contrast, the control of IAV of temperature to IAV of ET dramatically enhanced in north and central of the study area (35.14%) from 2000 to 2014 (Figure 9c), while the area dominated by IAV of radiation just decreased to 16.27% of the basin. However, there was a 39.46% of the area where IAV of ET was significantly regulated by none of the selected variables during this period.   (1983-1999 and 2000-2014, respectively). SDET, SDAT, SDSR, SDWD, SDPR and VC denote the intraannual standard deviation of monthly ET, temperature, radiation, wind speed and precipitation, and vegetation coverage, respectively. The color of each pixel indicates the dominant variable, of which the absolute value of the path coefficient is statistically significant (p < 0.05) and highest in the path analysis. NS indicates none of the selected variables was significantly related to SDET during the period.

Discussion
Improving understanding of changes in IAV of ET and the underlying drivers is an essential step for better modeling of water cycle in response to global change. Previous studies have suggested that the IAV of ET varied considerably over time with climatic variables [6][7][8]. However, both abiotic and biotic variables affect annual and seasonal ET [9,12,18], and ignoring the role of landscape characteristics (e.g., vegetation coverage) can introduce large biases in the prediction of IAV of ET. Our present work shows that the IAV of ET exhibited significant contrasting trends during the past decades across the Poyang Lake Basin, China (Table 1 and Figure 2). Moreover, the proposed variables (i.e., the climatic variables and vegetation coverage) could well capture the changes in ET variance ( Table 2 and Figures 5-7). Furthermore, during the different phases, the dominators of IAV of ET were different (Figure 9). . Dominators of the intra-annual variability of evapotranspiration (ET), measured by SD ET , over the Poyang Lake Basin during (a) the entire study period  and (b,c) the two phases (1983-1999 and 2000-2014, respectively). SD ET , SD AT , SD SR , SD WD , SD PR and VC denote the intraannual standard deviation of monthly ET, temperature, radiation, wind speed and precipitation, and vegetation coverage, respectively. The color of each pixel indicates the dominant variable, of which the absolute value of the path coefficient is statistically significant (p < 0.05) and highest in the path analysis. NS indicates none of the selected variables was significantly related to SD ET during the period.

Discussion
Improving understanding of changes in IAV of ET and the underlying drivers is an essential step for better modeling of water cycle in response to global change. Previous studies have suggested that the IAV of ET varied considerably over time with climatic variables [6][7][8]. However, both abiotic and biotic variables affect annual and seasonal ET [9,12,18], and ignoring the role of landscape characteristics (e.g., vegetation coverage) can introduce large biases in the prediction of IAV of ET. Our present work shows that the IAV of ET exhibited significant contrasting trends during the past decades across the Poyang Lake Basin, China (Table 1 and Figure 2). Moreover, the proposed variables (i.e., the climatic variables and vegetation coverage) could well capture the changes in ET variance (Table 2 and Figures 5-7). Furthermore, during the different phases, the dominators of IAV of ET were different (Figure 9).

Roles of the Climatic Variables
It has been well reported that climate change controls terrestrial ET at regional and global scales [25,47]. The findings of this study also support this view. The study area, the Poyang Lake Basin, belongs to a typical subtropical humid climate under energy-limited conditions [8,19]. Hence, energy, indicated by downward shortwave radiation, dominated IAV of ET in light of not only the significant correlation between SD ET and SD SR but also its control areas (Table 2 and Figure 8). Moreover, the turning points for both IAV of ET and radiation were found in 1999 by the piecewise regression analysis. Additionally, a synchronous decrease in ET and radiation in summer (June-August) was observed from 1983 to 1999. All of these manifested the role of radiation on the ET variability. Evaporation demand, which is mainly driven by temperature [47], is another dominator of ET in humid basins. An increase in air temperature would enhance vapor pressure deficit (VPD) and thereby potential evapotranspiration (ET p ). If soil moisture is sufficient, ET would be stimulated with an increasing ET p [48]. It has been reported that the Poyang Lake Basin has experienced an obvious warming since 1998 [49]. In our analysis, we found that SD AT also rose after 1998 because temperature increased in summer whereas decreased in winter. This was high consistent with the turning point of SD ET at the year of 1999. Among the selected variables, only SD AT and VC significantly changed after 1999. However, the correlation between IAV of ET and VC seemed to be severely weakened due to the slow-down of the revegetation (Figure 8). Hence, the increase in SD ET might be attributed to the enhanced fluctuation of seasonal temperature with asymmetric warming after 1999. Compared to energy and evaporation demand, the areas dominated by water supply (i.e., precipitation) were limited, mainly in the Ganjiang, Fuhe and Xiushui Basin. It may be associated with the relatively sufficient water supplement. As previously reported, ET variance is more sensitive to precipitation fluctuation in water-limited environments (e.g., arid climates) than other conditions (e.g., humid climates) [50,51]. Although wind speed plays an important role in evaporative demand, it slightly impacted IAV of ET over the study area, which is agreed to the previous study [12,19].

Roles of Vegetation Coverage
In addition to the climatic variables, increased controls of vegetation on terrestrial ET are observed over the entire world [24]. The Poyang Lake Basin experienced a dramatic reforestation from 1980s to 1990s, with a two-fold increase of the total vegetation coverage [18]. This restoration period was just consistent with that when SD ET significantly decreased, showing a tight coupling between VC and IAV of ET. After 2000, the growth of VC slowed down, and meanwhile SD ET started to exhibit an adverse trend across the basin. It indicates that a rapid restoration would promote the seasonal fluctuation of ET to become more stable (i.e., a decrease in SD ET ). Two potential causes accounted for the decline in SD ET with the increasing vegetation cover. On the one hand, an expansion of vegetation leaf area could diminish the amount of strong solar radiation on soil surface in summer, resulting in a decrease in soil evaporation [18,52]. On the other hand, high ET in a well vegetated area would take heat away and raise air relative humidity [53]. Hence, it could induce a decrease in air temperature and VPD to some extent, and further a decline of evaporation demand, especially in the humid regions [54]. That is, the cold island effect of afforestation would further weaken the terrestrial ET. Both of the pathways mentioned above would lead to lower summer ET and consequently minor IAV of ET. Compared to a sparse vegetation, restoration may produce a systematic shift between water cycle components [16,18]. Overall, an increase in VC could reduce the fluctuation of ET in a year cycle so that maintain the water content of soil. This is beneficial to enhance the stability of water cycle at the basin scale, as well as resistance of terrestrial ecosystem to climatic extremes. Notably, this positive effect of afforestation may just occur in humid basins, whereas it seems to be reversed over arid regions due to the strong evaporation demand [9]. Moreover, the decline in SD ET in response to revegetation might be observed in areas where VC dramatically increased. As shown in this study, despite a continuous but slight increase, SD ET seemed to be insensitive to VC after that was over 60%, implying a limited role on IAV of ET [12].

Uncertainties and Further Study
In our present work, a few limitations should be acknowledged. First, although both the AVHRR and MODIS ET data have been widely evaluated and applied, there are still uncertainties in their forcing data, retrieval algorithm and parametrized scheme [25]. One of the most noteworthy problems is the static land cover data adopted in the ET estimation, which were used to set the physiological parameters (e.g., potential stomatal conductance) and constraints on stomatal conductance (e.g., minimum air temperature) [20,23]. To this point, a previous study investigated the impacts of land cover change on the data quantity of AVHRR and MODIS ET during the past decades over the Poyang Lake Basin [18]. It showed that the AVHRR and MODIS ET data worked properly, because the conversions between the land cover types did not significantly alert the key parameters or marginally impacted to the variation in ET in the study area [20,55]. However, given land cover changes, the constant parametrization of land classification inevitably induced errors in the ET calculation. It is necessary to develop a reliable scheme in which land cover types could be well identified for ET inversion in future study. Second, the direct effect of each of the regulators on IAV of ET was examined by path analysis. Despite statistical independence among the four climatic variables over the most areas of the Poyang Lake Basin, the collinearity might be inevitable due to the temperate monsoon climate [56], which is another source of uncertainty in this study. Third, the intra-annual standard deviation of the climatic variables could not fully represent influences of extreme events (e.g., seasonal drought [57]), weakening the explanatory power of the abiotic regulator of ET. Forth, the result showed that there were lots of areas where all of the selected variables could not account for the change of SD ET during 2000-2014 ( Figure 8). This might be associated with the mountain topography at the upstream region of the sub-basins, which is needed to be further investigated. Finally, interferences of human activity, e.g., water management over the cropland of the study area, also leaded to uncertainties of the present results. Nevertheless, despite the uncertainties mentioned above, the findings of this work explicitly demonstrated that both the climatic variables and vegetation coverage play key roles in regulating IAV of ET, and the role of the variables varied in space and time. Hence, we suggest that the hydrological response and feedback should be interpreted by not only abiotic but also biotic variables (e.g., vegetation coverage) to avoid omissions of useful signals.

Conclusions
This study presented changes of intra-annual variability (IAV) of ET and the responses to climatic variables and vegetation coverage over the Poyang Lake Basin in China from 1983 to 2014. Our findings demonstrate that IAV of ET showed contrary trends across the study area during the past decades. It statistically significantly decreased with a linear slope of −0.52 mm/year before 2000, and then increased (slope = 0.06 mm/year) despite no significant. The changes in IAV of ET were largely attributed to the climatic variables and vegetation coverage by statistical analysis. Generally, IAV of solar radiation (a proxy of energy) and air temperature (a proxy of evaporation demand) dominated the changes of IAV of ET over 77.82% and 35.14% of the basin, respectively, at the two phases. Meanwhile, the increase in vegetation coverage through the rapid restoration significantly reduced IAV of ET across the study area (about 35% of the area). Overall, we argue that both abiotic and biotic variables should be taken into account in understanding of changes in IAV of ET at the basin scale.

Conflicts of Interest:
The authors declare no conflict of interest.