A Revised Temporal Scaling Method to Yield Better ET Estimates at a Regional Scale

This study presents a revised temporal scaling method based on a detection algorithm for the temporal stability of the evaporative fraction (EF) to estimate total daytime evapotranspiration (ET) at a regional scale. The study area is located in the Heihe River Basin, which is the second largest inland river basin in China. The remote sensing data and field observations used in this study were obtained from the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project. The half-hourly EF values (EFEC) calculated using meteorological observations from an eddy covariance (EC) system and an automatic meteorological station (AMS) represented the diurnal pattern of the EF across the majority of the study area. The remotely sensed instantaneous midday EF (EFASTER), which indicates the spatial distribution of the midday EF over the entire study area, was calculated from an Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image. The temporal stability of EFEC was examined using a detection algorithm. Intervals with inconsistent EFEC values were distinguished from those OPEN ACCESS Remote Sens. 2015, 7 6434 with consistent EFEC values; the total daytime ET (from 9:00 to 19:00) within these interval types was integrated separately. Validation of the total daytime ET at the satellite pixel scale was conducted using measurements from17 EC towers. Using the detection algorithm for the temporal stability of the EF and dynamic adjustment, the revised temporal scaling method resulted in a root-mean-square error (RMSE) of 0.54 (mm·d−1), a mean relative error (MRE) of 7.26% and a correlation coefficient (Corr.) of 0.81; all of these values were superior to those of the two other methods (i.e., the constant EF and variable EF methods). The revised method easily extends to other areas and exhibits a superior performance in flat and regularly-irrigated farmlands at the regional scale.


Introduction
Evapotranspiration (ET) plays an important role in hydrology, meteorology and agriculture by controlling the water cycle and energy transport within the biosphere, atmosphere and hydrosphere [1][2][3].Because the hydrological cycle is strongly affected by crop water consumption in semiarid agricultural regions, precise ET estimates are important for saving water through efficient irrigation practices [4,5].Remotely sensed data can provide instantaneous two-dimensional spatial information on ground surface properties, such as albedo, vegetation fraction, land surface temperature and net radiation (Rn).These properties are important variables for remotely sensed ET estimates [6][7][8].However, in many ET-related disciplines, total daily, weekly and monthly ET readings, rather than instantaneous values, are required at multiple scales.Therefore, to convert instantaneous spatial ET values to time-integrated values, temporal scaling is required.Current ET estimation models, such as the surface energy balance system (SEBS) [9], the surface energy balance algorithm for land (SEBAL) [10], and two-source models for evapotranspiration, such as TTME [11] and TSEBPS [12], primarily depend on observations during satellite overpasses.Several studies have been conducted to convert the instantaneous latent heat flux (λE) to the integrated ET, either using a numerical model that simulates the diurnal surface temperature under realistic conditions or using assumptions about the relationship between instantaneous ET and daily values [3].Sine function [13] and constant evaporative fraction (EF) [14] methods are commonly used for scaling instantaneous λE to total integrated ET [3,15].Jackson et al. [13] used a sine function to relate the ratio of instantaneous ET to daily ET to the diurnal trend of solar irradiance.Sugita and Brutsaert [14] proposed the constant EF method, which assumes that the EF is sufficiently constant during the daytime, to determine the regional daily ET.The sine function provides a good approximation of the change in diurnal solar irradiation from approximately sunrise to sunset on sunny days.However, this method has limited accuracy on cloudy days due to its empirical nature [15].The variability or conservation of the EF on a given individual day was affected by complicated combinations of factors, including weather conditions, soil moisture, topography, biophysical conditions, cloudiness and the advection of moisture and heat [16].Allen et al. [16] proposed a constant reference EF method to capture the impacts of advection and changing wind and humidity conditions in ET temporal scaling.Because it considers orographic effects, this method is more applicable over terrain areas.Tomás et al. [17] used a triangle method to estimate instantaneous EF during satellite overpass times over large semi-arid areas with uniform landscapes.Several studies based on in situ measurements and analyses that consist of land surface process modeling have indicated that the EF curve is concave upward and is relatively stable during the daytime under dry conditions [18][19][20][21][22]. Suleiman and Crago reported that the EF increases with the amount of vegetation, soil moisture, and air dryness [23].Hoedjes et al. [19] showed that the EF curve is markedly concave upward under wet conditions and developed a new heuristic parameterization based on routine daily meteorological observations and remotely sensed estimates of instantaneous EF to consider the dependence on soil and vegetation conditions.However, the EF values in the afternoon are unstable relative to those in the morning [22].The variable EF method [19] better approximates the true EF values under wet conditions in the morning, although the study did not address the large fluctuations in the afternoon.
We studied the Zhangye Oasis of the Heihe River Basin, which is the second largest inland river basin in the arid region of Northwest China [24].The main land use type in our study area is a typical irrigated farmland.The primary crop is seed corn.According to Yang [25], who conducted a study in the same area, the assumption that EF can be treated as constant during the daytime (so-called EF self-preservation) no longer holds over corn growing seasons because the ET diurnal pattern is primarily influenced by stomatal regulation.The constant EF method and the variable EF method have particular limitations in irrigated farmlands of arid regions because of the special land cover and climate conditions.A detection algorithm conducted in this revised ET temporal scaling method was applied to examine the temporal stability of the EF values, which were calculated from eddy covariance (EC) and automatic meteorological station (AMS) measurements.The primary objectives of this study were (1) to consider the temporal and spatial heterogeneity of the EF in the temporal scaling method; (2) to compare the results estimated using the revised method with those estimated using the other two temporal scaling methods (i.e., the constant EF method and the variable EF method); (3) to assess the applicability of the revised method to farmlands in arid regions; and (4) to obtain better regional ET products.

