SMAP Detects Soil Moisture Under Temperate Forest Canopies

Soil moisture dynamics in the presence of dense vegetation canopies are determinants of ecosystem function and biogeochemical cycles, but the capability of existing spaceborne sensors to support reliable and useful estimates is not known. New results from a recently initiated field experiment in the northeast United States show that the National Aeronautics and Space Administration (NASA) SMAP (Soil Moisture Active Passive) satellite is capable of retrieving soil moisture under temperate forest canopies. We present an analysis demonstrating that a parameterized emission model with the SMAP morning overpass brightness temperature resulted in a RMSD (root‐mean‐square difference) range of 0.047–0.057 m3/m3 and a Pearson correlation range of 0.75–0.85 depending on the experiment location and the SMAP polarization. The inversion approach included a minimal amount of ancillary data. This result demonstrates unequivocally that spaceborne L‐band radiometry is sensitive to soil moisture under temperate forest canopies, which has been uncertain because of lack of representative reference data

Sensitivity of TDS-1 GNSS-R reflectivity to soil moisture: global and regional differences and impact of different spatial scales

The potential of Global Navigation Satellite Systems-Reflectometry (GNSS-R) techniques to estimate land surface parameters such as soil moisture (SM) is experimentally studied using 2014-2017 global data from the UK TechDemoSat-1 (TDS-1) mission. The approach is based on the analysis of the sensitivity to SM of different observables extracted from the Delay Doppler Maps (DDM) computed by the Space GNSS Receiver-Remote Sensing Instrument (SGR-ReSI) instrument using the L1 (1575.42 MHz) left-hand circularly-polarized (LHCP) reflected signals emitted by the Global Positioning System (GPS) navigation satellites. The sensitivity of different GNSS-R observables to SM and its dependence on the incidence angle is analyzed. It is found that the sensitivity of the calibrated GNSS-R reflectivity to surface soil moisture is ~ 0.09 dB/% up to 30° incidence angle, and it decreases with increasing incidence angles, although differences are found depending on the spatial scale used for the ground-truth, and the region. The sensitivity to subsurface soil moisture has been also analyzed using a network of subsurface probes and hydrological models, apparently showing some dependence, but so far results are not conclusive. © 2018 by the authors.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

Uncertainty quantification of a radiative transfer model and a machine learning technique for use as observation operators in the assimilation of microwave observations into a land surface model to improve soil moisture and terrestrial snow

Soil moisture and terrestrial snow mass are two important hydrological states needed to accurately quantify terrestrial water storage and streamflow. Soil moisture and terrestrial snow mass can be measured using ground-based instrument networks, estimated using advanced land surface models, and retrieved via satellite imagery. However, each method has its own inherent sources of error and uncertainty. This leads to the application of data assimilation to obtain optimal estimates of soil moisture and snow mass. Before conducting data assimilation (DA) experiments, this dissertation explored the use of two different observation operators within a DA framework: a L-band radiative transfer model (RTM) for soil moisture and support vector machine (SVM) regression for soil terrestrial snow mass. First, L-band brightness temperature (Tb) estimated from the RTM after being calibrated against multi-angular SMOS Tb's showed good performance in both ascending and descending overpasses across North America except in regions with sub-grid scale lakes and dense forest. Detailed analysis of RTM-derived L-band Tb in terms of soil hydraulic parameters and vegetation types suggests the need for further improvement of RTM-derived Tb in regions with relatively large porosity, large wilting point, or grassland type vegetation. Secondly, a SVM regression technique was developed with explicit consideration of the first-order physics of photon scattering as a function of different training target sets, training window lengths, and delineation of snow wetness over snow-covered terrain. The overall results revealed that prediction accuracy of the SVM was strongly linked with the first-order physics of electromagnetic responses of different snow conditions. After careful evaluation of the observation operators, C-band backscatter observations over Western Colorado collected by Sentinel-1 were merged into an advanced land surface model using a SVM and a one-dimensional ensemble Kalman filter. In general, updated snow mass estimates using the Sentinel-1 DA framework showed modest improvements in comparison to ground-based measurements of snow water equivalent (SWE) and snow depth. These results motivate further application of the outlined assimilation schemes over larger regions in order to improve the characterization of the terrestrial hydrological cycle

The International Soil Moisture Network:Serving Earth system science for over a decade

