Comparison of Three Operative Models for Estimating the Surface Water Deficit using ASTER Reflective and Thermal Data

Three operative models with minimum input data requirements for estimating the partition of available surface energy into sensible and latent heat flux using ASTER data have been evaluated in a semiarid area in SE Spain. The non-evaporative fraction (NEF) is proposed as an indicator of the surface water deficit. The best results were achieved with NEF estimated using the “Simplified relationship” for unstable conditions (NEFSeguin) and with the S-SEBI (Simplified Surface Energy Balance Index) model corrected for atmospheric conditions (NEFS-SEBIt,) which both produced equivalent results. However, results with a third model, NEFCarlson, that estimates the exchange coefficient for sensible heat transfer from NDVI, were unrealistic for sites with scarce vegetation cover. These results are very promising for an operative monitoring of the surface water deficit, as validation with field data shows reasonable errors, within those reported in the literature (RMSE were 0.18 and 0.11 for the NEF, and 29.12 Wm-2 and 25.97 Wm-2 for sensible heat flux, with the Seguin and S-SEBIt models, respectively).


Introduction
The relationship between ecosystem latent heat (λE) and sensible heat (H) flux is critical to quantify the surface water deficit and to understand the hydrological cycle. The law of conservation of energy states that the available energy reaching a surface is dissipated as latent heat (λE) and/or sensible heat (H), the partition of which depends mostly on water availability. Factors affecting this relationship include longer-term interactions between biogeochemical cycles, disturbances and climate, and shorter-term interactions between plant physiology and the development of the atmospheric boundary layer [1].
Remote sensing is currently the only source of data providing frequent and spatially disaggregated estimates of radiometric temperature, albedo and surrogates of vegetation cover, variables that explain most of the partition of the available energy into H and λE. Therefore, model development using inputs from remote sensing data in the solar and thermal domain is a very active subject of research [2]. Still, there is always a trade-off between model parameterization requirements and operativity that has to be carefully evaluated in each case when selecting a methodology.
A widely used water deficit indicator is the evaporative fraction: EF=λE/(H+λE) [3,4]. However, in semiarid areas, λE and therefore the evaporative fraction (EF) of natural vegetation are very low during several days, with values sometimes lower than the accuracy levels of remote sensing models for estimating evapotranspiration (<1 mm for λE) [2,5].
We, therefore, propose using the non-evaporative (NEF) fraction to evaluate the surface water deficit, defined as: where Rn = net radiation; G = soil heat flux; and Rn-G = available energy.
In semiarid areas, the NEF should have a wider range of variability than the EF and a higher SNR (signal-to-noise ratio). It can be observed that the NEF is directly related to the Bowen ratio (β=H/λE). Introducing Equation (1) in the Bowen ratio gives Equation (2): The NEF is also similar to the Water Deficit Index (WDI) of [6] in which WDI=1-Ea/Eo, where Ea is actual evapotranspiration and Eo is potential evapotranspiration, using available energy for evapotranspiration (Rn-G) instead of the potential evapotranspiration, and latent heat (λE) which is equivalent to Ea.
The purpose of this study is to evaluate three simple operative models with minimum input data requirements for estimating the non-evaporative fraction (NEF) in southeast Spain. Two of the models are based on the "Simplified Relationship", estimating daily H and Rn [7] with two parameterizations for H, [8] and [9]. The third model is based on the S-SEBI (Simplified Surface Energy Balance Index) [4], which estimates the evaporative fraction (EF) directly from the relationship between albedo and surface temperature, modified for this study to account for atmospheric conditions. The study site, located in Almería (Spain) in the Mediterranean Basin, is characterized by escalating water demands for agriculture and tourism, leading to overexploitation of groundwater resources [10]. Given the strategic importance of recharge processes in the region's environmental and socioeconomic system, it is crucial to have more information about the spatial partition of surface energy.
ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data were used to calculate the non-evaporative fraction (NEF). ASTER is currently the only sensor collecting multispectral thermal infrared data at high spatial resolution, being very appropriate for model testing and direct ground comparison [11]. Furthermore, ASTER is on the Terra platform along with the MODIS sensor (Moderate Resolution Imaging Spectroradiometer), making it possible to scale up analyses from field to regional or global scales [12].
2. Description of the study region and data acquired