Data
The field observations and remote sensing images used in this study were derived from the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project.HiWATER is a comprehensive eco-hydrological experiment under the framework of the Heihe Plan and is based on the diverse needs of interdisciplinary research and existing observational infrastructures in the basin [26].The intensive observation period of the fieldwork began in May 2012 in the middle reaches of the river.

Field Observations
Within the Yingke and Daman irrigation districts of the Zhangye Oasis, an observation matrix that consists of 17 stations was built.Each station was equipped with an EC tower and an AMS.The land surface of the observation area is heterogeneous and is dominated by seed corn, corn planted with spring wheat, vegetables, orchards and residential areas (see Figure 1).The observations used in this study include humidity (sensor type: Avalon AV-14TH), 4 radiation components (total solar radiation, reflective shortwave radiation, land surface long wave radiation, and atmospheric long wave radiation; sensor type: Kipp & Zonen CNR 4, Eppley PIR & PSP), land surface temperature (sensor type: Avalon IRTC3), soil temperature (sensor type: Avalon AV-10T) and moisture (sensor type: Campbell CS616) profiles and turbulent fluxes (sensor type: CSAT3 + LI7500A) [27].
The raw high-frequency (10-Hz) data collected at the EC towers were successively processed through de-spiking [28], tilt correction [29], sonic virtual temperature corrections [30], time-lag calculations, frequency response corrections [31] and Webb, Pearman and Leuning (WPL) correlations [32].These processes were completed using the EdiRe software developed by Edinburgh University, UK.The half-hour-averaged λE and sensible heat flux (H) values were calculated from the post-processed EC data [33].The records with a friction velocity of less than 0.1 (m•s −1 ) were not used in the energy balance closure analysis or ET validation.The energy balance closure from 1 July to 15 August reached 90%.

Satellite Images
Meteorological data from the HiWATER middle reaches experiment had the highest integrity from June to August 2012.To maintain the soil moisture, seedlings were covered with plastic film.This plastic film can introduce uncertainty into the ET observations by blocking a portion of the vapor, which is not observable by the EC systems.After July, the corn grew up, the plastic film was broken, and the field ridges were closely spaced.The remaining plastic film residue in the soil did not have a large impact on the ET observations.Thus, the images used in our study were selected from July to August.During this period, three Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) L1B images acquired on 10 July, 2 August and 11 August of 2012 were used for estimating the ET from remote sensing data.The images were pre-processed through radiative, atmospheric and geometric corrections with Environment for Visualizing Images (ENVI).Atmospheric correction was conducted using ENVI Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH).First, the processed images were used to estimate instantaneous ground surface properties, including albedo, emissivity and vegetation fraction, at ASTER overpasses.Next, instantaneous surface heat fluxes, including λE, ground heat flux (G) and available energy (A), were estimated based on the ground surface properties.Then, the instantaneous EF and Bowen ratio during ASTER overpass times were calculated.The specific calculations are introduced in Section 3.4.

