Research Papers

Evaluation of a MODIS triangle-based evapotranspiration algorithm for semi-arid regions

[+] Author Affiliations
Jongyoun Kim

University of California, Department of Civil and Environmental Engineering, Los Angeles, California 90095

Terri S. Hogue

University of California, Department of Civil and Environmental Engineering, Los Angeles, California 90095

Colorado School of Mines, Civil and Environmental Engineering, 1500 Illinois Street, Golden, Colorado 80401

J. Appl. Remote Sens. 7(1), 073493 (Oct 28, 2013). doi:10.1117/1.JRS.7.073493
History: Received January 14, 2013; Revised September 11, 2013; Accepted October 1, 2013
Text Size: A A A

Open Access Open Access

Abstract.  The current study investigates a MODIS-based remote-sensing-based algorithm for the assessment of daily latent heat flux (LE), or evapotranspiration (ET), with a focus on semi-arid regions. The approach uses the triangle method, relying on remotely sensed inputs and a previously developed net radiation model. A major difference from previous studies is that we utilize only MODIS products for estimation of ET at the daily timestep for clear and cloudy days. The algorithm is evaluated at four flux tower locations in the San Pedro River basin in Arizona. The mean daily LE varies significantly between the sites, ranging from 144 to 179W/m2 in the riparian areas to 36 to 76W/m2 at the rangeland sites. Comparison of the flux towers shows good correlation and low root mean square error at the riparian sites (0.81 to 0.82 and 37 to 51W/m2, respectively) with slightly larger errors at the upland sites, where ET is strongly correlated to precipitation events. The model assumption of a linear variation in evaporative fraction across the triangular domain (LST/EVI space) results in more uncertainty under water-stressed conditions such as those found at the upland sites. Overall, the proposed MODIS-based algorithm provides reasonable estimates of riparian and upland plant water use and unique spatial and temporal information.

Figures in this Article

Evapotranspiration (ET) is a key parameter in numerous disciplines, including water resource management, agriculture, ecology, and climate change science. However, ET is arguably one of the most difficult hydrologic components to estimate, given the dependence on a range of climatological parameters (i.e., solar radiation, temperature, wind speed, vapor pressure, etc.) as well as soil physical properties, land cover, and the surrounding environment. Traditional ground-based measurements of land surface heat fluxes (e.g., Bowen ratio system and eddy covariance system) represent a relatively small area or footprint, and spatial interpolation is sometimes used to regionalize these values.1 A number of ET methods using meteorological and radiosonde observations are available, but the required data are only found typically over small areas of the global land surface. Also, utilization of reanalysis data has limitations due to its coarser scale.2,3 Alternatively, remotely sensed data have been recognized as a means to provide broader areal coverage, more frequent estimates, and moderate (or high) spatial resolution.35

A number of models with varying complexity have been developed to estimate ET based on combinations of remote-sensing observations, ground-based data, and/or high-resolution models. These include surface energy balance algorithm for land (SEBAL), mapping evapotranspiration at high resolution with internalized calibration (METRIC), atmosphere-land exchange inverse (ALEXI), and temperature-vegetation index triangle methods.514 Although significant progress has been made in the past 30 years on the use of satellite-based products to provide estimates of ET at daily to weekly time scales and higher spatial scales (1 m to 1 km), most efforts have been hindered by two key issues, including: (1) water fluxes and their associated stores must be estimated indirectly using algorithms that relate, for example, leaf area index, vegetation indices, and land surface temperature, to ET, and (2) key parameters, such as wind speed, air temperature, vapor pressure deficit, and vegetation height, are not available at regional scales or from remote-sensing platforms. In addition, deriving accurate spatial representation of ground-based measurements for large heterogeneous regions is problematic, especially in areas with limited observational networks. Residual methods, such as SEBAL and METRIC, minimize the requirements of surface measurements through unique internal calibration methods that retrieve air temperature and convert wind speed at blending height.68 However, these procedures (especially air temperature retrieval) also incur some degree of uncertainty. In order to overcome these problems, the temperature-vegetation index triangle method attempts to develop a regional parameterization of ET [evaporative fraction (EF)] through a relatively simple approach that solely uses satellite-derived surface parameters without complex parameterization and the need for extensive ground-based measurements. Further, the need for absolute accuracy from satellite-derived temperatures is avoided.10,11

Since its inception, several modifications of the triangle (or trapezoid) method have been developed and widely applied to study soil moisture, land use, and drought monitoring, using data from several satellite sensors, including MODIS, AVHRR, GOES, etc.12,1520 These studies established their own version of the triangle approach. For example, Jiang and Islam20 utilized a simplified Priestley–Taylor formulation to build up the triangle method and Jiang et al.18 introduced a trapezoid shape by taking into account correction parameters for vegetation transpiration that consider water stress in the root zone.18,20 Wang et al. combined day-night temperature differences by using MODIS Terra and Aqua products with NDVI in their triangle method.17 However, the idea behind these approaches is similar in that variations of temperatures for a given vegetation index are due to evaporative cooling effects and the proposed triangle methods are able to derive ET without auxiliary data or site-specific relationships.21,22

In the current study, we investigate a model for estimation of ET utilizing only MODIS satellite data and limited ancillary data. We combine the triangle method developed by Jiang and Islam12 and Wang et al.17 with thermal inertia information obtained from the MODIS sensor to estimate a regional EF. Additionally, MODIS remote-sensing data is used for estimating daily net radiation and ground heat flux, which in combination with the triangle-based EF, returns a daily actual evapotranspiration product. Specifically, the objectives of the current study are: (1) to evaluate the modified net radiation model from the Kim and Hogue23 framework over semi-arid regions, (2) to estimate the actual ET (LE) by integrating the Kim and Hogue23 net radiation model and additional MODIS products, (3) to investigate the ability or restriction of MODIS-based data to account for the spatial variation in fluxes for many landscapes, and (4) to evaluate the robustness of this parsimonious approach against ground-based flux tower observations representing various ecosystems in a semi-arid region.23

Initial validation of the proposed model is carried out at four flux tower sites in southern Arizona (Table 1; Fig. 1). We utilize two riparian sites and two upland sites to evaluate the performance in important biomes in the region. Specifically, we test our algorithm at a desert grassland (Kendall site), a savanna (Santa Rita site), a riparian mesquite (Charleston site), and a riparian grassland (Lewis Spring site). The climate is classified as semi-arid with average annual precipitation ranging between 354 and 458 mm, and 53% of the total rainfall occurring from July through September during the North American monsoon (NAM).2426 Precipitation during the NAM is characterized by localized short-duration, high-intensity convective thunderstorms, which closely tie to ecosystem response in this water-limited region.27,28 For the current study year (2005), the NAM was noted to have begun on July 18 and ended on September 10.29 One of the flux observation sites (Kendall) is located within the Walnut Gulch experimental watershed. This site is a typical of southwestern rangeland, where cattle graze on gentle hillslopes comprised mainly of C4 grasses with a few scattered shrubs.30 The Kendall site provides a set of soil temperature and moisture profile measurements, rain gauge, and runoff observations operated by the U.S. Department of Agriculture’s Agricultural Research Service (ARS). Another upland site is the Santa Rita site located in the Santa Rita Experimental Range, approximately 45 km south of Tucson, Arizona. Over the past several decades, the rangeland around Santa Rita mesquite site transitioned from a semi-desert grassland into a savanna.30 Vegetation consists of mesquite growing in a matrix of native and nonnative perennial grasses, subshrubs, and forbs. Soil at the upland site is a coarse-textured sandy loam.31 Beside the two upland sites (Kendall and Santa Rita), we also evaluate the model performance at two riparian sites (Charleston and Lewis Spring), which reside along the upper San Pedro Basin. The composition and amount of vegetation in the riparian area differ significantly from those in the terrestrial upland areas, reflecting the influence of water from the adjacent San Pedro River on soil moisture in the riparian zone. Vegetation at the Charleston site is predominantly dense woodland with a maximum vegetation height of 10 m. The Lewis Spring site is a giant sacaton grass.32 Using the model of Hsiegh et al., flux tower footprints ranged from 50 to 230 m for upland sites and 25 to 80 m for the riparian sites during daytime conditions.33

Table Grahic Jump Location
Table 1Location and characteristics of the four tower sites utilized for the study and single-pixel extraction.
Table Footer Note*Precipitation and reference crop evaporation (AZMET standard) are the annual averaged value for year 2005.
Graphic Jump LocationF1 :

