Snow stratigraphic heterogeneity within ground-based passive microwave radiometer footprints: implications for emission modeling

Two-dimensional measurements of snowpack properties (stratigraphic layering, density, grain size and temperature) were used as inputs to the multi-layer Helsinki University of Technology (HUT) microwave emission model at a centimeter-scale horizontal resolution, across a 4.5 m transect of ground-based passive microwave radiometer footprints near Churchill, Manitoba, Canada. Snowpack stratigraphy was complex (between six and eight layers) with only three layers extending continuously throughout the length of the transect. Distributions of one-dimensional simulations, accurately representing complex stratigraphic layering, were evaluated using measured brightness temperatures. Large biases (36 to 68 K) between simulated and measured brightness temperatures were minimized (-0.5 to 0.6 K), within measurement accuracy, through application of grain scaling factors (2.6 to 5.3) at different combinations of frequencies, polarizations and model extinction coefficients. Grain scaling factors compensated for uncertainty relating optical SSA to HUT effective grain size inputs and quantified relative differences in scattering and absorption properties of various extinction coefficients. The HUT model required accurate representation of ice lenses, particularly at horizontal polarization, and large grain scaling factors highlighted the need to consider microstructure beyond the size of individual grains. As variability of extinction coefficients was strongly influenced by the proportion of large (hoar) grains in a vertical profile, it is important to consider simulations from distributions of one-dimensional profiles rather than single profiles, especially in sub-Arctic snowpacks where stratigraphic variability can be high. Model sensitivity experiments suggested the level of error in field measurements and the new methodological framework used to apply them in a snow emission model were satisfactory. Layer amalgamation showed a three-layer representation of snowpack stratigraphy reduced the bias of a one-layer representation by about 50%

The influence of winter time boreal forest tree transmissivity on tree emission and passive microwave snow observations

Forest cover significantly attenuates natural upwelling ground microwave emission from seasonal terrestrial snow. This presents a major challenge for the accurate retrieval of snow from airborne or spaceborne passive microwave (PM) observations. Forest transmissivity is a key parameter describing tree emission because not only does it influence the proportion of sub-canopy upwelling microwave emission penetrating through the forest canopy, it also controls the forest thermal emission. Hence, it is a very important parameter for correcting the influence of forests on spaceborne or airborne observations of the Earth’s land surface. Under sub-zero temperatures, vegetation water content can be frozen influencing the microwave transmissivity of trees. Yet this phenomenon has not been verified through experimentation leaving significant uncertainty in tree emission modelling and spaceborne microwave observations. Therefore, a season-long experiment was designed to study this phenomenon. Ground-based radiometer observations of tree emission, spaceborne observations of forest emission, and model simulations of canopy emission were conducted during this experiment. Based on this experiment, the influence of physical temperature on tree transmissivity was verified, and a model developed to quantitatively describe this temperature-transmissivity relationship. An evaluation of this temperature-transmissivity relationship was conducted showing that both ground-based and spaceborne observations of tree emission are significantly influenced by this phenomenon. Furthermore, passive microwave spaceborne snow retrievals in forested regions are influenced by this phenomenon. Finally, an approach to reduce the influence of the temperature-transmissivity relationship on passive microwave spaceborne snow retrievals is demonstrated

New geophysics study programs at the University of Helsinki

Peer reviewe

Petrophysical properties of the Kylylahti Cu-Au-Zn sulphide mineralization and its host rocks

Non peer reviewe

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

The future of Earth observation in hydrology