Study region
The study region (Figure 1), located in the south-eastern Iberian Peninsula (Almería, Spain), comprises 3600 km 2 (36.95ºN, 2.58ºW). It is characterized by its heterogeneity, with altitudinal gradients ranging from sea level up to 2800 m (a.s.l.) in the Sierra Nevada Mountains. Precipitation and temperature regimes vary widely due to the orography [13]. Annual precipitation is the lowest in Tabernas lowlands, where it is less than 200 mm, while in the mountains it ranges from 400 mm to 700 mm, enough to sustain forest growth.
In the center of the study site, the karstic landscape of the Sierra de Gádor mountain range, covering 552 km 2 , consists of a series of thick carbonate rocks (limestones and dolomites), highly permeable and fractured with intercalated calcschists of low permeability underlain by impermeable metapelites [10]. The southern edge of this mountain range is the main source of recharge for the deep Triassic aquifers of the Campo de Dalías [10]. In general, the soils are very thin, rocky and very vulnerable to flash flooding and erosion, the most representative type being Lithic-Typic Haploxeroll according to the FAO classification [14]. The Sierra de Gádor mountains underwent intense, widespread deforestation during the 18 th and 19 th centuries, when the original oaks (Quercus ilex L. and Quercus faginea Lam.), olive trees (Olea europaea L.), poplars (Populus L. spp.) and strawberry-trees (Arbutus unedo Lam.) were cut down for ship construction, fuel and wood for mining activities. Other current pressures or disturbances in the Sierra de Gádor include construction, fire, agriculture and sheep grazing. Nowadays the area is dominated by a sparse shrubland, mostly Genista cinerea Vill, with rock outcrops, bare soil or some grasslands mainly composed of Festuca scariosa Lag. Around 73% of Sierra de Gádor presents this pattern with vegetation cover lower than the 50%.
Shrubland with a sparse cover of pine woodland (Pinus L. spp.) is the second largest natural land cover type (12% of the area). Only 1.5% of the land is covered by dense pine reforestation. 9% of the Sierra de Gádor is dedicated to agriculture (almond and olive trees) [15].
The rest of the study site includes part of the Sierra Nevada Natural Park, composed of pine forest with oak relicts and shrublands. To the northeast, there is the Tabernas Lowlands, an area of badlands with complex topography. Along the ephemeral Andarax River, which flows past the capital city of Almeria, there is a mosaic of citrus orchards and vineyards. One of the most salient features of the scene is the large area of plastic greenhouses spreading over more than 330 km 2 . This unique combination of land covers and uses makes this area an ideal pilot site for model testing. Two field research stations acquire data continuously in the study region, Llano de los Juanes in the Sierra de Gádor and Rambla Honda in the Tabernas Lowlands, shown by white arrows in Figure 1.   Figure 1). Latent and sensible heat flux are measured by an eddy covariance system using a three dimensional sonic anemometer CSTAT3 and a krypton hygrometer KH 2 0 (both from Campbell Scientific Inc., USA). Llano de los Juanes is a representative ~2 km 2 flat area located at an altitude of 1600 m in the high, well-developed karstic plain of the Sierra de Gádor. Fetch is sufficient for the vegetation height and sensors. Vegetation cover is 50-60% and consists mainly of patchy perennial dwarf shrubs (30-35%) dominated by Genista pumilla, Thymus serpylloides Bory and Hormathopylla spinosa L., and grasses (20-25%)

