Microwave Radiometry at Frequencies From 500 to 1400 MHz: An Emerging Technology for Earth Observations

icrowave radiometry has provided valuable spaceborne observations of Earth’s geophysical properties for decades. The recent SMOS, Aquarius, and SMAP satellites have demonstrated the value of measurements at 1400 MHz for observ- ing surface soil moisture, sea surface salinity, sea ice thickness, soil freeze/thaw state, and other geophysical variables. However, the information obtained is limited by penetration through the subsur- face at 1400 MHz and by a reduced sensitivity to surface salinity in cold or wind-roughened waters. Recent airborne experiments have shown the potential of brightness temperature measurements from 500–1400 MHz to address these limitations by enabling sensing of soil moisture and sea ice thickness to greater depths, sensing of temperature deep within ice sheets, improved sensing of sea salinity in cold waters, and enhanced sensitivity to soil moisture under veg- etation canopies. However, the absence of significant spectrum re- served for passive microwave measurements in the 500–1400 MHz band requires both an opportunistic sensing strategy and systems for reducing the impact of radio-frequency interference. Here, we summarize the potential advantages and applications of 500–1400 MHz microwave radiometry for Earth observation and review recent experiments and demonstrations of these concepts. We also describe the remaining questions and challenges to be addressed in advancing to future spaceborne operation of this technology along with recommendations for future research activities

Microwave emissions from snow

The radiation emitted from dry and wet snowpack in the microwave region (1 to 100 GHz) is discussed and related to ground observations. Results from theoretical model calculations match the brightness temperatures obtained by truck mounted, airborne and spaceborne microwave sensor systems. Snow wetness and internal layer structure complicate the snow parameter retrieval algorithm. Further understanding of electromagnetic interaction with snowpack may eventually provide a technique to probe the internal snow propertie

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

Atmospheric water parameters in mid-latitude cyclones observed by microwave radiometry and compared to model calculations

Existing and experimental algorithms for various parameters of atmospheric water content such as integrated water vapor, cloud water, precipitation, are used to examine the distribution of these quantities in mid latitude cyclones. The data was obtained from signals given by the special sensor microwave/imager (SSM/I) and compared with data from the nimbus scanning multichannel microwave radiometer (SMMR) for North Atlantic cyclones. The potential of microwave remote sensing for enhancing knowledge of the horizontal structure of these storms and to aid the development and testing of the cloud and precipitation aspects of limited area numerical models of cyclonic storms is investigated

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

Validation of SMOS sea ice thickness retrieval in the northern Baltic Sea

The Soil Moisture and Ocean Salinity (SMOS) mission observes brightness temperatures at a low microwave frequency of 1.4 GHz (L-band) with a daily coverage of the polar regions. L-band radiometry has been shown to provide information on the thickness of thin sea ice. Here, we apply a new emission model that has previously been used to investigate the impact of snow on thick Arctic sea ice. The model has not yet been used to retrieve ice thickness. In contrast to previous SMOS ice thickness retrievals, the new model allows us to include a snow layer in the brightness temperature simulations. Using ice thickness estimations from satellite thermal imagery, we simulate brightness temperatures during the ice growth season 2011 in the northern Baltic Sea. In both the simulations and the SMOS observations, brightness temperatures increase by more than 20 K, most likely due to an increase of ice thickness. Only if we include the snow in the model, the absolute values of the simulations and the observations agree well (mean deviations below 3.5 K). In a second comparison, we use high-resolution measurements of total ice thickness (sum of ice and snow thickness) from an electromagnetic (EM) sounding system to simulate brightness temperatures for 12 circular areas. While the SMOS observations and the simulations that use the EM modal ice thickness are highly correlated (r2=0.95), the simulated brightness temperatures are on average 12 K higher than observed by SMOS. This would correspond to an 8-cm overestimation of the modal ice thickness by the SMOS retrieval. In contrast, if the simulations take into account the shape of the EM ice thickness distributions (r2=0.87), the mean deviation between simulated and observed brightness temperatures is below 0.1 K

Towards long-term records of rain-on-snow events across the Arctic from satellite data

Rain-on-snow (ROS) events occur across many regions of the terrestrial Arctic in mid-winter. Snowpack properties are changing, and in extreme cases ice layers form which affect wildlife, vegetation and soils beyond the duration of the event. Specifically, satellite microwave observations have been shown to provide insight into known events. Only Ku-band radar (scatterometer) has been applied so far across the entire Arctic. Data availability at this frequency is limited, however. The utility of other frequencies from passive and active systems needs to be explored to develop a concept for long-term monitoring. The latter are of specific interest as they can be potentially provided at higher spatial resolution. Radar records have been shown to capture the associated snow structure change based on time-series analyses. This approach is also applicable when data gaps exist and has capabilities to evaluate the impact severity of events. Active as well as passive microwave sensors can also detect wet snow at the timing of an ROS event if an acquisition is available. The wet snow retrieval methodology is, however, rather mature compared to the identification of snow structure change since ambiguous scattering behaviour needs consideration. C-band radar is of special interest due to good data availability including a range of nominal spatial resolutions (10 m–12.5 km). Scatterometer and SAR (synthetic aperture radar) data have therefore been investigated. The temperature dependence of C-band backscatter at VV (V – vertical) polarization observable down to −40 ◦C is identified as a major issue for ROS retrieval but can be addressed by a combination with a passive microwave wet snow indicator (demonstrated for Metop ASCAT – Advanced Scatterometer – and SMOS – Soil Moisture and Ocean Salinity). Results were compared to in situ observations (snowpit records, caribou migration data) and Ku-band products. Ice crusts were found in the snowpack after detected events (overall accuracy 82 %). The more crusts (events) there are, the higher the winter season backscatter increase at C-band will be. ROS events captured on the Yamal and Seward peninsulas have had severe impacts on reindeer and caribou, respectively, due to ice crust formation. SAR specifically from Sentinel-1 is promising regarding ice layer identification at better spatial details for all available polarizations. The fusion of multiple types of microwave satellite observations is suggested for the creation of a climate data record, but the consideration of performance differences due to spatial and temporal cover, as well as microwave frequency, is crucial. Retrieval is most robust in the tundra biome, where results are comparable between sensors. Records can be used to identify extremes and to apply the results for impact studies at regional scale