In just the past 5 years, the field of Earth observation has progressed beyond the offerings of conventional space-agency-based platforms to include a plethora of sensing opportunities afforded by CubeSats, unmanned aerial vehicles (UAVs), and smartphone technologies that are being embraced by both for-profit companies and individual researchers. Over the previous decades, space agency efforts have brought forth well-known and immensely useful satellites such as the Landsat series and the Gravity Research and Climate Experiment (GRACE) system, with costs typically of the order of 1 billion dollars per satellite and with concept-to-launch timelines of the order of 2 decades (for new missions). More recently, the proliferation of smart-phones has helped to miniaturize sensors and energy requirements, facilitating advances in the use of CubeSats that can be launched by the dozens, while providing ultra-high (3-5 m) resolution sensing of the Earth on a daily basis. Start-up companies that did not exist a decade ago now operate more satellites in orbit than any space agency, and at costs that are a mere fraction of traditional satellite missions. With these advances come new space-borne measurements, such as real-time high-definition video for tracking air pollution, storm-cell development, flood propagation, precipitation monitoring, or even for constructing digital surfaces using structure-from-motion techniques. Closer to the surface, measurements from small unmanned drones and tethered balloons have mapped snow depths, floods, and estimated evaporation at sub-metre resolutions, pushing back on spatio-temporal constraints and delivering new process insights. At ground level, precipitation has been measured using signal attenuation between antennae mounted on cell phone towers, while the proliferation of mobile devices has enabled citizen scientists to catalogue photos of environmental conditions, estimate daily average temperatures from battery state, and sense other hydrologically important variables such as channel depths using commercially available wireless devices. Global internet access is being pursued via high-altitude balloons, solar planes, and hundreds of planned satellite launches, providing a means to exploit the "internet of things" as an entirely new measurement domain. Such global access will enable real-time collection of data from billions of smartphones or from remote research platforms. This future will produce petabytes of data that can only be accessed via cloud storage and will require new analytical approaches to interpret. The extent to which today's hydrologic models can usefully ingest such massive data volumes is unclear. Nor is it clear whether this deluge of data will be usefully exploited, either because the measurements are superfluous, inconsistent, not accurate enough, or simply because we lack the capacity to process and analyse them. What is apparent is that the tools and techniques afforded by this array of novel and game-changing sensing platforms present our community with a unique opportunity to develop new insights that advance fundamental aspects of the hydrological sciences. To accomplish this will require more than just an application of the technology: in some cases, it will demand a radical rethink on how we utilize and exploit these new observing systems

Analyse de la modélisation de l'émission multi-fréquences micro-onde des sols et de la neige, incluant les croutes de glace à l'aide du modèle Microwave Emission Model of Layered Snowpacks (MEMLS).