Study area in Arizona with a Landsat scene (path 35 and row 38; false-color display of the 30 m Landsat5 VNIR imagery) including the four flux towers (Kendall, Charleston, Lewis Spring, and Santa Rita) marked as triangles.

In Situ Measurements

Study flux towers consist of a net radiometer instrument and several soil heat flow transducers with soil temperature/moisture sensors for estimating the available energy. Sensible and latent heat fluxes are measured by eddy covariance (consisting of a sonic anemometer and a fast-response open-path water vapor/CO2 sensor) at all study sites. Soil heat fluxes are measured with soil heat flux plates and are corrected for heat storage above the plate. Thus, relevant components of the surface energy balance, including measurements of net radiation, soil heat flux, sensible and latent heat fluxes, are provided at all study flux towers. In addition, ancillary meteorological observations (air temperature, wind speed, vapor, atmospheric pressure, precipitation, etc.) are available at the four study sites. Measurements at the grassland and shrubland upland sites (i.e., Kendall and Santa Rita) are made at 2 to 3 m above the local terrain. Measurements at the riparian sites (i.e., Charleston and Lewis Spring) are made at 1.5 to 2.5 m above mean canopy height. All sites have data recorded at 30-min intervals. When the daily average is calculated, days with less than 75% of data available are excluded.

Satellite Observations (MODIS)

