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

Synergistic optical and microwave remote sensing approaches for soil moisture mapping at high resolution

Aplicat embargament des de la data de defensa fins al dia 1 d'octubre de 2022Soil moisture is an essential climate variable that plays a crucial role linking the Earth’s water, energy, and carbon cycles. It is responsible for the water exchange between the Earth’s surface and the atmosphere, and provides key information about soil evaporation, plant transpiration, and the allocation of precipitation into runoff, surface flow and infiltration. Therefore, an accurate estimation of soil moisture is needed to enhance our current climate and meteorological forecasting skills, and to improve our current understanding of the hydrological cycle and its extremes (e.g., droughts and floods). L-band Microwave passive and active sensors have been used during the last decades to estimate soil moisture, since there is a strong relationship between this variable and the soil dielectric properties. Currently, there are two operational L-band missions specifically devoted to globally measure soil moisture: the ESA’s Soil Moisture and the Ocean Salinity (SMOS), launched in November 2009; and the NASA’s Soil Moisture Active Passive (SMAP), launched in January 2015. The spatial resolution of the SMOS and SMAP radiometers, in the order of tens of kilometers (~40 km), is adequate for global applications. However, to fulfill the needs of a growing number of applications at local or regional scale, higher spatial detail (< 1 km) is required. To bridge this gap and improve the spatial resolution of the soil moisture maps, a variety of spatial enhancement or spatial (sub-pixel) disaggregation approaches have been proposed. This Ph.D. Thesis focuses on the study of the Earth’s surface soil moisture from remotely sensed observations. This work includes the implementation of several soil moisture retrieval techniques and the development, implementation, validation and comparison of different spatial enhancement or downscaling techniques, applied at local, regional, and continental scale. To meet these objectives, synergies between several active/passive microwave sensors (SMOS, SMAP and Sentinel-1) and optical/thermal sensors (MODIS) have been explored. The results are presented as follows: - Spatially consistent downscaling approach for SMOS using an adaptive moving window A passive microwave/optical downscaling algorithm for SMOS is proposed to obtain fine-scale soil moisture maps (1 km) from the native resolution (~40 km) of the instrument. This algorithm introduces the concept of a shape-adaptive window as a central improvement of the disaggregation technique presented by Piles et al. (2014), allowing its application at continental scales. - Assessment of multi-scale SMOS and SMAP soil moisture products across the Iberian Peninsula The temporal and spatial characteristics of SMOS and SMAP soil moisture products at coarse- and fine-scales are assessed in order to learn about their distinct features and the rationale behind them, tracing back to the physical assumptions they are based upon. - Impact of incidence angle diversity on soil moisture retrievals at coarse and fine scales An incidence angle (32.5°, 42.5° and 52.5°)-adaptive calibration of radiative transfer effective parameters single scattering albedo and soil roughness has been carried out, highlighting the importance of such parameterization to accurately estimate soil moisture at coarse-resolution. Then, these parameterizations are used to examine the potential application of a physically-based active-passive downscaling approach to upcoming microwave missions, namely CIMR, ROSE-L and Sentinel-1 Next Generation. Soil moisture maps obtained for the Iberian Peninsula at the three different angles, and at coarse and fine scales are inter-compared using in situ measurements and model data as benchmarks.La humedad del suelo es una variable climática esencial que juega un papel crucial en la relación de los ciclos del agua, la energía y el carbono de la Tierra. Es responsable del intercambio de agua entre la superficie de la Tierra y la atmósfera, y proporciona información crucial sobre la evaporación del suelo, la transpiración de las plantas y la distribución de la precipitación en escorrentía, flujo superficial e infiltración. Por lo tanto, es necesaria una estimación precisa de la humedad del suelo para mejorar las predicciones climáticas y meteorológicas, y comprender mejor el ciclo hidrológico y sus extremos (v.g., sequías e inundaciones). Los sensores pasivos y activos en banda L se han usado durante las últimas décadas para estimar la humedad del suelo debido a la relación directa que existe entre esta variable y las propiedades dieléctricas del suelo. Actualmente, hay dos misiones operativas en banda L específicamente dedicadas a medir la humedad del suelo a escala global: la misión Soil Moisture and Ocean Salinity (SMOS) de la ESA, lanzada en noviembre de 2009; y la misión Soil Moisture Active Passive (SMAP) de la NASA, lanzada en enero de 2015. La resolución espacial de los radiómetros SMOS y SMAP, del orden de unas decenas de kilómetros (~40 km), es adecuada para aplicaciones a escala global. Sin embargo, para satisfacer las necesidades de un número creciente de aplicaciones a escala local o regional, se requiere más detalle espacial (<1 km). Para solventar esta limitación y mejorar la resolución espacial de los mapas de humedad, se han propuesto diferentes técnicas de mejora o desagregación espacial. Esta Tesis se centra en el estudio de la humedad de la superficie terrestre a partir de datos obtenidos a través de teledetección. Este trabajo incluye la implementación de distintos algoritmos de recuperación de la humedad del suelo y el desarrollo, implementación, validación y comparación de distintas técnicas de desagregación, aplicadas a escala local, regional y continental. Para cumplir estos objetivos, se han explorado sinergias entre diferentes sensores de microondas activos/pasivos (SMOS, SMAP y Sentinel-1) y sensores ópticos/térmicos. Los resultados se presentan de la siguiente manera: - Técnica de desagregación espacialmente consistente, basada en una ventana móvil adaptativa, aplicada a los datos SMOS Se propone un algoritmo de desagregación del píxel basado en datos obtenidos de medidas radiométricas de microondas en banda L y datos ópticos, para mejorar la resolución espacial de los mapas de humedad del suelo desde la resolución nativa del instrumento (~40 km) hasta resoluciones de 1 km. El algoritmo introduce el concepto de una ventana de contorno adaptativo, como mejora principal sobre la técnica de desagregación presentada en Piles et al. (2014), permitiendo su implementación a escala continental. - Análisis multiescalar de productos de humedad del suelo SMAP y SMOS sobre la Península Ibérica Se han evaluado las características temporales y espaciales de distintos productos de humedad del suelo SMOS y SMAP, a baja y a alta resolución, para conocer sus características distintivas y comprender las razones de sus diferencias. Para ello, ha sido necesario rastrear los supuestos físicos en los que se basan. - Impacto del ángulo de incidencia en la recuperación de la humedad del suelo a baja y a alta resolución Se ha llevado a cabo una calibración adaptada al ángulo de incidencia (32.5°, 42.5° y 52.5°) de los parámetros efectivos, albedo de dispersión simple y rugosidad del suelo, descritos en el modelo de transferencia radiativa � − �, incidiendo en la importancia de esta parametrización para estimar la humedad del suelo de forma precisa a baja resolución. El resultado de las mismas se ha utilizado para estudiar la potencial aplicación de un algoritmo activo/pasivo de desagregación basado en la física para las próximas misiones de microondas, llamadas CIMR, ROSE-L y Sentinel-1 Next Generation. Los mapas de humedad recuperados a los tres ángulos de incidencia, tanto a baja como a alta resolución, se han obtenido para la Península Ibérica y se han comparado entre ellos usando como referencia mediciones de humedad in situ.Postprint (published version

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

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

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

