Thickness retrieval and emissivity modeling of thin sea ice at L-band for SMOS satellite observations

In this study we have developed an empirical retrieval for thickness of young and first-year ice during the freeze up period for the L-band passive microwave radiometer Microwave Imaging Radiometer with Aperture Synthesis (MIRAS) on the Soil Moisture and Ocean Salinity (SMOS) satellite. The retrieval is based on intensity and polarization difference using the incidence angle range of 40° to 50° and is validated using data from airborne EM-Bird, Moderate-resolution Imaging Spectroradiometer (MODIS) thermal imagery, and self consistency checks for ice thicknesses up to 50 cm with an error of 30 % on average. In addition, we modeled the microwave emission for Arctic first-year ice using the sea ice version of the Microwave Emission Model of Layered Snowpacks (MEMLS). The sea ice conditions used as input for MEMLS were generated using a thermodynamic energy balance model (based on the Crocus model) driven by reanalysis data from European Centre for Medium-Range Weather Forecasts (ECMWF). From unexpected features in the modeled microwave emission and disagreements with the empirically trained SMOS retrieval several shortcomings of the energy balance model and MEMLS were identified and corrected. The corrections include a treatment of mismatch of layer definition between the energy balance model and MEMLS, an adaptation of the reflection coefficient for lossy media in MEMLS, and several smaller corrections. For comparison, two simple models ignoring volume scattering, one incoherent and one coherent, were set up and were found to be able to reproduce the results of the more complex MEMLS model on average. With the simple models, the effects of thin coherent layers, the snow cover, the interface roughness and three different dielectric mixture models for sea ice were explored. It was found that the choice of the mixture model is essential for the relation of sea ice thickness to brightness temperatures in L-band, suggesting sea ice thickness sensitivities from few centimeters to several meters for salinity conditions of the global oceans. The interface properties, especially at the sea ice bottom, were found to be a major uncertainty source when modeling the microwave emission of thin sea ice. In addition, the variability in snow depth, the interface roughness, and the ice surface salinity and temperature were found to have a similar influence on the resulting brightness temperatures, with a strong effect on horizontally (up to 30 K) and weak effect on vertically polarized radiation (up to 10 K) for temperatures below 260 K. A model for simulating coherent microwave emission for thickness distributions of ice and snow was prepared to overcome weaknesses from the single thickness coherent and incoherent models. Comparison to the incoherent model showed that for realistic snow depth distributions obtained from Operation IceBridge (OIB) coherence effects can change the brightness temperatures on the scale of a SMOS footprint up to 10 K in horizontal polarization. These findings suggest that the retrieval for the thickness of thin sea ice with satellite based L-band sensors yield higher uncertainties than expected from earlier studies

Remote Sensing of Environmental Changes in Cold Regions

This Special Issue gathers papers reporting recent advances in the remote sensing of cold regions. It includes contributions presenting improvements in modeling microwave emissions from snow, assessment of satellite-based sea ice concentration products, satellite monitoring of ice jam and glacier lake outburst floods, satellite mapping of snow depth and soil freeze/thaw states, near-nadir interferometric imaging of surface water bodies, and remote sensing-based assessment of high arctic lake environment and vegetation recovery from wildfire disturbances in Alaska. A comprehensive review is presented to summarize the achievements, challenges, and opportunities of cold land remote sensing

Estimation de l'humidité du sol à haute résolution spatio-temporelle : une nouvelle approche basée sur la synergie des observations micro-ondes actives/passives et optiques/thermiques

