Influence of Surface Roughness Spatial Variability and Temporal Dynamics on the Retrieval of Soil Moisture from SAR Observations

Radar-based surface soil moisture retrieval has been subject of intense research during the last decades. However, several difficulties hamper the operational estimation of soil moisture based on currently available spaceborne sensors. The main difficulty experienced so far results from the strong influence of other surface characteristics, mainly roughness, on the backscattering coefficient, which hinders the soil moisture inversion. This is especially true for single configuration observations where the solution to the surface backscattering problem is ill-posed. Over agricultural areas cultivated with winter cereal crops, roughness can be assumed to remain constant along the growing cycle allowing the use of simplified approaches that facilitate the estimation of the moisture content of soils. However, the field scale spatial variability and temporal variations of roughness can introduce errors in the estimation of soil moisture that are difficult to evaluate. The objective of this study is to assess the impact of roughness spatial variability and roughness temporal variations on the retrieval of soil moisture from radar observations. A series of laser profilometer measurements were performed over several fields in an experimental watershed from September 2004 to March 2005. The influence of the observed roughness variability and its temporal variations on the retrieval of soil moisture is studied using simulations performed with the Integral Equation Model, considering different sensor configurations. Results show that both field scale roughness spatial variability and its temporal variations are aspects that need to be taken into account, since they can introduce large errors on the retrieved soil moisture values

Coupled modelling of land surface microwave interactions using ENVISAT ASAR data

In the last decades microwave remote sensing has proven its capability to provide valuable information about the land surface. New sensor generations as e.g. ENVISAT ASAR are capable to provide frequent imagery with an high information content. To make use of these multiple imaging capabilities, sophisticated parameter inversion and assimilation strategies have to be applied. A profound understanding of the microwave interactions at the land surface is therefore essential. The objective of the presented work is the analysis and quantitative description of the backscattering processes of vegetated areas by means of microwave backscattering models. The effect of changing imaging geometries is investigated and models for the description of bare soil and vegetation backscattering are developed. Spatially distributed model parameterisation is realized by synergistic coupling of the microwave scattering models with a physically based land surface process model. This enables the simulation of realistic SAR images, based on bioand geophysical parameters. The adequate preprocessing of the datasets is crucial for quantitative image analysis. A stringent preprocessing and sophisticated terrain geocoding and correction procedure is therefore suggested. It corrects the geometric and radiometric distortions of the image products and is taken as the basis for further analysis steps. A problem in recently available microwave backscattering models is the inadequate parameterisation of the surface roughness. It is shown, that the use of classical roughness descriptors, as the rms height and autocorrelation length, will lead to ambiguous model parameterisations. A new two parameter bare soil backscattering model is therefore recommended to overcome this drawback. It is derived from theoretical electromagnetic model simulations. The new bare soil surface scattering model allows for the accurate description of the bare soil backscattering coefficients. A new surface roughness parameter is introduced in this context, capable to describe the surface roughness components, affecting the backscattering coefficient. It is shown, that this parameter can be directly related to the intrinsic fractal properties of the surface. Spatially distributed information about the surface roughness is needed to derive land surface parameters from SAR imagery. An algorithm for the derivation of the new surface roughness parameter is therefore suggested. It is shown, that it can be derived directly from multitemporal SAR imagery. Starting from that point, the bare soil backscattering model is used to assess the vegetation influence on the signal. By comparison of the residuals between measured backscattering coefficients and those predicted by the bare soil backscattering model, the vegetation influence on the signal can be quantified. Significant difference between cereals (wheat and triticale) and maize is observed in this context. It is shown, that the vegetation influence on the signal can be directly derived from alternating polarisation data for cereal fields. It is dependant on plant biophysical variables as vegetation biomass and water content. The backscattering behaviour of a maize stand is significantly different from that of other cereals, due to its completely different density and shape of the plants. A dihedral corner reflection between the soil and the stalk is identified as the major source of backscattering from the vegetation. A semiempirical maize backscattering model is suggested to quantify the influences of the canopy over the vegetation period. Thus, the different scattering contributions of the soil and vegetation components are successfully separated. The combination of the bare soil and vegetation backscattering models allows for the accurate prediction of the backscattering coefficient for a wide range of surface conditions and variable incidence angles. To enable the spatially distributed simulation of the SAR backscattering coefficient, an interface to a process oriented land surface model is established, which provides the necessary input variables for the backscattering model. Using this synergistic, coupled modelling approach, a realistic simulation of SAR images becomes possible based on land surface model output variables. It is shown, that this coupled modelling approach leads to promising and accurate estimates of the backscattering coefficients. The remaining residuals between simulated and measured backscatter values are analysed to identify the sources of uncertainty in the model. A detailed field based analysis of the simulation results revealed that imprecise soil moisture predictions by the land surface model are a major source of uncertainty, which can be related to imprecise soil texture distribution and soil hydrological properties. The sensitivity of the backscattering coefficient to the soil moisture content of the upper soil layer can be used to generate soil moisture maps from SAR imagery. An algorithm for the inversion of soil moisture from the upper soil layer is suggested and validated. It makes use of initial soil moisture values, provided by the land surface process model. Soil moisture values are inverted by means of the coupled land surface backscattering model. The retrieved soil moisture results have an RMSE of 3.5 Vol %, which is comparable to the measurement accuracy of the reference field data. The developed models allow for the accurate prediction of the SAR backscattering coefficient. The various soil and vegetation scattering contributions can be separated. The direct interface to a physically based land surface process model allows for the spatially distributed modelling of the backscattering coefficient and the direct assimilation of remote sensing data into a land surface process model. The developed models allow for the derivation of static and dynamic landsurface parameters, as e.g. surface roughness, soil texture, soil moisture and biomass from remote sensing data and their assimilation in process models. They are therefore reliable tools, which can be used for sophisticated practice oriented problem solutions in manifold manner in the earth and environmental sciences

