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

Remote Sensing of Complex Permittivity and Penetration Depth of Soils Using P-Band SAR Polarimetry

A P-band SAR moisture estimation method is introduced for complex soil permittivity and penetration depth estimation using fully polarimetric P-band SAR signals. This method combines eigen- and model-based decomposition techniques for separation of the total backscattering signal into three scattering components (soil, dihedral, and volume). The incorporation of a soil scattering model allows for the first time the estimation of complex soil permittivity and permittivity-based penetration depth. The proposed method needs no prior assumptions on land cover characteristics and is applicable to a variety of vegetation types. The technique is demonstrated for airborne P-band SAR measurements acquired during the AirMOSS campaign (2012–2015). The estimated complex permittivity agrees well with climate and soil conditions at different monitoring sites. Based on frequency and permittivity, P-band penetration depths vary from 5 cm to 35 cm. This value range is in accordance with previous studies in the literature. Comparison of the results is challenging due to the sparsity of vertical soil in situ sampling. It was found that the disagreement between in situ measurements and SAR-based estimates originates from the discrepancy between the in situ measuring depth of the top-soil layer (0–5 cm) and the median penetration depth of the P-band waves (24.5–27 cm)

ALOS-2 L-band SAR backscatter data improves the estimation and temporal transferability of wildfire effects on soil properties under different post-fire vegetation responses

Remote sensing techniques are of particular interest for monitoring wildfire effects on soil properties, which may be highly context-dependent in large and heterogeneous burned landscapes. Despite the physical sense of synthetic aperture radar (SAR) backscatter data for characterizing soil spatial variability in burned areas, this approach remains completely unexplored. This study aimed to evaluate the performance of SAR backscatter data in C-band (Sentinel-1) and L-band (ALOS-2) for monitoring fire effects on soil organic carbon and nutrients (total nitrogen and available phosphorous) at short term in a heterogeneous Mediterranean landscape mosaic made of shrublands and forests that was affected by a large wildfire. The ability of SAR backscatter coefficients and several band transformations of both sensors for retrieving soil properties measured in the field in immediate post-fire situation (one month after fire) was tested through a model averaging approach. The temporal transferability of SAR-based models from one month to one year after wildfire was also evaluated, which allowed to assess short-term changes in soil properties at large scale as a function of pre-fire plant community type. The retrieval of soil properties in immediate post-fire conditions featured a higher overall fit and predictive capacity from ALOS-2 L-band SAR backscatter data than from Sentinel-1 C-band SAR data, with the absence of noticeable under and overestimation effects. The transferability of the ALOS-2 based model to one year after wildfire exhibited similar performance to that of the model calibration scenario (immediate post-fire conditions). Soil organic carbon and available phosphorous content was significantly higher one year after wildfire than immediately after the fire disturbance. Conversely, the short-term change in soil total nitrogen was ecosystem-dependent. Our results support the applicability of L-band SAR backscatter data for monitoring short-term variability of fire effects on soil properties, reducing data gathering costs within large and heterogeneous burned landscapesS

Soil moisture and water stage estimation using precipitation radar

In south-western United States, soil moisture data is important for drought studies in the region which is experiencing a drought for many years, whereas in South Florida, water stage data is required by hydrologists to monitor the hydrological flow in wetlands. Soil moisture data and water stage data are not sufficiently available due to sparse monitoring stations. Installation of dense measuring stations over an extended area is costly and labor intensive. Therefore, there is a need to develop an alternative method of measuring soil moisture and water stage. Microwave remote sensing has proven to be a useful tool in the measurement of various surface variables from space. This research explores the capability of microwave remote sensing to measure soil moisture and water stage on the earth from space. Tropical Rainfall Measuring Mission Precipitation Radar (TRMMPR) provides the Ku -band backscatter measurements that are used to measure soil moisture and water stage. Models that relate soil moisture and water stage to TRMMPR backscatter (σ°) are developed. The dependence of σ° on the dielectrical and physical characteristics of the land surface is used as the basis of this research. The soil moisture content affects σ° by changing the dielectric constant of the surface whereas the vegetation density affects σ° by changing the physical characteristics of the surface. Vegetation density in the model is represented by Normalized Difference Vegetation Index (NDVI). Dependence of σ° on partial submergence of vegetation in inundated areas is used to measure water stage in wetlands of South Florida. The effects of the exposed vegetation above the water surface on the model are assessed by comparing two cases of model run3 (a) that includes NDVI in the model, and (b) that does not include NDVI in the model. Eleven years of data is used in this research where 75% of the data is used for calibration of the model and 25% of the data is used for validation. The estimated values of soil moisture and water stage are compared to the observed values and the performance of the models is assessed by calculating correlation coefficients, calculating root mean square errors, and plotting non-exceedance probability plots for the absolute error between observed and modeled values. The soil moisture and water stage models work reasonably well and are able to estimate soil moisture and water stage with low errors. The soil moisture model works better in low vegetated areas because low vegetation allows the incident radiation to penetrate through the canopy cover and provide measurements from underlying surfaces. The water stage model works better in shrublands where there are no tree trunks and the model has an immediate impact from the vegetation canopy. This research provides an alternate way of measurement of soil moisture and water stage using remote sensing

