Suivi des ressources en eau par télédétection multi-capteur: désagrégation de données spatiales et modélisation descendante des processus hydrologiques

Certaines variables extraites par télédétection sont potentiellement très utiles en hydrologie. C'est le cas de l'humidité du sol en surface qui contrôle la partition des précipitations en évaporation du sol, infiltration et ruissellement de surface. C'est aussi le cas de la température de surface qui lorsque l’énergie n’est pas limitante est une signature de l’évapotranspiration. Pourtant la résolution spatiale à laquelle ces données sont disponibles depuis l’espace n'est pas toujours compatible avec les échelles d'application. Dans ce contexte, la désagrégation de données apparaît comme un moyen d'améliorer la résolution spatiale des observations disponibles et de spatialiser les processus hydrologiques à des échelles multiples à l’aide de modèles descendants, c’est-à-dire basés sur les données spatiales. Au cours de cet exposé, je présenterai 1) des méthodes de désagrégation des données d’humidité du sol et de température de surface, 2) des modèles descendants de l’évaporation du sol et de l’évapotranspiration des surfaces et 3) une généralisation de ces approches permettant à long terme une spatialisation d’autres données spatiales et des processus hydrologiques associés

Coupled land surface and radiative transfer models for the analysis of passive microwave satellite observations

Soil moisture is one of the key variables controlling the water and energy exchanges between Earth’s surface and the atmosphere. Therefore, remote sensing based soil moisture information has potential applications in many disciplines. Besides numerical weather forecasting and climate research these include agriculture and hydrologic applications like flood and drought forecasting. The first satellite specifically designed to deliver operational soil moisture products, SMOS (Soil Moisture and Ocean Salinity), was launched 2009 by the European Space Agency (ESA). SMOS is a passive microwave radiometer working in the L-band of the microwave domain, corresponding to a frequency of roughly 1.4 GHz and relies on a new concept. The microwave radiation emitted by the Earth’s surface is measured as brightness temperatures in several look angles. A radiative transfer model is used in an inversion algorithm to retrieve soil moisture and vegetation optical depth, a measure for the vegetation attenuation of the soil’s microwave emission. For the application of passive microwave remote sensing products a proper validation and uncertainty assessment is essential. As these sensors have typical spatial resolutions in the order of 40 – 50 km, a validation that relies solely on ground measurements is costly and labour intensive. Here, environmental modelling can make a valuable contribution. Therefore the present thesis concentrates on the question which contribution coupled land surface and radiative transfer models can make to the validation and analysis of passive microwave remote sensing products. The objective is to study whether it is possible to explain known problems in the SMOS soil moisture products and to identify potential approaches to improve the data quality. The land surface model PROMET (PRocesses Of Mass and Energy Transfer) and the radiative transfer model L-MEB (L-band microwave emission of the Biosphere) are coupled to simulate land surface states, e.g. temperatures and soil moisture, and the resulting microwave emission. L-MEB is also used in the SMOS soil moisture processor to retrieve soil moisture and vegetation optical depth simultaneously from the measured microwave emission. The study area of this work is the Upper Danube Catchment, located mostly in Southern Germany. Since model validation is essential if model data are to be used as reference, both models are validated on different spatial scales with measurements. The uncertainties of the models are quantified. The root mean squared error between modelled and measured soil moisture at several measuring stations on the point scale is 0.065 m3/m3. On the SMOS scale it is 0.039 m3/m3. The correlation coefficient on the point scale is 0.84. As it is essential for the soil moisture retrieval from passive microwave data that the radiative transfer modelling works under local conditions, the coupled models are used to assess the radiative transfer modelling with L-MEB on the local and SMOS scales in the Upper Danube Catchment. In doing so, the emission characteristics of rape are described for the first time and the soil moisture retrieval abilities of L-MEB are assessed with a newly developed LMEB parameterization. The results show that the radiative transfer modelling works well under most conditions in the study area. The root mean squared error between modelled and airborne measured brightness temperatures on the SMOS scale is less than 6 – 9 K for the different look angles. The coupled models are used to analyse SMOS brightness temperatures and vegetation optical depth data in the Upper Danube Catchment in Southern Germany. Since the SMOS soil moisture products are degraded in Southern Germany and in different other parts of the world these analyses are used to narrow down possible reasons for this. The thorough analysis of SMOS brightness temperatures for the year 2011 reveals that the quality of the measurements is degraded like in the SMOS soil moisture product. This points towards radio frequency interference problems (RFI), that are known, but have not yet been studied thoroughly. This is consistent with the characteristics of the problems observed in the SMOS soil moisture products. In addition to that it is observed that the brightness temperatures in the lower look angles are less reliable. This finding could be used to improve the brightness temperature filtering before the soil moisture retrieval. An analysis of SMOS optical depth data in 2011 reveals that this parameter does not contain valuable information about vegetation. Instead, an unexpected correlation with SMOS soil moisture is found. This points towards problems with the SMOS soil moisture retrieval, possibly under the influence of RFI. The present thesis demonstrates that coupled land surface and radiative transfer models can make a valuable contribution to the validation and analysis of passive microwave remote sensing products. The unique approach of this work incorporates modelling with a high spatial and temporal resolution on different scales. This makes detailed process studies on the local scale as well as analyses of satellite data on the SMOS scale possible. This could be exploited for the validation of future satellite missions, e.g. SMAP (Soil Moisture Active and Passive) which is currently being prepared by NASA (National Aeronautics and Space Administration). Since RFI seems to have a considerable influence on the SMOS data due to the gained insights and the quality of the SMOS products is very good in other parts of the world, the RFI containment and mitigation efforts carried out since the launch of SMOS should be continued