A model to estimate daily albedo from remote sensing data : accuracy assessment of MODIS MCD43 product

L’albedo superficial és un paràmetre físic que afecta al clima de la Terra i, a més, suposa una de les majors incerteses radiatives en l’actual modelització climàtica. Aquest paràmetre és molt variable tant a nivell espacial com temporal degut als canvis en les propietats de les superfícies i als canvis en les condicions d’il•luminació. En conseqüència, es requereix una resolució temporal diària de l’albedo per a realitzar estudis climàtics. L’augment de la resolució espacial dels models climàtics fa necessari l’estudi de les seues característiques espacials a nivell global. La teledetecció proporciona l’única opció pràctica de proporcionar dades d’albedo a nivell global amb alta qualitat i alta resolució tant espacial com temporal. El sensor MODerate Resolution Imaging Spectroradiometer (MODIS) a bord dels satèl•lits Terra i Aqua presenta unes característiques adequades per a l’estimació d’aquest paràmetre. En el present treball realitzem diversos estudis buscant les possibles fonts d’incertesa del producte oficial d’albedo de MODIS (MCD43). A més, presentem un model que millora la resolució temporal d’aquest paràmetre.Surface albedo is a critical land physical parameter affecting the earth’s climate and is among the main radiative uncertainties in current climate modelling. This parameter is highly variable in space and time, both as a result of changes in surface properties and as a function of changes in the illumination conditions. Consequently, an albedo daily temporal resolution is required for climate studies. The increasing spatial resolution of modern climate models makes it necessary to examine its spatial features. Satellite remote sensing provides the only practical way of producing high-quality global albedo data sets with high spatial and temporal resolutions. The MODerate Resolution Imaging Spectroradiometer (MODIS) sensor on board the Terra and Aqua satellites presents the required sampling characteristics in order to derive the this parameter. In this PhD we develop several studies looking for the improvement of the official MODIS albedo product (MCD43) accuracy. Moreover, we present a model that improves the temporal resolution of this parameter

Harmonization of remote sensing land surface products : correction of clear-sky bias and characterization of directional effects