Assimilation of SMOS Retrievals in the Land Information System

The Soil Moisture and Ocean Salinity (SMOS) satellite provides retrievals of soil moisture in the upper 5 cm with a 30-50 km resolution and a mission accuracy requirement of 0.04 cm(sub 3 cm(sub -3). These observations can be used to improve land surface model soil moisture states through data assimilation. In this paper, SMOS soil moisture retrievals are assimilated into the Noah land surface model via an Ensemble Kalman Filter within the NASA Land Information System. Bias correction is implemented using Cumulative Distribution Function (CDF) matching, with points aggregated by either land cover or soil type to reduce sampling error in generating the CDFs. An experiment was run for the warm season of 2011 to test SMOS data assimilation and to compare assimilation methods. Verification of soil moisture analyses in the 0-10 cm upper layer and root zone (0-1 m) was conducted using in situ measurements from several observing networks in the central and southeastern United States. This experiment showed that SMOS data assimilation significantly increased the anomaly correlation of Noah soil moisture with station measurements from 0.45 to 0.57 in the 0-10 cm layer. Time series at specific stations demonstrate the ability of SMOS DA to increase the dynamic range of soil moisture in a manner consistent with station measurements. Among the bias correction methods, the correction based on soil type performed best at bias reduction but also reduced correlations. The vegetation-based correction did not produce any significant differences compared to using a simple uniform correction curve

L-Band Vegetation optical depth and effective scattering albedo estimation from SMAP

Over land the vegetation canopy affects the microwave brightness temperature by emission, scattering and attenuation of surface soil emission. Attenuation, as represented by vegetation optical depth (VOD), is a potentially useful ecological indicator. The NASA Soil Moisture Active Passive (SMAP) mission carries significant potential for VOD estimates because of its radio frequency interference mitigation efforts and because the L-band signal penetrates deeper into the vegetation canopy than the higher frequency bands used for many previous VOD retrievals. In this study, we apply the multi-temporal dual-channel retrieval algorithm (MT-DCA) to derive global VOD, soil moisture, and effective scattering albedo estimates from SMAP Backus-Gilbert enhanced brightness temperatures posted on a 9 km grid and with three day revisit time. SMAP VOD values from the MT-DCA follow expected global distributions and are shown to be highly correlated with canopy height. They are also broadly similar in magnitude (though not always in seasonal amplitude) to European Space Agency Soil Moisture and Ocean Salinity (SMOS) VOD. The SMOS VOD values are based on angular brightness temperature information while the SMAP measurements are at a constant incidence angle, requiring an alternate approach to VOD retrieval presented in this study. Globally, albedo values tend to be high over regions with heterogeneous land cover types. The estimated effective scattering albedo values are generally higher than those used in previous soil moisture estimation algorithms and linked to biome classifications. MT-DCA retrievals of soil moisture show only small random differences with soil moisture retrievals from the Baseline SMAP algorithm, which uses a prior estimate of VOD based on land cover and optical data. However, significant biases exist between the two datasets. The soil moisture biases follow the pattern of differences between the MT-DCA retrieved and Baseline-assigned VOD values

Microwave Indices from Active and Passive Sensors for Remote Sensing Applications

Past research has comprehensively assessed the capabilities of satellite sensors operating at microwave frequencies, both active (SAR, scatterometers) and passive (radiometers), for the remote sensing of Earth’s surface. Besides brightness temperature and backscattering coefficient, microwave indices, defined as a combination of data collected at different frequencies and polarizations, revealed a good sensitivity to hydrological cycle parameters such as surface soil moisture, vegetation water content, and snow depth and its water equivalent. The differences between microwave backscattering and emission at more frequencies and polarizations have been well established in relation to these parameters, enabling operational retrieval algorithms based on microwave indices to be developed. This Special Issue aims at providing an overview of microwave signal capabilities in estimating the main land parameters of the hydrological cycle, e.g., soil moisture, vegetation water content, and snow water equivalent, on both local and global scales, with a particular focus on the applications of microwave indices

Downscaling GLDAS Soil Moisture Data in East Asia through Fusion of Multi-Sensors by Optimizing Modified Regression Trees

Soilmoisture is a key part of Earth&apos;s climate systems, including agricultural and hydrological cycles. Soil moisture data from satellite and numerical models is typically provided at a global scale with coarse spatial resolution, which is not enough for local and regional applications. In this study, a soil moisture downscaling model was developed using satellite-derived variables targeting Global Land Data Assimilation System (GLDAS) soil moisture as a reference dataset in East Asia based on the optimization of a modified regression tree. A total of six variables, Advanced Microwave Scanning Radiometer 2 (AMSR2) and Advanced SCATterometer (ASCAT) soil moisture products, Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM), and MODerate resolution Imaging Spectroradiometer (MODIS) products, including Land Surface Temperature, Normalized Difference Vegetation Index, and land cover, were used as input variables. The optimization was conducted through a pruning approach for operational use, and finally 59 rules were extracted based on root mean square errors (RMSEs) and correlation coefficients (r). The developed downscaling model showed a good modeling performance (r = 0.79, RMSE = 0.056 m(3)center dot m(3), and slope = 0.74). The 1 km downscaled soil moisture showed similar time series patterns with both GLDAS and ground soil moisture and good correlation with ground soil moisture (average r = 0.47, average RMSD = 0.038 m(3)center dot m(3)) at 14 ground stations. The spatial distribution of 1 km downscaled soil moisture reflected seasonal and regional characteristics well, although the model did not result in good performance over a few areas such as Southern China due to very high cloud cover rates. The results of this study are expected to be helpful in operational use to monitor soil moisture throughout East Asia since the downscaling model produces daily high resolution (1 km) real time soil moisture with a low computational demand. This study yielded a promising result to operationally produce daily high resolution soil moisture data from multiple satellite sources, although there are yet several limitations. In future research, more variables including Global Precipitation Measurement (GPM) precipitation, Soil Moisture Active Passive (SMAP) soil moisture, and other vegetation indices will be integrated to improve the performance of the proposed soil moisture downscaling model.ope

