Skip to main content
Estimating evaporation is important when manag- ing water resources and cultivating crops. Evaporation can be estimated using land surface heat flux models and remotely sensed land surface temperatures (LST), which have recently become... more
Estimating evaporation is important when manag-
ing water resources and cultivating crops. Evaporation can be
estimated using land surface heat flux models and remotely
sensed land surface temperatures (LST), which have recently
become obtainable in very high resolution using lightweight
thermal cameras and Unmanned Aerial Vehicles (UAVs). In
this study a thermal camera was mounted on a UAV and ap-
plied into the field of heat fluxes and hydrology by concate-
nating thermal images into mosaics of LST and using these
as input for the two-source energy balance (TSEB) modelling
scheme. Thermal images are obtained with a fixed-wing
UAV overflying a barley field in western Denmark during
the growing season of 2014 and a spatial resolution of 0.20 m
is obtained in final LST mosaics. Two models are used: the
original TSEB model (TSEB-PT) and a dual-temperature-
difference (DTD) model. In contrast to the TSEB-PT model,
the DTD model accounts for the bias that is likely present
in remotely sensed LST. TSEB-PT and DTD have already
been well tested, however only during sunny weather condi-
tions and with satellite images serving as thermal input. The
aim of this study is to assess whether a lightweight thermal
camera mounted on a UAV is able to provide data of suf-
ficient quality to constitute as model input and thus attain
accurate and high spatial and temporal resolution surface en-
ergy heat fluxes, with special focus on latent heat flux (evapo-
ration). Furthermore, this study evaluates the performance of
the TSEB scheme during cloudy and overcast weather con-
ditions, which is feasible due to the low data retrieval alti-
tude (due to low UAV flying altitude) compared to satellite
thermal data that are only available during clear-sky condi-
tions. TSEB-PT and DTD fluxes are compared and validated
against eddy covariance measurements and the comparison
shows that both TSEB-PT and DTD simulations are in good
agreement with eddy covariance measurements, with DTD
obtaining the best results. The DTD model provides results
comparable to studies estimating evaporation with similar
experimental setups, but with LST retrieved from satellites
instead of a UAV. Further, systematic irrigation patterns on
the barley field provide confidence in the veracity of the spa-
tially distributed evaporation revealed by model output maps.
Lastly, this study outlines and discusses the thermal UAV im-
age processing that results in mosaics suited for model input.
This study shows that the UAV platform and the lightweight
thermal camera provide high spatial and temporal resolution
data valid for model input and for other potential applications
requiring high-resolution and consistent LST.
Research Interests:
El objetivo de esta investigación es predecir la probabilidad de ocurrencia de incendios por causa natural (rayo) para la España peninsular mediante el ajuste de un modelo logístico. El modelo utiliza como variable dependiente los partes... more
El objetivo de esta investigación es predecir la probabilidad de ocurrencia de incendios por causa natural (rayo) para la España peninsular mediante el ajuste de un modelo logístico. El modelo utiliza como variable dependiente los partes de incendios por causa de rayo y como variables independientes los índices meteorológicos (FFMC, DMC, 1h y 10h), las descargas de rayos, la topografía (pendiente, altitud y exposiciones), las coberturas forestales y los tipos de clima. Los resultados muestran que el DMC medio, los rayos totales, el clima oceánico y mediterráneo de montaña son las variables más explicativas en el proceso. El modelo ofrece una precisión global aceptable (64,12%). Esto demuestra que las condiciones climáticas, el grado de humedad del combustible y la actividad de tormentas están relacionadas con la probabilidad de ocurrencia de incendios.
... MIKE BASIN model was used to calculate water balance in Hong-Thai Binh river basin, whereMIKE NAM model was also used for inflow calculation in this basin. ... 3. Application of MIKE BASIN model to calculate the water balance in... more
... MIKE BASIN model was used to calculate water balance in Hong-Thai Binh river basin, whereMIKE NAM model was also used for inflow calculation in this basin. ... 3. Application of MIKE BASIN model to calculate the water balance in Hong-Thai Binh river basin 3.1. ...
Evapotranspiration is one of the key components of the water balance. Knowing its value is crucial for managing water resources, especially in areas with water scarcity. Despite that its determination at regional scales is a difficult... more
Evapotranspiration is one of the key components of the water balance. Knowing its value is crucial for managing water resources, especially in areas with water scarcity. Despite that its determination at regional scales is a difficult task, since there is not a direct way for its measurement; remote sensing has proved to be a powerful tool to solve that issue. In this work the triangle method, based on the relation of the surface temperature with a vegetation index, was applied on Landsat5‐TM data over a heterogeneous patch of 5x5 km in the Henares river basin – part of the Tagus river basin in central Spain. On‐site measurements were carried out during the growing season of an experimental wheat field in 2010, including sensible heat flux (H) estimated with a scintillometer, soil heat flux (G), net radiation (Rn) and other typical meteorological data, which were used to obtain actual evapotranspiration (ET) rates at field scale, as a residual of the surface energy balance. A mask r...
ABSTRACT This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution and sharing with... more
ABSTRACT This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution and sharing with colleagues. Other uses, including reproduction and distribution, or selling or licensing copies, or posting to personal, institutional or third party websites are prohibited. In most cases authors are permitted to post their version of the article (e.g. in Word or Tex form) to their personal website or institutional repository. Authors requiring further information regarding Elsevier's archiving and manuscript policies are encouraged to visit:
ABSTRACT The REFLEX 2012 campaign was initiated as part of a training course on the organization of an airborne campaign to support advancement of the understanding of land-atmosphere interaction processes. This article describes the... more
ABSTRACT The REFLEX 2012 campaign was initiated as part of a training course on the organization of an airborne campaign to support advancement of the understanding of land-atmosphere interaction processes. This article describes the campaign, its objectives and observations, remote as well as in situ. The observations took place at the experimental Las Tiesas farm in an agricultural area in the south of Spain. During the period of ten days, measurements were made to capture the main processes controlling the local and regional land-atmosphere exchanges. Apart from multi-temporal, multi-directional and multi-spatial space-borne and airborne observations, measurements of the local meteorology, energy fluxes, soil temperature profiles, soil moisture profiles, surface temperature, canopy structure as well as leaf-level measurements were carried out. Additional thermo-dynamical monitoring took place at selected sites. After presenting the different types of measurements, some examples are given to illustrate the potential of the observations made.
ABSTRACT This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution and sharing with... more
ABSTRACT This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution and sharing with colleagues. Other uses, including reproduction and distribution, or selling or licensing copies, or posting to personal, institutional or third party websites are prohibited. In most cases authors are permitted to post their version of the article (e.g. in Word or Tex form) to their personal website or institutional repository. Authors requiring further information regarding Elsevier's archiving and manuscript policies are encouraged to visit:
Fire risk indices are useful tools for fire prevention actions by fire managers. A fire ignition is either the result of lightning or human activities. In European Mediterranean countries most forest fires are due to human activities.... more
Fire risk indices are useful tools for fire prevention actions by fire managers. A fire ignition is either the result of lightning or human activities. In European Mediterranean countries most forest fires are due to human activities. However, lightning is still an important fire ignition source in some regions. Integration of lightning and human fire occurrence probability into fire risk indices would be necessary to have a complete picture of the causal agents and their relative importance in fire occurrence. We present two methods for the integration of lightning and human fire occurrence probability models at 1 × 1 km grid cell resolution in two regions of Spain: Madrid, which presents a high fire incidence due to human activities; and Aragón, one of the most affected regions in Spain by lightning-fires. For validation, independent fire ignition points were used to compute the Receiver Operating Characteristic (ROC)-Area Under de Curve (AUC) and the Mahalanobis Distance. Results in Madrid are satisfactory for the human fire occurrence probability model (AUC∼0.7) but less suitable for the lightning and the integrated models. In Aragón the fit for the human model is reasonable (AUC∼0.7) whereas for the integration methods is practically useless (AUC∼0.58).
Land Surface Temperature (LST) is one of the key inputs for Soil-Vegetation-Atmosphere transfer modeling in terrestrial ecosystems. In the frame of BIOSPEC (Linking spectral information at different spatial scales with biophysical... more
Land Surface Temperature (LST) is one of the key inputs for Soil-Vegetation-Atmosphere transfer modeling in terrestrial ecosystems. In the frame of BIOSPEC (Linking spectral information at different spatial scales with biophysical parameters of Mediterranean vegetation in the context of global change) and FLUXPEC (Monitoring changes in water and carbon fluxes from remote and proximal sensing in Mediterranean “dehesa” ecosystem) projects LST retrieved from Landsat data is required to integrate ground-based observations of energy, water, and carbon fluxes with multi-scale remotely-sensed data and assess water and carbon balance in ecologically fragile heterogeneous ecosystem of Mediterranean wooded grassland (dehesa). Thus, three methods based on the Radiative Transfer Equation were used to extract LST from a series of 2009–2011 Landsat-5 TM images to assess the applicability for temperature input generation to a Landsat-MODIS LST integration. When compared to surface temperatures simulated using MODerate resolution atmospheric TRANsmission 5 (MODTRAN 5) with atmospheric profiles inputs (LSTref), values from Single-Channel (SC) algorithm are the closest (root-mean-square deviation (RMSD) = 0.50 °C); procedure based on the online Radiative Transfer Equation Atmospheric Correction Parameters Calculator (RTE-ACPC) shows RMSD = 0.85 °C; Mono-Window algorithm (MW) presents the highest RMSD (2.34 °C) with systematical LST underestimation (bias = 1.81 °C). Differences between Landsat-retrieved LST and MODIS LST are in the range of 2 to 4 °C and can be explained mainly by differences in observation geometry, emissivity, and time mismatch between Landsat and MODIS overpasses. There is a seasonal bias in Landsat-MODIS LST differences due to greater variations in surface emissivity and thermal contrasts between landcover components.
The temperature-based two-source model (TSM) of Norman et al. (1995) has not been properly evaluated under the water stress conditions that are typical in natural Mediterranean drylands. In such areas, the asynchrony between precipitation... more
The temperature-based two-source model (TSM) of Norman et al. (1995) has not been properly evaluated under the water stress conditions that are typical in natural Mediterranean drylands. In such areas, the asynchrony between precipitation and energy supply strongly reduces evapotranspiration, E (or latent heat flux, LE, if expressed in energy terms), making sensible heat flux (H) the dominant turbulent heat flux. In this study, we present a detailed analysis of the main environmental factors affecting the TSM effectiveness under such challenging conditions. The accuracy of the TSM, evaluated via errors in 15-min H estimates, was shown to have a diurnal variation. Accuracy was clearly reduced for solar elevation angles lower than 25° and during marginal hours of daytime, before 10 am and after 3 pm. The surface to air temperature difference (TR − Ta) and the wind speed were the two environmental factors showing the strongest effect on the TSM accuracy. In contrast with results observed in other ecosystems, in this Mediterranean tussock grassland the TSM accuracy was not clearly reduced by cloudiness and it was improved under highly stressed vegetation conditions. The parallel resistances scheme of the TSM (TSMP) showed overall lower errors and a lower tendency to underestimate at high H values, but the series resistances scheme of the TSM (TSMS) increased the model accuracy under some specific circumstances such as low energy supply and atmospheric neutral conditions.

Finally, two extrapolation methods to obtain daytime (Rn > 55 W m−2) turbulent fluxes from the 15-min estimates of TSM were compared: (i) assuming the self-preservation of the evaporative and the non-evaporative fraction (EF and NEF method) and (ii) averaging the total daytime instantaneous fluxes (Averaging method). Despite the assumption of daytime self-preservation of EF and NEF was showed consistent, this method retrieved less accurate daytime estimates of H, and E than the Averaging method as a result of inaccuracies affecting estimates of EF and NEF from the TSM at our site. Moreover, better daytime estimates of H and E were obtained when using instantaneous fluxes from the TSMP than from the TSMS. Thus, reliable daytime estimates of H were obtained from the TSMP in a Mediterranean dryland, with mean errors of 20% and high correlations (R2 = 0.85). However, daytime E was strongly overestimated (125%) using the TSM by both methods, although a good correlation with eddy covariance measurements was found (R2 = 0.84).
The Dual Temperature Difference (DTD) model, introduced by Norman et al. (2000), uses a two source energy balance modelling scheme driven by remotely sensed observations of diurnal changes in land surface temperature (LST) to estimate... more
The Dual Temperature Difference (DTD) model, introduced by Norman et al. (2000), uses a two source energy balance modelling scheme driven by remotely sensed observations of diurnal changes in land surface temperature (LST) to estimate surface energy fluxes. By using a time-differential temperature measurement as input, the approach reduces model sensitivity to errors in absolute temperature retrieval. The original formulation of the DTD required an early morning LST observation (approximately 1 h after sunrise) when surface fluxes are minimal, limiting application to data provided by geostationary satellites at sub-hourly temporal resolution. The DTD model has been applied primarily during the active growth phase of agricultural crops and rangeland vegetation grasses, and has not been rigorously evaluated during senescence or in forested ecosystems. In this paper we present modifications to the DTD model that enable applications using thermal observations from polar orbiting satellites, such as Terra and Aqua, with day and night overpass times over the area of interest. This allows the application of the DTD model in high latitude regions where large viewing angles preclude the use of geostationary satellites, and also exploits the higher spatial resolution provided by polar orbiting satellites. A method for estimating nocturnal surface fluxes and a scheme for estimating the fraction of green vegetation are developed and evaluated. Modification for green vegetation fraction leads to significantly improved estimation of the heat fluxes from the vegetation canopy during senescence and in forests. When the modified DTD model is run with LST measurements acquired with the Moderate Resolution Imaging Spectroradiometer (MODIS) on board the Terra and Aqua satellites, generally satisfactory agreement with field measurements is obtained for a number of ecosystems in Denmark and the United States. Finally, regional maps of energy fluxes are produced for the Danish Hydrological ObsErvatory (HOBE) in western Denmark, indicating realistic patterns based on land use.
A two-source model based on surface temperature is tested in a Mediterranean dryland.Reliable estimates of sensible heat flux were achieved but not for latent heat flux.The series resistance network was the best for disaggregating soil... more
A two-source model based on surface temperature is tested in a Mediterranean dryland.Reliable estimates of sensible heat flux were achieved but not for latent heat flux.The series resistance network was the best for disaggregating soil and canopy fluxes.The iterative procedure included in the TSM was proven effective in drylands.A two-source model (TSM) for surface energy balance, considering explicitly soil and vegetation components, was tested under water stress conditions. The TSM evaluated estimates the sensible heat flux (H) using the surface-air thermal gradient and the latent heat flux (LE) as a residual from the surface energy balance equation. The analysis was performed in a semiarid Mediterranean tussock grassland in southeast Spain, where H is the dominant flux and LE rates are low, challenging conditions under which the TSM has not been validated before. We evaluated two different resistance schemes: series and parallel; as well as the iterative algorithm included in the TSM to disaggregate the soil-surface composite temperature into its separate components. Continuous field measurements of composite soil–vegetation surface temperature (TR) and bare soil temperature (Ts) from thermal infrared sensors were used for model testing along with canopy temperature estimates (T′c), derived from TR and Ts.Comparisons with Eddy covariance and field data showed that the TSM produced reliable estimates of net radiation (Rn) and H fluxes, with errors of ~ 30% and ~ 10%, respectively, but not for LE, with errors ~ 90%. Despite of lower errors (~ 10%) in estimating H using parallel resistance, the series scheme was more robust showing slightly higher correlations (r2 = 0.78–0.80 vs. r2 = 0.75–0.77) and allowing a better disaggregation of soil and canopy fluxes. Differences between model runs using the iterative algorithm to disaggregate TR and the simplified version that uses separate inputs of Ts and T′c were minor. This demonstrates the robustness of the iterative procedure to disaggregate a composite soil–vegetation temperature into separate soil and vegetation components in semiarid environments with good prospects for image applications.