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 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

Assessment of Multi-Scale SMOS and SMAP Soil Moisture Products across the Iberian Peninsula

In the last decade, technological advances led to the launch of two satellite missions dedicated to measure the Earth's surface soil moisture (SSM): the ESA's Soil Moisture and Ocean Salinity (SMOS) launched in 2009, and the NASA's Soil Moisture Active Passive (SMAP) launched in 2015. The two satellites have an L-band microwave radiometer on-board to measure the Earth's surface emission. These measurements (brightness temperatures TB) are then used to generate global maps of SSM every three days with a spatial resolution of about 30-40 km and a target accuracy of 0.04 m3/m3. To meet local applications needs, different approaches have been proposed to spatially disaggregate SMOS and SMAP TB or their SSM products. They rely on synergies between multi-sensor observations and are built upon different physical assumptions. In this study, temporal and spatial characteristics of six operational SSM products derived from SMOS and SMAP are assessed in order to diagnose their distinct features, and the rationale behind them. The study is focused on the Iberian Peninsula and covers the period from April 2015 to December 2017. A temporal inter-comparison analysis is carried out using in situ SSM data from the Soil Moisture Measurements Station Network of the University of Salamanca (REMEDHUS) to evaluate the impact of the spatial scale of the different products (1, 3, 9, 25, and 36 km), and their correspondence in terms of temporal dynamics. A spatial analysis is conducted for the whole Iberian Peninsula with emphasis on the added-value that the enhanced resolution products provide based on the microwave-optical (SMOS/ERA5/MODIS) or the active-passive microwave (SMAP/Sentinel-1) sensor fusion. Our results show overall agreement among time series of the products regardless their spatial scale when compared to in situ measurements. Still, higher spatial resolutions would be needed to capture local features such as small irrigated areas that are not dominant at the 1-km pixel scale. The degree to which spatial features are resolved by the enhanced resolution products depend on the multi-sensor synergies employed (at TB or soil moisture level), and on the nature of the fine-scale information used. The largest disparities between these products occur in forested areas, which may be related to the reduced sensitivity of high-resolution active microwave and optical data to soil properties under dense vegetation. Keywords: soil moisture; moisture variability; temporal dynamics; moisture patterns; spatial disaggregation; Soil Moisture Active Passive (SMAP); Soil Moisure and Ocean Salinity (SMOS); REMEDHUSSobre la continuidad de las misiones satelitales debanda L. Nuevos paradigmas en productos y aplicaciones, grant numbers ESP2017-89463-C3-2-R (UPC part) andESP2017-89463-C3-1-R (ICM part)Unidad de Excelencia María de Maeztu MDM-2016-060

Downscaling SMAP Soil Moisture Data Using MODIS Data

Soil moisture level is an important index in studying environmental changes. High resolution soil moisture data is in high demand for agricultural and weather forecasting purpose. Current daily large-scale soil moisture projects fail to provide sufficient resolution for medium or small region research. To acquire high-resolution soil moisture data, different kinds of methods are put into practice, including multivariate statistical regression, weight aggregation and so on. In this research, SMAP (Soil Moisture Active Passive) level 3 data with 36-km resolution are successfully downscaled by MODIS (Moderate Resolution Imaging Spectroradiometer) 1-km LST (Land Surface Temperature) product, NDVI (Difference Vegetation Index) product, SRTM (Shuttle Radar Topography Mission) DEM (Digital Elevation Model), and TWI (Topographic Wetness Index). Three regression models are built based on these supplemental indexes correlated with the SMAP retrieval. All downscaled results are validated with SMAPVEX15 field data. The research aims to establish and validate the multivariate regression method for downscaling low-resolution remote sensing image (such as SMAP) with local field observations. Based on the validation results, the research suggests the regression models have a decent fit. The downscaled soil moisture data indicating the method is applicable to small region research

The SMAP and Copernicus Sentinel 1A/B Microwave Active-Passive High Resolution Surface Soil Moisture Product and Its Applications

SMAP project released a new enhanced high-resolution (3km and 1 km) soil moisture active-passive product. This product is obtained by combining the SMAP radiometer data and the Sentinel-1A and -1B Synthetic Aperture Radar (SAR) data. The approach used for this product draws heavily from the heritage SMAP active-passive algorithm. Modifications in the SMAP active-passive algorithm are done to accommodate the Copernicus Program's Sentinel-1A and -1B multi-angular C-band SAR data. Assessment of the SMAP and Sentinel active-passive algorithm has been conducted and results show feasibility of estimating surface soil moisture at high-resolution in regions with low vegetation density (~< 3 kg/sq.m). A new version of this product is released to public in May 2018. This high resolution (3 km and 1 km) soil moisture product with reasonable accuracy of 0.05 m3/m3 is useful for agriculture, flood mapping, watershed/rangeland management, and ecological/hydrological applications

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

Long-term and high-resolution global time series of brightness temperature from copula-based fusion of SMAP enhanced and SMOS data

Long and consistent soil moisture time series at adequate spatial resolution are key to foster the application of soil moisture observations and remotely-sensed products in climate and numerical weather prediction models. The two L-band soil moisture satellite missions SMAP (Soil Moisture Active Passive) and SMOS (Soil Moisture and Ocean Salinity) are able to provide soil moisture estimates on global scales and in kilometer accuracy. However, the SMOS data record has an appropriate length of 7.5 years since late 2009, but with a coarse resolution of 25km only. In contrast, a spatially-enhanced SMAP product is available at a higher resolution of 9 km, but for a shorter time period (since March 2015 only). Being the fundamental observable from passive microwave sensors, reliable brightness temperatures (Tbs) are a mandatory precondition for satellite-based soil moisture products. We therefore develop, evaluate and apply a copula-based data fusion approach for combining SMAP Enhanced (SMAP_E) and SMOS brightness Temperature (Tb) data. The approach exploits both linear and non-linear dependencies between the two satellite-based Tb products and allows one to generate conditional SMAP_E-like random samples during the pre-SMAP period. Our resulting global Copula-combined SMOS-SMAP_E (CoSMOP) Tbs are statistically consistent with SMAP_E brightness temperatures, have a spatial resolution of 9km and cover the period from 2010 to 2018. A comparison with Service Soil Climate Analysis Network (SCAN)-sites over the Contiguous United States (CONUS) domain shows that the approach successfully reduces the average RMSE of the original SMOS data by 15%. At certain locations, improvements of 40% and more can be observed. Moreover, the median NSE can be enhanced from zero to almost 0.5. Hence, CoSMOP, which will be made freely available to the public, provides a first step towards a global, long-term, high-resolution and multi-sensor brightness temperature product, and thereby, also soil moisture