Les capteurs micro-ondes passifs SMOS et SMAP fournissent des données d'humidité du sol (SM) à une résolution d'environ 40 km avec un intervalle de 2 à 3 jours à l' échelle mondiale et une profondeur de détection de 0 à 5 cm. Ces données sont très pertinentes pour les applications cli- matiques et météorologiques. Cependant, pour les applications à échelle régionales (l'hydrologie) ou locales (l'agriculture), des données de SM à une haute résolution spatiale (typiquement 100 m ou plus fine) seraient nécessaires. Les données collectées par les capteurs optiques/thermiques et les radars peuvent fournir des indicateurs de SM à haute résolution spatiale, mais ces deux approches alternatives ont des limites. En particulier, les données optiques/thermiques ne sont pas disponibles sous les nuages et sous les couverts végétaux. Quant aux données radar, elles sont sensibles à la rugosité du sol et à la structure de la végétation, qui sont tous deux difficiles à caractériser depuis l'espace. De plus, la résolution temporelle de ces données est d'environ 6 jours. Dans ce contexte, la ligne directrice de la thèse est de proposer une nouvelle approche qui combine pour la première fois des capteurs passifs micro-ondes, optiques/thermiques et actifs micro-ondes (radar) pour estimer SM sur de grandes étendues à une résolution de 100 m chaque jour. Notre hypothèse est d'abord de nous appuyer sur une méthode de désagrégation existante (DISPATCH) des données SMOS/SMAP pour atteindre la résolution cible obtenue par les radars. A l'origine, DISPATCH est basé sur l'efficacité d' évaporation du sol (SEE) estimée sur des pixels partiellement végétalisés à partir de données optiques/thermiques (généralement MODIS) de température de surface et de couverture végétale à résolution de 1 km. Les données désagrégées de SM sont ensuite combinées avec une méthode d'inversion de SM basée sur les données radar afin d'exploiter les capacités de détection des radars Sentinel-1. Enfin, les capacités de l'assimilation des donnés satellitaires de SM dans un modèle de bilan hydrique du sol sont évaluées en termes de prédiction de SM à une résolution de 100 m et à une échelle temporelle quotidienne.Dans une première étape, l'algorithme DISPATCH est amélioré par rapport à sa version actuelle, principalement 1) en étendant son applicabilité aux pixels optiques entièrement végétalisés en utilisant l'indice de sécheresse de la végétation basé sur la température et un produit de couverture végétale amélioré, et 2) en augmentant la résolution de désagrégation de 1 km à 100 m en utilisant les données optiques/thermiques de Landsat (en plus de MODIS). Le produit de SM désagrégé à la résolution de 100 m est validé avec des mesures in situ collectées sur des zones irriguées au Maroc, indiquant une corrélation spatiale quotidienne variant de 0,5 à 0,9. Dans un deuxième étape, un nouvel algorithme est construit en développant une synergie entre les données DISPATCH et radar à 100 m de résolution. En pratique, le produit SM issu de DISPATCH les jours de ciel clair est d'abord utilisé pour calibrer un modèle de transfert radiatif radar en mode direct. Ensuite, le modèle de transfert radiatif radar ainsi calibré est utilisé en mode inverse pour estimer SM à la résolution spatio-temporelle de Sentinel-1. Sur les sites de validation, les résultats indiquent une corrélation entre les mesures satellitaires et in situ, de l'ordre de 0,66 à 0,81 pour un indice de végétation inférieur à 0,6. Dans une troisième et dernière étape, une méthode d'assimilation optimale est utilisée pour interpoler dans le temps les données de SM à la résolution de 100 m. La dynamique du produit SM dérivé de l'assimilation de SM DISPATCH à 100 m de résolution est cohérente avec les événements d'irrigation. Cette approche peut être facilement appliquée sur de grandes zones, en considérant que toutes les données (télédétection et météorologique) requises en entrée sont disponibles à l' échelle globale.SMOS and SMAP passive microwave sensors provide soil moisture (SM) data at 40 km resolution every 2-3 days globally, with a 0-5 cm sensing depth relevant for climatic and meteorological applications. However, SM data would be required at a higher (typically 100 m or finer) spatial resolution for many other regional (hydrology) or local (agriculture) applications. Optical/thermal and radar sensors can be used for retrieving SM proxies at such high spatial resolution, but both techniques have limitations. In particular, optical/thermal data are not available under clouds and under plant canopies. Moreover, radar data are sensitive to soil roughness and vegetation structure, which are challenging to characterize from outer space, and have a repeat cycle of at least six days, limiting the observations' temporal frequency. In this context, the leading principle of the thesis is to propose a new approach that combines passive microwave, optical/thermal, and active microwave (radar) sensors for the first time to retrieve SM data at 100 m resolution on a daily temporal scale. Our assumption is first to rely on an existing disaggregation method (DISPATCH) of SMOS/SMAP SM data to meet the target resolution achieved by radars. DISPATCH is originally based on the soil evaporative efficiency (SEE) retrieved over partially vegetated pixels from 1 km resolution optical/thermal (typically MODIS) surface temperature and vegetation cover data. The disaggregated SM data is then combined with a radar-based SM retrieval method to exploit the sensing capabilities of the Sentinel-1 radars. Finally, the efficacy of the assimilation of satellite-based SM data in a soil water balance model is assessed in terms of SM predictions at the 100 m resolution and daily temporal scale. As a first step, the DISPATCH algorithm is improved from its current version by mainly 1) extending its applicability to fully vegetated optical pixels using the temperature vegetation dryness index and an enhanced vegetation cover product, and 2) increasing the targeted downscaling resolution from 1 km to 100 m using Landsat (in addition to MODIS) optical/thermal data. The 100 m resolution disaggregated SM product is validated with in situ measurements collected over irrigated areas in Morocco, showing a daily spatial correlation in the range of 0.5-0.9. As a second step, a new algorithm is built on a synergy between DISPATCH and radar 100 m resolution data. In practice, the DISPATCH SM product available on clear sky days is first used to calibrate a radar radiative transfer model in the direct mode. Then the calibrated radar radia- tive transfer model is used in the inverse mode to estimate SM at the spatio-temporal resolution of Sentinel-1. Results indicate a positive correlation between satellite and in situ measurements in the range of 0.66 to 0.81 for a vegetation index lower than 0.6. As a third and final step, an optimal assimilation method is used to interpolate 100 m resolution SM data in time. The assimilation exercise is undertaken over irrigated crop fields in Spain. The analyzed SM product derived from the assimilation of 100 m resolution DISPATCH SM is consistent with irrigation events. This approach can be readily applied over large areas, given that all the required input (remote sensing and meteorological) data are available globally