An artificial neural network approach for soil moisture retrieval using passive microwave data

Soil moisture is a key variable that defines land surface-atmosphere (boundary layer) interactions, by contributing directly to the surface energy and water balance. Soil moisture values derived from remote sensing platforms only accounts for the near surface soil layers, generally the top 5cm. Passive microwave data at L-band (1.4 GHz, 21cm wavelength) measurements are shown to be a very effective observation for surface soil moisture retrieval. The first space-borne L-band mission dedicated to observing soil moisture, the European Space Agency's (ESA) Soil Moisture and Ocean Salinity (SMOS) mission, was launched on 2nd November 2009.Artificial Neural Network (ANN) methods have been used to empirically ascertain the complex statistical relationship between soil moisture and brightness temperature in the presence of vegetation cover. The current problem faced by this method is its inability to predict soil moisture values that are 'out-of-range' of the training data.In this research, an optimization model is developed for the Backpropagation Neural Network model. This optimization model utilizes the combination of the mean and standard deviation of the soil moisture values, together with the prediction process at different pre-determined, equal size regions to cope with the spatial and temporal variation of soil moisture values. This optimized model coupled with an ANN of optimum architecture, in terms of inputs and the number of neurons in the hidden layers, is developed to predict scale-to-scale and downscaling of soil moisture values. The dependency on the accuracy of the mean and standard deviation values of soil moisture data is also studied in this research by simulating the soil moisture values using a multiple regression model. This model obtains very encouraging results for these research problems.The data used to develop and evaluate the model in this research has been obtained from the National Airborne Field Experiments in 2005

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

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