Constant EF Method (cEF)
The cEF method assumes that the daily EF remains relatively constant despite particular daily variations.This 'self-preservation' of the EF has allowed for the use of midday EF values for converting instantaneous remotely sensed λE to total daily values [18,19,34].From ASTER L1B images, the instantaneous midday EF (EFASTER) was estimated during ASTER overpass times as follows: Then, the EFASTER values were multiplied by the total daytime Rn (from 9:00 to 19:00) minus the total daytime G (from 9:00 to 19:00) to estimate the total daytime ET (from 9:00 to 19:00).Remotely sensed λEASTER values were estimated using a revised surface resistance model [35].GASTER values were estimated using the same parameterization proposed in SEBS [9].Other remotely sensed land surface variables and half-hourly Rni and Gi were the same in the three methods.The specific calculations are introduced in Section 3.4.

Variable EF Method (vEF)
The 'self-preservation' assumption of the EF is not valid under wet conditions [18,19].During the daytime, the EF of bare soil is more stable than that of regularly-irrigated farmlands whose EF curve is concave upward.To distinguish farmlands under wet conditions from bare soil under dry conditions, Hoedjes et al. [19] used a dynamic coefficient to adjust the EF values, which better approximate the true value under wet conditions in the morning [15].The adjustment factor (r) can be defined as follows: where sim i EF is the simulated EF of the i th time step (half hour), which is a function of the incoming solar radiation (S ↓ , W•m −2 ) and the relative humidity (RH, %): The Bowen ratio during an ASTER overpass (βASTER) was estimated as follows: Under dry conditions, βASTER exceeds 1.5; thus, EFASTER was used for temporal scaling.When βASTER did not exceed 1.5, EFASTER was optimized by the ri value of the corresponding time step: (5)