Coupled land surface and radiative transfer models for the analysis of passive microwave satellite observations

Soil moisture is one of the key variables controlling the water and energy exchanges between Earth’s surface and the atmosphere. Therefore, remote sensing based soil moisture information has potential applications in many disciplines. Besides numerical weather forecasting and climate research these include agriculture and hydrologic applications like flood and drought forecasting. The first satellite specifically designed to deliver operational soil moisture products, SMOS (Soil Moisture and Ocean Salinity), was launched 2009 by the European Space Agency (ESA). SMOS is a passive microwave radiometer working in the L-band of the microwave domain, corresponding to a frequency of roughly 1.4 GHz and relies on a new concept. The microwave radiation emitted by the Earth’s surface is measured as brightness temperatures in several look angles. A radiative transfer model is used in an inversion algorithm to retrieve soil moisture and vegetation optical depth, a measure for the vegetation attenuation of the soil’s microwave emission. For the application of passive microwave remote sensing products a proper validation and uncertainty assessment is essential. As these sensors have typical spatial resolutions in the order of 40 – 50 km, a validation that relies solely on ground measurements is costly and labour intensive. Here, environmental modelling can make a valuable contribution. Therefore the present thesis concentrates on the question which contribution coupled land surface and radiative transfer models can make to the validation and analysis of passive microwave remote sensing products. The objective is to study whether it is possible to explain known problems in the SMOS soil moisture products and to identify potential approaches to improve the data quality. The land surface model PROMET (PRocesses Of Mass and Energy Transfer) and the radiative transfer model L-MEB (L-band microwave emission of the Biosphere) are coupled to simulate land surface states, e.g. temperatures and soil moisture, and the resulting microwave emission. L-MEB is also used in the SMOS soil moisture processor to retrieve soil moisture and vegetation optical depth simultaneously from the measured microwave emission. The study area of this work is the Upper Danube Catchment, located mostly in Southern Germany. Since model validation is essential if model data are to be used as reference, both models are validated on different spatial scales with measurements. The uncertainties of the models are quantified. The root mean squared error between modelled and measured soil moisture at several measuring stations on the point scale is 0.065 m3/m3. On the SMOS scale it is 0.039 m3/m3. The correlation coefficient on the point scale is 0.84. As it is essential for the soil moisture retrieval from passive microwave data that the radiative transfer modelling works under local conditions, the coupled models are used to assess the radiative transfer modelling with L-MEB on the local and SMOS scales in the Upper Danube Catchment. In doing so, the emission characteristics of rape are described for the first time and the soil moisture retrieval abilities of L-MEB are assessed with a newly developed LMEB parameterization. The results show that the radiative transfer modelling works well under most conditions in the study area. The root mean squared error between modelled and airborne measured brightness temperatures on the SMOS scale is less than 6 – 9 K for the different look angles. The coupled models are used to analyse SMOS brightness temperatures and vegetation optical depth data in the Upper Danube Catchment in Southern Germany. Since the SMOS soil moisture products are degraded in Southern Germany and in different other parts of the world these analyses are used to narrow down possible reasons for this. The thorough analysis of SMOS brightness temperatures for the year 2011 reveals that the quality of the measurements is degraded like in the SMOS soil moisture product. This points towards radio frequency interference problems (RFI), that are known, but have not yet been studied thoroughly. This is consistent with the characteristics of the problems observed in the SMOS soil moisture products. In addition to that it is observed that the brightness temperatures in the lower look angles are less reliable. This finding could be used to improve the brightness temperature filtering before the soil moisture retrieval. An analysis of SMOS optical depth data in 2011 reveals that this parameter does not contain valuable information about vegetation. Instead, an unexpected correlation with SMOS soil moisture is found. This points towards problems with the SMOS soil moisture retrieval, possibly under the influence of RFI. The present thesis demonstrates that coupled land surface and radiative transfer models can make a valuable contribution to the validation and analysis of passive microwave remote sensing products. The unique approach of this work incorporates modelling with a high spatial and temporal resolution on different scales. This makes detailed process studies on the local scale as well as analyses of satellite data on the SMOS scale possible. This could be exploited for the validation of future satellite missions, e.g. SMAP (Soil Moisture Active and Passive) which is currently being prepared by NASA (National Aeronautics and Space Administration). Since RFI seems to have a considerable influence on the SMOS data due to the gained insights and the quality of the SMOS products is very good in other parts of the world, the RFI containment and mitigation efforts carried out since the launch of SMOS should be continued