Tese de doutoramento, Ciências Geofísicas e da Geoinformação (Deteção Remota), Universidade de Lisboa, Faculdade de Ciências, 2018Land surface temperature (LST) is the mean radiative skin temperature of an area of land resulting from the mean energy balance at the surface. LST is an important climatological variable and a diagnostic parameter of land surface conditions, since it is the primary variable determining the upward thermal radiation and one of the main controllers of sensible and latent heat fluxes between the surface and the atmosphere. The reliable and long-term estimation of LST is therefore highly relevant for a wide range of applications, including, amongst others: (i) land surface model validation and monitoring; (ii) data assimilation; (iii) hydrological applications; and (iv) climate monitoring. Remote sensing constitutes the most effective method to observe LST over large areas and on a regular basis. Satellite LST products generally rely on measurements in the thermal infrared (IR) atmospheric window, i.e., within the 8-13 micrometer range. Beside the relatively weak atmospheric attenuation under clear sky conditions, this band includes the peak of the Earth’s spectral radiance, considering surface temperature of the order of 300K (leading to maximum emission at approximately 9.6 micrometer, according to Wien’s Displacement Law). The estimation of LST from remote sensing instruments operating in the IR is being routinely performed for nearly 3 decades. Nevertheless, there is still a long list of open issues, some of them to be addressed in this PhD thesis. First, the viewing position of the different remote sensing platforms may lead to variability of the retrieved surface temperature that depends on the surface heterogeneity of the pixel – dominant land cover, orography. This effect introduces significant discrepancies among LST estimations from different sensors, overlapping in space and time, that are not related to uncertainties in the methodologies or input data used. Furthermore, these directional effects deviate LST products from an ideally defined LST, which should correspond to the ensemble directional radiometric temperature of all surface elements within the FOV. In this thesis, a geometric model is presented that allows the upscaling of in situ measurements to the any viewing configuration. This model allowed generating a synthetic database of directional LST that was used consistently to evaluate different parametric models of directional LST. Ultimately, a methodology is proposed that allows the operational use of such parametric models to correct angular effects on the retrieved LST. Second, the use of infrared data limits the retrieval of LST to clear sky conditions, since clouds “close” the atmospheric window. This effect introduces a clear-sky bias in IR LST datasets that is difficult to quantify since it varies in space and time. In addition, the cloud clearing requirement severely limits the space-time sampling of IR measurements. Passive microwave (MW) measurements are much less affected by clouds than IR observations. LST estimates can in principle be derived from MW measurements, regardless of the cloud conditions. However, retrieving LST from MW and matching those estimations with IR-derived values is challenging and there have been only a few attempts so far. In this thesis, a methodology is presented to retrieve LST from passive MW observations. The MW LST dataset is examined comprehensively against in situ measurements and multiple IR LST products. Finally, the MW LST data is used to assess the spatial-temporal patterns of the clear-sky bias at global scale.Fundação para a Ciência e a Tecnologia, SFRH/BD/9646

Analysis of surface soil moisture patterns in an agricultural landscape utilizing measurements and ecohydrological modeling

Soil moisture and its distribution in space and time plays a decisive role in terrestrial water and energy cycles. It controls the partitioning of precipitation into infiltration and runoff as well as the partitioning of solar radiation into latent and sensible heat flux. Therefore it has a strong impact on numerous processes, e.g., controlling floods, crop yield, erosion, and climate processes. Soil moisture, and surface soil moisture in particular, is highly variable in space and time and its spatial and temporal patterns in an agricultural landscape are affected by multiple natural (precipitation, soil, etc.) and agricultural (soil management, fertilization etc.) parameters. Against this background, the current study investigates the spatial and temporal patterns of surface soil moisture in an agricultural landscape, to determine the dominant parameters and the underlying processes controlling these patterns. The study was conducted on different spatial scales, from the field scale to the whole catchment scale of the river Rur (2364 km2) in Western Germany, because observed patterns are intrinsically connected to the scale on which they are observed. For the investigation three different approaches were used: Analysis based on A) Field measurements, B) Radar remote sensing, and C) Ecohydrological modeling. Extensive field measurements were carried out in a small arable land and grassland test site, measuring surface soil moisture, plant parameters, meteorological parameters, and soil parameters. These measurements were used to analyze the small scale (field scale) patterns of surface soil moisture and for the validation of the two other methods. Since large scale investigations based on field measurements are generally not feasible, surface soil moisture maps from radar remote sensing and ecohydrological modeling were used to analyze large scale patterns of surface soil moisture and their scaling properties. Precipitation, vegetation patterns, topography and soil properties were found to be the dominant parameters for soil moisture patterns in an agriculturally used landscape. Precipitation can be assumed to be homogeneous on the small scale, but can be very heterogeneous on the large scale at the same time. Evapotranspiration causes high small scale variability, especially during the growing season. If analyzed on coarser resolutions, this small scale pattern is smoothed out. Topography is a source of small scale patterns only on wet surface soil moisture states, because of the lateral redistribution of water during or shortly after precipitation events. Soils have a major influence on the variability of surface soil moisture on all scales, due to the large heterogeneity of soil properties within a given soil type (small scale) and between different soil types (large scale). Altogether, the variability of surface soil moisture increases with an increasing size of the investigation area and with an increasing resolution within the investigation area. During the course of the year surface soil moisture variability and its scaling properties are being influenced by different parameters with temporally varying intensities. During the growing season, at the time of high small scale variability of evapotranspiration, the variability of surface soil moisture is high and decreases much stronger with decreasing spatial resolution of the investigation, than during times outside the growing season. In the beginning and towards the end of the year (outside the growing season, when the soil is wet) the patterns and their scaling properties are mainly determined by soil properties. Precipitation events generally superimpose their large scale patterns for a short period of time and diminish the small scale variability induced by evapotranspiration. This thesis improves the knowledge about surface soil moisture patterns in agriculturally used areas and their underlying processes. The results of the scaling analysis indicate the potential to use vegetation and precipitation parameters for downscaling purposes. Understanding the subscale soil moisture heterogeneity is, for example, particularly relevant to better utilize coarse scale soil moisture data derived from SMOS (Soil Moisture and Ocean Salinity) or the upcoming SMAP (Soil Moisture Active Passive) satellite measurements