The Moderate Resolution Imaging Spectroradiometer (MODIS) sensor, with 36 spectral bands (20 reflective solar and 16 thermal emissive bands), provides unprecedented information regarding vegetation and surface energy, which are critical for developing a remotely sensed ET model.34 The proposed MODIS-based ET model requires a total of 15 variables obtained from MODIS atmospheric and land surface products for the calculation of net radiation and ground heat flux, as well as vegetation index (VI) and land surface temperature (LST) products for the development of the EF (Table 2). The required Aqua MODIS data are extracted from the NASA EOS Gateway (https://wist.echo.nasa.gov/api/) and the Land Processes Distributed Active Archive Center (LP DAAC) (http://modis.gsfc.nasa.gov) with a standard hierarchical data format (HDF). The albedo data used in the current study is a Terra and Aqua combined product (MCD43A1). For simplicity, we utilize a solar zenith angle equal to local solar noon and an optical depth of 0.2 as default values in calculating the actual (blue-sky) albedo based on a known black-and-white sky albedo (http://daac.ornl.gov/MODIS/MODIS-menu/MCD43.html). Daily (MYD11A1) and 8-day composite (MYD11A2) LST products were obtained at a 1-km spatial resolution. Both daytime and nighttime products from MYD11A2 are acquired in order to calculate the difference of daytime and nighttime of LST. The VI was obtained using 16-day composite datasets with a 250-m resolution from both Terra (MOD13Q1) and Aqua (MYD13Q1). The MODIS composite technique (constrained-view angle-maximum value composite) includes determination of the two greatest VI values for each pixel per 16-day composite interval. The VI value observed with the nearest to nadir view is then selected for inclusion in the composite product.35 An 8-day phasing in the production of the 16-day composites between Terra and Aqua allows the generation of a combined 8-day time series of VI data. We note that the compositing process may result in phenological changes going undetected, particularly where localized short-duration, high-intensity convective thunderstorms are characterized.36

Table Grahic Jump Location
Table 2Summary of MODIS products used in this study with relevant spatial and temporal information.

Since enhanced vegetation index (EVI) is optimized to improve the vegetation signal and reduce soil background influence, we utilize EVI instead of normalized difference vegetation index (NDVI) in the current study.37 We note that the selection of VI may affect the shape of the final EVI-LST triangle; however, the EVI-calculated triangle method has been shown to perform as well as, or better than, than NDVI-calculated triangle method in most cases.38

Estimation of Available Energy (Rn−G)

To estimate the daily available energy (RnG), which in combination with EF returns the daily actual ET, we propose a satellite-based stand-alone methodology.

The net radiation for clear sky conditions is estimated according to Display Formula

Rnclear=(1A)Rsclear+RlclearRlclear,(1)
where A is the surface albedo, Rsclear is the downward shortwave radiation, Rlclear is the downward longwave radiation, and Rlclear is the upward longwave radiation for clear sky condition. The solar radiation scheme known as the PS parametric model is used to estimate the downward shortwave radiation for clear skies.39 This scheme accounts for many factors, including gaseous absorption, Rayleigh scattering, and scattering and absorption by ozone and aerosols. Further details of this scheme can be found in Kim and Hogue.23

Upward longwave radiation is expressed using the Stefan–Boltzmann equation: Display Formula

Rlclear=εsσTs4,(2)
where εs is surface emissivity, σ is the Stefan–Boltzmann constant (5.67×108W/m2K4), and Ts is the surface temperature.

A parameterization scheme derived by Brutsaert40 is chosen to estimate downwelling longwave radiation. This methodology is chosen based on the comparison among eight longwave parameterization schemes.4047Display Formula

Rlclear=εaσTa4,(3)
where εa=1.723(e0/Ta)1/7.

In the Brutsaert scheme, air emissivity (εa) is determined by water vapor pressure (e0) and air temperature (Ta). Water vapor pressure is obtained from MYD05, and air temperature is retrieved through Eqs. (6) and (7).

Under cloudy conditions, we use a formulation of downward shortwave radiation proposed by Bisht and Bras.48Display Formula

Rscloudy=Rsclear[(1N)+Neτ/cosθ],(4)
where N is the cloud cover fraction, τ is cloud optical thickness from MYD06, and θ is solar zenith angle from MYD03. In the current study, Rsclear is the theoretical clear-day shortwave radiation, which is generated by a simple interpolation scheme between adjacent clear days.

The downward longwave radiation for cloudy pixels is estimated as proposed by Bisht and Bras.48Display Formula

Rlcloudy=εaσTa4+(1εa)εcσTc4,(5)
where εc is cloud emissivity, and Tc is the cloud temperature from MYD06. To retrieve air temperature (Ta), Bisht and Bras48 proposed a method using the offsets between air temperature from the MYD07 data and surface temperature from the MYD06 data. However, we calculated the interpolated ratio (Ratioint) between air temperature from the MYD07 data and surface temperature from the MYD06 data using inverse distance weighting [Eq. (6)] and back-calculated air temperature using Eq. (7). Display Formula
Ratioint=Ta_MYD07/Ts_06(6)
Display Formula
Ta=Ratioint(Ts_06),(7)
where Ts_06 is surface temperature from MYD06.

To estimate the upward longwave radiation for cloudy-sky pixels, we obtain the surface temperature (Ts_06) from MYD06 to substitute for the MYD11 data product, and surface emissivity (εs_A2) is obtained from MYD11A2, which is an 8-day composite product. Display Formula

Rlcloudy=εs_A2σTs_064.(8)

Using Eqs. (1)–(8), we can now calculate net radiation estimates for all sky conditions at the satellite overpass (instantaneous) time. Finally, we convert the net radiation estimates from the instantaneous time to a daily average net radiation estimate based on a sinusoidal function, assuming that net radiation begins to rise (become positive) at sunrise and declines at sunset.23,49,50Display Formula

Rn¯=Rni×2πsin[(titsunrisetsunsettsunrise)π],(9)
where Rn¯ and Rni are daily average and instantaneous net radiation in W/m2, respectively. tsunrise, tsunset, and ti are sunrise and sunset times obtained from the U.S. Naval Observatory for the satellite over-passing time.

Soil heat flux (G) is estimated from the net radiation estimates and the MODIS EVI, as proposed in the Kim and Hogue23 methodology. However, due to the relatively poor performance of the ground heat flux scheme noted in Kim and Hogue23 at the semi-arid study site, we incorporate a modified ground heat flux scheme into our previous algorithm for application in the Arizona region. Given that empirical schemes have strong sensor dependency, we updated the previously formulated Moran and Jackson2 relationship used in our prior work. We also evaluate the performance of various MODIS-based VI (NDVI, VI, and EVI) in approximating flux tower observations of ground heat flux for the study sites. Huete et al. demonstrated that EVI is optimized to improve the vegetation signal and reduce the soil background influences in the semi-arid regions.37 EVI also provided the best performance in our study, with optimized parameters of 0.22 and 1.4, among the vegetation indices tested. The revised ground heat flux scheme is developed as follows: Display Formula

G¯=0.22·exp(1.4EVI)·Rn¯,(10)
where Rn¯ is the net radiation, and G¯ is the ground heat flux at daily average in W/m2.

Estimation of Evaporative Fraction

The EF is defined as the ratio of LE to available energy.15Display Formula

EF=LERnG.(11)

We utilize a model developed by Wang et al., which derives EF through interpretation of the Priestley–Taylor parameter α using the trapezoidal distribution of a ΔT (aqua day-night temperature difference)—NDVI spatial variation.17 Again, in the current study, we incorporate EVI instead of NDVI in our model. We also chose to utilize the temporal variation of Ts, rather than Ts itself, to avoid a significant absolute bias from the MODIS LST product as noted in previous studies.49,51 EF is formulated as Display Formula

EF=αΔΔ+γ,(12)
where γ is the psychrometric constant (unit: hPa/K), and Δ is the slope of saturated vapor pressure at the air temperature (unit: hPa/K), which can be calculated as Display Formula
Δ=26297.77(Ta29.65)2exp[17.67(Ta273.15)Ta29.65].(13)

Due to the instability of air temperature retrieval from the MOD07 product and no significant difference (less than 5%) between the use of air and surface temperatures reported by Wang et al., surface temperature is used to estimate Δ instead of air temperature.17

At each pixel, ET consists of surface evaporation from a bare soil and transpiration from vegetation. For a given EVI, variations of surface temperature occur primarily as a result of stressed surface soil water content. At ΔT=0, on the cold edge, bare soil and vegetation will evaporate at the potential rate, implying that ET is not limited by soil moisture. As ΔT increases, evaporation from bare soil decreases and become negligible on the warm edge.

Based on the above interpretation of the (ΔT and EVI) space, we can derive Display Formula

ΔTiΔTminΔTmaxΔTmin=(αmaxαmin)(αiαmin)αmaxαmin.(14)
Display Formula
αi=ΔTmaxΔTiΔTmaxΔTmin(αmaxαmin)+αmin.(15)

From Jiang and Islam,12αmax(=1.26) and αmin can be calculated as12Display Formula

αmin=1.26fveg=1.26EVIiEVIminEVImaxEVImin.(16)

By substituting αi from Eq. (15) into Eq. (12), we obtain EF for any point within the boundary of the triangle domain as Display Formula

EF=ΔΔ+γ[ΔTmaxΔTiΔTmaxΔTmin(αmaxαmin)+αmin].(17)

Equation (17) can also be rewritten as Display Formula

EF=1.26ΔΔ+γ[ΔTmaxΔTiΔTmaxΔTmin(1fveg)+fveg].(18)

Since a range of LST and vegetation indices are required to correctly represent the triangle boundary, we employ an 8-day composite LST (MYD11A2; the average LST of all cloud-free data in the compositing window) which minimizes cloud obstacles in determining ΔTmax and ΔTmin. However, ΔTi is derived from the daily LST (MYD11A1) to represent the LST at the satellite overpass. Thus, Eq. (18) can ultimately be expressed as Display Formula

EF=1.26ΔΔ+γ[ΔTA2maxΔTA1iΔTA2maxΔTA2min(1fveg)+fveg],(19)
where ΔTA2max and ΔTA2min are the parameters determined by a schematic plot of ΔT [MYD11A2 (8-day composite)]—EVI (8-day composite), and ΔTA1i is the day and night temperature difference obtained from MYD11A1 (daily). Display Formula
LE¯=EF(Rn¯G¯).(20)

Finally, daily LE is derived from Eq. (20), given estimated net radiation (Rn) [Eq. (9)], ground heat flux (G) [Eq. (10)] at a daily basis and evaporative fraction [Eq. (19)] with the assumption of constant EF during a day.52

We implemented the above algorithm for the entire 2005 year and based on the relevant MODIS remote-sensing products (Table 2). We evaluate various algorithm products on an instantaneous, daily or 8-day basis. In the first part of the results (Sec. 5.1), we compare MODIS-based net radiation estimates against the net radiation observations at the flux tower sites. Second, we examine the performance of the triangle method to derive the EF. Finally, we evaluate the ET results against the latent heat flux observation from the four flux tower sites.

Validation of Net Radiation

The initial net radiation model developed in the Kim and Hogue23 methodology incorporates the cloud product from MOD08 (1×1deg) using regional fitting coefficients for downward shortwave radiation under cloudy conditions; however, this study uses the cloud product from MOD06 (1×1km and 5×5km) and does not require regional calibration. Bisht and Bras48 discussed this method in detail and validated their algorithm against a range of flux towers throughout the United States. Although a semi-arid site (Desert Rock in Nevada) is included in the Bisht and Bras48 validation, the single semi-arid/dry site was not sufficient to provide strong conclusions on the performance of the net radiation model for semi-arid regions. Thus, before application of the proposed ET model, net radiation estimates derived in the current study are further evaluated against observed net radiation from the flux tower sites.

First, we present results of the comparison between estimated instantaneous downward shortwave radiation for all sky conditions against ground-measured instantaneous solar radiation. For clear sky conditions, the shortwave radiation is estimated using the PS model from Kim and Hogue,23 but for shortwave radiation under cloudy sky conditions, we calculate shortwave radiation using the Eq. (4) suggested by Bisht and Bras,48 along with a regional parameterized scheme from the Kim and Hogue23 methodology. In the comparison during clear sky days, shortwave radiation estimates have bias errors that range from 1 to 5%, and correlations are on the range of 0.94 to 0.98 across all study sites (Fig. 2). Although the Bisht and Bras48 shortwave radiation scheme has significant advantages in that it is independent from ground-based measurements, the scatter plot between the estimated and observed quantities (Fig. 2) indicates that the result from the Bisht and Bras48 model tends to be more scattered than the shortwave radiation estimates derived from Kim and Hogue23 across all study sites. Bias errors in instantaneous shortwave radiation obtained from the Kim and Hogue23 method range between 12 and 17% of the mean value, while errors from the Bisht and Bras48 method were slightly higher, between 17 and 22% across the study sites. Root mean square errors (RMSEs) reported in Kim and Hogue23 were also lower than those reported in Bisht and Bras.48 Our results are also similar to a previous study, in which GOES-based solar insolation estimates have errors ranging from 5 to 10% and 15 to 30% for clear and all sky conditions, respectively.53

Graphic Jump LocationF2 :

Comparison of downward shortwave radiation observations and the MODIS-derived shortwave radiation using Bisht and Bras,48 referred as BB10, and Kim and Hogue,23 referred as KH08, at four flux tower sites in Arizona. Both the correlation coefficient (R) and root mean square error (RMSE) are shown.

The downward longwave radiation scheme,41 as mentioned in Kim and Hogue,23 showed larger uncertainty in semi-arid regions.23 Therefore, we compared several existing algorithms that use more readily available meteorological observations such as air temperature and humidity.54 Although these simpler algorithms may have more uncertainty relative to the more complex methods, they may be more useful for a variety of disciplinary satellite applications. Eight longwave parameterization schemes4047 were assessed at the Santa Rita site, which is the only site that provides longwave radiation observations.4047 The Brutsaert,40 Garratt,45 and Keding46 algorithms show the best overall performance in terms of correlation and RMSE (Fig. 3). Among them, the Brutsaert40 algorithm showed the best performance with a correlation of 0.94 and an RMSE of 30W/m2. Due to the limited observed longwave radiation data available, we are not able to perform a more comprehensive assessment of downscaling algorithms across diverse biomes including deserts, semi-arid rangelands, and riparian corridors. However, based on the results noted above, we incorporate the Brutsaert40 longwave radiation scheme into our net radiation model for our study area.

Graphic Jump LocationF3 :

Comparison of tested downward longwave radiation schemes at the Santa Rita site.4047 Both the correlation coefficient (R) and RMSE are shown.

Comparison between the observed net radiation and the modeled net radiation at satellite overpass time [Fig. 4(a)] tends to be scattered due to the cumulative effect of bias in both the shortwave and longwave radiation schemes. However, this bias was mostly related to the underestimation in the shortwave radiation. Although not shown here, the MODIS-derived upward longwave radiation is slightly overestimated (up to 52W/m2 for the Santa Rita site) due to the positive bias of the LST 1-km product MODIS (MYD11). Previous studies note that the generalized split-window LST algorithm used in the 1-km MODIS LST retrieval overestimates LST in semi-arid and arid regions.51,55,56 The overall RMSE of net radiation derived from the MODIS products ranges from 153 to 234W/m2. The largest bias is observed at the Charleston site, where the modeled net radiation is, on average, 197W/m2 underestimated compared with the observations during the MODIS overpass time. Derived daily scatter plots for each site from the sinusoidal model are shown in Fig. 4(b). Results from the daily net radiation comparison are almost similar or slightly improved relative to the instantaneous net radiation comparisons in terms of %Bias with 32 and 20% at the Lewis Spring and Kendall sites, respectively, in the daily comparison. The RMSE error at Kendall is 69W/m2 and 97 to 116W/m2 for the remaining flux tower sites. Correlation coefficients from 0.65 to 0.69 are found at the daily scale. Due to the lack of detailed spatial information around the flux towers, it is difficult to explicitly identify sources of errors and uncertainties. However, we note that significant uncertainty may result from the obvious scale differences between the MODIS-based value and in situ data. Several model inputs (i.e., MOD06 cloud fraction and surface temperature and MOD11 surface temperature) are relatively coarse (1 to 5 km), which ultimately adds uncertainty and reduces accuracy when evaluating at the flux tower scale. In addition, topographic effects and uncertainties from instrumental measurements and remote-sensing retrieval are also a factor.57 Relative to results from the semi-arid site (Desert Rock, Nevada) in the Bisht and Bras48 study, which reported an RMSE of 41W/m2 and a correlation (R) of 0.88, results at our study sites show relatively higher errors. However, a recent study by Tang et al. noted an RMSE of 57 and 84W/m2 and a correlation of 0.50 and 0.34 for MODIS-based net radiation model over the Kendall and Audubon Research Ranch flux tower sites during satellite overpass time in clear days.57

Graphic Jump LocationF4 :

Comparison of net radiation observations and the MODIS-derived net radiation (a) at satellite overpassing time and (b) daily estimate. Correlation coefficient (R), RMSE, bias, and mean of observation (mean) are also presented.

Sensitivity of Evaporative Fraction Model

Before validation of the final ET product, we evaluate the sensitivity of the triangle relationship (i.e., LST and EVI responses). Since the trapezoidal domain may be different during the growing (mid-May to early September) or senescent seasons, we separate results and explore triangle development during the two distinct periods. Figure 5(a) shows the relationship between the EF and the difference between daytime and nighttime LST from the Aqua platform, and Fig. 5(b) shows the correlation between the EF and EVI at individual sites. The observed EF is calculated from observed (flux tower) net radiation, ground heat flux, and latent heat flux through the relationship EF=LE/(RnG). First, correlations between EF versus LST and EF versus EVI are very different during the growing and nongrowing seasons. For all sites, the correlations are significantly higher during the growing season than during the nongrowing season, due to the fact that EF is closely related to root-zone soil moisture during conditions of extreme soil water deficit (nongrowing season), which is difficult to represent only through either LST or EVI variables. But, in general, LST appears better correlated with EF than EVI during the nongrowing season, because the influence of soil moisture is more dominant when vegetation is small. However, during the growing season, the correlation of both LST and EVI with EF are reasonably high, ranging in magnitude from 0.63 to 0.86 and 0.50 to 0.86, respectively, indicating that LST and EVI can be used to parameterize EF. More importantly, it is clear that the slopes in the scatter plots (Fig. 5) are different for the upland grass (Kendall), savanna (Santa Rita), and riparian sites (i.e., Charleston and Lewis Spring), indicating that the influence of LST and EVI relative to EF varies from site to site. Given these results, we note that a single-triangular LST-EVI relationship may not be satisfactory for all of our study sites. Since riparian vegetation is distinctly different from the upland vegetation (density and species and better access to water compared with upland vegetation), the dominant factors affecting ET patterns are different. Ultimately, precipitation may be a more appropriate variable for predicting ET at these dry upland sites, as noted in Nagler et al.58

Graphic Jump LocationF5 :

Regression equations for (a) evaporative fraction (EF) versus land surface temperature (LST) and (b) evaporative fraction (EF) versus EVI at the four flux tower sites. Data points are MODIS-derived 8-day composite LST (MYD11A2) and 8-day composite EVI (MOD13Q1 and MYD13Q1). Correlation coefficient (R) and regression line (gray solid and dashed line) are presented.

We evaluate the EF derived from the triangle method against observed (flux tower) EF [Fig. 6(a)]. The EF comparison demonstrates that the performance of the proposed triangle model in riparian areas is reasonably good, but the comparison is not as good at the semi-arid rangelands or upland sites. Particularly, the lowest correlation occurs for the Kendall site, which has the smallest vegetation cover and limited water supply, which ultimately results in extremely low EF even during growing season. Comparison at the Savanna mesquite (Santa Rita) site shows a relatively good correlation (0.65) during the growing season. At the riparian sites, the correlation varies from 0.62 to 0.80 and mean absolute error (MAE) varies from 0.10 to 0.15 (unitless) during the growing season. However, the triangle-derived EF during the nongrowing season is generally overestimated at all sites. Similarly, results reported by Tang et al. showed that the triangle method significantly overestimates EF when compared with observed EF from in situ measurements from May to September during 2004 to 2007 at two flux tower sites in Arizona.57 Another study by Wang et al. showed that the correlation coefficients of EF ranged from 0.081 to 0.89 at a range of sites over the southern Great Plains with data collected from May to October in 2004.17 They also presented that EF is nonuniform regardless of the vegetation uniformity having low correlations when the mean soil moisture content is low.

Graphic Jump LocationF6 :

(a) Scatter plot of EF for the 8-day composite and (b) the latent heat (LE) at the daily scale at four flux sites. The correlation coefficient (R), RMSE, mean absolute error (MAE), and mean of observation (mean) are also presented for growing and nongrowing seasons.

Validation of Evapotranspiration

Based on the net radiation estimates (Sec. 5.1) and the triangle-derived EF (Sec. 5.2), we calculate the actual ET (LE) and validate against the ground-based LE at the four study sites [Fig. 6(b)]. Evaluating our MODIS-derived net radiation results against observed net radiation from flux towers for all sites, we hypothesize that the MODIS-derived net radiation contributes significant uncertainty or bias in the ET model. Compared with the result from the EF [Fig. 6(a)], the estimated daily LE is slightly improved for all sites, likely due to the increased number of days used in the comparison. Particularly, Lewis Spring has the highest correlation (0.82) and low RMSE (36.67W/m2). Results from Charleston also show similar values (0.81, 51W/m2 and 55W/m2 for correlation, RMSE, and MAE, respectively). Overall, the LE rates for the riparian areas are consistently higher and show lower variability, reflecting the typical features of riparian LE and the influence of a consistent supply of soil water. However, we note that underestimation of LE in riparian sites during the growing season may be related to the inclusion of LE from dry, nonriparian portions within the MODIS pixels. Unlike the riparian sites, the LE at the upland sites tends to be more strongly tied to precipitation events. During the summer monsoon season, the variability of LE is especially dependent on precipitation. We also observe that the LE is moderately correlated with precipitation (correlations of 0.58 and 0.63 for Kendall and Santa Rita sites, respectively). In addition, both of the upland sites (Santa Rita and Kendall) are sparsely vegetated and contain more bare soil exposure. Hence, the patchy distribution of vegetation at the upland sites likely causes a higher level of noise due to the scale difference between the satellite-derived VI (250m) and the flux tower (footprint).

It is important to note the additional degree of uncertainty in this analysis due to the accuracy of eddy covariance systems and our comparison to these estimates. From previous work in this region,59 LE measurements from eddy covariance tower are noted to be systematically underestimated due to the lack of energy balance closure.30 In the current study, we applied a Bowen ratio correction method (BR method), in which the available energy is repartitioned into H and LE by conserving their ratio to enforce closure in the energy components.57 We observed that the correct LE is consistently greater than the observed LE (MAE: 6 to 7W/m2 for Kendall and Santa Rita). These differences were much larger at both riparian sites showing about RMSE of 30W/m2 and MAE of 19 to 22W/m2, which is about 18 to 26% of the mean. Our results are similar to previous studies,30 noting that LE is systematically underestimated by 17 to 27% from the lack of energy balance closure at these sites (i.e., Kendall and Santa Rita in Arizona).59 Also, due to the long, narrow strips of vegetation along the riparian corridor, flux measurements at the riparian sites may not meet the required fetch (i.e., upwind area), which can be a significant source of error.60 Performance of the proposed triangle model for an 8-day time series is highlighted in Fig. 7. The results are significantly improved over the daily results. The correlation is 0.72 on average for all four sites, and the RMSE is also decreased. The multiple relationships between ET, LST, and EVI are complex in semi-arid regions. However, the various relationships may be categorized depending on the land surface type (upland versus riparian). Water use efficiency is also variable depending on plant species, as seen at the grass-dominated Kendall site and the Mesquite Santa Rita site.

Graphic Jump LocationF7 :

Eight-day time series (a) and scatter plots (b) of observed LE and the MODIS-derived latent heat (LE) at four flux tower sites along with precipitation (mm/day). The correlation (R), RMSE, MAE, and mean of observation (mean) are shown.

In the current study, a remote-sensing based scheme, based only on MODIS products, is developed for the assessment of daily ET distribution for all sky conditions at regional scales with a focus on performance in semi-arid regions. The approach is based on the triangle method, which relies on two remotely sensed inputs (EVI and LST) and a previously developed net radiation (Rn) model which uses a range of MODIS products. A major improvement over previous studies is that we attempt to utilize only MODIS satellite products for a simple and direct estimate of actual ET using limited ancillary data. Ground-based data is only used for regional coefficients developed in the previous PET model from Kim and Hogue.12,17,23 The proposed model shows the potential for operational monitoring of ET, since the algorithm is not associated with ground-based meteorological data, and all required datasets can be obtained from freely available MODIS platforms. Therefore, we advocate that the model is useful for a range of applications, especially where in situ data are not readily available.

The proposed approach was initially tested in southern Arizona. Flux tower data was available in the riparian corridor of the upper San Pedro River Basin (i.e., Charleston and Lewis Spring Sacaton) in the rangeland Walnut Gulch watershed (Kendall) and in a natural upland region (Santa Rita). The net radiation model used in this study was shown to systematically underestimate the surface net radiation, resulting in underestimation as large as 122W/m2 during satellite overpass. These distinct discrepancies might result from the relatively unsatisfactory performance of shortwave radiation scheme using Bisht and Bras.48

The observed LE rates are significantly different between the riparian and rangeland upland sites with Charleston and Lewis Spring having approximately 144 to 179W/m2 of LE, while stands of shrubs on Santa Rita and grass on Kendall have much lower ET rates of 36 to 76W/m2 during the growing season (mid May to early September). Although the proposed ET model provides reasonable estimates of LE in the riparian sites, it is obvious that heterogeneity within the MODIS pixel (including both riparian and nonriparian) contributes to the underestimation of ET during the growing season. At the rangeland sites, LE estimation is slightly worse. Our initial results demonstrate some limitations of the proposed triangle method. Since the model presumes a linear variation in EF across the triangular domain of LST/EVI space, the assumption may result in larger errors under water-stressed conditions. Because deeper root zone soil moisture is a controlling factor in severely water stress areas, the triangular relationship between LST and EVI may not represent this soil moisture that is available to plants. Further, the relationship between EVI and root-zone soil moisture is dependent on the vegetation species and local climatic patterns. Alternatively, a modified model by Stisen et al., where decomposition is regarded as nonlinear, may provide improved estimation of EF by allowing some degree of water stress and heterogeneity within the pixel (i.e., mixture of vegetation and bare soil).15

We note that there are still significant challenges for large-scale, regional ET modeling in semi-arid regions, and ongoing work is aimed at reducing this uncertainty. However, results from the current study indicate that MODIS-based algorithms can be utilized to track temporal and spatial changes in LE across heterogeneous, semi-arid domains. The independence from ancillary data and near real-time applicability makes the proposed MODIS-based ET model well suited as a near real-time operational model, which if we advocate can lead to improved water management in water-stressed regions.

Financial support for this work was partially supported by a NASA Earth System Science (ESS) Fellowship (#NNX07AO53H) as well as a NSF CAREER award (#EAR0846662).

Andersen  F. H., “Hydrological modelling in a semi-arid area using remote sensing data,” Doctoral Thesis, Department of Geography and Geology, University of Copenhagen, Copenhagen (2008).
Moran  M. S., Jackson  R. D., “Assessing the spatial distribution of evapotranspiration using remote sensed inputs,” J. Environ. Qual.. 20, , 725 –737 (1991), CrossRef. 0047-2425 
Nishida  K. et al., “An operational remote sensing algorithm of land surface,” J. Geophys. Res.. 108, (D9 ), 4270  (2003), CrossRef. 0148-0227 
Li  Z. L. et al., “A review of current methodologies for regional evapotranspiration estimation from remotely sensed data,” Sensors. 9, , 3801 –3853 (2009), CrossRef. 0746-9462 
Kalma  J. D., Jupp  D. L. B., “Estimating evaporation from pasture using infrared thermometry: evaluation of a one-layer resistance model,” Agr. Forest Meteorol.. 51, , 223 –246 (1990), CrossRef. 0168-1923 
Bastiaanssen  W. G. M. et al., “A remote sensing surface energy balance algorithm for Land (SEBAL), Part 1: formulation,” J. Hydrol.. 212–213, , 198 –212 (1998), CrossRef. 0022-1694 
Bastiaanssen  W. G. M., “SEBAL-based sensible and latent heat fluxes in the irrigated Gediz Basin, Turkey,” J. Hydrol.. 229, (1–2 ), 87 –100 (2000), CrossRef. 0022-1694 
Allen  R. G., Tasumi  M., Trezza  R., “Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC) model,” ASCE J. Irrig. Drain. Eng.. 133, (4 ), 380 –394 (2007), CrossRef. 0733-9437 
Norman  J. M., Kustas  W. B., Humes  K. S., “Two-source approach for estimating soil and vegetation energy fluxes in observations of directional radiometric surface temperature,” Agr. Forest Meteorol.. 77, , 263 –293 (1995), CrossRef. 0168-1923 
Gillies  R. R., Carlson  T. N., “Thermal remote sensing of surface soil water content with partial vegetation cower for incorporation into mesoscale prediction models,” J. Appl. Meteorol.. 34, , 745 –756 (1995), CrossRef. 0894-8763 
Nishida  K. et al., “Development of an evapotranspiration index from Aqua/MODIS for monitoring surface moisture status,” IEEE Trans. Geosci. Remote Sens.. 41, (2 ), 493 –501 (2003), CrossRef. 0196-2892 
Jiang  L., Islam  S., “Estimation of surface evaporation map over southern Great Plains using remote sensing data,” Water Resour. Res.. 37, (2 ), 329 –340 (2001), CrossRef. 0043-1397 
Long  D., Singh  V. P., “A two-source trapezoid model for evapotranspiration (TTME) from satellite imagery,” Remote Sens. Environ.. 121, , 370 –388 (2012), CrossRef. 0034-4257 
Long  D., Singh  V. P., Scanlon  B. R., “Deriving theoretical boundaries to address scale dependencies of triangle models for evapotranspiration estimation,” J. Geophys. Res.—Atmos.. 117, , D07188  (2012), CrossRef. 0148-0227 
Stisen  S. et al., “Combining the triangle method with thermal inertia to estimate regional evapotranspiration—applied to MSG-SEVIRI data in the Senegal River basin,” Remote Sens. Environ.. 112, , 1242 –1255 (2008), CrossRef. 0034-4257 
Venturini  V. et al., “Comparison of EFs estimated from AVHRR and MODIS sensors over South Florida,” Remote Sens. Environ.. 93, , 77 –86 (2004), CrossRef. 0034-4257 
Wang  K., Li  Z., Cribb  M., “Estimation of evaporative fraction from a combination of day and night land surface temperatures and NDVI: a new method to determine the Priestley–Taylor parameter,” Remote Sens. Environ.. 102, , 293 –305 (2006), CrossRef. 0034-4257 
Jiang  L., Islam  S., Carlson  T., “Towards the uncertainty of surface fluxes estimation and implications for using a simplified approach with remote sensing data,” Can. J. Remote Sens.. 30, (5 ), 769 –787 (2004), CrossRef. 0703-8992 
Jiang  L., Islam  S., “An intercomparison of regional latent heat flux estimation using remote sensing data,” Int. J. Remote Sens.. 24, (11 ), 2221 –2236 (2003), CrossRef. 0143-1161 
Jiang  L., Islam  S., “A methodology for estimation of surface evapotranspiration overlarge areas using remote sensing observations,” Geophys. Res. Lett.. 26, (17 ), 2773 –2776 (1999), CrossRef. 0094-8276 
Tang  R., Li  Z. L., Tang  B., “An application of the Ts-VI triangle method with enhanced edges determination for evapotranspiration estimation from MODIS data in arid and semi arid regions: implementation and validation,” Remote Sens. Environ.. 114, , 540 –551 (2010), CrossRef. 0034-4257 
Xiangsong  Z. et al., “Estimating evapotranspiration using triangle method with topographic correction from MODIS data in Taihu Basin, China,” in  Int. Conf. Wireless Communications, Networking and Mobile Computing , pp. 1 –5 (2010).
Kim  J., Hogue  T. S., “Evaluation of a MODIS-based potential evapotranspiration product at the point-scale,” J. Hydrometeorol.. 9, , 444 –460 (2008), CrossRef. 1525-755X 
Polyakov  V. et al., “Long-term runoff and sediment yields from small semi-arid watersheds in southern Arizona,” Water Resour. Res.. 46, , W09512  (2010), CrossRef. 0043-1397 
Higgins  R. W., Mo  K. C., Yao  Y., “Interannual variability of the U.S. summer precipitation regime with emphasis on the southwestern monsoon,” J. Clim.. 11, , 2582 –2606 (1998), CrossRef. 0894-8755 
Zhu  C., Lettenmaier  D., Cavazos  T., “Role of antecedent land surface conditions on North American Monsoon Rainfall variability,” J. Clim.. 18, , 3104 –3121 (2005), CrossRef. 0894-8755 
Nagler  P. L. et al., “Predicting riparian evapotranspiration from MODIS vegetation indices and meteorological data,” Remote Sens. Environ.. 94, (1 ), 17 –30 (2005), CrossRef. 0034-4257 
Emmerich  W. E., “Carbon dioxide fluxes in a semiarid environment with high carbonate soils,” Agr. Forest Meteorol.. 116, , 91 –102 (2003), CrossRef. 0168-1923 
National Weather Service Forecast Office “Rainfall across Southeast Arizona during 2005 Monsoon,” http://www.wrh.noaa.gov/twc/monsoon/season/2005monsoon.php#graph ( January 2012).
Scott  R. L. et al., “Carbon dioxide exchange in a semidesert grassland through drought induced vegetation change,” J. Geophys. Res.. 115, , G03026  (2010). CrossRef. 0148-0227 
Potts  D. L. et al., “Sensitivity of mesquite shrubland CO2 exchange to precipitation on contrasting landscape settings,” Ecology. 89, , 2900 –2910 (2008), CrossRef. 0012-9658 
Scott  R. et al., “Inter-annual and seasonal variation in fluxes of water and carbon dioxide from a riparian woodland ecosystem,” Agr. Forest Meteorol.. 122, , 65 –84 (2004), CrossRef. 0168-1923 
Hsiegh  C. I., Katul  G. G., Chi  T. W., “An approximate analytical model for footprint estimation of scalar fluxes in thermally stratified atmospheric flows,” Adv. Water Resour.. 23, , 765 –772 (2000), CrossRef. 0309-1708 
Justice  C. O. et al., “The MODIS fire products,” Remote Sens. Environ.. 83, , 244 –262 (2002), CrossRef. 0034-4257 
Van Leeuwen  W. J. D., Huete  A., Laing  T. W., “MODIS vegetation index compositing approach: A prototype with AVHRR data,” Remote Sens. Environ.. 69, (3 ), 264 –280 (1999).CrossRef
Wessels  K. J., Bachoo  A. K., Archibald  S., “Influence of composite period and date of observation on phenological metrics extracted from MODIS data,” in  33rd Int. Symp. Remote Sensing of Environment ,  ISRSE ,  Stresa, Italy  (2009).
Huete  A. et al., “Overview of the radiometric and biophysical performance of the MODIS vegetation indices,” Remote Sens. Environ.. 83, , 195 –213 (2002), CrossRef. 0034-4257 
Yang  X. et al., “Modified triangle method to estimate soil moisture status with moderate resolution imaging spectroradiometer (MODIS) products,” in  ISPRS Congress Beijing 2008 , pp. 555 –560 (2008).
Paulescu  M., Schlett  Z., “Performance assessment of global solar irradiation models under Romanian climate,” Renew. Energy. 29, , 767 –777 (2004), CrossRef. 0960-1481 
Brutsaert  W., “On a derivable formula for long-wave radiation from clear skies,” Water Resour. Res.. 11, , 742 –744 (1975), CrossRef. 0043-1397 
Prata  A. J., “A new long-wave formula for estimating downward clear-sky radiation at the surface,” Quart. J. R. Meteorol. Soc.. 122, , 1127 –1151 (1996), CrossRef. 0035-9009 
Dilley  A. C., O’Brien  D. M., “Estimating downward clear sky long-wave irradiance at the surface from screen temperature and precipitable water,” Quart. J. R. Meteorol. Soc.. 124, , 1391 –1401 (1998), CrossRef. 0035-9009 
Angstrom  A., “A study of the radiation of the atmosphere,” Smithson. Misc. Collect.. 65, , 1 –159 (1918). 0096-8749 
Idso  S. B., “A set of equations for full spectrum and 8 to 14 mm and 10.5 to 12.5 mm thermal radiation from cloudless skies,” Water Resour. Res.. 17, , 295 –304 (1981), CrossRef. 0043-1397 
Garratt  J. A., “Extreme maximum land surface temperatures,” J. Appl. Meteorol. Climatol.. 31, , 1096 –1105 (1992), CrossRef. 1558-8424 
Keding  I., “Klimatologische Untersuchung über die atmosphärische Gegenstrahlung und Vergleich von Berechnungsverfahren anhand langjähriger Messungen im Oberrheintal,” Ber. Dtsch. Wetterd. 178, , 72  (1989).
Satterlund  D. R., “An improved equation for estimating longwave radiation from the atmosphere,” Water Resour. Res.. 15, , 1649 –1650 (1979), CrossRef. 0043-1397 
Bisht  G., Bras  R. L., “Estimation of net radiation from the MODIS data under all sky conditions—Southern Great Plains case study,” Remote Sens. Environ.. 114, (7 ), 1522 –1534 (2010), CrossRef. 0034-4257 
Bisht  G. et al., “Estimation of the net radiation using MODIS (Moderate Resolution Imaging Spectroradiometer) data for clear sky days,” Remote Sens. Environ.. 97, , 52 –67 (2005), CrossRef. 0034-4257 
Liu  W. et al., “Estimation of the land surface instantaneous net radiation and its diurnal cycle integrating multi-source remote sensing data under clear sky,” in  IEEE International Goescience and Remote Sensing Symposium (IGARSS) , pp. 2109 –2112 (2011).
Wan  Z. et al., “Validation of the land surface temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data,” Remote Sens. Environ.. 83, , 163 –180 (2002), CrossRef. 0034-4257 
Allen  R. G. et al., “Satellite-based ET estimation in agriculture using SEBAL and METRIC,” Hydrol. Process.. 25, (26 ), 4011 –4027 (2011), CrossRef. 1099-1085 
Pinker  R. T., Frouin  R., Li  Z., “A review of satellite methods to derive surface shortwave irradiance,” Remote Sens. Environ.. 51, , 105 –124 (1995), CrossRef. 0034-4257 
Flerchinger  G. N. et al., “Comparison of algorithms for incoming atmospheric longwave radiation,” Water Resour. Res.. 45, , W03423  (2009), CrossRef. 0043-1397 
Wan  Z. et al., “Quality assessment and validation of the modis global land surface temperature,” Int. J. Remote Sens.. 25, , 261 –274 (2004), CrossRef. 0143-1161 
Hulley  G. C., Hook  S. J., “The North American ASTER Land Surface Emissivity Database (NAALSED) version 2.0,” Remote Sens. Environ.. 113, (9 ),1967 –1975 (2009), CrossRef. 0034-4257 
Tang  R., Li  Z. L., Chen  K. S., “Validating MODIS-derived land surface evapotranspiration with in situ measurements at two Ameriflux sites in a semiarid region,” J. Geophys. Res.. 116, , D04106  (2011), CrossRef. 0148-0227 
Nagler  P. L. et al., “Relationship between evapotranspiration and precipitation pulses in a semiarid rangeland estimated by moisture flux towers and MODIS vegetation indices,” J. Arid Environ.. 70, , 443 –462 (2007), CrossRef. 0140-1963 
Scott  R., “Using watershed water balance to evaluate the accuracy of eddy covariance evaporation measurements for three semiarid ecosystems,” Agr. Forest Meteorol.. 150, (2 ), 219 –225 (2010), CrossRef. 0168-1923 
Goodrich  D. C., Scott  R., Qi  J., “Seasonal estimates of riparian evapotranspiration using remote and in situ measurements,” Agr. Forest Meteorol.. 105, , 281 –309 (2000), CrossRef. 0168-1923 

Jongyoun Kim received her MS degree in atmospheric and oceanic sciences and her MS and PhD degrees in civil and environmental engineering from the University of California, Los Angeles (UCLA). She is currently a project scientist with the Department of Botany and Plant Sciences, UC Riverside. She specializes in the development and application of remote sensing products within hydrologic modeling systems and in the evaluation of land–atmosphere interactions at a range of temporal and spatial scales.

Terri S. Hogue is an associate professor in the Department of Civil and Environmental Engineering at the Colorado School of Mines. She received her BS from the University of Wisconsin-Eau Claire, and MS and PhD degrees from the Department of Hydrology and Water Resources at the University of Arizona. Her research centers on understanding hydrologic and land surface processes, with much of her work focused in semi-arid regions. Projects include wildfire impacts, urbanization and ecosystem dynamics, and hydrologic response to climate change.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Jongyoun Kim and Terri S. Hogue
"Evaluation of a MODIS triangle-based evapotranspiration algorithm for semi-arid regions", J. Appl. Remote Sens. 7(1), 073493 (Oct 28, 2013). ; http://dx.doi.org/10.1117/1.JRS.7.073493


Figures

Graphic Jump LocationF1 :

Study area in Arizona with a Landsat scene (path 35 and row 38; false-color display of the 30 m Landsat5 VNIR imagery) including the four flux towers (Kendall, Charleston, Lewis Spring, and Santa Rita) marked as triangles.

Graphic Jump LocationF3 :

Comparison of tested downward longwave radiation schemes at the Santa Rita site.4047 Both the correlation coefficient (R) and RMSE are shown.

Graphic Jump LocationF2 :

Comparison of downward shortwave radiation observations and the MODIS-derived shortwave radiation using Bisht and Bras,48 referred as BB10, and Kim and Hogue,23 referred as KH08, at four flux tower sites in Arizona. Both the correlation coefficient (R) and root mean square error (RMSE) are shown.

Graphic Jump LocationF4 :

Comparison of net radiation observations and the MODIS-derived net radiation (a) at satellite overpassing time and (b) daily estimate. Correlation coefficient (R), RMSE, bias, and mean of observation (mean) are also presented.

Graphic Jump LocationF6 :

(a) Scatter plot of EF for the 8-day composite and (b) the latent heat (LE) at the daily scale at four flux sites. The correlation coefficient (R), RMSE, mean absolute error (MAE), and mean of observation (mean) are also presented for growing and nongrowing seasons.

Graphic Jump LocationF7 :

Eight-day time series (a) and scatter plots (b) of observed LE and the MODIS-derived latent heat (LE) at four flux tower sites along with precipitation (mm/day). The correlation (R), RMSE, MAE, and mean of observation (mean) are shown.

Graphic Jump LocationF5 :

Regression equations for (a) evaporative fraction (EF) versus land surface temperature (LST) and (b) evaporative fraction (EF) versus EVI at the four flux tower sites. Data points are MODIS-derived 8-day composite LST (MYD11A2) and 8-day composite EVI (MOD13Q1 and MYD13Q1). Correlation coefficient (R) and regression line (gray solid and dashed line) are presented.

Tables

Table Grahic Jump Location
Table 1Location and characteristics of the four tower sites utilized for the study and single-pixel extraction.
Table Footer Note*Precipitation and reference crop evaporation (AZMET standard) are the annual averaged value for year 2005.
Table Grahic Jump Location
Table 2Summary of MODIS products used in this study with relevant spatial and temporal information.

References

Andersen  F. H., “Hydrological modelling in a semi-arid area using remote sensing data,” Doctoral Thesis, Department of Geography and Geology, University of Copenhagen, Copenhagen (2008).
Moran  M. S., Jackson  R. D., “Assessing the spatial distribution of evapotranspiration using remote sensed inputs,” J. Environ. Qual.. 20, , 725 –737 (1991), CrossRef. 0047-2425 
Nishida  K. et al., “An operational remote sensing algorithm of land surface,” J. Geophys. Res.. 108, (D9 ), 4270  (2003), CrossRef. 0148-0227 
Li  Z. L. et al., “A review of current methodologies for regional evapotranspiration estimation from remotely sensed data,” Sensors. 9, , 3801 –3853 (2009), CrossRef. 0746-9462 
Kalma  J. D., Jupp  D. L. B., “Estimating evaporation from pasture using infrared thermometry: evaluation of a one-layer resistance model,” Agr. Forest Meteorol.. 51, , 223 –246 (1990), CrossRef. 0168-1923 
Bastiaanssen  W. G. M. et al., “A remote sensing surface energy balance algorithm for Land (SEBAL), Part 1: formulation,” J. Hydrol.. 212–213, , 198 –212 (1998), CrossRef. 0022-1694 
Bastiaanssen  W. G. M., “SEBAL-based sensible and latent heat fluxes in the irrigated Gediz Basin, Turkey,” J. Hydrol.. 229, (1–2 ), 87 –100 (2000), CrossRef. 0022-1694 
Allen  R. G., Tasumi  M., Trezza  R., “Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC) model,” ASCE J. Irrig. Drain. Eng.. 133, (4 ), 380 –394 (2007), CrossRef. 0733-9437 
Norman  J. M., Kustas  W. B., Humes  K. S., “Two-source approach for estimating soil and vegetation energy fluxes in observations of directional radiometric surface temperature,” Agr. Forest Meteorol.. 77, , 263 –293 (1995), CrossRef. 0168-1923 
Gillies  R. R., Carlson  T. N., “Thermal remote sensing of surface soil water content with partial vegetation cower for incorporation into mesoscale prediction models,” J. Appl. Meteorol.. 34, , 745 –756 (1995), CrossRef. 0894-8763 
Nishida  K. et al., “Development of an evapotranspiration index from Aqua/MODIS for monitoring surface moisture status,” IEEE Trans. Geosci. Remote Sens.. 41, (2 ), 493 –501 (2003), CrossRef. 0196-2892 
Jiang  L., Islam  S., “Estimation of surface evaporation map over southern Great Plains using remote sensing data,” Water Resour. Res.. 37, (2 ), 329 –340 (2001), CrossRef. 0043-1397 
Long  D., Singh  V. P., “A two-source trapezoid model for evapotranspiration (TTME) from satellite imagery,” Remote Sens. Environ.. 121, , 370 –388 (2012), CrossRef. 0034-4257 
Long  D., Singh  V. P., Scanlon  B. R., “Deriving theoretical boundaries to address scale dependencies of triangle models for evapotranspiration estimation,” J. Geophys. Res.—Atmos.. 117, , D07188  (2012), CrossRef. 0148-0227 
Stisen  S. et al., “Combining the triangle method with thermal inertia to estimate regional evapotranspiration—applied to MSG-SEVIRI data in the Senegal River basin,” Remote Sens. Environ.. 112, , 1242 –1255 (2008), CrossRef. 0034-4257 
Venturini  V. et al., “Comparison of EFs estimated from AVHRR and MODIS sensors over South Florida,” Remote Sens. Environ.. 93, , 77 –86 (2004), CrossRef. 0034-4257 
Wang  K., Li  Z., Cribb  M., “Estimation of evaporative fraction from a combination of day and night land surface temperatures and NDVI: a new method to determine the Priestley–Taylor parameter,” Remote Sens. Environ.. 102, , 293 –305 (2006), CrossRef. 0034-4257 
Jiang  L., Islam  S., Carlson  T., “Towards the uncertainty of surface fluxes estimation and implications for using a simplified approach with remote sensing data,” Can. J. Remote Sens.. 30, (5 ), 769 –787 (2004), CrossRef. 0703-8992 
Jiang  L., Islam  S., “An intercomparison of regional latent heat flux estimation using remote sensing data,” Int. J. Remote Sens.. 24, (11 ), 2221 –2236 (2003), CrossRef. 0143-1161 
Jiang  L., Islam  S., “A methodology for estimation of surface evapotranspiration overlarge areas using remote sensing observations,” Geophys. Res. Lett.. 26, (17 ), 2773 –2776 (1999), CrossRef. 0094-8276 
Tang  R., Li  Z. L., Tang  B., “An application of the Ts-VI triangle method with enhanced edges determination for evapotranspiration estimation from MODIS data in arid and semi arid regions: implementation and validation,” Remote Sens. Environ.. 114, , 540 –551 (2010), CrossRef. 0034-4257 
Xiangsong  Z. et al., “Estimating evapotranspiration using triangle method with topographic correction from MODIS data in Taihu Basin, China,” in  Int. Conf. Wireless Communications, Networking and Mobile Computing , pp. 1 –5 (2010).
Kim  J., Hogue  T. S., “Evaluation of a MODIS-based potential evapotranspiration product at the point-scale,” J. Hydrometeorol.. 9, , 444 –460 (2008), CrossRef. 1525-755X 
Polyakov  V. et al., “Long-term runoff and sediment yields from small semi-arid watersheds in southern Arizona,” Water Resour. Res.. 46, , W09512  (2010), CrossRef. 0043-1397 
Higgins  R. W., Mo  K. C., Yao  Y., “Interannual variability of the U.S. summer precipitation regime with emphasis on the southwestern monsoon,” J. Clim.. 11, , 2582 –2606 (1998), CrossRef. 0894-8755 
Zhu  C., Lettenmaier  D., Cavazos  T., “Role of antecedent land surface conditions on North American Monsoon Rainfall variability,” J. Clim.. 18, , 3104 –3121 (2005), CrossRef. 0894-8755 
Nagler  P. L. et al., “Predicting riparian evapotranspiration from MODIS vegetation indices and meteorological data,” Remote Sens. Environ.. 94, (1 ), 17 –30 (2005), CrossRef. 0034-4257 
Emmerich  W. E., “Carbon dioxide fluxes in a semiarid environment with high carbonate soils,” Agr. Forest Meteorol.. 116, , 91 –102 (2003), CrossRef. 0168-1923 
National Weather Service Forecast Office “Rainfall across Southeast Arizona during 2005 Monsoon,” http://www.wrh.noaa.gov/twc/monsoon/season/2005monsoon.php#graph ( January 2012).
Scott  R. L. et al., “Carbon dioxide exchange in a semidesert grassland through drought induced vegetation change,” J. Geophys. Res.. 115, , G03026  (2010). CrossRef. 0148-0227 
Potts  D. L. et al., “Sensitivity of mesquite shrubland CO2 exchange to precipitation on contrasting landscape settings,” Ecology. 89, , 2900 –2910 (2008), CrossRef. 0012-9658 
Scott  R. et al., “Inter-annual and seasonal variation in fluxes of water and carbon dioxide from a riparian woodland ecosystem,” Agr. Forest Meteorol.. 122, , 65 –84 (2004), CrossRef. 0168-1923 
Hsiegh  C. I., Katul  G. G., Chi  T. W., “An approximate analytical model for footprint estimation of scalar fluxes in thermally stratified atmospheric flows,” Adv. Water Resour.. 23, , 765 –772 (2000), CrossRef. 0309-1708 
Justice  C. O. et al., “The MODIS fire products,” Remote Sens. Environ.. 83, , 244 –262 (2002), CrossRef. 0034-4257 
Van Leeuwen  W. J. D., Huete  A., Laing  T. W., “MODIS vegetation index compositing approach: A prototype with AVHRR data,” Remote Sens. Environ.. 69, (3 ), 264 –280 (1999).CrossRef
Wessels  K. J., Bachoo  A. K., Archibald  S., “Influence of composite period and date of observation on phenological metrics extracted from MODIS data,” in  33rd Int. Symp. Remote Sensing of Environment ,  ISRSE ,  Stresa, Italy  (2009).
Huete  A. et al., “Overview of the radiometric and biophysical performance of the MODIS vegetation indices,” Remote Sens. Environ.. 83, , 195 –213 (2002), CrossRef. 0034-4257 
Yang  X. et al., “Modified triangle method to estimate soil moisture status with moderate resolution imaging spectroradiometer (MODIS) products,” in  ISPRS Congress Beijing 2008 , pp. 555 –560 (2008).
Paulescu  M., Schlett  Z., “Performance assessment of global solar irradiation models under Romanian climate,” Renew. Energy. 29, , 767 –777 (2004), CrossRef. 0960-1481 
Brutsaert  W., “On a derivable formula for long-wave radiation from clear skies,” Water Resour. Res.. 11, , 742 –744 (1975), CrossRef. 0043-1397 
Prata  A. J., “A new long-wave formula for estimating downward clear-sky radiation at the surface,” Quart. J. R. Meteorol. Soc.. 122, , 1127 –1151 (1996), CrossRef. 0035-9009 
Dilley  A. C., O’Brien  D. M., “Estimating downward clear sky long-wave irradiance at the surface from screen temperature and precipitable water,” Quart. J. R. Meteorol. Soc.. 124, , 1391 –1401 (1998), CrossRef. 0035-9009 
Angstrom  A., “A study of the radiation of the atmosphere,” Smithson. Misc. Collect.. 65, , 1 –159 (1918). 0096-8749 
Idso  S. B., “A set of equations for full spectrum and 8 to 14 mm and 10.5 to 12.5 mm thermal radiation from cloudless skies,” Water Resour. Res.. 17, , 295 –304 (1981), CrossRef. 0043-1397 
Garratt  J. A., “Extreme maximum land surface temperatures,” J. Appl. Meteorol. Climatol.. 31, , 1096 –1105 (1992), CrossRef. 1558-8424 
Keding  I., “Klimatologische Untersuchung über die atmosphärische Gegenstrahlung und Vergleich von Berechnungsverfahren anhand langjähriger Messungen im Oberrheintal,” Ber. Dtsch. Wetterd. 178, , 72  (1989).
Satterlund  D. R., “An improved equation for estimating longwave radiation from the atmosphere,” Water Resour. Res.. 15, , 1649 –1650 (1979), CrossRef. 0043-1397 
Bisht  G., Bras  R. L., “Estimation of net radiation from the MODIS data under all sky conditions—Southern Great Plains case study,” Remote Sens. Environ.. 114, (7 ), 1522 –1534 (2010), CrossRef. 0034-4257 
Bisht  G. et al., “Estimation of the net radiation using MODIS (Moderate Resolution Imaging Spectroradiometer) data for clear sky days,” Remote Sens. Environ.. 97, , 52 –67 (2005), CrossRef. 0034-4257 
Liu  W. et al., “Estimation of the land surface instantaneous net radiation and its diurnal cycle integrating multi-source remote sensing data under clear sky,” in  IEEE International Goescience and Remote Sensing Symposium (IGARSS) , pp. 2109 –2112 (2011).
Wan  Z. et al., “Validation of the land surface temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data,” Remote Sens. Environ.. 83, , 163 –180 (2002), CrossRef. 0034-4257 
Allen  R. G. et al., “Satellite-based ET estimation in agriculture using SEBAL and METRIC,” Hydrol. Process.. 25, (26 ), 4011 –4027 (2011), CrossRef. 1099-1085 
Pinker  R. T., Frouin  R., Li  Z., “A review of satellite methods to derive surface shortwave irradiance,” Remote Sens. Environ.. 51, , 105 –124 (1995), CrossRef. 0034-4257 
Flerchinger  G. N. et al., “Comparison of algorithms for incoming atmospheric longwave radiation,” Water Resour. Res.. 45, , W03423  (2009), CrossRef. 0043-1397 
Wan  Z. et al., “Quality assessment and validation of the modis global land surface temperature,” Int. J. Remote Sens.. 25, , 261 –274 (2004), CrossRef. 0143-1161 
Hulley  G. C., Hook  S. J., “The North American ASTER Land Surface Emissivity Database (NAALSED) version 2.0,” Remote Sens. Environ.. 113, (9 ),1967 –1975 (2009), CrossRef. 0034-4257 
Tang  R., Li  Z. L., Chen  K. S., “Validating MODIS-derived land surface evapotranspiration with in situ measurements at two Ameriflux sites in a semiarid region,” J. Geophys. Res.. 116, , D04106  (2011), CrossRef. 0148-0227 
Nagler  P. L. et al., “Relationship between evapotranspiration and precipitation pulses in a semiarid rangeland estimated by moisture flux towers and MODIS vegetation indices,” J. Arid Environ.. 70, , 443 –462 (2007), CrossRef. 0140-1963 
Scott  R., “Using watershed water balance to evaluate the accuracy of eddy covariance evaporation measurements for three semiarid ecosystems,” Agr. Forest Meteorol.. 150, (2 ), 219 –225 (2010), CrossRef. 0168-1923 
Goodrich  D. C., Scott  R., Qi  J., “Seasonal estimates of riparian evapotranspiration using remote and in situ measurements,” Agr. Forest Meteorol.. 105, , 281 –309 (2000), CrossRef. 0168-1923 

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.