On the Statistical and Scaling Properties of Observed and Simulated Soil Moisture

abstract: Soil moisture (θ) is a fundamental variable controlling the exchange of water and energy at the land surface. As a result, the characterization of the statistical properties of θ across multiple scales is essential for many applications including flood prediction, drought monitoring, and weather forecasting. Empirical evidences have demonstrated the existence of emergent relationships and scale invariance properties in θ fields collected from the ground and airborne sensors during intensive field campaigns, mostly in natural landscapes. This dissertation advances the characterization of these relations and statistical properties of θ by (1) analyzing the role of irrigation, and (2) investigating how these properties change in time and across different landscape conditions through θ outputs of a distributed hydrologic model. First, θ observations from two field campaigns in Australia are used to explore how the presence of irrigated fields modifies the spatial distribution of θ and the associated scale invariance properties. Results reveal that the impact of irrigation is larger in drier regions or conditions, where irrigation creates a drastic contrast with the surrounding areas. Second, a physically-based distributed hydrologic model is applied in a regional basin in northern Mexico to generate hyperresolution θ fields, which are useful to conduct analyses in regions and times where θ has not been monitored. For this aim, strategies are proposed to address data, model validation, and computational challenges associated with hyperresolution hydrologic simulations. Third, analyses are carried out to investigate whether the hyperresolution simulated θ fields reproduce the statistical and scaling properties observed from the ground or remote sensors. Results confirm that (i) the relations between spatial mean and standard deviation of θ derived from the model outputs are very similar to those observed in other areas, and (ii) simulated θ fields exhibit the scale invariance properties that are consistent with those analyzed from aircraft-derived estimates. The simulated θ fields are then used to explore the influence of physical controls on the statistical properties, finding that soil properties significantly affect spatial variability and multifractality. The knowledge acquired through this dissertation provides insights on θ statistical properties in regions and landscape conditions that were never investigated before; supports the refinement of the calibration of multifractal downscaling models; and contributes to the improvement of hyperresolution hydrologic modeling.Dissertation/ThesisDoctoral Dissertation Civil, Environmental and Sustainable Engineering 201

Estimation de l'humidité du sol à haute résolution spatio-temporelle : une nouvelle approche basée sur la synergie des observations micro-ondes actives/passives et optiques/thermiques