Résumé : L'étude du couvert nival est essentielle afin de mieux comprendre les processus climatiques et hydrologiques. De plus, avec les changements climatiques observés dans l'hémisphère nord, des événements de dégel-regel ou de pluie hivernale sont de plus en plus courants et produisent des croutes de glace dans le couvert nival affectant les moeurs des communautés arctiques en plus de menacer la survie de la faune arctique. La télédétection micro-ondes passives (MOP) démontre un grand potentiel de caractérisation du couvert nival. Toutefois, a fin de bien comprendre les mesures satellitaires, une modélisation adéquate du signal est nécessaire. L'objectif principal de cette thèse est d'analyser le transfert radiatif (TR) MOP des sols, de la neige et de la glace a fin de mieux caractériser les propriétés géophysiques du couvert nival par télédétection. De plus, un indice de détection des croutes de glace par télédétection MOP a été développé. Pour ce faire, le modèle Microwave Emission Model of Layered Snowpacks (MEMLS) a été étudié et calibré afin de minimiser les erreurs des températures de brillance simulées en présences de croutes de glace. La première amélioration faite à la modélisation du TR MOP de la neige a été la caractérisation de la taille des grains de neige. Deux nouveaux instruments, utilisant la réflectance dans le proche infrarouge, ont été développés afin de mesurer la surface spécifique de la neige (SSA). Il a été démontré que la SSA est un paramètre plus précis et plus objectif pour caractériser la taille des grains de neige. Les deux instruments ont démontré une incertitude de 10% sur la mesure de la SSA. De plus, la SSA a été calibré pour la modélisation MOP a n de minimiser l'erreur sur la modélisation de la température de brillance. Il a été démontré qu'un facteur multiplicatif [phi] = 1.3 appliqué au paramètre de taille des grains de neige dans MEMLS, paramètre dérivé de la SSA, est nécessaire afin de minimiser l'erreur des simulations. La deuxième amélioration apportée à la modélisation du TR MOP a été l'estimation de l'émission du sol. Des mesures radiométriques MOP in-situ ainsi que des profils de températures de sols organiques arctiques gelés ont été acquis et caractérisés a fin de simuler l'émission MOP de ces sols. Des constantes diélectriques effectives à 10.7, 19 et 37 GHz ainsi qu'une rugosité de surface effective des sols ont été déterminés pour simuler l'émission des sols. Une erreur quadratique moyenne (RMSE) de 4.65 K entre les simulations et les mesures MOP a été obtenue. Suite à la calibration du TR MOP du sol et de la neige, un module de TR de la glace a été implémenté dans MEMLS. Avec ce nouveau module, il a été possible de démontré que l'approximation de Born améliorée, déjà implémenté dans MEMLS, pouvait être utilisé pour simuler des croutes de glace pure à condition que la couche de glace soit caractérisée par une densité de 917 kg m[indice supérieur _3] et une taille des grains de neige de 0 mm. Il a aussi été démontré que, pour des sites caractérisés par des croutes de glace, les températures de brillances simulées des couverts de neige avec des croutes de glace ayant les propriétés mesurées in-situ (RMSE=11.3 K), avaient une erreur similaire aux températures de brillances simulées des couverts de neige pour des sites n'ayant pas de croutes de glace (RMSE=11.5 K). Avec le modèle MEMLS validé pour la simulation du TR MOP du sol, de la neige et de la glace, un indice de détection des croutes de glace par télédétection MOP a été développé. Il a été démontré que le ratio de polarisation (PR) était très affecté par la présence de croutes de glace dans le couvert de neige. Avec des simulations des PR à 10.7, 19 et 37 GHz sur des sites mesurés à Churchill (Manitoba, Canada), il a été possible de déterminer des seuils entre la moyenne hivernale des PR et les valeurs des PR mesurés indiquant la présence de croutes de glace. Ces seuils ont été appliqués sur une série temporelle de PR de 33 hivers d'un pixel du Nunavik (Québec, Canada) où les conditions de sols étaient similaires à ceux observés à Churchill. Plusieurs croutes de glace ont été détectées depuis 1995 et les mêmes événements entre 2002 et 2009 que (Roy, 2014) ont été détectés. Avec une validation in-situ, il serait possible de confirmer ces événements de croutes de glace mais (Roy, 2014) a démontré que ces événements ne pouvaient être expliqués que par la présence de croutes de glace dans le couvert de neige. Ces mêmes seuils sur les PR ont été appliqués sur un pixel de l'Île Banks (Territoires du Nord-Ouest, Canada). L'événement répertorié par (Grenfell et Putkonen, 2008) a été détecté. Plusieurs autres événements de croutes de glace ont été détectés dans les années 1990 et 2000 avec ces seuils. Tous ces événements ont suivi une période où les températures de l'air étaient près ou supérieures au point de congélation et sont rapidement retombées sous le point de congélation. Les températures de l'air peuvent être utilisées pour confirmer la possibilité de présence de croutes de glace mais seul la validation in-situ peut définitivement confirmer la présence de ces croutes.Abstract : Snow cover studies are essential to better understand climatic and hydrologic processes. With recent climate change observed in the northern hemisphere, more frequent rain-on-snow and meltrefreeze events have been reported, which affect the habits of the northern comunities and the survival of arctique wildlife. Passive microwave remote sensing has proven to be a great tool to characterize the state of snow cover. Nonetheless, proper modeling of the microwave signal is needed in order to understand how the parameters of the snowpack affect the measured signal. The main objective of this study is to analyze the soil, snow and ice radiative transfer in order to better characterize snow cover properties and develop an ice lens detection index with satellite passive microwave brightness temperatures. To do so, the passive microwave radiative transfer modeling of the Microwave Emission Model of Layered Snowpacks (MEMLS) was improved in order to minimize the errors on the brightness temperature simulations in the presence of ice lenses. The first improvement to passive microwave radiative transfer modeling of snow made was the snow grain size parameterization. Two new instruments, based on short wave infrared reflectance to measure the snow specific surface area (SSA) were developed. This parameter was shown to be a more accurate and objective to characterize snow grain size. The instruments showed an uncertainty of 10% to measure the SSA of snow. Also, the SSA of snow was calibrated for passive microwave modeling in order to reduce the errors on the simulated brightness temperatures. It was showed that a correction factor of φ = 1.3 needed to be applied to the grain size parameter of MEMLS, obtain through the SSA measurements, to minimize the simulation error. The second improvement to passive microwave radiative transfer modeling was the estimation of passive microwave soil emission. In-situ microwave measurements and physical temperature profiles of frozen organic arctic soils were acquired and characterized to improve the modeling of the soil emission. Effective permittivities at 10.7, 19 and 37 GHz and effective surface roughness were determined for this type of soil and the soil brightness temperature simulations were obtain with a minimal root mean square error (RMSE) of 4.65K. With the snow grain size and soil contributions to the emitted brightness temperature optimized, it was then possible to implement a passive microwave radiative transfer module of ice into MEMLS. With this module, it was possible to demonstrate that the improved Born approximation already implemented in MEMLS was equivalent to simulating a pure ice lens when the density of the layer was set to 917 kg m−3 and the grain size to 0 mm. This study also showed that by simulating ice lenses within the snow with there measured properties, the RMSE of the simulations (RMSE= 11.3 K) was similar to the RMSE for simulations of snowpacks where no ice lenses were measured (only snow, RMSE= 11.5 K). With the validated MEMLS model for snowpacks with ice lenses, an ice index was created. It is shown here that the polarization ratio (PR) was strongly affected by the presence of ice lenses within the snowpack. With simulations of the PR at 10.7, 19 and 37 GHz from measured snowpack properties in Chucrhill (Manitoba, Canada), thresholds between the measured PR and the mean winter PR were determined to detect the presence of ice within the snowpack. These thresholds were applied to a timeseries of nearly 34 years for a pixel in Nunavik (Quebec, Canada) where the soil surface is similar to that of the Churchill site. Many ice lenses are detected since 1995 with these thresholds and the same events as Roy (2014) were detected. With in-situ validation, it would be possible to confirm the precision of these thresholds but Roy (2014) showed that these events can not be explained by anything else than the presence of an ice layer within the snowpack. The same thresholds were applied to a pixel on Banks island (North-West Territories, Canada). The 2003 event that was reported by Grenfell et Putkonen (2008) was detected by the thresholds. Other events in the years 1990 and 2000’s were detected with these thresholds. These events all follow periods where the air temperature were warm and were followed by a quick drop in air temperature which could be used to validate the presence of ice layer within the snowpack. Nonetheless, without in-situ validation, these events can not be confirmed

