Remote Sensing Daily Evaporative Fraction Parameterization Scheme Driven by Day–night Differences in Surface Parameters: Improvement and Validation

In a previous study, a daily evaporative fraction (EF) parameterization scheme was derived based on day–night differences in surface temperature, air temperature, and net radiation. Considering the advantage that incoming solar radiation can be readily retrieved from remotely sensed data in comparison with surface net radiation, this study simplified the daily EF parameterization scheme using incoming solar radiation as an input. Daily EF estimates from the simplified scheme were nearly equivalent to the results from the original scheme. In situ measurements from six Ameriflux sites with different land covers were used to validate the new simplified EF parameterization scheme. Results showed that daily EF estimates for clear skies were consistent with the in situ EF corrected OPEN ACCESS 4370 by the residual energy method, showing a coefficient of determination of 0.586 and a root mean square error of 0.152. Similar results were also obtained for partly clear sky conditions. The non-closure of the measured energy and heat fluxes and the uncertainty in determining fractional vegetation cover were likely to cause discrepancies in estimated daily EF and measured counterparts. The daily EF estimates of different land covers indicate that the constant coefficients in the simplified EF parameterization scheme are not strongly site-specific.


Introduction
Studies on surface energy and water budgets require accurate evapotranspiration (ET) information at large spatial scales [1,2].Satellite-based data can be used for the ET estimation [3,4].Physical ET models based on remotely sensed data (e.g., one-source models and two-source models [5][6][7][8][9][10][11]) require ground-based measurements.This hinders the application of these models particularly to ungauged areas with few ground-based measurements.Therefore, it is advantageous to develop ET models that completely depend on surface parameters/variables retrieved from remote sensing.Surface parameters/variables that can be extracted from remotely sensed data mainly include surface temperature, surface emissivity, albedo, vegetation indices, leaf area index, soil moisture, net radiation, solar radiation, and air temperature [12][13][14][15][16][17][18][19][20][21][22].Some simple methods can estimate surface energy fluxes by directly using these remotely sensed surface parameters/variables (e.g., empirical equations and the triangle feature space [23][24][25][26][27][28]).However, the dependence of site and spatial domain is the main disadvantage of these methods [29,30].In addition, these methods are mainly driven by remotely sensed surface variables from a one-time observation.
The temporal variation in remotely sensed parameters/variables contains important information on surface fluxes, and it has been shown to be insensitive to retrieval error from remote sensing [31][32][33][34][35][36][37].Based on surface energy balance, the bulk heat transfer equation for sensible heat flux, and the assumption of fairly invariant evaporative fraction (EF) during the daytime, Lu et al. (2013) [38] derived a new parameterization scheme of directly estimating daily EF, which is parameterized as a function of temporal variations in surface temperature, air temperature, and net radiation.The input requirements can be easily met by remotely sensed data from the MODIS sensor and sensors onboard geostationary meteorological satellites.The parameterization scheme was applied to the Yucheng station in North China and encouraging results were obtained.Considering the relatively easy estimation of incoming solar radiation from remotely sensed data, this study will simplify the daily EF parameterization scheme using incoming solar radiation as an input.Details about the simplification will be presented in Section 2. To further clarify the potential ability of the simplified EF parameterization scheme and examine whether it is site-specific, ground-based measurements from Ameriflux sites with different vegetation types will be used to validate the EF estimates.Data and processing will be described in Section 3. Results and the uncertainties for EF estimates are analyzed in Section 4, and summaries and conclusions will be provided in Section 5.

Method
On the basis of the surface energy balance equation, surface latent heat flux (LE) at any time t in a day is expressed as follows: where R n is the net radiation, G is the soil heat flux, H is the sensible heat flux, and Q is the surface available energy, which is equal to R n − G and all variables are expressed in units of W/m 2 .The derivative of Equation ( 1) with respect to t is written as follows: H is theoretically determined by aerodynamic temperature (T aero ).Because T aero is not easily measured, surface radiative temperature (T s ) is frequently used to replace T aero for H estimates.By defining radiometric heat conductance (P rad ) corresponding to T s [39,40], H is expressed as follows: where T a is the air temperature (K).During a short study period of dt, provided that P rad depending on atmosphere and surface properties does not change considerably, the variation of H with time t is written as follows: From a study by Carlson et al. (1995) [23], P rad almost linearly varies with fractional vegetation coverage (f c ). G is frequently assumed to be proportional to R n , depending on f c [5].On the basis of the assumptions of relatively invariant EF during the daytime and low evaporation at night [41,42], daily average EF can be substituted by an instantaneous EF during the daytime.As a result, daily EF is parameterized as [38]: The variable dt in above equations is not required for a long period.It is generally universal that T s , T a , and R n on cloud-free days vary linearly with time during the morning [31,33].Therefore, the slope of the linear relationship can be used in Equation (6).Additionally, the day-night differences in T s , T a , and R n can be also applied because the day-night differences in surface parameters are almost positively proportioned to their variations during the morning, i.e.,