Revised Variable EF Method (vEFr)
The adjusted EF curve, which was simulated using the vEF method, displayed a concave upward shape that better approximates the true value under wet conditions in the morning.However, the EF values in the afternoon (after 14:00) are unstable relative to those in the morning [22].The vEF method does not address the large afternoon fluctuations.To address the high temporal heterogeneity of the EF in the afternoon under wet conditions, a detection algorithm for the temporal stability of the EF was conducted in the revised temporal scaling method for irrigated fields.The half-hourly EF (EF EC ) values calculated from the meteorological observations from tower AMS/EC12were used as references to represent the EF diurnal pattern across most of the study area: where the subscript i is the EF time step (half hour), and A is the available energy, which is equal to the Rn observed by the net radiometer on tower AMS12 minus G.In addition, the half-hourly λEi values were obtained from tower EC12.The half-hourly Gi values were calculated using soil moisture and temperature profiles observed by tower AMS12 based on a temperature prediction-correction method suggested by Yang and Wang [36].The logic behind the detection algorithm for determining the temporal stability of the EF is shown in Figure 2. The detection algorithm used a moving window to search for a continuous time interval with relatively stable EF EC values.When the absolute value of the difference between EC i EF and the mean value of the stable EF EC values was lower than a particular threshold, EC i EF was considered relatively stable.However, determining the mean value of the stable EF EC values and the threshold was problematic.The atmosphere is generally more stable in the morning than in the afternoon.Thus, the EF values in the morning (from approximately 9:00 to 12:00) are steady relative to those in the afternoon [22].A moving window (2.5 h) was used to scan the EF EC values from 9:00 to 14:00.The mean values (u) of the EF EC and their standard deviations (σ) within each window were calculated.The window with the minimum σ (σ_min) was considered as a continuous time interval with less unstable EF values.The mean value for the window with the minimum σ (uσ_min) was considered as the mean of the stable EF EC values, and the σ_min provided the threshold.Then, the moving window scanned all EF EC values from 9:00 to 19:00 to identify the stable EF EC .Each Therefore, the total daytime ET (from 9:00 to 19:00) was calculated as follows: where t represents a half hour (in seconds); L is the latent heat of vaporization (2.45 × 10 6 J•kg −1 ); and ρw is the density of water (1000 kg•m −3 ).Within the interval that EF EC remains stable, the midday EFASTER was optimized via r in Equation ( 5) under wet conditions (βASTER below 1.5) while keeping the original value under dry conditions (βASTER exceeding 1.5).When EFi EC exhibited large variations, EFASTER was replaced by EFi EC in the corresponding time step during the temporal scaling process.
Although EF EC values are point measurements that ignore the spatial heterogeneity of the EF, the error induced by replacing EFASTER with EFi EC should be smaller than that produce by ignoring the temporal heterogeneity for conditions in which EF greatly fluctuates.This point is discussed in Section 5.2.The input data required in these three methods are shown in Table 1, and the logic behind the vEFr method is shown in Figure 3.A validation of the land surface ET at the satellite pixel scale was conducted using measurements collected from 17 EC systems.An innovative validation method for remotely sensed ET using EC measurements [37] was used in our study.The observations from the EC system represent a weighted average over a range that primarily contributes to the measurement, namely the source area, which depends on the measurement height, wind direction/velocity, atmospheric stability and the underlying surface conditions.It has been suggested that all pixels within the EC system source area should be considered as validation pixels [35,37].By determining the footprint based on the model of Kormann and Meixner [38], the source area and the footprint-relative weights (relative contribution) during the passage times of the three images can be acquired.The remotely sensed ET estimates were overlapped using the relative contribution of EC systems.The weighted average remotely sensed ET (ETweighted) was computed as follows: where n is the number of pixels within the source area; wp and ETp are the relative weight and ET estimate, respectively, of each pixel within the source area.

Remotely Sensed Estimation for Land Surface Variables
The land surface variables of albedo and emissivity cannot be obtained at times other than during the ASTER overpass.Thus, the albedo and emissivity are assumed to be constant during the daytime.
The assumption that the albedo is constant throughout the day may be unrealistic [39].However, this assumption is not critical for accurately retrieving the available energy [19].The surface albedo (α) was estimated from the reflectance in bands 1, 3, 5, 6, 8 and 9 of the ASTER L1B images [40] The land surface emissivity was calculated using the normalized difference vegetation index (NDVI) threshold method [41,42] as follows: Where εv and εs are the emissivity of the full canopy (0.986) and bare soil (0.97215), respectively [42]; Rv and Rs are the temperature ratios of the full canopy and bare soil, respectively, which are estimated according to the fractional canopy coverage (fc) [42]: and dε is estimated as follows [41]: Here, fc is calculated using the NDVI [43]: where NDVImax and NDVImin are 0.7 and 0.05, respectively.The land surface temperature (LST) in the i th time step for one pixel is composed of the surface temperatures of the bare soil (Ts) and full canopies (Tv) in the same time step [44].These two limiting cases were spatially interpolated using fc: where the subscript i represents the i th time step (half hour).
According to the principles of microclimatology, the LST varies sinusoidally during the day.The Ts and Tv curves were observed by infrared radiometers installed on the AMS towers located in bare soil and full canopies (see Figure 4).The surface temperatures of the canopy and bare soil in the i th time step were simulated as follows: where and av are the mean and amplitude of the surface temperature curves of the canopy, respectively; and as are the mean and amplitude of the surface temperature curves of the bare soil, respectively; t* is the time (in hours on a 24-hour basis) after sunrise; and ω is the cycle of the surface temperature curves (which should be equivalent to 24 h).However, a half cycle is less than 12 h because the overnight land surface temperature does not have a sinusoidal curve.In addition, the Tv and Ts estimates match the observations during the daytime better when the half cycle is set to 10 h.φ is the initial phase of the temperature curve and is related to the local latitude (l) and the solar declination (sd), which is defined as follow: Tower AMS/EC04 was placed on a bare field next to a settlement, and tower AMS/EC12 was placed on farmland with densely planted seed corn.The Tv and Ts curves were fitted using the observations of infrared radiometers installed on towers AMS04 and AMS12, respectively.The parameters and root-mean-square error (RMSE) of the Tv and Ts estimates are shown in Table 2.  Half-hourly Rni values were calculated using instantaneously remotely sensed α and ε, half-hourly LSTi values that were spatially interpolated using fc and half-hourly observations of total solar radiation (RSi) and atmospheric long wave radiation (RSi): Then, half-hourly Gi values were estimated using the same parameterization proposed in SEBS [9]: The ratios between the soil heat flux and the net radiation for the full vegetation canopy ( c  ) and bare soil ( s  ) were 0.05 [45] and 0.315 [46], respectively.

