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

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

Désagrégation de l'humidité du sol issue des produits satellitaires micro-ondes passives et exploration de son utilisation pour l'amélioration de la modélisation et la prévision hydrologique

De plus en plus de produits satellitaires en micro-ondes passives sont disponibles. Cependant, leur large résolution spatiale (25-50 km) n’en font pas un outil adéquat pour des applications hydrologiques à une échelle locale telles que la modélisation et la prévision hydrologiques. Dans de nombreuses études, une désagrégation d’échelle de l’humidité du sol des produits satellites micro-ondes est faite puis validée avec des mesures in-situ. Toutefois, l’utilisation de ces données issues d’une désagrégation d’échelle n’a pas encore été pleinement étudiée pour des applications en hydrologie. Ainsi, l’objectif de cette thèse est de proposer une méthode de désagrégation d’échelle de l’humidité du sol issue de données satellitaires en micro-ondes passives (Satellite Passive Microwave Active and Passive - SMAP) à différentes résolutions spatiales afin d’évaluer leur apport sur l’amélioration potentielle des modélisations et prévisions hydrologiques. À partir d’un modèle de forêt aléatoire, une désagrégation d’échelle de l’humidité du sol de SMAP l’amène de 36-km de résolution initialement à des produits finaux à 9-, 3- et 1-km de résolution. Les prédicteurs utilisés sont à haute résolution spatiale et de sources différentes telles que Sentinel-1A, MODIS et SRTM. L'humidité du sol issue de cette désagrégation d’échelle est ensuite assimilée dans un modèle hydrologique distribué à base physique pour tenter d’améliorer les sorties de débit. Ces expériences sont menées sur les bassins versants des rivières Susquehanna (de grande taille) et Upper-Susquehanna (en comparaison de petite taille), tous deux situés aux États-Unis. De plus, le modèle assimile aussi des données d’humidité du sol en profondeur issue d’une extrapolation verticale des données SMAP. Par ailleurs, les données d’humidité du sol SMAP et les mesures in-situ sont combinées par la technique de fusion conditionnelle. Ce produit de fusion SMAP/in-situ est assimilé dans le modèle hydrologique pour tenter d’améliorer la prévision hydrologique sur le bassin versant Au Saumon situé au Québec. Les résultats montrent que l'utilisation de l’humidité du sol à fine résolution spatiale issue de la désagrégation d’échelle améliore la représentation de la variabilité spatiale de l’humidité du sol. En effet, le produit à 1- km de résolution fournit plus de détails que les produits à 3- et 9-km ou que le produit SMAP de base à 36-km de résolution. De même, l’utilisation du produit de fusion SMAP/ in-situ améliore la qualité et la représentation spatiale de l’humidité du sol. Sur le bassin versant Susquehanna, la modélisation hydrologique s’améliore avec l’assimilation du produit de désagrégation d’échelle à 9-km, sans avoir recours à des résolutions plus fines. En revanche, sur le bassin versant Upper-Susquehanna, c’est le produit avec la résolution spatiale la plus fine à 1- km qui offre les meilleurs résultats de modélisation hydrologique. L’assimilation de l’humidité du sol en profondeur issue de l’extrapolation verticale des données SMAP n’améliore que peu la qualité du modèle hydrologique. Par contre, l’assimilation du produit de fusion SMAP/in-situ sur le bassin versant Au Saumon améliore la qualité de la prévision du débit, même si celle-ci n’est pas très significative.Abstract: The availability of satellite passive microwave soil moisture is increasing, yet its spatial resolution (i.e., 25-50 km) is too coarse to use for local scale hydrological applications such as streamflow simulation and forecasting. Many studies have attempted to downscale satellite passive microwave soil moisture products for their validation with in-situ soil moisture measurements. However, their use for hydrological applications has not yet been fully explored. Thus, the objective of this thesis is to downscale the satellite passive microwave soil moisture (i.e., Satellite Microwave Active and Passive - SMAP) to a range of spatial resolutions and explore its value in improving streamflow simulation and forecasting. The random forest machine learning technique was used to downscale the SMAP soil moisture from 36-km to 9-, 3- and 1-km spatial resolutions. A combination of host of high-resolution predictors derived from different sources including Sentinel-1A, MODIS and SRTM were used for downscaling. The downscaled SMAP soil moisture was then assimilated into a physically-based distributed hydrological model for improving streamflow simulation for Susquehanna (larger in size) and Upper Susquehanna (relatively smaller in size) watersheds, located in the United States. In addition, the vertically extrapolated SMAP soil moisture was assimilated into the model. On the other hand, the SMAP and in-situ soil moisture were merged using the conditional merging technique and the merged SMAP/in-situ soil moisture was then assimilated into the model to improve streamflow forecast over the au Saumon watershed. The results show that the downscaling improved the spatial variability of soil moisture. Indeed, the 1-km downscaled SMAP soil moisture presented a higher spatial detail of soil moisture than the 3-, 9- or original resolution (36-km) SMAP product. Similarly, the merging of SMAP and in-situ soil moisture improved the accuracy as well as spatial representation soil moisture. Interestingly, the assimilation of the 9-km downscaled SMAP soil moisture significantly improved the accuracy of streamflow simulation for the Susquehanna watershed without the need of going to higher spatial resolution, whereas for the Upper Susquehanna watershed the 1-km downscaled SMAP showed better results than the coarser resolutions. The assimilation of vertically extrapolated SMAP soil moisture only slightly further improved the accuracy of the streamflow simulation. On the other hand, the assimilation of merged SMAP/in-situ soil moisture for the au Saumon watershed improved the accuracy of streamflow forecast, yet the improvement was not that significant. Overall, this study demonstrated the potential of satellite passive microwave soil moisture for streamflow simulation and forecasting

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

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