(
) where ΔT s , ΔT a , and ΔR n are the day-night differences in T s , T a , and R n or their variation during the morning.Considering the application for remotely sensed data from a geostationary meteorological satellite or a polar orbit satellite with a MODIS sensor, inputs for equation (7) [38].
Incoming solar radiation (R g ) is the main driving source of surface energy [43].Figure 1 shows the relationships between R n and R g based on ALEX-simulated data and in situ measurements.It can be inferred that R n is generally positively proportional to R g with a coefficient of determination (R 2 ) larger than 0.98.Therefore, ΔR n in Equation ( 7) can be further replaced by ΔR g but with different coefficients.Using ΔR g as an input variable of Equation ( 7) can obviate the calculation for longwave radiation.Furthermore, estimating R g from remotely sensed data is easier than estimating R n .Another purpose for using ΔR g is that there is no incoming solar radiation at night, and only a one-time measurement of R g in a day can meet the input requirements of Equation (7).Constant coefficients A, B, and C in Equation ( 7) corresponding to the input variable ΔR g are obtained using the ALEX-simulated data under different surface and atmospheric conditions.Daily EF is consequently expressed as follows: where ΔT s , ΔT a , and ΔR g are the differences of T s , T a , and R g between 1:30 P.M. and 1:30 A.M, respectively.With known f c , ΔT s , ΔT a , and ΔR g , daily EF can be calculated by Equation (8).When daily EF estimates with simulated data as inputs are compared with ALEX-simulated values, R 2 is 0.716 and a root mean square error (RMSE) is 0.129 (see Figure 2a).This result is similar to the estimation with ΔR n as an input, which obtained an R 2 of 0.719 and an RMSE of 0.130 [38].As shown in Figure 2b, daily EF estimates from two input schemes are nearly the same.Therefore, the substitute for ΔR n with ΔR g seems to be reasonable for daily EF estimates.

Study Sites
The input variables of Equation ( 8) can be obtained from remotely sensed data.Because there are numerous uncertainties in the processes of retrieving surface parameters/variables [44,45], the ground-based measurements from six Ameriflux sites with different land covers are used to validate the performance of the EF parameterization scheme of Equation ( 8).Among these stations, five stations (US-Ne2, US-Seg, US-SRC, US-Me6, and US-UMB) are located in the United States, and one (BR-Sa3) is located in Brazil.Land use and land cover at the stations include croplands (CRO), grasslands (GRA), open shrublands (OSH), evergreen needleleaf forests (ENF), deciduous broadleaf forests (DBF), and evergreen broadleaf forests (EBF), respectively (see Figure 3).The surface elevation varies from 234 to 1622 m, and the soil textures include deep silty clay loams, loamy sand, sandy loam, sandy lake deposits, and yellow dystrophic latosol.Basic information about these stations is listed in Table 1.