Rambla Honda research site
The Rambla Honda research site is located in a dry valley near Tabernas, Almería, Spain (37º8'N, 2º22'W, 630 m altitude). For a detailed description of the site, see [17]. The valley has been abandoned for several decades and the only activity is now restricted to small-scale sheepherding.
Three perennial species dominate the landscape, Retama sphaerocarpa (L.) Boiss shrubs across the valley floor, Stipa tenacissima L. tussocks on the steep sides of the valley and Anthyllis cytisoides L. shrubs on alluvial fans between the two. The valley floor has deep loamy soils overlying mica schist bedrock. The average annual rainfall is 220 mm with a dry season from June to September.
Experiments related to hydrology and erosion [18,19], surface energy balance and evapotranspiration [5,[20][21][22][23][24][25] and vegetation ecology [26][27][28] among others have been performed at the site during the last decade. The eddy covariance system, currently located at Llano de los Juanes, was acquiring data in Rambla Honda between February 2002 until July 2003. Evapotranspiration was also measured by the BREB (Bowen Ratio Energy Balance) and transpiration with sap flow measurements. This allowed a detailed SVAT (Soil-Vegetation-Atmosphere Transfer) model to be calibrated for Retama sphaerocarpa with a 10% error in daily evapotranspiration compared to the BREB. This SVAT model [20] is an extension of [29] model which combined the [30] and [31] approaches to account for energy partition in sparse vegetation by explicitly separating contributions from soil and vegetation. The parameterization was further extended to the other two dominant species in the area, Anthyllis cytisoides, and Stipa tenacissima with errors (as % of the mean ±std) of 5.6±1.8 and 7.3±0.07, respectively, compared to the BREB. Therefore, in the three species, errors were within the uncertainty of the validation system. Because model calibration was performed for a time span representative of the variability encountered in surface and climate variables at longer time scales, the model can now be run, as the input variables required for parameterization are still being acquired at the site.

Remote sensing and spatial data
ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data on July 18, 2004, July 9, 2004 and June 19, 2005 at 11.00 UTC were acquired for this study. ASTER is an imaging instrument on-board of Terra, a satellite launched in December 1999 as part of the NASA's Earth Observing System (EOS). ASTER is a cooperative effort between NASA, Japan's Ministry of Economy, Trade and Industry (METI) and its Earth Remote Sensing Data Analysis Center (ERSDAC) [32]. ASTER scans a swath of 60 km on the ground every 16 days with a swath angle of ±2.4. The sensor has nine reflective bands and five bands in the thermal infrared (TIR) region.
The ASTER products used in our research included surface reflectance (2AST07) with a spatial resolution of 15 m (VNIR) and 30 m (SWIR), and kinetic temperature at 90 m (2AST0) with a surface temperature absolute precision of 1-4 K. No incidences have been reported for these scenes. The three images did not cover exactly the same area due to ASTER's off-nadir sensor pointing capability. A digital elevation model (DEM) from the USGS (United States Geological Survey) with a 30 m resolution and a digital 0.5-m pixel orthophoto (from the Andalusian Regional Government) were used at different stages of the study.
Half-hourly air temperatures (ºC) at the time of the satellite overpass (11.00 UTC) were acquired from meteorological stations for validation purposes. Ten or eleven stations were available for each image depending on scene coverage. Seven of the stations belong to the EEZA (Estación Experimental de Zonas Áridas), and the rest to the Andalusian Regional Government, (Red de Información Agroclimática de Andalucía).

Estimating the non-evaporative fraction from the "Simplified Relationship"
Daily NEF (non-evaporative fraction) was estimated from ASTER and ancillary data using the ratio between daily sensible heat (H) and net radiation (Rn): H/Rnd. Soil heat flux (G) can be considered negligible at a daily scale compared to the other components of the surface energy balance [2,8], as shown in Equation (3): where EF is the evaporative fraction, λE is latent heat flux (Wm -2 ), H is the sensible heat flux (Wm -2 ) and Rnd is the daily net radiation (Wm -2 ).

Daily net radiation (Rnd)
Daily net radiation (Rnd) was calculated as the balance between incoming (↓) and outgoing fluxes (↑) of shortwave (Rs) and longwave (Rl) radiation. By agreement, incoming fluxes are positive and outgoing negative. Net radiation is the sum of net shortwave (Rns) and net longwave radiation (Rnl) [2].
First, Rni, instantaneous net radiation at the time of image acquisition, was calculated by estimating its four components: The shortwave net radiation using remote sensing data was calculated as in Equation (5): Rs↓ is the incoming solar radiation or incoming shortwave radiation, which was estimated at the time of the satellite overpass (11.00 UTC) using a solar radiation model [33] accounting for elevation, aspect, latitude and longitude, solar geometry, atmospheric transmissivity, and the influence of the surrounding topography. α is the broadband surface albedo estimated according to [34] as follows: Where ρ i is reflectance at the surface for the band indicated by the subscript i. Reflectance was acquired from ASTER product 2AST07.
The longwave energy components are related to surface and atmospheric temperatures by the Stefan-Boltzmann Law. The outgoing longwave radiation was calculated at the time of image acquisition as in Equation (7): Where σ is the Stefan-Boltzmann constant, (5.67x10 -8 W m -2 ), T s is surface temperature (K), and ε s is broadband emissivity for the surface, estimated based on the logarithmic relationship to NDVI as proposed by [35]: Radiometric surface temperature, T s , was acquired directly from the ASTER kinetic temperature product retrieved by the TES (Temperature Emissivity Separation) algorithm.
An empirical function was used for the incoming longwave radiation Rl↓ [36].
where T air is the air temperature and c and d are constants (0.261 and 7.77 10 -4 K 2 , respectively). Daily net radiation (Rnd) (Wm -2 ) was calculated from Rni by assuming Rnd/Rni≈0.3 ±0.03 at midday as proposed by [8].

Sensible heat flux (H)
The sensible heat flux (H) can be estimated by the turbulent transport from the surface to the lower atmosphere based on surface layer similarity of mean temperature and wind speed profiles using the resistance formula: where Ts is the land surface temperature and Tair is the air temperature, both at the time of image acquisition; r h is a turbulent exchange coefficient dependent on wind speed, aerodynamic roughness length, roughness length for heat transfer and Monin-Obukov length [37] , and ρ and Cp are air density and specific heat at constant pressure, respectively. Due to differences between radiometric surface temperature and aerodynamic surface temperature, another resistance term Kb -1 (excess resistance) was added to the denominator of Equation (10) to account for the difference between roughness length for heat (z o ) and momentum (z oh ): Kb -1 =ln(z o /z oh ). However, in practice, it is hard to get large-scale spatial estimates for all the variables in the resistance terms of H, therefore more operational parameterizations have been proposed. One of the most widely used approaches to solving the surface energy balance that explicitly calculates H and Rnd is the so-called "Simplified Relationship" [7,8,38] which states that daily evapotranspiration (λE d ) can be estimated from the difference between daily net radiation (Rnd) and daily sensible heat flux (H), by estimating H from the difference between instantaneous surface (Ts) and air temperatures (T air ) near midday, as in Equation (11): The "Simplified Relationship" has been verified empirically and theoretically [8,[39][40][41]. B can be defined as a mean exchange coefficient of sensible heat transfer. According to this relationship, the surface-atmosphere temperature gradient at midday, related to the instantaneous sensible heat flux at midday by means of B, can be considered representative of the influence of daily H in the energy balance by assuming that the evaporative fraction is constant throughout the day [3,8,42]. Considering this and equation (10), B could be estimated as in (12): [8] proposed two values for B as a first approximation, 0.25 mmK -1 day -1 for stable atmospheric conditions (T s -T a <0) and 0.18 mmK -1 day -1 for unstable conditions (T s -T air > 0). At our study site at the time of image acquisition, unstable conditions tend to be prevalent [20].
Another operative approach for estimating H that also builds upon the simplified relationship was proposed by [9]. They showed that the main factor affecting resistance to heat transfer, and therefore B, was vegetation cover, and established a linear relationship between fractional cover and B, and between fractional cover and an exponent n close to 1 affecting the (T s -T air ) term. To estimate fractional cover NDVI was rescaled between values associated with sites where fractional cover is 0 (bare soil sites) and sites with vegetation cover =1, associated with dense forest sites. In our study site NDVI from associated to those two extremes was: 0.16 ± 0.012 (mean ± standard deviation) for bare soil sites and 0.68 ± 0.20 from complete vegetation cover. Mean values from bare soil and complete vegetation were taken to calculate B and n.
Hereinafter we will refer to NEF and H as calculated using the [8] model for stable conditions as NEF Seguin and H Seguin and the one by [9] as NEF Carlson. and H Carlson. Air temperature (T air ) Air temperature (T air ) is used to estimate sensible heat flux and net radiation. To avoid relying on meteorological information, T air was estimated from the images using the NDVI-Ts triangle as proposed by [9] in an approach similar to [43] and [43,44]. The apex of the NDVI-Ts space (high NDVI and low temperature) should correspond to pixels with high NDVI located at the wet edge of the triangle, and can be assumed to be at T air . T s at the apex is found by locating minimum surface temperature areas in the scene. Those with the highest NDVI corresponding to forest patches are identified, and the average T s for that selected region is calculated. Due to the altitudinal gradients in the study region, T air must be corrected using the pixels at the region forming the apex as a reference altitude. Then positive corrections for altitude can be made for pixels below the baseline and viceversa for pixels above it, considering a lapse rate of 6.5ºC per 1000 m. This works better than considering a single T air for the whole area, by assuming constant meteorological conditions at the blending height [9] or the dry and wet pixel method used by [3], evaluated in preliminary tests (results not shown), and which is probably better suited for flat areas.

Estimating the non-evaporative fraction from S-SEBI (Simplified-Surface Energy Balance Index)
Another method of estimating the NEF, other than explicitly estimating surface energy balance variables, was derived from the S-SEBI (Simplified Surface Energy Balance Index) [4]. S-SEBI directly estimates the evaporative fraction, EF, on a pixel basis by means of a contextual relationship between surface temperature (T s ) and albedo. It has been applied to estimate the evaporative fraction for crops and natural vegetation at different spatial scales [45][46][47] and assumes that atmospheric conditions remain relatively constant across the study region, and requires enough wet and dry pixels in the scene for hydrological contrast. Figure 2. Scheme of the S-SEBI model (adapted from [4]) showing the relationship between surface albedo (α) and surface temperature (Ts). "A" and B" correspond to the evaporation-controlled and radiation-controlled domains, respectively. Ts obs is the observed temperature, Ts λE is the temperature in the evaporation-controlled domain and Ts H is the temperature in the radiationcontrolled domain.
Boundary conditions for the model can be extracted from the scatter plot of Ts and albedo. At low reflectance, Ts is almost constant with increasing albedo (Line "A" in Figure 2). This would be the case of saturated water surfaces or irrigated areas where all available energy is used for evapotranspiration (EF≈1; with maximum evapotranspiration, λE max ). When albedo increases, Ts also increases because of reduced evapotranspiration at the expense of increased sensible heat flux. This is also within the "evaporation-controlled domain". When albedo increases beyond a certain level, there is an inflection point, and Ts begins to decrease with albedo (line "B" in Figure 2). At this point soil moisture is so low that no evaporation can take place and all the available energy is dissipated as sensible heat flux (EF≈0 with maximum sensible heat flux, H max ). However, the increase in albedo produces a decrease in shortwave net radiation reducing Ts. This is the "radiation controlled domain". Rescaling each observed Ts according to these boundary conditions, "A" and "B" in Figure 2 allows to obtain EF as expressed in Equation (14) (14) where Ts obs is the observed temperature, Ts λE is the temperature at the lower boundary function or the evaporation-controlled domain and Ts H is the temperature at the upper boundary function or radiationcontrolled domain.
We modified the S-SEBI to account for variation in atmospheric conditions across the study site. Therefore, in this study, surface temperature (Ts) was standardized by T air using the following equation to calculate NEF: Where DT obs is the observed difference between surface and air temperature, DT λE is the difference between surface and air temperature at the lower boundary function or evaporation-controlled domain, and DT H is the difference between surface and air temperature at the upper boundary function or radiation-controlled domain.
Both the upper and lower boundary functions were calculated by quantile regression [48] from the scatter plot of albedo vs. Ts for the study region, using the 5% and the 95% quantiles, respectively. H was estimated by multiplying NEF by Rnd, assuming daily G to be zero, for comparison with H Seguin and H S -SEBIt .
The NEF and H calculated by this approach, in which the original S-SEBI (Simplified-Surface Energy Balance Index) formulation for calculating the EF (evaporative fraction) was modified by including T air , are hereinafter referred to as NEF S-SEBIt and H S-SEBIt .

Validation of the non-evaporative fraction model results
It is extremely complicated to validate surface energy fluxes estimated from remote sensing data due to limitated availability of measured surface fluxes for several surface types over large scales. In addition, field measurements and remote sensing footprints are not always comparable. In this paper we propose two validation procedures (a) Comparison of representative semi-arid vegetation surface types for which these field surface fluxes were available: Retama, Anthyllis and Stipa, and a mixture of shrubs and forbs (b) Evaluation of NEF estimated at sites of known vegetation and land cover. The eddy covariance system located at Llano de los Juanes research site has been acquiring data over the shrublands, being the fluxes representative of an area of approximately 2 km 2 since Fall 2003. At the Rambla Honda research site, outputs from the detailed SVAT model parameterized by [20] and [49] for Retama, Anthyllis and Stipa, with errors of < 10%, were used as field references. For comparison between the image and modeled data, patches of those species in Rambla Honda were selected based on field visits and the aerial photo (0.5 m). The Rambla Honda field site was only present in the July 18 th 2004 scene.
In addition, we used the wetland named "Cañada de las Norias" located in the greenhouse area as a validation site. For validation purposes we considered a field value for daily H=0, and therefore also for NEF=0, similarly to [3] and [4]. The wetland comprises 135 ha with a maximum depth of 2 m. The riparian vegetation is composed of Phragmites australis, Tamarix canariensis, and Tamarix africana, the latter also appears within the water table, and in the shallower parts Typha domingensis and Scirpus litoralis. Within the wetland macroalgi from Entermorpha and Cladophora genus, indicative of high eutrofization, tend to replace aquatic macrophytes [50]. Solids and algi increase water turbidity and reduce the effective penetration of solar radiation in the water column, which reduces the water storage term at a daily scale (G) [51] that becomes almost negligible in the case of vegetated wetlands [52].
In shallow lakes such as "Cañada de las Norias", heat accumulates during the day and supplies evaporative heat loss at night. Thus, although hourly changes in the storage term (G) can be high, on a daily basis (24 hr) G is smaller. For instance, in a similar wetland of the same area hourly G fits in summer a sinusoidal curve with midday peaks of around 50 Wm -2 [51,53]. In the Daimiel wetlands in Central Spain, daily energy storage (24 hr) on open water presents oscillations of ±39.9 Wm -2 [54]. In a semiarid shallow wetland in Nebraska, G for open water ranged between -76 and 60 Wm -2 with those extremes associated with rainy or cloudy days [55].
Regarding the sensible heat flux (H) on a daily basis its contribution to the wetland energy balance is minimum. In a semi-arid wetland daily H presented values ranging between ±20 Wm -2 . In a wetland nearby "Cañada de las Norias", H contribution to heat budget was found to be less than 2% [53].
Because NEF Seguin and NEF Carlson models assume G=0, which is correct over land surfaces, modeled results over the lake were corrected just for validation considering NEF=H/(Rn-G) instead of NEF=H/Rn by assuming a daily G value in the wetland at the most between ±50 Wm -2 (around ±23% of Rn)-.
Estimated means of H, Rnd and NEF from each patch in the image and observed daily means from the eddy covariance or the SVAT model were compared in terms of R 2 , RMSE (Root Mean Square Error), MAE (Mean Absolute Error), p, slope and intercept of the linear regression.

Comparison of models estimating the non-evaporative fraction (NEF)
Spatial patterns with the NEF S-SEBIt and NEF Seguin were observed to be more similar to each other on the three dates than with the NEF Carlson model, although NEF S-SEBIt always yielded higher values. Across the study site, the lowest NEF with NEF S-SEBIt and NEF Seguin corresponded to water surfaces (sea and lakes), and high-altitude mountain forests. The highest NEF were located in the Tabernas lowlands which is plausible at this time of the year. However, NEF Carlson values in the Tabernas lowlands where vegetation cover is scarce, are unrealistically low. Figure 3 shows an example of results for NEF (non-evaporative fraction) on July18, 2004. Linear regression of NEF on the three dates using all the pixels in each scene, including sea water, and the greenhouse area, shows (Table 1) very good correspondence between NEF Seguin and NEF S-SEBIt. but not between NEF Carlson and the other two models.  When model performance is examined more in detail for certain surface types on each date as a first approach to validating our results (Figure 4), it is observed that the most coherent patterns in NEF are provided by NEF Seguin and NEF S-SEBIt with a high NEF for dry, bare soil sites and low NEF in mountain forest sites. The NEF Carlson seems reasonable for forest types, with values similar to NEF Seguin. However, over bare soil surfaces, NEF Carlson greatly underestimates NEF yielding even lower levels than forests (e.g., the limestone quarry would be evapotranspiring at the same rate as water according to NEF Carlson ). The [9] approach has been successfully used in Argentina at sites with NDVI over 0.5 [56], however, in semiarid sites with low NDVI, estimating the B exchange coefficient for sensible heat transfer based solely on NDVI produces unrealistic patterns. It is therefore preferable to consider a constant value of B as at least a first approximation. In the Sierra de Gádor, B does not vary much: 0. It is important to mention that NEF S-SEBIt and NEF Seguin are very similar, but an offset equivalent to the NEF for water (lake at the middle of the greenhouse area) is observed for NEF S-SEBIt suggesting a low wet edge, which would be easy to recalibrate. The fact that R 2 >0.96 and the slope between NEF S-SEBIt and NEF Seguin is close to 1 on the three dates is very important, as these two models are calculated with very different approaches and NEF S-SEBIt requires fewer input variables than NEF Seguin . Figure 4. NEF (non-evaporative fraction) over selected surface types calculated by the NEF Seguin , NEF Carlson and NEF S-SEBIt models. The first set of surfaces corresponds to undisturbed sites.

NEF S-SEBI t NDVI
Sierra N and Sierra S are pine forests on the northern and southern slopes of the Sierra Nevada Mountains, three densities of oak relicts correspond to: oaks (dense), oaks (sparse) and oaks. Pines correspond old reforested sites. The second set is for disturbed sites: greenhouses (Greenh), a strong burn scar (burnt), a limestone quarry (quarry), almond orchards (almond), and an abandoned mine (mine). The third set comprises miscellaneous sites: the Tabernas badlands, a lake, a golf coarse (golf), irrigated citrus orchards along the Andarax River (orchards); Ll. Juanes is the Llano de los Juanes, eph.river is the ephemeral Andarax River and R.Honda is the Rambla Honda research site.

Field validation of the non-evaporative fraction models
Air temperature Air temperature was required to estimate longwave incoming radiation, sensible heat flux and the S-SEBI t evaporative fraction ( Table 2).
The overall adjustment is good (< 2ºC), but T air estimates are subject to local errors. Altitude is not the only factor affecting T air , but using this approach has the advantage of not having to use meteorological station data. Also, any systematic error in T s retrieval will propagate in T air . These errors should therefore be compensated when calculating T s -T air differences for estimating sensible heat flux. In our case, this approach yields better results than the [3] wet and dry pixel method. Table 2. Air temperature estimates at the study site. Mean Absolute Error (MAE) is the average absolute difference in residuals between estimated and measured air temperature at the meteorological stations. MAE after adjustment is the temperature corrected for a lapse rate of 6.5ºC per 1000 m. Reference altitude is the altitude of the pixels the apex temperature was taken from. Daily net radiation (Rnd) Results from estimating the Rnd using ASTER data show very good agreement with field data (Table 3), with an overall error < 5% of Rnd and close to the 1:1 line. Reported Rnd accuracy of the net radiometer is around ±10% (NR-lite by Kipp&Zonen), which is the case for some remote sensing studies [2]. Other remote sensing studies have reported an Rnd RMSE of 20-30 Wm -2 [57][58][59] or 40-50 Wm -2 [47,60,61] depending on data, model and surface type used for validation.
Emissivity values estimated in this work are reasonable according to reported values for soil and vegetation [35,62] The emissivity values used in the van de Griend and Owe [35,62] model were 0.91 for bare soil, between 0.94-0.97 for shrubs and grasslands with partial cover depending on vegetation type, and 0.98-0.99 for complete vegetation cover. At the Rambla Honda field site measured emissivity is highly variable, between 0.94-0.97 for soil and 0.94-0.97 for vegetation ( [22]). Given the pixel size (90 m) the mix of soil and vegetation makes it very difficult to get pure pixels. In this area, estimated emissivity using Aster data ranged between a minimum of 0.91 up to a maximum of 0.95 depending on vegetation type and cover.
In the whole study region, maximum emissivity correspond to sites with complete cover such as forests, golf course, and irrigated orchards (0.97-0.98 depending on date). Estimated emissivity decreases with vegetation cover reaching values of 0.92-0.94 at the shrublands in Sierra de Gádor. The lowest emissivity values correspond to almond orchards dominated by bare soil signature (0.94-0.92), limestone quarry and mining areas (0.90-0.92). Table 3. Daily Net Radiation Rnd (Wm -2 ) of Retama, Anthyllis, Stipa, shrubs, and bare soil at the Llano de los Juanes and Rambla Honda research sites. "Field" is the Rnd field measurement, and ASTER is the Rnd estimated using ASTER and ancillary data. AE is the absolute value of (Rnd field -Rnd Aster ). The % Error is calculated as (Rnd field -Rnd Aster )100/ Rnd field. For an overall error evaluation, the MAE (mean absolute error), the average AE, the mean % error, (% Error), R 2 (Pearson correlation coefficient), p (probability), slope and regression intercept between field and ASTER results were calculated. The non-evaporative fraction (NEF) and the sensible heat flux H were validated using Eddy Covariance data at the Llano de los Juanes research site and a SVAT model previously calibrated with the Eddy Covariance and Bowen Ratio systems at the Rambla Honda research site (Figure 1). A lake located at the greenhouse area was considered to evaporate at potential rate (H~0).

Date
Validation results are similar for NEF and H (Wm -2 ) ( Tables 4 and 5). H Seguin and therefore, NEF Seguin , provide the best overall performance being closer to the 1:1 line. Although the R 2 is higher for H Carlson and therefore, so is NEF Carlson , their RMSE is still the highest of the three models.
At the Llano de los Juanes research site, ASTER results underestimate H compared to Eddy Covariance measurements using H Seguin and H Carlson around 30%. In addition to the simplicity of the modeling approaches, there is an error propagation from input data. Thus, although reported errors in T s are within acceptable quality levels (< 4 K) they contribute to final error combined with the error in T air estimates (< 2 K) and in the aerodynamic resistance that, in this case, is probably too high.
We should also be aware that the Eddy Covariance and Bowen Ratio techniques are subject to error. Uncertainty levels in the Eddy Covariance are around 20% [63] and 10% in the Bowen Ratio technique [64,65]. Moreover, in semi-arid areas with sparse vegetation cover, error in energy fluxes tends to be even higher, around 25% [25].
Daily H estimates from remote sensing models usually contribute the highest uncertainty to the surface energy balance, with errors at a daily scale of around 20-30% or 1 mm day -1 , equivalent to ~29 Wm -2 [2]. In our case, the RMSEs for H Seguin and H S-SEBIt , are below that threshold, but not for H Carlson, Individual errors range between 3% and 30% for the worst cases (Llano de los Juanes and H Carlson ) being the lowest in general for H S-SEBIt .
In general, reported range of errors in H varies widely depending on surface type, image data, time average period, and model used. [66] consider an error of around 50 Wm -2 as acceptable for instantaneous H and 23 Wm -2 for daily H. In the literature, best case errors for instantaneous (halfhourly) and daily fluxes are around 10-22 Wm -2 [2,57] and can reach up to 50%, even using sophisticated models when parameterization is not good. It is more complicated to get good estimates over heterogeneous semiarid areas than over agricultural or humid sites [67]. [68] obtained an RMSE of 43.35 Wm -2 for instantaneous H, while [69] found errors of around 40% in the Great Basin Desert, and [70] obtained an RMSE=47 Wm -2 for instantaneous H using dual angle observations over a semiarid grassland in Mexico.
Our results for the NEF (non-evaporative fraction) are within errors reported for evaporative fraction at a daily scale from the SEBAL model, with a more complex parameterization, with RMSEs in the daily evaporative fraction (EF=1-NEF) between 0.10-0.20 [71]. [72] obtained an RMSE for daily EF of 0.13, and [45] between 0.09-0.05 using S-SEBIt over European forests compared to Euroflux data.
Results from field validation confirm previous results about model performance for selected surface types as shown in Figures 3 and 4. Given the simplicity of the models used in this study, our results are reasonable, being the errors within ranges reported in the literature.  Table 5. Field validation of the NEF (non-evaporative fraction), NEF=H/(H+λE) NEF Seguin , NEF Carlson and NEF S-SEBIt models. AE is the absolute error (absolute difference between model and field observations). When using the lake for validation, two cases have been considered G lake =50 Wm -2 and G lake =-50 Wm -2 . For an overall error evaluation, the MAE (mean absolute error), which is the average AE, the R 2 (Pearson correlation coefficient), p (probability), slope and regression intercept between field and ASTER data were calculated (n=9 observations). For NEF Seguin and NEF Carlson overall error estimation has been performed twice, one with the dataset including the lake when G lake =-50 Wm -2 and another with the dataset including the lake when G lake = -50 Wm -2 , the latter in parenthesis the

Conclusions
Three operative models for estimating the non-evaporative fraction (NEF) as an indicator of the surface water deficit, were evaluated in a semiarid area of southeast Spain. Of the three models evaluated, the NEF computed with the Simplified Relationship for unstable conditions [8] was found to be equivalent to NEF estimated from S-SEBIt (Simplified-Surface Energy Balance Index ) [4] modified in our study to account for varying atmospheric conditions (R 2 >0.97, p<0.0001). This is important, as the two algorithms are very different. One explicitly solves the variables in the surface energy balance, while the other uses a contextual relationship between albedo and surface temperature requiring fewer input variables. On the other hand, the spatial patterns obtained with Carlson's model [9] were different from the other two models.
Comparison to field data showed that net radiation (Rn) modeled with ASTER data produces very good results (R 2 =0.91, p<0.0001), with an overall error within reported instrumental accuracy (< 5%). Validation of sensible heat flux (H) and NEF showed the most coherent spatial patterns with Seguin's approach [8] followed by the S-SEBI t model. Of the three models, Carlson´s approach had the highest RMSE. In semiarid sites with low NDVI, estimating the B exchange coefficient for sensible heat transfer based solely on NDVI produces unrealistic patterns, making it preferable to use a constant value of B, at least as a first approximation.
The errors found for NEF Seguin and NEF S-SEBIt are within those reported in the literature (RMSE for NEF are ∼ 0.18 and 0.13 and for H are 29.12 Wm -2 and 25.97 Wm -2 for Seguin's and S-SEBI t respectively). Given the simplicity of the models used, these results are very promising for operative monitoring of the surface water deficit and the partition of surface energy into sensible and latent heat flux.