Coupled land surface and radiative transfer models for the analysis of passive microwave satellite observations

Soil moisture is one of the key variables controlling the water and energy exchanges between Earth’s surface and the atmosphere. Therefore, remote sensing based soil moisture information has potential applications in many disciplines. Besides numerical weather forecasting and climate research these include agriculture and hydrologic applications like flood and drought forecasting. The first satellite specifically designed to deliver operational soil moisture products, SMOS (Soil Moisture and Ocean Salinity), was launched 2009 by the European Space Agency (ESA). SMOS is a passive microwave radiometer working in the L-band of the microwave domain, corresponding to a frequency of roughly 1.4 GHz and relies on a new concept. The microwave radiation emitted by the Earth’s surface is measured as brightness temperatures in several look angles. A radiative transfer model is used in an inversion algorithm to retrieve soil moisture and vegetation optical depth, a measure for the vegetation attenuation of the soil’s microwave emission. For the application of passive microwave remote sensing products a proper validation and uncertainty assessment is essential. As these sensors have typical spatial resolutions in the order of 40 – 50 km, a validation that relies solely on ground measurements is costly and labour intensive. Here, environmental modelling can make a valuable contribution. Therefore the present thesis concentrates on the question which contribution coupled land surface and radiative transfer models can make to the validation and analysis of passive microwave remote sensing products. The objective is to study whether it is possible to explain known problems in the SMOS soil moisture products and to identify potential approaches to improve the data quality. The land surface model PROMET (PRocesses Of Mass and Energy Transfer) and the radiative transfer model L-MEB (L-band microwave emission of the Biosphere) are coupled to simulate land surface states, e.g. temperatures and soil moisture, and the resulting microwave emission. L-MEB is also used in the SMOS soil moisture processor to retrieve soil moisture and vegetation optical depth simultaneously from the measured microwave emission. The study area of this work is the Upper Danube Catchment, located mostly in Southern Germany. Since model validation is essential if model data are to be used as reference, both models are validated on different spatial scales with measurements. The uncertainties of the models are quantified. The root mean squared error between modelled and measured soil moisture at several measuring stations on the point scale is 0.065 m3/m3. On the SMOS scale it is 0.039 m3/m3. The correlation coefficient on the point scale is 0.84. As it is essential for the soil moisture retrieval from passive microwave data that the radiative transfer modelling works under local conditions, the coupled models are used to assess the radiative transfer modelling with L-MEB on the local and SMOS scales in the Upper Danube Catchment. In doing so, the emission characteristics of rape are described for the first time and the soil moisture retrieval abilities of L-MEB are assessed with a newly developed LMEB parameterization. The results show that the radiative transfer modelling works well under most conditions in the study area. The root mean squared error between modelled and airborne measured brightness temperatures on the SMOS scale is less than 6 – 9 K for the different look angles. The coupled models are used to analyse SMOS brightness temperatures and vegetation optical depth data in the Upper Danube Catchment in Southern Germany. Since the SMOS soil moisture products are degraded in Southern Germany and in different other parts of the world these analyses are used to narrow down possible reasons for this. The thorough analysis of SMOS brightness temperatures for the year 2011 reveals that the quality of the measurements is degraded like in the SMOS soil moisture product. This points towards radio frequency interference problems (RFI), that are known, but have not yet been studied thoroughly. This is consistent with the characteristics of the problems observed in the SMOS soil moisture products. In addition to that it is observed that the brightness temperatures in the lower look angles are less reliable. This finding could be used to improve the brightness temperature filtering before the soil moisture retrieval. An analysis of SMOS optical depth data in 2011 reveals that this parameter does not contain valuable information about vegetation. Instead, an unexpected correlation with SMOS soil moisture is found. This points towards problems with the SMOS soil moisture retrieval, possibly under the influence of RFI. The present thesis demonstrates that coupled land surface and radiative transfer models can make a valuable contribution to the validation and analysis of passive microwave remote sensing products. The unique approach of this work incorporates modelling with a high spatial and temporal resolution on different scales. This makes detailed process studies on the local scale as well as analyses of satellite data on the SMOS scale possible. This could be exploited for the validation of future satellite missions, e.g. SMAP (Soil Moisture Active and Passive) which is currently being prepared by NASA (National Aeronautics and Space Administration). Since RFI seems to have a considerable influence on the SMOS data due to the gained insights and the quality of the SMOS products is very good in other parts of the world, the RFI containment and mitigation efforts carried out since the launch of SMOS should be continued