In situ Measurements
In situ measurements (meteorological variables, radiation data, and fluxes data from the eddy covariance (EC)) are collected from the Ameriflux network website [46] for the US-Ne2 station in 2011, the US-Seg station in 2010, the US-SRC station in 2011, the US-Me6 station in 2012, the US-UMB station in 2012, and the BR-Sa3 station in 2003.Collected T a and R g are directly used to drive the EF parameterization scheme of Equation (8).R n , H, LE, and G are used to calculate the EC-measured daily EF and then to evaluate the EF estimates.The daily EF is calculated as the ratio of daily averaged LE to daily averaged R n , neglecting the effect of G on EF at the daily scale.T s is not directly measured at these sites.Incoming and outgoing longwave radiation (L d and L u ) was used to calculate T s by the equation of with an assumption of 0.98 of surface emissivity (ε) in this study (σ is the Stefan-Boltzmann constant).Averaged values of the main meteorological variables of precipitation, air temperature, and relative humidity in one year were calculated based on the collected data (see Table 2).Annual precipitation at six sites varied from 182 to 1461 mm, annual mean air temperature ranged from 8 to 25 °C , and annual mean relative humidity ranged from 33% to 79%.There were missing measurements because of sensor malfunction, power outages, unfavorable weather, and other unknown reasons.Gaps in the missing data at the US-Ne2, US-Seg, and US-Me6 stations were filled by adopting the approach that combines measurement, interpolation, and empirical data synthesis [47], whereas the missing flux data at the US-SRC, US-UMB, and BR-Sa3 stations were not filled.Days with missing flux data were not considered in this study.Data at the US-Seg, US-SRC, US-Me6, and BR-Sa3 stations were averaged at a 30-min interval, whereas data at the Us-Ne2 and US-UMB stations were within 1 h.The EF parameterization scheme in Section 2 was derived under the assumption that EF was relatively invariant during the daytime, which most frequently occurs on cloud-free days with high radiation and humidity [41].An automatic procedure based on Tang et al. (2013) [48] was developed to select available data for days in which EF was relatively invariant.Automated techniques were designed considering that (1) the maximum of R g in a day occurs in the midday between 11:00 and 13:00; (2) R g increases monotonically from sunrise to midday and decreases monotonically from midday to sunset; (3) days with daily (24 h) averaged R g < 100 W/m 2 because of low solar radiation, which may occur in winter or be caused by full cloud cover during a day, are excluded; (4) days with daily mean T a < 0 due to low solar radiation and with ΔT s and ΔT a less than 0 that were probably influenced by nighttime cloud are excluded.Days having a measured daily EF < 0 and > 1 were also not utilized.A total of 408 cloud-free days were selected at these stations including 95 days in 2011 at the US-Ne2 station, 71 days in 2010 at the US-Seg station, 118 days in 2011 at the US-SRC station, 79 days in 2012 at the US-Me6 station, and 45 days in 2012 at the US-UMB station (see filled circles in Figure 4).Cloud-free days were not observed at the BR-Sa3 station in 2003 because of strong convective weather in the tropical region.Because remote sensing observations represent a snapshot in time, it is not easy to judge whether an entire day is clear using remotely sensed data.Therefore, the performance of EF parameterization for partly clear days, i.e., no cloud cover during the afternoon here, will also be discussed in this study.For the selection of partly clear days, the automatic procedure for judging clear days is equally applicable but not to regard the monotonic increase from sunrise to midday.A total of 77 partly clear days (marked by open circles in Figure 4) were selected including 12, 10, 10, 15, 14, 16 partly clear days at the US-Ne2 station, at the US-Seg station, at the US-SRC station, at the US-Me6 station, at the US-UMB station, and at the BR-Sa3 station, respectively.
The energy balance non-closure in EC-based measurements was frequently reported in many previous studies [49].As shown in Figure 5