The SAR Handbook: Comprehensive Methodologies for Forest Monitoring and Biomass Estimation

This Synthetic Aperture Radar (SAR) handbook of applied methods for forest monitoring and biomass estimation has been developed by SERVIR in collaboration with SilvaCarbon to address pressing needs in the development of operational forest monitoring services. Despite the existence of SAR technology with all-weather capability for over 30 years, the applied use of this technology for operational purposes has proven difficult. This handbook seeks to provide understandable, easy-to-assimilate technical material to remote sensing specialists that may not have expertise on SAR but are interested in leveraging SAR technology in the forestry sector

Sensitivity of L-band vegetation optical depth to carbon stocks in tropical forests: a comparison to higher frequencies and optical indices

Supplementary data to this article can be found online at https://doi.org/10.1016/j.rse.2019.111303.Monitoring vegetation carbon in tropical regions is essential to the global carbon assessment and to evaluate the actions oriented to the reduction of forest degradation. Mainly, satellite optical vegetation indices and LiDAR data have been used to this purpose. These two techniques are limited by cloud cover and are sensitive only to the top of vegetation. In addition, the vegetation attenuation to the soil microwave emission, represented by the vegetation optical depth (VOD), has been applied for biomass estimation using frequencies ranging from 4 to 30¿GHz (C- to K-bands). Atmosphere is transparent to microwaves and their sensitivity to canopy layers depends on the frequency, with lower frequencies having greater penetration depths. In this regard, L-band VOD (1.4¿GHz) is expected to enhance the ability to estimate carbon stocks. This study compares the sensitivity of different VOD products (from L, C, and X-bands) and an optical vegetation index (EVI) to the above-ground carbon density (ACD). It quantifies the contribution of ACD and forest cover proportion to the VOD/EVI signals. The study is conducted in Peru, southern Colombia and Panama, where ACD maps have been derived from airborne LiDAR. Results confirm the enhanced sensitivity of L-band VOD to ACD when compared to higher frequency bands, and show that the sensitivity of all VOD bands decreases in the densest forests. ACD explains 34% and forest cover 30% of L-band VOD variance, and these proportions gradually decrease for EVI, C-, and X-band VOD, respectively. Results are consistent through different categories of altitude and carbon density. This pattern is found in most of the studied regions and in flooded forests. Results also show that C-, X-band VOD and EVI provide complementary information to L-band VOD, especially in flooded forests and in mountains, indicating that synergistic approaches could lead to improved retrievals in these regions. Although the assessment of vegetation carbon in the densest forests requires further research, results from this study support the use of new L-band VOD estimates for mapping the carbon of tropical forests.Peer ReviewedPostprint (author's final draft