Results and Validation
According to the results of the detection algorithm for the temporal stability of the EF, the EF EC values remained stable over particular daytime intervals and fluctuated before sunset and after sunrise (see Table 3).The EF values (see Figure 5) remained relatively stable from 10:00 to 15:00 with some variability during the early morning and late afternoon based on the surface and atmospheric conditions.The soil moisture on 10 July was higher than that on the other two days because most of the fields had recently been irrigated.Compared with the other two days, the EF on 10 July was less stable.However, the soil moisture on 11 August was the lowest of the three days.Only the field near towerEC14 had been irrigated after 31 July.Thus, the EF from tower EC14 was higher than that from the other stations.The solid azure line in Figure 5 shows the optimized EF values calculated using the vEF method.The optimized EF values were approximately the true values for the corresponding EF EC values in the morning.However, the vEF method did not simulate the significant increase in the EF after 15:00.Among these three ET temporal scaling methods, the cEF method resulted in the lowest ET estimates (see Figure 6).The mean ET values from the cEF method on 10 July, 2 August and 11 August of 2012 were 4.89, 4.01 and 2.89 (mm•d −1 ), respectively.The transpiration on 11 August was the lowest of the three days because the corn tassels had been cut off.The distributions of the three ET results were similar.Villages, areas with regular shapes, roads, and straight lines in the figures displayed low ET values.
From the validation results (see Tables 4-6), the estimates from the vEFr method with a RMSE of (0.54 mm•d −1 ), a mean relative error (MRE) of 7.26% and a correlation coefficient (Corr.) of 0.81 were most similar to the EC measurements compared with the other two temporal scaling methods (see Figure 7 and Table 7).Relative to the vEFr method, the cEF method underestimated the total daytime ET.Because the EF curve is typically concave upward, the use of a constant EF value for an ASTER overpass underestimates the total daytime ET.The vEF method accounted for the spatial heterogeneity of the EF and minimized this disadvantage of the cEF method to some extent.However, the vEF method did not capture the EF fluctuation in the early morning or in the late afternoon, which affected the results of the vEF method relative to the vEFr method.The estimates from both the cEF and vEF methods have lower Corr.with the EC measurements than the vEFr method.The three methods underestimated the ET in areas with high fc (see Tables 4-6), such as the fields of EC12, EC13 and EC14, compared with areas having low fc.This difference is discussed in Section 5.4.

Discussion
The variability or conservation of the EF on a given day is affected by a complex combination of factors, such as the weather, soil moisture, topography, biophysical conditions, cloud cover, moisture advection and heat advection.Thus, EF varies depending on various underlying surface and meteorological conditions.Four particular problems caused errors in the ET estimation during the temporal scaling process.