Monitoring crops water needs at high spatio-temporal resolution by synergy of optical/thermal and radar observations

L'optimisation de la gestion de l'eau en agriculture est essentielle dans les zones semi-arides afin de préserver les ressources en eau qui sont déjà faibles et erratiques dues à des actions humaines et au changement climatique. Cette thèse vise à utiliser la synergie des observations de télédétection multispectrales (données radar, optiques et thermiques) pour un suivi à haute résolution spatio-temporelle des besoins en eau des cultures. Dans ce contexte, différentes approches utilisant divers capteurs (Landsat-7/8, Sentinel-1 et MODIS) ont été developpées pour apporter une information sur l'humidité du sol (SM) et le stress hydrique des cultures à une échelle spatio-temporelle pertinente pour la gestion de l'irrigation. Ce travail va parfaitement dans le sens des objectifs du projet REC "Root zone soil moisture Estimates at the daily and agricultural parcel scales for Crop irrigation management and water use impact: a multi-sensor remote sensing approach" (http://rec.isardsat.com/) qui visent à estimer l'humidité du sol dans la zone racinaire (RZSM) afin d'optimiser la gestion de l'eau d'irrigation. Des approches innovantes et prometteuses sont mises en place pour estimer l'évapotranspiration (ET), RZSM, la température de surface du sol (LST) et le stress hydrique de la végétation à travers des indices de SM dérivés des observations multispectrales à haute résolution spatio-temporelle. Les méthodologies proposées reposent sur des méthodes basées sur l'imagerie, la modélisation du transfert radiatif et la modélisation du bilan hydrique et d'énergie et sont appliquées dans une région à climat semi-aride (centre du Maroc). Dans le cadre de ma thèse, trois axes ont été explorés. Dans le premier axe, un indice de RZSM dérivé de LST-Landsat est utilisé pour estimer l'ET sur des parcelles de blé et des sols nus. L'estimation par modélisation de ET a été explorée en utilisant l'équation de Penman-monteith modifiée obtenue en introduisant une relation empirique simple entre la résistance de surface (rc) et l'indice de RZSM. Ce dernier est estimé à partir de la température de surface (LST) dérivée de Landsat, combinée avec les températures extrêmes (en conditions humides et sèches) simulée par un modèle de bilan d'énergie de surface piloté par le forçage météorologique et la fraction de couverture végétale dérivée de Landsat. La méthode utilisée est calibrée et validée sur deux parcelles de blé situées dans la même zone près de Marrakech au Maroc. Dans l'axe suivant, une méthode permettant de récupérer la SM de la surface (0-5 cm) à une résolution spatiale et temporelle élevée est développée à partir d'une synergie entre données radar (Sentinel-1) et thermique (Landsat) et en utilisant un modèle de bilan d'énergie du sol. L'approche développée a été validée sur des parcelles agricoles en sol nu et elle donne une estimation précise de la SM avec une différence quadratique moyenne en comparant à la SM in situ, égale à 0,03 m3 m-3. Dans le dernier axe, une nouvelle méthode est développée pour désagréger la MODIS LST de 1 km à 100 m de résolution en intégrant le SM proche de la surface dérivée des données radar Sentinel-1 et l'indice de végétation optique dérivé des observations Landsat. Le nouvel algorithme, qui inclut la rétrodiffusion S-1 en tant qu'entrée dans la désagrégation, produit des résultats plus stables et robustes au cours de l'année sélectionnée. Dont, 3,35 °C était le RMSE le plus bas et 0,75 le coefficient de corrélation le plus élevé évalués en utilisant le nouvel algorithme.Optimizing water management in agriculture is essential over semi-arid areas in order to preserve water resources which are already low and erratic due to human actions and climate change. This thesis aims to use the synergy of multispectral remote sensing observations (radar, optical and thermal data) for high spatio-temporal resolution monitoring of crops water needs. In this context, different approaches using various sensors (Landsat-7/8, Sentinel-1 and MODIS) have been developed to provide information on the crop Soil Moisture (SM) and water stress at a spatio-temporal scale relevant to irrigation management. This work fits well the REC "Root zone soil moisture Estimates at the daily and agricultural parcel scales for Crop irrigation management and water use impact: a multi-sensor remote sensing approach" (http://rec.isardsat.com/) project objectives, which aim to estimate the Root Zone Soil Moisture (RZSM) for optimizing the management of irrigation water. Innovative and promising approaches are set up to estimate evapotranspiration (ET), RZSM, land surface temperature (LST) and vegetation water stress through SM indices derived from multispectral observations with high spatio-temporal resolution. The proposed methodologies rely on image-based methods, radiative transfer modelling and water and energy balance modelling and are applied in a semi-arid climate region (central Morocco). In the frame of my PhD thesis, three axes have been investigated. In the first axis, a Landsat LST-derived RZSM index is used to estimate the ET over wheat parcels and bare soil. The ET modelling estimation is explored using a modified Penman-Monteith equation obtained by introducing a simple empirical relationship between surface resistance (rc) and a RZSM index. The later is estimated from Landsat-derived land surface temperature (LST) combined with the LST endmembers (in wet and dry conditions) simulated by a surface energy balance model driven by meteorological forcing and Landsat-derived fractional vegetation cover. The investigated method is calibrated and validated over two wheat parcels located in the same area near Marrakech City in Morocco. In the next axis, a method to retrieve near surface (0-5 cm) SM at high spatial and temporal resolution is developed from a synergy between radar (Sentinel-1) and thermal (Landsat) data and by using a soil energy balance model. The developed approach is validated over bare soil agricultural fields and gives an accurate estimates of near surface SM with a root mean square difference compared to in situ SM equal to 0.03 m3 m-3. In the final axis a new method is developed to disaggregate the 1 km resolution MODIS LST at 100 m resolution by integrating the near surface SM derived from Sentinel-1 radar data and the optical-vegetation index derived from Landsat observations. The new algorithm including the S-1 backscatter as input to the disaggregation, produces more stable and robust results during the selected year. Where, 3.35 °C and 0.75 were the lowest RMSE and the highest correlation coefficient assessed using the new algorithm

Research reports: 1991 NASA/ASEE Summer Faculty Fellowship Program

The basic objectives of the programs, which are in the 28th year of operation nationally, are: (1) to further the professional knowledge of qualified engineering and science faculty members; (2) to stimulate an exchange of ideas between participants and NASA; (3) to enrich and refresh the research and teaching activities of the participants' institutions; and (4) to contribute to the research objectives of the NASA Centers. The faculty fellows spent 10 weeks at MSFC engaged in a research project compatible with their interests and background and worked in collaboration with a NASA/MSFC colleague. This is a compilation of their research reports for summer 1991

Summaries of the Sixth Annual JPL Airborne Earth Science Workshop

This publication contains the summaries for the Sixth Annual JPL Airborne Earth Science Workshop, held in Pasadena, California, on March 4-8, 1996. The main workshop is divided into two smaller workshops as follows: (1) The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) workshop, on March 4-6. The summaries for this workshop appear in Volume 1; (2) The Airborne Synthetic Aperture Radar (AIRSAR) workshop, on March 6-8. The summaries for this workshop appear in Volume 2