Vincent Rivalland

An Analytical Model of Evaporation Efficiency for Unsaturated Soil Surfaces with an Arbitrary Thickness

Analytical expressions of evaporative efficiency over bare soil (defined as the ratio of actual to potential soil evaporation) have been limited to soil layers with a fixed depth and/or to specific atmospheric conditions. To fill the gap, a new analytical model is developed for arbitrary soil thicknesses and varying boundary layer conditions. The soil evaporative efficiency is written [0.5 – 0.5 cos(πθL/ θmax)]^P with θL being the water content in the soil layer of thickness L, θmax the soil moisture at saturation and P a function of L and potential soil evaporation. This formulation predicts soil evaporative efficiency in both energy-driven and moisture-driven conditions, which correspond to P 0.5 respectively. For P = 0.5, an equilibrium state is identified when retention forces in the soil compensate the evaporative demand above the soil surface. The approach is applied to in situ measurements of actual evaporation, potential evaporation and soil moisture at five different depths (5, 10, 30 and 60/100 cm) collected in summer at two sites in southwestern France. It is found that (i) soil evaporative efficiency cannot be considered as a function of soil moisture only, since it also depends on potential evaporation, (ii) retention forces in the soil increase in reaction to an increase of potential evaporation and (iii) the model is able to accurately predict soil evaporation process for soil layers with an arbitrary thickness up to 100 cm. This new model representation is expected to facilitate the coupling of land surface models with multi-sensor (multi-sensing-depth) remote sensing data.

On the Use of Hydrological Models and Satellite Data to Study the Water Budget of River Basins Affected by Human Activities: Examples from the Garonne Basin of France

Natural and anthropogenic forcing factors and their changes significantly impact water resources in many river basins around the world. Information on such changes can be derived from fine scale in situ and satellite observations, used in combination with hydrological models. The latter need to account for hydrological changes caused by human activities to correctly estimate the actual water resource. In this study, we consider the catchment area of the Garonne river (in France) to investigate the capabilities of space-based observations and up-to-date hydrological modeling in estimating water resources of a river basin modified by human activities and a changing climate. Using the ISBA–MODCOU and SWAT hydrological models, we find that the water resources of the Garonne basin display a negative climate trend since 1960. The snow component of the two models is validated using the moderate-resolution imaging spectroradiometer snow cover extent climatology. Crop sowing dates based on remote sensing studies are also considered in the validation procedure. Use of this dataset improves the simulated evapotranspiration and river discharge amounts when compared to conventional data. Finally, we investigate the benefit of using the MAELIA multi-agent model that accounts for a realistic agricultural and management scenario. Among other results, we find that changes in crop systems have significant impacts on water uptake for agriculture. This work constitutes a basis for the construction of a future modeling framework of the sociological and hydrological system of the Garonne river region.

A Software Tool for Atmospheric Correction and Surface Temperature Estimation of Landsat Infrared Thermal Data