In 2009, the International Soil Moisture Network (ISMN) was initiated as a community effort, funded by the European Space Agency, to serve as a centralised data hosting facility for globally available in situ soil moisture measurements (Dorigo et al., 2011b, a). The ISMN brings together in situ soil moisture measurements collected and freely shared by a multitude of organisations, harmonises them in terms of units and sampling rates, applies advanced quality control, and stores them in a database. Users can freely retrieve the data from this database through an online web portal (https://ismn.earth/en/, last access: 28 October 2021). Meanwhile, the ISMN has evolved into the primary in situ soil moisture reference database worldwide, as evidenced by more than 3000 active users and over 1000 scientific publications referencing the data sets provided by the network. As of July 2021, the ISMN now contains the data of 71 networks and 2842 stations located all over the globe, with a time period spanning from 1952 to the present. The number of networks and stations covered by the ISMN is still growing, and approximately 70 % of the data sets contained in the database continue to be updated on a regular or irregular basis. The main scope of this paper is to inform readers about the evolution of the ISMN over the past decade, including a description of network and data set updates and quality control procedures. A comprehensive review of the existing literature making use of ISMN data is also provided in order to identify current limitations in functionality and data usage and to shape priorities for the next decade of operations of this unique community-based data repository

The Evaluation of the Additive Information Contained in New Ecohydrological Measurements

The development of ecohydrological frameworks and theories under the ongoing global climate crisis depends on the development of new and advanced ecohydrological measurements. Currently, numerous of datasets have been collected at plot and ecosystem levels to understand the complex interactions of along the soil, plant, and atmosphere continuum. The development of new measurements costs considerable effort, time, and funding. Therefore, it is important to quantify if new measurements contain useful information for ecohydrological studies, and how the information encoded in the new measurements are transferred to understand and predict key environmental variables. In this dissertation, we show that predictions in vegetation dynamics can be improved by adding observations from advanced satellite and in situ sensor systems. First, we evaluated new satellite soil moisture and vegetation optical depth measurements by Soil Moisture Passive Active (SMAP). We find that there is more opportunity for SMAP soil moisture and vegetation observations to be useful in locations where the daily vegetation climatology cannot adequately reflect observed vegetation dynamics. We also find that the uncertainties SMAP dual channel algorithm (DCA) are largely contributed by uncertainty of the algorithm’s inputs (horizontal and vertical polarized brightness temperature), while the informational uncertainty of the SMAP DCA model itself is more related to the retrieval quality of soil moisture. The informational redundancy and synergy of the two brightness temperature measurements from SMAP are tightly related to the SMAP soil moisture retrieval quality. Finally, we apply a similar informational analysis to new isotopic measurements at the National Ecological Observation Network (NEON). We find that majority of the information from the isotope measurements collected by NEON is unique, which cannot be obtained by other meteorological variables. Carbon isotope (δ13C) provides more additional information about LH in arid locations, while the water isotope (δ2H) provides more additional information about LH at locations with higher aridity, lower mean annual temperature, and lower mean site elevation. These studies show that informational analysis is useful to evaluated how additional information is encoded in new ecohydrological measurements

Advanced methods for earth observation data synergy for geophysical parameter retrieval

The first part of the thesis focuses on the analysis of relevant factors to estimate the response time between satellite-based and in-situ soil moisture (SM) using a Dynamic Time Warping (DTW). DTW was applied to the SMOS L4 SM, and was compared to in-situ root-zone SM in the REMEDHUS network in Western Spain. The method was customized to control the evolution of time lag during wetting and drying conditions. Climate factors in combination with crop growing seasons were studied to reveal SM-related processes. The heterogeneity of land use was analyzed using high-resolution images of NDVI from Sentinel-2 to provide information about the level of spatial representativity of SMOS data to each in-situ station. The comparison of long-term precipitation records and potential evapotranspiration allowed estimation of SM seasons describing different SM conditions depending on climate and soil properties. The second part of the thesis focuses on data-driven methods for sea ice segmentation and parameter retrieval. A Bayesian framework is employed to segment sets of multi-source satellite data. The Bayesian unsupervised learning algorithm allows to investigate the ‘hidden link’ between multiple data. The statistical properties are accounted for by a Gaussian Mixture Model, and the spatial interactions are reflected using Hidden Markov Random Fields. The algorithm segments spatial data into a number of classes, which are represented as a latent field in physical space and as clusters in feature space. In a first application, a two-step probabilistic approach based on Expectation-Maximization and the Bayesian segmentation algorithm was used to segment SAR images to discriminate surface water from sea ice types. Information on surface roughness is contained in the radar backscattering images which can be - in principle - used to detect melt ponds and to estimate high-resolution sea ice concentration (SIC). In a second study, the algorithm was applied to multi-incidence angle TB data from the SMOS L1C product to harness the its sensitivity to thin ice. The spatial patterns clearly discriminate well-determined areas of open water, old sea ice and a transition zone, which is sensitive to thin sea ice thickness (SIT) and SIC. In a third application, SMOS and the AMSR2 data are used to examine the joint effect of CIMR-like observations. The information contained in the low-frequency channels allows to reveal ranges of thin sea ice, and thicker ice can be determined from the relationship between the high-frequency channels and changing conditions as the sea ice ages. The proposed approach is suitable for merging large data sets and provides metrics for class analysis, and to make informed choices about integrating data from future missions into sea ice products. A regression neural network approach was investigated with the goal to infer SIT using TB data from the Flexible Microwave Payload 2 (FMPL-2) of the FSSCat mission. Two models - covering thin ice up to 0.6m and the full-range of SIT - were trained on Arctic data using ground truth data derived from the SMOS and Cryosat-2. This work demonstrates that moderate-cost CubeSat missions can provide valuable data for applications in Earth observation.La primera parte de la tesis se centra en el análisis de los factores relevantes para estimar el tiempo de respuesta entre la humedad del suelo (SM) basada en el satélite y la in-situ, utilizando una deformación temporal dinámica (DTW). El DTW se aplicó al SMOS L4 SM, y se comparó con la SM in-situ en la red REMEDHUS en el oeste de España. El método se adaptó para controlar la evolución del desfase temporal durante diferentes condiciones de humedad y secado. Se estudiaron los factores climáticos en combinación con los períodos de crecimiento de los cultivos para revelar los procesos relacionados con la SM. La heterogeneidad del uso del suelo se analizó utilizando imágenes de alta resolución de NDVI de Sentinel-2 para proporcionar información sobre el nivel de representatividad espacial de los datos de SMOS a cada estación in situ. La comparación de los patrones de precipitación a largo plazo y la evapotranspiración potencial permitió estimar las estaciones de SM que describen diferentes condiciones de SM en función del clima y las propiedades del suelo. La segunda parte de esta tesis se centra en métodos dirigidos por datos para la segmentación del hielo marino y la obtención de parámetros. Se emplea un método de inferencia bayesiano para segmentar conjuntos de datos satelitales de múltiples fuentes. El algoritmo de aprendizaje bayesiano no supervisado permite investigar el “vínculo oculto” entre múltiples datos. Las propiedades estadísticas se contabilizan mediante un modelo de mezcla gaussiana, y las interacciones espaciales se reflejan mediante campos aleatorios ocultos de Markov. El algoritmo segmenta los datos espaciales en una serie de clases, que se representan como un campo latente en el espacio físico y como clústeres en el espacio de las variables. En una primera aplicación, se utilizó un enfoque probabilístico de dos pasos basado en la maximización de expectativas y el algoritmo de segmentación bayesiano para segmentar imágenes SAR con el objetivo de discriminar el agua superficial de los tipos de hielo marino. La información sobre la rugosidad de la superficie está contenida en las imágenes de backscattering del radar, que puede utilizarse -en principio- para detectar estanques de deshielo y estimar la concentración de hielo marino (SIC) de alta resolución. En un segundo estudio, el algoritmo se aplicó a los datos TB de múltiples ángulos de incidencia del producto SMOS L1C para aprovechar su sensibilidad al hielo fino. Los patrones espaciales discriminan claramente áreas bien determinadas de aguas abiertas, hielo marino viejo y una zona de transición, que es sensible al espesor del hielo marino fino (SIT) y al SIC. En una tercera aplicación, se utilizan los datos de SMOS y de AMSR2 para examinar el efecto conjunto de las observaciones tipo CIMR. La información contenida en los canales de baja frecuencia permite revelar rangos de hielo marino delgado, y el hielo más grueso puede determinarse a partir de la relación entre los canales de alta frecuencia y las condiciones cambiantes a medida que el hielo marino envejece. El enfoque propuesto es adecuado para fusionar grandes conjuntos de datos y proporciona métricas para el análisis de clases, y para tomar decisiones informadas sobre la integración de datos de futuras misiones en los productos de hielo marino. Se investigó un enfoque de red neuronal de regresión con el objetivo de inferir el SIT utilizando datos de TB de la carga útil de microondas flexible 2 (FMPL-2) de la misión FSSCat. Se entrenaron dos modelos - que cubren el hielo fino hasta 0.6 m y el rango completo del SIT - con datos del Ártico utilizando datos de “ground truth” derivados del SMOS y del Cryosat-2. Este trabajo demuestra que las misiones CubeSat de coste moderado pueden proporcionar datos valiosos para aplicaciones de observación de la Tierra.Postprint (published version

Tropical Peatland Hydrology Simulated With a Global Land Surface Model

Tropical peatlands are among the most carbon-dense ecosystems on Earth, and their water storage dynamics strongly control these carbon stocks. The hydrological functioning of tropical peatlands differs from that of northern peatlands, which has not yet been accounted for in global land surface models (LSMs). Here, we integrated tropical peat-specific hydrology modules into a global LSM for the first time, by utilizing the peatland-specific model structure adaptation (PEATCLSM) of the NASA Catchment Land Surface Model (CLSM). We developed literature-based parameter sets for natural (PEATCLSM(Trop,Nat)) and drained (PEATCLSM(Trop,Drain)) tropical peatlands. Simulations with PEATCLSM(Trop,Nat) were compared against those with the default CLSM version and the northern version of PEATCLSM (PEATCLSM(North,Nat)) with tropical vegetation input. All simulations were forced with global meteorological reanalysis input data for the major tropical peatland regions in Central and South America, the Congo Basin, and Southeast Asia. The evaluation against a unique and extensive data set of in situ water level and eddy covariance-derived evapotranspiration showed an overall improvement in bias and correlation compared to the default CLSM version. Over Southeast Asia, an additional simulation with PEATCLSM(Trop,Drain) was run to address the large fraction of drained tropical peatlands in this region. PEATCLSM(Trop,Drain) outperformed CLSM, PEATCLSM(North,Nat), and PEATCLSM(Trop,Nat) over drained sites. Despite the overall improvements of PEATCLSM(Trop,Nat) over CLSM, there are strong differences in performance between the three study regions. We attribute these performance differences to regional differences in accuracy of meteorological forcing data, and differences in peatland hydrologic response that are not yet captured by our model.Peer reviewe

Hillslope-scale soil moisture estimation with a physically-based ecohydrology model and L-band microwave remote sensing observations from space

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Civil and Environmental Engineering, 2009.Includes bibliographical references (p. 469-488).Soil moisture is a critical hydrosphere state variable that links the global water, energy, and carbon cycles. Knowledge of soil moisture at scales of individual hillslopes (10's to 100's of meters) is critical to advancing applications such as landslide prediction, rainfall-runoff modeling, and wildland fire fuel load assessment. This thesis develops a data assimilation framework that employs the ensemble Kalman Filter (EnKF) to estimate the spatial distribution of soil moisture at hillslope scales by combining uncertain model estimates with noisy active and passive L-band microwave observations. Uncertainty in the modeled soil moisture state is estimated through Monte Carlo simulations with an existing spatially distributed ecohydrology model. Application of the EnKF to estimate hillslope-scale soil moisture in a watershed critically depends on: (1) identification of factors contributing to uncertainty in soil moisture, (2) adequate representation of the sources of uncertainty in soil moisture, and (3) formulation of an observing system to estimate the geophysically observable quantities based on the modeled soil moisture. Uncertainty in the modeled soil moisture distribution arises principally from uncertainty in the hydrometeorological forcings and imperfect knowledge of the soil parameters required as input to the model. Three stochastic models are used in combination to simulate uncertain hourly hydrometeorological forcings for the model. Soil parameter sets are generated using a stochastic approach that samples low probability but potentially high consequence parameter values and preserves correlation among the parameters. The observing system recognizes the role of the model in organizing the factors effecting emission and reflection of L-band microwave energy and emphasizes the role of topography in determining the satellite viewing geometry at hillslope scales.(cont.) Experiments in which true soil moisture conditions were simulated by the model and used to produce synthetic observations at spatial scales significantly coarser than the model resolution reveal that sequential assimilation of observations improves the hillslope-scale near-surface moisture estimate. Results suggest that the data assimilation framework is an effective means of disaggregating coarse-scale observations according to the model physics represented by the ecohydrology model. The thesis concludes with a discussion of contributions, implications, and future directions of this work.by Alejandro Nicolas Flores.Ph.D