Energy and Water Cycles in the Third Pole

As the most prominent and complicated terrain on the globe, the Tibetan Plateau (TP) is often called the “Roof of the World”, “Third Pole” or “Asian Water Tower”. The energy and water cycles in the Third Pole have great impacts on the atmospheric circulation, Asian monsoon system and global climate change. On the other hand, the TP and the surrounding higher elevation area are also experiencing evident and rapid environmental changes under the background of global warming. As the headwater area of major rivers in Asia, the TP’s environmental changes—such as glacial retreat, snow melting, lake expanding and permafrost degradation—pose potential long-term threats to water resources of the local and surrounding regions. To promote quantitative understanding of energy and water cycles of the TP, several field campaigns, including GAME/Tibet, CAMP/Tibet and TORP, have been carried out. A large amount of data have been collected to gain a better understanding of the atmospheric boundary layer structure, turbulent heat fluxes and their coupling with atmospheric circulation and hydrological processes. The focus of this reprint is to present recent advances in quantifying land–atmosphere interactions, the water cycle and its components, energy balance components, climate change and hydrological feedbacks by in situ measurements, remote sensing or numerical modelling approaches in the “Third Pole” region

Measurements and modelling of seasonal snow characteristics for interpreting passive microwave observations

Information on snow water equivalent (SWE) of seasonal snow is used for various purposes, including longterm climate monitoring and river discharge forecasting. Global monitoring of SWE is made feasible through remote sensing. Currently, passive microwave observations are utilized for SWE retrievals. The main challenges in the interpretation of microwave observations include the spatial variability of snow characteristics and the inaccurate characterization of snow microstructure in retrieval algorithms. Even a minor variability in snow microstructure has a notable impact to microwave emission from the snowpack. This thesis work aims to improve snow microstructure modelling and measurement methods, and understanding the influence of snow microstructure to passive microwave observations, in order to enable a more accurate SWE estimation from remote sensing observations. The thesis work applies two types of models: physical snow models and radiative transfer models that simulate microwave emission. The physical snow models use meteorological driving data to simulate physical snow characteristics, such as SWE and snow microstructure. Models are used for different purposes such as hydrological simulations and avalanche forecasting. On the other hand, microwave emission models use physical snow characteristics for predicting microwave emission from a snowpack. Microwave emission models are applied for the interpretation of spaceborne passive microwave remote sensing observations, for example. In this study, physical snow model simulations and microwave emission model simulations are compared with field observations to investigate problems in characterizing snow for microwave emission models. An extensive set of manual field measurements of snow characteristics is used for the comparisons. The measurements are collected from taiga snow in Sodankylä, northern Finland. The representativeness of the measurements is defined by investigating the spatial and temporal variability of snow characteristics. The work includes studies on microwave emission modelling from natural snowpacks and from excavated snow slabs. Radiometric observations of microwave emission from natural snowpacks are compared with simulations from three microwave emission models coupled with three physical snow models. Additionally, homogenous snow samples are excavated from the natural snowpack during the Arctic Snow Microstructure Experiment, and the incident snow characteristics and microwave emission characteristics are measured with an experimental set-up developed for this study. Predictions of two microwave emission models are compared with the radiometric observations of collected snow samples. The results indicate that none of the model configurations can accurately simulate the microwave emission from natural snowpack or snow samples. The results also suggest that the characterization of microstructure in the applied microwave emission models is not adequate