Mapping Soil Moisture from Remotely Sensed and In-situ Data with Statistical Methods

Soil moisture is an important factor for accurate prediction of agricultural productivity and rainfall runoff with hydrological models. Remote sensing satellites such as Soil Moisture Active Passive (SMAP) offer synoptic views of soil moisture distribution at a regional-to-global scale. To use the soil moisture product from these satellites, however, requires a downscaling of the data from an usually large instantaneous field of view (i.e. 36 km) to the watershed analysis scales ranging from 30 m to 1 km. In addition, validation of the soil moisture products using the ground station observations without an upscaling treatment would lead to cross-level fallacy. In the literature of geographical analysis, scale is one of the top research concens because of the needs for multi-source geospatial data fusion. This dissertation research introduced a multi-level soil moisture data assimilation and processing methodology framework based on spatial information theories. The research contains three sections: downscaling using machine learning and geographically weighted regression, upscaling ground network observation to calibrate satellite data, and spatial and temporal multi-scale data assimilation using spatio-temporal interpolation. (1) Soil moisture downscaling In the first section, a downscaling method is designed using 1-km geospatial data to obtain subpixel soil moisture from the 9-km soil moisture product of the SMAP satellite. The geospatial data includes normalized difference vegetation index (NDVI), land surface temperature (LST), gross primary productivity (GPP), topographical moisture index (TMI), with all resampled to 1-km resolution. The machine learning algorithm – random forest was used to create a prediction model of the soil moisture at a 1-km resolution. The 1-km soil moisture product was compared with the ground samples from the West Texas Mesonet (WTM) station data. The residual was then interpolated to compensate the unpredicted variability of the model. The entire process was based on the concept of regression kriging- where the regression was done by the random forest model. Results show that the downscaling approach was able to achieve better accuracy than the current statistical downscaling methods. (2) Station network data upscaling The Texas Soil Observation Network (TxSON) network was designed to test the feasibility of upscaling the in-situ data to match the scale of the SMAP data. I advanced the upscaling method by using the Voronoi polygons and block kriging with a Gaussian kernel aggregation. The upscaling algorithm was calibrated using different spatial aggregation parameters, such as the fishnet cell size and Gaussian kernel standard deviation. The use of the kriging can significantly reduce the spatial autocorrelation among the TxSON stations because of its declustering ability. The result proved the new upscaling method was better than the traditional ones. (3) Multi-scale data fusion in a spatio-temporal framework None of the current works for soil moisture statistical downscaling honors time and space equally. It is important, however, that the soil moisture products are consistent in both domains. In this section, the space-time kriging model for soil moisture downscaling and upscaling computation framework designed in the last two sections is implemented to create a spatio-temporal integrated solution to soil moisture multi-scale mapping. The present work has its novelty in using spatial statistics to reconcile the scale difference from satellite data and ground observations, and therefore proposes new theories and solutions for dealing with the modifiable areal unit problem (MAUP) incurred in soil moisture mapping from satellite and ground stations