Impact of day/night time land surface temperature in soil moisture disaggregation algorithms

Since its launch in 2009, the ESA’s SMOS mission is providing global soil moisture (SM) maps at ~40 km, using the first L-band microwave radiometer on space. Its spatial resolution meets the needs of global applications, but prevents the use of the data in regional or local applications, which require higher spatial resolutions (~1-10 km). SM disaggregation algorithms based generally on the land surface temperature (LST) and vegetation indices have been developed to bridge this gap. This study analyzes the SM-LST relationship at a variety of LST acquisition times and its influence on SM disaggregation algorithms. Two years of in situ and satellite data over the central part of the river Duero basin and the Iberian Peninsula are used. In situ results show a strong anticorrelation of SM to daily maximum LST (R˜-0.5 to -0.8). This is confirmed with SMOS SM and MODIS LST Terra/Aqua at day time-overpasses (R˜-0.4 to -0.7). Better statistics are obtained when using MODIS LST day (R˜0.55 to 0.85; ubRMSD˜0.04 to 0.06 m3 /m3 ) than LST night (R˜0.45 to 0.80; ubRMSD˜0.04 to 0.07 m3 /m3 ) in the SM disaggregation. An averaged ensemble of day and night MODIS LST Terra/Aqua disaggregated SM estimates also leads to robust statistics (R˜0.55 to 0.85; ubRMSD˜0.04 to 0.07 m3 /m3 ) with a coverage improvement of ~10-20 %.Peer ReviewedPostprint (published version

An Initial Assessment of a SMAP Soil Moisture Disaggregation Scheme Using TIR Surface Evaporation Data over the Continental United States

The Soil Moisture Active Passive (SMAP) mission is dedicated toward global soil moisture mapping. Typically, an L-band microwave radiometer has spatial resolution on the order of 36-40 km, which is too coarse for many specific hydro-meteorological and agricultural applications. With the failure of the SMAP active radar within three months of becoming operational, an intermediate (9-km) and finer (3-km) scale soil moisture product solely from the SMAP mission is no longer possible. Therefore, the focus of this study is a disaggregation of the 36-km resolution SMAP passive-only surface soil moisture (SSM) using the Soil Evaporative Efficiency (SEE) approach to spatial scales of 3-km and 9-km. The SEE was computed using thermal-infrared (TIR) estimation of surface evaporation over Continental U.S. (CONUS). The disaggregation results were compared with the 3 months of SMAP-Active (SMAP-A) and Active/Passive (AP) products, while comparisons with SMAP-Enhanced (SMAP-E), SMAP-Passive (SMAP-P), as well as with more than 180 Soil Climate Analysis Network (SCAN) stations across CONUS were performed for a 19 month period. At the 9-km spatial scale, the TIR-Downscaled data correlated strongly with the SMAP-E SSM both spatially (r = 0.90) and temporally (r = 0.87). In comparison with SCAN observations, overall correlations of 0.49 and 0.47; bias of 0.022 and 0.019 and unbiased RMSD of 0.105 and 0.100 were found for SMAP-E and TIR-Downscaled SSM across the Continental U.S., respectively. At 3-km scale, TIR-Downscaled and SMAP-A had a mean temporal correlation of only 0.27. In terms of gain statistics, the highest percentage of SCAN sites with positive gains (>55%) was observed with the TIR-Downscaled SSM at 9-km. Overall, the TIR-based downscaled SSM showed strong correspondence with SMAP-E; compared to SCAN, and overall both SMAP-E and TIR-Downscaled performed similarly, however, gain statistics show that TIR-Downscaled SSM slightly outperformed SMAP-E

SMOS based high resolution soil moisture estimates for Desert locust preventive management

This paper presents the first attempt to include soil moisture information from remote sensing in the tools available to desert locust managers. The soil moisture requirements were first assessed with the users. The main objectives of this paper are: i) to describe and validate the algorithms used to produce a soil moisture dataset at 1 km resolution relevant to desert locust management based on DisPATCh methodology applied to SMOS and ii) the development of an innovative approach to derive high-resolution (100 m) soil moisture products from Sentinel-1 in synergy with SMOS data. For the purpose of soil moisture validation, 4 soil moisture stations where installed in desert areas (one in each user country). The soil moisture 1 km product was thoroughly validated and its accuracy is amongst the best available soil moisture products. Current comparison with in-situ soil moisture stations shows good values of correlation (R>0.7R>0.7) and low RMSE (below 0.04 m3 m−3). The low number of acquisitions on wet dates has limited the development of the soil moisture 100 m product over the Users Areas. The Soil Moisture product at 1 km will be integrated into the national and global Desert Locust early warning systems in national locust centres and at DLIS-FAO, respectively

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