Land surface temperature (LST) is an important variable involved in the Earths surface energy and water budgets and a key component in many aspects of environmental research. The Landsat program, jointly carried out by NASA and the USGS, has been recording thermal infrared data for the past 40 years. Nevertheless, LST data products for Landsat remain unavailable. The atmospheric correction (AC) method commonly used for mono-window Landsat thermal data requires detailed information concerning the vertical structure (temperature, pressure) and the composition (water vapor, ozone) of the atmosphere. For a given coordinate, this information is generally obtained through either radio-sounding or atmospheric model simulations and is passed to the radiative transfer model (RTM) to estimate the local atmospheric correction parameters. Although this approach yields accurate LST data, results are relevant only near this given coordinate. To meet the scientific communitys demand for high-resolution LST maps, we developed a new software tool dedicated to processing Landsat thermal data. The proposed tool improves on the commonly-used AC algorithm by incorporating spatial variations occurring in the Earths atmosphere composition. The ERA-Interim dataset (ECMWFmeteorological organization) was used to retrieve vertical atmospheric conditions, which are available at a global scale with a resolution of 0.125 degrees and a temporal resolution of 6 h. A temporal and spatial linear interpolation of meteorological variables was performed to match the acquisition dates and coordinates of the Landsat images. The atmospheric correction parameters were then estimated on the basis of this reconstructed atmospheric grid using the commercial RTMsoftware MODTRAN. The needed surface emissivity was derived from the common vegetation index NDVI, obtained from the red and near-infrared (NIR) bands of the same Landsat image. This permitted an estimation of LST for the entire image without degradation in resolution. The software tool, named LANDARTs, which stands for Landsat automatic retrieval of surface temperatures, is fully automatic and coded in the programming language Python. In the present paper, LANDARTs was used for the local and spatial validation of surface temperature obtained from a Landsat dataset covering two climatically contrasting zones: southwestern France and central Tunisia. Long-term datasets of in situ surface temperature measurements for both locations were compared to corresponding Landsat LST data. This temporal comparison yielded RMSE values ranging from 1.84 • C–2.55 • C. Landsat surface temperature data obtained with LANDARTs were then spatially compared using the ASTER data products of kinetic surface temperatures (AST08) for both geographical zones. This comparison yielded a satisfactory RMSE of about 2.55 • C. Finally, a sensitivity analysis for the effect of spatial validation on the LST correction process showed a variability of up to 2 • C for an entire Landsat image, confirming that the proposed spatial approach improved the accuracy of Landsat LST estimations.

Evaluation and aggregation properties of thermal infra-red-based evapotranspiration algorithms from 100 m to the km scale over a semi-arid irrigated agricultural area

Evapotranspiration (ET) estimates are particularly needed for monitoring the available water of arid lands. Remote sensing data offer the ideal spatial and temporal coverage needed by irrigation water management institutions to deal with increasing pressure on available water. Low spatial resolution (LR) products present strong advantages. They cover larger zones and are acquired more frequently than high spatial resolution (HR) products. Current sensors such as Moderate-Resolution Imaging Spectroradiometer (MODIS) offer a long record history. However, validation of ET products at LR remains a difficult task. In this context, the objective of this study is to evaluate scaling properties of ET fluxes obtained at high and low resolution by two commonly used Energy Balance models, the Surface Energy Balance System (SEBS) and the Two-Source Energy Balance model (TSEB). Both are forced by local meteorological observations and remote sensing data in Visible, Near Infra-Red and Thermal Infra-Red spectral domains. Remotely sensed data stem from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and MODIS sensors, respectively, resampled at 100 m and 1000 m resolutions. The study zone is a square area of 4 by 4 km2 located in a semi-arid irrigated agricultural zone in the northwest of Mexico. Wheat is the dominant crop, followed by maize and vegetables. The HR ASTER dataset includes seven dates between the 30 December 2007 and 13 May 2008 and the LR MODIS products were retrieved for the same overpasses. ET retrievals from HR ASTER products provided reference ET maps at LR once linearly aggregated at the km scale. The quality of this retrieval was assessed using eddy covariance data at seven locations within the 4 by 4 km2 square. To investigate the impact of input aggregation, we first compared to the reference dataset all fluxes obtained by running TSEB and SEBS models using ASTER reflectances and radiances previously aggregated at the km scale. Second, we compared to the same reference dataset all fluxes obtained with SEBS and TSEB models using MODIS data. LR fluxes obtained by both models driven by aggregated ASTER input data compared well with the reference simulations and illustrated the relatively good accuracy achieved using aggregated inputs (relative bias of about 3.5% for SEBS and decreased to less than 1% for TSEB). Results also showed that MODIS ET estimates compared well with the reference simulation (relative bias was down to about 2% for SEBS and 3% for TSEB). Discrepancies were mainly related to fraction cover mapping for TSEB and to surface roughness length mapping for SEBS. This was consistent with the sensitivity analysis of those parameters previously published. To improve accuracy from LR estimates obtained using the 1 km surface temperature product provided by MODIS, we tested three statistical and one deterministic aggregation rules for the most sensible input parameter, the surface roughness length. The harmonic and geometric averages appeared to be the most accurate.