Fractional Vegetation Cover f c
Apart from the variables ΔT s , ΔT a , and ΔR g , f c is also required as an input for the EF parameterization scheme presented as Equation (8).Because of the lack of in situ measurements for vegetation cover, MOD13Q1 and MYD13Q1 products in the spatial resolution of 250 m, which provide the normalized difference vegetation index (NDVI) every 16 days, are combined to calculate f c every 8 days by linear mixture modeling [50], and computed using Equation ( 9 where NDVI min and NDVI max are the NDVI for bare soil and dense vegetation at full coverage, respectively.NDVI min and NDVI max are not only determined by surface characteristics but also are dependent on remotely sensed image data.For bare soil, NDVI tends to vary between −0.1 and 0.2 [13], but in this study, NDVI min is given as 0. For dense vegetation, NDVI depends on vegetation type, age, and leaf water content.Carlson et al. (1990) [51] mentioned that the upper limit of NDVI for dense vegetation is in the range of 0.5-0.8,whereas Prihodko and Goward (1997) [52] used the observation of reflectance and transmittance in the red and near-infrared wavelengths from the FIFE information system to compute the NDVI of an infinitely thick canopy for each species, and obtained an average NDVI of 0.86 for the full canopy.Determination of the value of NDVI max for different surface characteristics is not easy.NDVI max is set as 0.86 for all sites here.Because the scattering coefficients for the leaves and the soil, the leaf-area index, the leaf angle distribution, and the angle of incident radiation affect the NDVI, the vegetation cover problem becomes more complex [13].The errors in EF estimation caused by the uncertainty of f c will be discussed in Section 4.2.

Daily EF Estimates for Clear Skies and Partly Clear Skies
For 408 cloud-free days, with known ΔT s , ΔT a , and ΔR g from in situ measurements and f c from remotely sensed data, daily EF was estimated using Equation (8).Compared with EC-based measurements, daily EF was overestimated by 0.212 (see black triangles in Figure 6a).The overestimation was mainly caused by the underestimation of LE from EC measurements.Two methods were frequently used to correct EC-measured LE, i.e., the residual energy (RE) method and the Bowen ratio (BR) method [49].The LE corrected by the RE method is equal to R n − G − H with the assumption of relatively accurate measurements for H.The BR method is to assume that the non-closure of energy balance is not only caused by EC-measured LE but also by EC-measured H.The imbalance of energy is partitioned into H and LE according to the Bowen ratio, i.e., The corrected LE is then used to calculate EC-measured daily EF.When the estimated EF from Equation ( 8) is compared with the BR/RE-corrected EF, the overestimation in EF is reduced (see green circles and blue symbols in Figure 6a, respectively).Daily EF estimates are closer to the RE-corrected EF, showing an RMSE of 0.152 and an R 2 of 0.586.For selected 77 partly clear days, daily EF estimates from Equation ( 8) are also closer to the RE-corrected daily EF with an R 2 of 0.549 and an RMSE of 0.171 (see Figure 6b).The result indicates that the EF parameterization scheme of Equation ( 8) can also be applied to estimate daily EF for partly clear skies, which mainly occur during the afternoon.

Uncertainty in EF Estimates Caused by f c
Uncertainty in determining input variable f c is also a source of error for daily EF estimates.Because f c is from the 16-day MODIS NDVI products at spatial resolution of 250 m, the different temporal and spatial scales between in situ measurements and the MODIS sensor could result in discrepancies between estimated daily EF and measured values.Sensitivity and error analysis for f c made by [38] concluded that 0.2 variations in f c can lead to approximately 0.1 mean variations in daily EF estimates.Different schemes for calculating f c and the assignments for NDVI max and NDVI min can produce different daily EF estimates (see Figure 7 and Table 4).f c from a non-linear model relating to NDVI [53], i.e., the square of equation ( 6), is less than that from the linear model with an RMSE of 0.183 (see filled triangles in Figure 7a).The underestimation in f c will cause the overestimation of daily EF according to the mathematical formulation of equation (8).Therefore, when f c from the non-linear model is used as an input for equation (8), daily EF estimates are greater than the results from the linear model with an RMSE of ~0.1 (see filled triangles in Figure 7b).The larger bias is on the lower range of EF values, most of which are related to the condition of low f c .When the NDVI max and NDVI min in the linear model of Equation ( 6) are given as 0.66 and 0.2, respectively [51,53], f c is less than the values calculated by 0 and 0.86 as inputs for low f c .The opposite is true for high f c (see open circles in Figure 7a).For daily EF estimates, the RMSE caused by different assignments for NDVI max and NDVI min is also ~0.1 (see open circles in Figure 7b).Those open circles deviating from the 1:1 line in Figure 7b correspond to low f c conditions, which are attributed to underestimation of f c .The larger bias also occurs at the lower EF condition.Therefore, daily EF estimation from Equation ( 8) is more sensitive to low f c , particularly for a low ET surface.The implicit assumption in Equation ( 8) that EF is relatively invariant during the daytime may be valid at this situation.Sensitivity of EF estimates to low ET surfaces is similar to the findings of Long et al. (2014) [54], where uncertainties among the ET products increased with aridity.A low ET surface is generally related to low soil water content.Further study is required to determine if soil water content can be used as an input parameter for Equation ( 8), and to find the best approach for parameterizing the EF equation.

Performance of EF Parameterization Scheme for Different Land Covers
The six stations used in this study correspond to six different land covers.Daily EF estimates for clear days and partly clear days were both considered in this section.The estimated EF for each station is compared with the RE-corrected EF, which is shown in Figure 8.Most of the absolute biases of daily EF estimates from Equation (8) for different vegetation types are less than 0.1, and the RMSEs are not more than 0.2.The estimated daily EF for CRO and OSH surfaces are correlated more strongly to the RE-corrected EF with the R 2 of 0.640 and 0.650, respectively.For those surfaces covered by GRA, ENF, DBF, and EBF vegetation types, most results appear to be consistent with the RE-corrected EF, but some scatter points and relatively concentrated distribution of daily EF lead to the low R 2 .It seems that surface elevation and soil texture have no obvious contribution to the errors of EF estimates, i.e., the errors of EF estimates do not strongly depend on surface elevation and soil texture.Errors of daily EF estimates, daily averaged R g , daily total precipitation, f c , and ΔT s − ΔT a of six stations are displayed in Figures 9-14 to further evaluate the reason for discrepancies.(1) For the surface covered by CRO vegetation types, which is at the US-Ne2 station, a majority of the EF estimates are satisfied (see Figure 8a), and the errors vary from −0.2 to 0.2 (see Figure 9).Several large errors (>0.3) for EF estimates occur on DOY 75, 79, 128, 129, 293, 294, 359, and 362 (marked by red circles in Figure 9).The overestimation on DOY 75 may be mainly caused by the small value of ΔT s − ΔT a , which is mostly from measurement errors or is influenced by high wind speed (wind speed in the afternoon at approximately 13:30 on DOY 75 is ~9 m/s).The error on DOY 79 could be attributed to the effect of clouds, and obvious low averaged R g was observed on the day.The overestimation on DOY 128, 129, 293, and 294 may be due to the underestimation in f c .On DOY 359 and 362, low R g that occurred in winter may be the main reason leading to the negative biases of EF estimates.This is because under the condition of low solar radiation, the implicit assumption of relatively invariant EF in Equation ( 8) is not well met.In general, under the condition of strong solar radiation, sufficient soil water content, dense vegetation cover, and low ΔT s − ΔT a , errors in EF estimates are relatively small.(2) At the US-Seg station covered by GRA vegetation type, the climate is semi-arid, and the annual precipitation of the region is ~230 mm, with a major portion (>50%) occurring in a summer rainy season, July to September [55].In 2010, the total precipitation was only 182 mm (see Table 2).Therefore, EF from measurements at the US-Seg station was generally low, even less than 0.6 (see Figure 8b).Most of the EF estimates from Equation ( 8) generally fall into this range but with a slight overestimation.According to Equation ( 8), the underestimation in f c would lead to an overestimation in EF.Therefore, low f c of ~0.2 in 2010 (may be underestimated) results in the overestimation of EF at the US-Seg station.On DOY 6, 336, and 358, larger errors may be caused by low solar radiation in winter (see Figure 10).The relatively concentrated distribution of daily EF at the US-Seg station in 2010, approximately from 0.2 to 0.6, may lead to the low R 2 of 0.148 in Figure 8b.(3) At the US-SRC station, most of daily EF values lie in the range of 0.1 to 0.6 (see Figure 8c).This may be related to low precipitation of ~229 mm and low vegetation cover of ~0.3 in 2011 (see Table 2 and Figure 11, respectively).Although the R 2 between the estimated EF and measured values is high, there is an overestimation when measured EF is low (See Figure 8c).As analyzed in Section 4.2, Equation ( 8) is more sensitive to low f c for daily EF estimates, particularly for a low ET surface.It is difficult to accurately determine f c particularly at the sparse vegetation surface covered by OSH types.Therefore, the potential underestimation in f c may lead to the overestimation of daily EF at the US-SRC station.The obvious underestimation on DOY 5 was primarily because of the effect of clouds.The positive biases on DOY 93 and DOY 149 were from the relative low values of ΔT s − ΔT a probably due to high wind speed (wind speeds at approximately 13:30 on DOY 93 and DOY 149 were ~6 m/s and 8 m/s, respectively).The energy balance non-closure in EC-based measurements at this station with the EBR of 0.79 (see Table 3) may have led to discrepancies between the estimated EF and measured values.
(4) At the US-Me6 station covered by evergreen needleleaf forest, the simplified EF parameterization scheme gave the worst results with an RMSE of 0.192 (see Figure 8d), and most of the EF was underestimated (see Figure 12).This may be related to the local climate having warm dry summers and cool wet winters.Most precipitation (~97% of annual rainfall) falls between October and June.The climatic condition of the warm dry summer and cool wet winter may not easily meet the assumption of Equation ( 8) because EF during the daytime is relatively invariant when there is high radiation and humidity.The underestimation of EF on DOY 70 and 123 and the overestimation on DOY 139 and 171 may result from clouds, whereas the overestimation on DOY 138 and 161 was evidently from low values of ΔT s − ΔT a due to high wind speed of ~6 m/s.There was no precipitation on DOY 189, 192, 193, and 195, i.e., there was insufficient water in the soil, resulting in less ET on these days.The low EF from Equation (8) seems more reasonable; therefore, large negative biases on DOY 189, 192, 193, and 195 were attributed to the uncertainty in measurements or the errors from the RE correction method.As shown in Table 3, energy balance non-closure is most noticeable for EC-based measurements that occurred at this station.(5) For the surface covered with deciduous broadleaf forest, daily EF mainly varies from 0.5 to 1 (see Figure 8e).This is related to the temperature and humidity climate at the US-UMB station with a mean temperature of 9 °C and total precipitation of 609 mm in 2012 (see Table 2).The overestimation on DOY 70 may be caused by low solar radiation and underestimation in f c (see Figure 13).During the summer days from June to September, the surface is completely covered by broadleaf vegetation, and soil moisture is ample.Under this condition, high ET value is more reasonable.Therefore, the large overestimation of EF on DOY 245, 259, 271, 273, and 274 may be from the uncertainty in measurements or errors for energy imbalance correction.(6) The rainy tropical climate results in high EF of ~0.8 at the BR-Sa3 station covered by evergreen broadleaf forest (EBF) during one year (see Figure 8f), being coincident with the estimated values close to 0.8.The concentrated distribution of EF during one year leads to low R 2 .Although the data displayed in Figure 14 occurred on partly clear days, there is not an obvious large bias for daily EF estimates at this station.The abundant rainfall and dense vegetation led to the low values of ΔT s − ΔT a and reduced the effects of clouds on EF estimates.

Conclusions
Lu et al. (2013) [38] developed a method for estimating daily evaporative fraction (EF) derived by day-night differences in surface temperature (T s ), air temperature (T a ), and net radiation (R n ).Because incoming solar radiation (R g ) can be readily estimated from remotely sensed data and there is a strong relationship between R g and R n , this study simplified the equation for estimating daily EF using R g .The accuracy of the simplified EF parameterization scheme with ΔR g as an input is nearly equivalent to that with ΔR n as an input.
When in situ measurements from Ameriflux sites are used to drive the simplified parameterization scheme, daily EF estimates for clear days are closer to the RE-corrected daily EF with a coefficient of determination of 0.586 and a root mean square error of 0.152.It is not an easy thing to judge the sky condition using several remote sensing observations.However, the developed method can be used to estimate daily EF for partly clear days with the accuracy similar to that for clear sky condition.Errors from the correction for energy imbalance in EC-based measurements and the uncertainty in calculating fractional vegetation cover (f c ) are the primary factors that lead to discrepancies between estimated and measured daily EF.Daily EF estimates from the simplified parameterization scheme are more sensitive to low f c conditions with low ET surfaces.Analyses of daily EF estimates at six vegetation surfaces indicate that the constant coefficients in the simplified EF parameterization scheme are not strongly site-specific, and could be applied to different land covers with reasonable accuracy if surface vegetation cover can be accurately determined.EF estimates are more accurate for regions that are characteristically warm and humid with strong solar radiation and adequate soil water content.
The input requirements of relatively few variables that can be easily obtained from remotely sensed data, i.e., the day-night difference in surface and air temperatures, can reduce uncertainties in EF estimates due to retrieval error from remote sensing.Only by using the parameters/variables of f c , ΔT s , ΔT a , and ΔR g , can Equation (8) estimate daily EF, which can be further used to partition the surface available energy into sensible and latent heat fluxes.Although the accuracy of the EF parameterization scheme is acceptable, there is still room for improving the model, particularly for low ET surfaces with low f c .Further study should be directed toward improving daily EF estimates by introducing other retrieval variables from remote sensing, such as soil moisture.

Figure 1 .
Figure 1.Relationships between surface net radiation (R n ) and incoming solar radiation (R g ) based on ALEX-simulated data with different surface characteristics for three representative clear days and in situ measurements at the US-SRC station for clear skies in 2011 (see Section 3).Atmosphere forcing driven the ALEX model was measured at the Yucheng station on 24 April, 29 May, and 21 July 2010, and more details can be found in the Appendix of [38].

Figure 2 .
Figure 2. Comparisons of daily evaporative fraction (EF) estimates from Equation (8) with (a) ALEX-simulated values and (b) the results using ΔR n .

Figure 3 .
Figure 3. Distribution of sites for validation of EF parameterization scheme (The background is the land use type of 2012 classified according to the University of Maryland scheme, which is extracted from the Land Cover Type Climate Modeling Grid product MCD12C1 with 0.05 degree spatial resolution.).

Figure 4 .
Figure 4. Distribution of clear days and partly clear days in one year.
, daily average H + LE are all less than daily average R n − G for clear and partly clear skies at the six sites.The mean energy balance ratio (EBR) is ~0.83, and the RMSE is ~25 W/m 2 .For each site, the daily-scale statistics of energy balance closure of EC-based measurements including EBR, R 2 , RMSE, and BIAS are listed in Table3.The most serious energy balance non-closure in EC-based measurements occurs at the US-Me6 station with a mean EBR of 0.78, and the RMSE between H + LE and R n − G is 35.2W/m 2 .At the BR-Sa3 station, the EC-measured daily average H + LE is closer to R n − G with an EBR of 0.93.In general, the inconsistent source areas between radiation flux (R n − G) and turbulent flux (H + LE), and neglecting some terms in surface energy balance equation such as storage heat in forest vegetation, the energy due to photosynthesis and advection, may be the main reasons for non-closure in the energy balance from EC-based measurements.

Figure 5 .
Figure 5. Comparisons of daily average H + LE from EC-based measurements with daily average R n − G for clear and partly clear skies at six Ameriflux sites.

Figure 6 .
Figure 6.Comparisons of daily EF estimates with values from EC measurements, the Bowen ratio (BR) correction method, and the residual energy (RE) correction method for (a) clear skies and (b) partly clear skies.

Figure 7 .
Figure 7. Comparisons of (a) f c from different schemes and (b) daily EF estimates corresponding to different f c (f c_non-linear represents the results from a non-linear model; f c_0.2/0.66 denotes the results for NDVI min = 0.2 and NDVI max = 0.66).

Table 4 .
Results of f c from different schemes and corresponding daily EF estimates compared with those values from the linear model (f c_non-linear represents the results from a non-linear model; f c_0.2/0.66 denotes the results for NDVI min = 0.2 and NDVI max = 0.66).

Figure 9 .
Figure 9. Error in EF estimates, daily averaged incoming solar radiation (R g ), fractional vegetation cover (f c ), and ΔT s − ΔT a at the US-Ne2 station for all selected clear days and partly clear days, and daily precipitation (P) for all days in one year.

Figure 10 .
Figure10.Error in EF estimates, daily averaged incoming solar radiation (R g ), fractional vegetation cover (f c ), and ΔT s − ΔT a at the US-Seg station for all selected clear days and partly clear days, and daily precipitation (P) for all days in one year.

Figure 11 .
Figure 11.Error in EF estimates, daily averaged incoming solar radiation (R g ), fractional vegetation cover (f c ), and ΔT s − ΔT a at the US-SRC station for all selected clear days and partly clear days, and daily precipitation (P) for all days in one year.

Figure 12 .
Figure 12.Error in EF estimates, daily averaged incoming solar radiation (R g ), fractional vegetation cover (f c ), and ΔT s − ΔT a at the US-Me6 station for all selected clear days and partly clear days, and daily precipitation (P) for all days in one year.

Figure 13 .Figure 14 .
Figure 13.Error in EF estimates, daily averaged incoming solar radiation (R g ), fractional vegetation cover (f c ), and ΔT s − ΔT a at the US-UMB station for all selected clear days and partly clear days, and daily precipitation (P) for all days in one year are based on five different schemes including (1) change rate during the morning, (2) differences between 1:30 P.M. (MODIS/Aqua daytime overpass time) and 1:30 A.M. (MODIS/Aqua nighttime overpass time), (3) differences between 10:30 A.M. (MODIS/Terra daytime overpass time) and 10:30 P.M.

Table 2 .
Summary of in situ measurements at six stations (  is the annual mean air temperature, °C ; RH is the annual mean relative humidity, %.)

Table 3 .
Statistics of energy balance closure of EC-based measurements at the daily scale (EBR is the energy balance ratio, which is the mean value of the ratio of daily average H+LE to daily average R n − G )