Les capteurs micro-ondes passifs SMOS et SMAP fournissent des données d'humidité du sol (SM) à une résolution d'environ 40 km avec un intervalle de 2 à 3 jours à l' échelle mondiale et une profondeur de détection de 0 à 5 cm. Ces données sont très pertinentes pour les applications cli- matiques et météorologiques. Cependant, pour les applications à échelle régionales (l'hydrologie) ou locales (l'agriculture), des données de SM à une haute résolution spatiale (typiquement 100 m ou plus fine) seraient nécessaires. Les données collectées par les capteurs optiques/thermiques et les radars peuvent fournir des indicateurs de SM à haute résolution spatiale, mais ces deux approches alternatives ont des limites. En particulier, les données optiques/thermiques ne sont pas disponibles sous les nuages et sous les couverts végétaux. Quant aux données radar, elles sont sensibles à la rugosité du sol et à la structure de la végétation, qui sont tous deux difficiles à caractériser depuis l'espace. De plus, la résolution temporelle de ces données est d'environ 6 jours. Dans ce contexte, la ligne directrice de la thèse est de proposer une nouvelle approche qui combine pour la première fois des capteurs passifs micro-ondes, optiques/thermiques et actifs micro-ondes (radar) pour estimer SM sur de grandes étendues à une résolution de 100 m chaque jour. Notre hypothèse est d'abord de nous appuyer sur une méthode de désagrégation existante (DISPATCH) des données SMOS/SMAP pour atteindre la résolution cible obtenue par les radars. A l'origine, DISPATCH est basé sur l'efficacité d' évaporation du sol (SEE) estimée sur des pixels partiellement végétalisés à partir de données optiques/thermiques (généralement MODIS) de température de surface et de couverture végétale à résolution de 1 km. Les données désagrégées de SM sont ensuite combinées avec une méthode d'inversion de SM basée sur les données radar afin d'exploiter les capacités de détection des radars Sentinel-1. Enfin, les capacités de l'assimilation des donnés satellitaires de SM dans un modèle de bilan hydrique du sol sont évaluées en termes de prédiction de SM à une résolution de 100 m et à une échelle temporelle quotidienne.Dans une première étape, l'algorithme DISPATCH est amélioré par rapport à sa version actuelle, principalement 1) en étendant son applicabilité aux pixels optiques entièrement végétalisés en utilisant l'indice de sécheresse de la végétation basé sur la température et un produit de couverture végétale amélioré, et 2) en augmentant la résolution de désagrégation de 1 km à 100 m en utilisant les données optiques/thermiques de Landsat (en plus de MODIS). Le produit de SM désagrégé à la résolution de 100 m est validé avec des mesures in situ collectées sur des zones irriguées au Maroc, indiquant une corrélation spatiale quotidienne variant de 0,5 à 0,9. Dans un deuxième étape, un nouvel algorithme est construit en développant une synergie entre les données DISPATCH et radar à 100 m de résolution. En pratique, le produit SM issu de DISPATCH les jours de ciel clair est d'abord utilisé pour calibrer un modèle de transfert radiatif radar en mode direct. Ensuite, le modèle de transfert radiatif radar ainsi calibré est utilisé en mode inverse pour estimer SM à la résolution spatio-temporelle de Sentinel-1. Sur les sites de validation, les résultats indiquent une corrélation entre les mesures satellitaires et in situ, de l'ordre de 0,66 à 0,81 pour un indice de végétation inférieur à 0,6. Dans une troisième et dernière étape, une méthode d'assimilation optimale est utilisée pour interpoler dans le temps les données de SM à la résolution de 100 m. La dynamique du produit SM dérivé de l'assimilation de SM DISPATCH à 100 m de résolution est cohérente avec les événements d'irrigation. Cette approche peut être facilement appliquée sur de grandes zones, en considérant que toutes les données (télédétection et météorologique) requises en entrée sont disponibles à l' échelle globale.SMOS and SMAP passive microwave sensors provide soil moisture (SM) data at 40 km resolution every 2-3 days globally, with a 0-5 cm sensing depth relevant for climatic and meteorological applications. However, SM data would be required at a higher (typically 100 m or finer) spatial resolution for many other regional (hydrology) or local (agriculture) applications. Optical/thermal and radar sensors can be used for retrieving SM proxies at such high spatial resolution, but both techniques have limitations. In particular, optical/thermal data are not available under clouds and under plant canopies. Moreover, radar data are sensitive to soil roughness and vegetation structure, which are challenging to characterize from outer space, and have a repeat cycle of at least six days, limiting the observations' temporal frequency. In this context, the leading principle of the thesis is to propose a new approach that combines passive microwave, optical/thermal, and active microwave (radar) sensors for the first time to retrieve SM data at 100 m resolution on a daily temporal scale. Our assumption is first to rely on an existing disaggregation method (DISPATCH) of SMOS/SMAP SM data to meet the target resolution achieved by radars. DISPATCH is originally based on the soil evaporative efficiency (SEE) retrieved over partially vegetated pixels from 1 km resolution optical/thermal (typically MODIS) surface temperature and vegetation cover data. The disaggregated SM data is then combined with a radar-based SM retrieval method to exploit the sensing capabilities of the Sentinel-1 radars. Finally, the efficacy of the assimilation of satellite-based SM data in a soil water balance model is assessed in terms of SM predictions at the 100 m resolution and daily temporal scale. As a first step, the DISPATCH algorithm is improved from its current version by mainly 1) extending its applicability to fully vegetated optical pixels using the temperature vegetation dryness index and an enhanced vegetation cover product, and 2) increasing the targeted downscaling resolution from 1 km to 100 m using Landsat (in addition to MODIS) optical/thermal data. The 100 m resolution disaggregated SM product is validated with in situ measurements collected over irrigated areas in Morocco, showing a daily spatial correlation in the range of 0.5-0.9. As a second step, a new algorithm is built on a synergy between DISPATCH and radar 100 m resolution data. In practice, the DISPATCH SM product available on clear sky days is first used to calibrate a radar radiative transfer model in the direct mode. Then the calibrated radar radia- tive transfer model is used in the inverse mode to estimate SM at the spatio-temporal resolution of Sentinel-1. Results indicate a positive correlation between satellite and in situ measurements in the range of 0.66 to 0.81 for a vegetation index lower than 0.6. As a third and final step, an optimal assimilation method is used to interpolate 100 m resolution SM data in time. The assimilation exercise is undertaken over irrigated crop fields in Spain. The analyzed SM product derived from the assimilation of 100 m resolution DISPATCH SM is consistent with irrigation events. This approach can be readily applied over large areas, given that all the required input (remote sensing and meteorological) data are available globally