Spatial Heterogeneity of the EF
Two primary factors caused the spatial heterogeneity of the EF in our study area.First, the EF of farmlands differs from that of bare soil.The EF curves for regularly-irrigated farmlands are clearly concave upward, whereas those of bare soil during the daytime are more stable.The vEFr method uses Hoedjes's adjustment factor on farmlands under wet conditions when EFi EC is classified as stable.Second, particular abnormal meteorological and hydrological conditions in small areas or heterogeneous properties of the shallow subsurface also caused the spatial heterogeneity of the EF.The soil moisture conditions were not uniform because the fields could not all be irrigated simultaneously.Thus, the EF values were highly spatially heterogeneous when particular fields were irrigated but others were not.In our study, the midday EF value for each pixel used in the temporal scaling was calculated based on the remotely sensed EFASTER.However, the adjustment factor for the EF optimization and EF EC for the temporal stability detection of the EF were calculated using in situ observations at the point scale.The second type of spatial heterogeneity was not considered in our study.The underestimation near tower EC03 on 10 July and 11 August represents a good example of this variation.These fields were irrigated two or three days before the ASTER overpass, i.e., later than the irrigation of other fields; thus, their soil moisture was much higher than that of the other fields during the ASTER overpass.In addition, heterogeneous crop types can also lead to spatial heterogeneity in the EF.The underestimations near towers EC01, EC16 and EC17 are good examples of this variation.This spatial heterogeneity was addressed to an extent in this study.Only the first type of spatial heterogeneity was addressed using a dynamic coefficient that slightly adjusts the EFASTER value of farmlands.However, the second type of spatial heterogeneity requires extensive meteorological observations; the instruments required to make these observations are expensive.Thus, it was difficult to apply such methodology beyond the HiWATER context.

Temporal Heterogeneity of the EF
Solar radiation increases after 7:00 and decreases after14:00; however, ET increases after 8:00 and decreases after 15:00 (see Figure 8).Thus, the EF produces a concave-upward curve with a linear decrease in the early morning and a sharp increase in the late afternoon, and the ET can be underestimated if a constant midday EF value during the satellite overpass is used in the temporal scaling process.The cEF method clearly underestimated the ET relative to the other two methods (see Table 7).The ETWatch model, which uses the cEF method for temporal scaling, has been shown to have a RMSE of 0.99 (mm•d −1 ) [37,47].This problem was addressed by the vEF method to some extent, although the large fluctuations in the EF in the late afternoon were not identified with this method.If the temporal heterogeneity of the EF is not accounted for by temporal scaling, then the ET can be underestimated by 0.9 (mm•d −1 ) on a summer day.This amount is equivalent to approximately16% of the total daytime ET in farmland in arid regions on a summer day.Using our detection algorithm for the temporal stability of the EF, the ET temporal scaling method yielded more accurate values (30 m resolution), with an RMSE of 0.54 (mm•d −1 ), MRE of 7.26% and Corr. of 0.81, which are improvements over the other two methods.EF EC should represent the properties of the EF across most of the area because the spatial heterogeneity caused by particularly abnormal small-scale meteorological and hydrological states was not considered.Thus, the ET should be estimated several days after irrigation, when the soil moisture becomes homogeneous across the area.

The Uncertainty in Input
Instrumentation errors and the estimation errors for the instantaneous surface heat fluxes were analyzed in this study.The uncertainty in the EC measurements was approximately 16% [48].The observation error of the net radiation (CNR 1) was approximately 20 (W•m −2 ) [49].These errors were considered when calculating the EF EC values that were used as reference EF values in the detection algorithm for the temporal stability of the EF.The remotely sensed instantaneous Rn, G and λE values had errors of 2-3 (W•m −2 ), 9-12 (W•m −2 ), and 7-11 (W•m −2 ), respectively [35].These errors were considered in the calculation of EFASTER and βASTER, which were used as the midday constant EF value and the determination index for wet conditions, respectively.

The Saturation of fc
fc is another factor that affects the ET estimates.After July, seed corn in the Zhangye Oasis usually reaches 1.6 m in height.The crop leaves are large and overlapping, and the field ridges are closely spaced.Transpiration significantly contributes to the total ET at this time, and the percentage of transpiration to the total ET reaches 85.69% ± 5.86% [50].According to the equation for fc (Equation ( 12)) [42], when an NDVI value surpasses the maximum value, fc becomes 100%.However, when fc becomes saturated (close to 100%), the LAI continues to increase (see Figure 9).fc cannot account for the effects of stacked leaves.Moreover, fc was not only used in the estimation of the instantaneous emissivity [41,42], G [9], and λE [35] at the ASTER overpasses but also used in the derivation of the LST.This saturation effect of fc on the ET estimation requires further analysis.The ET should be underestimated in densely planted farmlands, where as EFASTER should be underestimated during an ASTER overpass based on these errors.Consequently, the three methods underestimated the ET in the fields with high fc, particularly in the fields near towers EC06, EC12, EC13 and EC14.The ET is underestimated by 0.3 (mm•d −1 ) when fc is nearly 100%.