Multi-temporal evaluation of soil moisture and land surface temperature dynamics using in situ and satellite observations

Soil moisture (SM) is an important component of the Earth’s surface water balance and by extension the energy balance, regulating the land surface temperature (LST) and evapotranspiration (ET). Nowadays, there are two missions dedicated to monitoring the Earth’s surface SM using L-band radiometers: ESA’s Soil Moisture and Ocean Salinity (SMOS) and NASA’s Soil Moisture Active Passive (SMAP). LST is remotely sensed using thermal infrared (TIR) sensors on-board satellites, such as NASA’s Terra/Aqua MODIS or ESA & EUMETSAT’s MSG SEVIRI. This study provides an assessment of SM and LST dynamics at daily and seasonal scales, using 4 years (2011–2014) of in situ and satellite observations over the central part of the river Duero basin in Spain. Specifically, the agreement of instantaneous SM with a variety of LST-derived parameters is analyzed to better understand the fundamental link of the SM–LST relationship through ET and thermal inertia. Ground-based SM and LST measurements from the REMEDHUS network are compared to SMOS SM and MODIS LST spaceborne observations. ET is obtained from the HidroMORE regional hydrological model. At the daily scale, a strong anticorrelation is observed between in situ SM and maximum LST (R ˜ -0.6 to -0.8), and between SMOS SM and MODIS LST Terra/Aqua day (R ˜ - 0.7). At the seasonal scale, results show a stronger anticorrelation in autumn, spring and summer (in situ R ˜ -0.5 to -0.7; satellite R ˜ -0.4 to -0.7) indicating SM–LST coupling, than in winter (in situ R ˜ +0.3; satellite R ˜ -0.3) indicating SM–LST decoupling. These different behaviors evidence changes from water-limited to energy-limited moisture flux across seasons, which are confirmed by the observed ET evolution. In water-limited periods, SM is extracted from the soil through ET until critical SM is reached. A method to estimate the soil critical SM is proposed. For REMEDHUS, the critical SM is estimated to be ~0.12 m3/m3 , stable over the study period and consistent between in situ and satellite observations. A better understanding of the SM–LST link could not only help improving the representation of LST in current hydrological and climate prediction models, but also refining SM retrieval or microwave-optical disaggregation algorithms, related to ET and vegetation status.Peer ReviewedPostprint (published version

Disaggregation of SMOS soil moisture over West Africa using the Temperature and Vegetation Dryness Index based on SEVIRI land surface parameters

The overarching objective of this study was to produce a disaggregated SMOS Soil Moisture (SM) product using land surface parameters from a geostationary satellite in a region covering a diverse range of ecosystem types. SEVIRI data at 15 minute temporal resolution were used to derive the Temperature and Vegetation Dryness Index (TVDI) that served as SM proxy within the disaggregation process. West Africa (3 N, 26 W; 28 N, 26 E) was selected as a case study as it presents both an important North-South climate gradient and a diverse range of ecosystem types. The main challenge was to set up a methodology applicable over a large area that overcomes the constraints of SMOS (low spatial resolution) and TVDI (requires similar atmospheric forcing and triangular shape formed when plotting morning rise temperature versus fraction of vegetation cover) in order to produce a 0.05 degree resolution disaggregated SMOS SM product at sub-continental scale. Consistent cloud cover appeared as one of the main constraints for deriving TVDI, especially during the rainy season and in the southern parts of the region and a large adjustment window (105x105 SEVIRI pixels) was therefore deemed necessary. Both the original and the disaggregated SMOS SM products described well the seasonal dynamics observed at six locations of in situ observations. However, there was an overestimation in both products for sites in the humid southern regions; most likely caused by the presence of forest. Both TVDI and the associated disaggregated SM product was found to be highly sensitive to algorithm input parameters; especially of conditions of high fraction of vegetation cover. Additionally, seasonal dynamics in TVDI did not follow the seasonal patters of SM. Still, its spatial heterogeneity was found to be a good proxy for disaggregating SMOS SM data; main river networks and spatial patterns of SM extremes (i.e. droughts and floods) not seen in the original SMOS SM product were revealed in the disaggregated SM product for a test case of July-September 2012. The disaggregation methodology thereby successfully increased the spatial resolution of SMOS SM, with potential application for local drought/flood monitoring of importance for the livelihood of the population of West Africa

A review of spatial downscaling of satellite remotely sensed soil moisture

Satellite remote sensing technology has been widely used to estimate surface soil moisture. Numerous efforts have been devoted to develop global soil moisture products. However, these global soil moisture products, normally retrieved from microwave remote sensing data, are typically not suitable for regional hydrological and agricultural applications such as irrigation management and flood predictions, due to their coarse spatial resolution. Therefore, various downscaling methods have been proposed to improve the coarse resolution soil moisture products. The purpose of this paper is to review existing methods for downscaling satellite remotely sensed soil moisture. These methods are assessed and compared in terms of their advantages and limitations. This review also provides the accuracy level of these methods based on published validation studies. In the final part, problems and future trends associated with these methods are analyzed

Disaggregation of SMOS soil moisture to 100m resolution using MODIS optical/thermal and sentinel-1 radar data: evaluation over a bare soil site in morocco

The 40 km resolution SMOS (Soil Moisture and Ocean Salinity) soil moisture, previously disaggregated at a 1 km resolution using the DISPATCH (DISaggregation based on Physical And Theoretical scale CHange) method based on MODIS optical/thermal data, is further disaggregated to 100 m resolution using Sentinel-1 backscattering coefficient (σ°). For this purpose, three distinct radar-based disaggregation methods are tested by linking the spatio-temporal variability of σ° and soil moisture data at the 1 km and 100 m resolution. The three methods are: (1) the weight method, which estimates soil moisture at 100 m resolution at a certain time as a function of σ° ratio (100 m to 1 km resolution) and the 1 km DISPATCH products of the same time; (2) the regression method which estimates soil moisture as a function of σ° where the regression parameters (e.g., intercept and slope) vary in space and time; and (3) the Cumulative Distribution Function (CDF) method, which estimates 100 m resolution soil moisture from the cumulative probability of 100 m resolution backscatter and the maximum to minimum 1 km resolution (DISPATCH) soil moisture difference. In each case, disaggregation results are evaluated against in situ measurements collected between 1 January 2016 and 11 October 2016 over a bare soil site in central Morocco. The determination coefficient (R2) between 1 km resolution DISPATCH and localized in situ soil moisture is 0.31. The regression and CDF methods have marginal effect on improving the DISPATCH accuracy at the station scale with a R2 between remotely sensed and in situ soil moisture of 0.29 and 0.34, respectively. By contrast, the weight method significantly improves the correlation between remotely sensed and in situ soil moisture with a R2 of 0.52. Likewise, the soil moisture estimates show low root mean square difference with in situ measurements (RMSD = 0.032 m3 m−3).This work is a contribution to the REC project funded by the European Commission Horizon 2020 Programme for Research and Innovation (H2020) in the context of the Marie Skłodowska-Curie Research and Innovation Staff Exchange (RISE) action under grant agreement no: 645642. In addition, this work has been partially funded by a public grant of Ministerio de Economía y Competitividad (DI-14-06587) and AGAUR-Generalitat de Catalunya (DI-2015-058)

Impact of day/night time land surface temperature in soil moisture disaggregation algorithms

18 pages, 5 figures, 1 tableSince its launch in 2009, the ESA’s SMOS mission is providing global soil moisture (SM) maps at ~40 km, using the first L-band microwave radiometer on space. Its spatial resolution meets the needs of global applications, but prevents the use of the data in regional or local applications, which require higher spatial resolutions (~1-10 km). SM disaggregation algorithms based generally on the land surface temperature (LST) and vegetation indices have been developed to bridge this gap. This study analyzes the SM-LST relationship at a variety of LST acquisition times and its influence on SM disaggregation algorithms. Two years of in situ and satellite data over the central part of the river Duero basin and the Iberian Peninsula are used. In situ results show a strong anticorrelation of SM to daily maximum LST (R≈0.5 to -0.8). This is confirmed with SMOS SM and MODIS LST Terra/Aqua at day time-overpasses (R≈-0.4 to -0.7). Better statistics are obtained when using MODIS LST day (R≈0.55 to 0.85; ubRMSD≈0.04 to 0.06 m/m) than LST night (R≈0.45 to 0.80; ubRMSD≈0.04 to 0.07 m/m) in the SM disaggregation. An averaged ensemble of day and night MODIS LST Terra/Aqua disaggregated SM estimates also leads to robust statistics (R≈0.55 to 0.85; ubRMSD≈0.04 to 0.07 m/m) with a coverage improvement of~10-20 %This work was supported by the Spanish Ministry of Economy and Competitiveness, through a Formación Personal Investigador (FPI) grant BES-2011-043322, the project PROMISES: ESP2015-67549-C3, ERDF (European Regional Development Fund) and the BBVA foundationPeer Reviewe

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