The Transferability and Limitations of vEFr
Humidity, four radiation components, LST, and half-hourly λE and G are required for the vEFr method.According to the input data of the vEFr method, one ordinary AMS and one automatic infrared thermometer are required to use the vEFr method in other areas.If conditions permit, it is more preferable to utilize one or more EC systems.The reason is that the eddy covariance method [51] is regarded as the most reliable method for calculating half-hourly λE values [48], which are used to calculate the reference EF EC values used in the detection algorithm for the temporal stability of the EF and to validate the final ET estimate.If the EC system is unavailable under limited conditions, the vertical gradient method [52] using two layers of AMS observations is an alternative to the eddy covariance method [51] for calculating half-hourly λE values and the total daytime ET.
The vEFr method has two limitations in application.First, the orographic effect was not considered in the vEFr method, so the impacts of advection, changing wind and humidity conditions could not be captured in ET temporal scaling.Thus, the vEFr method should be applied in flat farmlands.Second, the EF EC values that were used as reference EF values in the detection algorithm for the temporal stability of the EF was a point observation, so spatial heterogeneity in EF caused by heterogeneous soil moisture conditions and crop types would bring about larger errors in estimation.Thus, the vEFr method was more applicable for regularly-irrigated farmlands with relatively homogeneous crops at the regional scale.In addition, due to the second limitation of the vEFr method, remotely sensed images with a spatial resolution of about 30 m (i.e., Landsat 8 OLI/TIRS image) were more preferable to implement the vEFr method.

Conclusions
In regularly-irrigated farmland, the instability of the EF in the afternoon is the main cause of the ET underestimation by the cEF and vEF temporal scaling methods.The revised vEFr temporal scaling method yielded better ET estimates because it considers the temporal heterogeneity of the EF in the afternoon.The vEFr method easily extends to other areas and exhibits a superior performance in flat and regularly-irrigated farmlands at the regional scale.

Figure 1 .
Figure 1.Land cover and the meteorological instrumentation configuration in the study area.

Figure 2 .
Figure 2. Flowchart of the detection algorithm.

Figure 3 .
Figure 3. Flowchart of the revised temporal scaling method.

Figure 5 .
Figure 5.The EF values calculated from EC and AMS observations on (a) 10 July, (b) 2 August, and (c) 11 August of 2012.Note: the EF values in azure were optimized using the vEF method for EC05 on 10th July and EC08 on the other two days.

Figure 7 .
Figure 7.Comparison of the EC measurements and ETdt based on the estimates from the (a) cEF, (b) vEF and (c) vEFr methods.

Figure 8 .
Figure 8.Comparison of λE from observations of tower EC12 and the total solar radiation from AMS12 on (a) 10 July, (b) 2 August and (c) 11 August of 2012.

Figure 9 .
Figure 9.Comparison of the LAI and fc at 17 AMS/EC stations on (a) 10 July, (b) 2 August and (c) 11 August of 2012.

Table 1 .
Input for the three methods.
Note: The remotely sensed EF ASTER , β ASTER , land surface emissivity, albedo and f c used in the three temporal scaling methods were estimated from ASTER L1B images.The time period for the total solar radiation, atmospheric long wave radiation, adjustment factor, and EF EC is 9:00 to 19:00.

Table 2 .
Regression parameters of the Ts and Tv curves and the root-mean-square errors (RMSEs).

Table 3 .
Stability of EF EC on 10 July, 2 August, and 11 August of 2012.

Table 7 .
Validation statistics for the three ET temporal scaling methods.