Characterizing tundra snow sub-pixel variability to improve brightness temperature estimation in satellite SWE retrievals

Topography and vegetation play a major role in sub-pixel variability of Arctic snowpack properties, but are not considered in current passive microwave (PMW) satellite SWE retrievals. Simulation of sub-pixel variability of snow properties is also problematic when downscaling snow and climate models. In this study, we simplified observed variability of snowpack properties (depth, density, microstructure) in a two-layer model with mean values and distributions of two multi-year tundra dataset so they could be incorporated in SWE retrieval schemes. Spatial variation of snow depth was parametrized by a log-normal distribution with mean (μsd) values and coefficients of variation (CVsd). Snow depth variability (CVsd) was found to increase as a function of the area measured by a Remotely Piloted Aircraft System (RPAS). Distributions of snow specific surface area (SSA) and density were found for the wind slab (WS) and depth hoar (DH) layers. The mean depth hoar fraction (DHF) was found to be higher in Trail Valley Creek (TVC) than in Cambridge Bay (CB) where TVC is at a lower latitude with a sub-arctic shrub tundra compared to CB which is a graminoid tundra. DHF were fitted with a gaussian process and predicted from snow depth. Simulations of brightness temperatures using the Snow Microwave Radiative Transfer (SMRT) model incorporating snow depth and DHF variation were evaluated with measurements from the Special Sensor Microwave/Imager and Sounder (SSMIS) sensor. Variation in snow depth (CVsd) is proposed as an effective parameter to account for sub-pixel variability in PMW emission, improving simulation by 8K. SMRT simulations using a CVsd of 0.9 best matched CVsd observations from spatial datasets for areas > 3 km2, which is comparable to the 3.125 km pixel size of the Equal-Area Scalable Earth (EASE) grid 2.0 enhanced resolution at 37 GHz

Snow Properties Retrieval Using Passive Microwave Observations

Seasonal snow cover, the second-largest component of the cryosphere, is crucial in controlling the climate system, through its important role in modifying Earth’s albedo. The temporal variability of snow extent and its physical properties in the seasonal cycle also make up a significant element to the cryospheric energy balance. Thus, seasonal snowcover should be monitored not only for its climatological impacts but also for its rolein the surface-water supply, ground-water recharge, and its insolation properties at local scales. Snowpack physical properties strongly influence the emissions from the substratum, making feasible snow property retrieval by means of the surface brightness temperature observed by passive microwave sensors. Depending on the observing spatial resolution, the time series records of daily snow coverage and a snowpacks most-critical properties such as the snow depth and snow water equivalent (SWE) could be helpful in applications ranging from modeling snow variations in a small catchment to global climatologic studies. However, the challenge of including spaceborne snow water equivalent (SWE) products in operational hydrological and hydroclimate modeling applications is very demanding with limited uptake by these systems. Various causes have been attributed to this lack of up-take but most stem from insufficient SWE accuracy. The root causes of this challenge includes the coarse spatial resolution of passive microwave (PM) observations that observe highly aggregated snowpack properties at the spaceborne scale, and inadequacies during the retrieval process that are caused by uncertainties with the forward emission modeling of snow and challenges to find robust parameterizations of the models. While the spatial resolution problem is largely in the realm of engineering design and constrained by physical restrictions, a better understanding of the whole range of retrieval methodologies can provide the clarity needed to move the thinking forward in this important field. Following a review on snow depth and SWE retrieval methods using passive microwave remote sensing observations, this research employs a forward emission model to simulate snowpacks emission and compare the results to the PM airborne observations. Airborne radiometer observations coordinated with ground-based in-situ snow measurements were acquired in the Canadian high Arctic near Eureka, NT, in April 2011. The observed brightness temperatures (Tb) at 37 GHz from typical moderate density dry snow in mid-latitudes decreases with increasing snow water equivalent (SWE) due to the volume scattering of the ground emissions by the overlying snow. At a certain point, however, as SWE increases, the emission from the snowpack offsets the scattering of the sub-nivean emission. In tundra snow, the Tb slope reversal occurs at shallower snow thicknesses. While it has been postulated that the inflection point in the seasonal time series of observed Tb V 37 GHz of tundra snow is controlled by the formation of a thick wind slab layer, the simulation of this effect has yet to be confirmed. Therefore, the Dense Media Radiative Transfer Theory forMulti Layered (DMRT-ML) snowpack is used to predict the passive microwave response from airborne observations over shallow, dense, slab-layered tundra snow. The DMRT-ML was parameterized with the in-situ snow measurements using a two-layer snowpack and run in two configurations: a depth hoar and a wind slab dominated pack. Snow depth retrieval from passive microwave observations without a-priori information is a highly underdetermined system. An accurate estimate of snow depth necessitates a-priori information of snowpack properties, such as grain size, density, physical temperature and stratigraphy, and, very importantly, a minimization of this a prior information requirement. In previous studies, a Bayesian Algorithm for Snow Water Equivalent (SWE) Estimation (BASE) have been developed, which uses the Monte Carlo Markov Chain (MCMC) method to estimate SWE for taiga and alpine snow from 4-frequency ground-based radiometer Tb. In our study, BASE is used in tundra snow for datasets of 464 footprints inthe Eureka region coupled with airborne passive microwave observations—the same fieldstudy that forward modelling was evaluated. The algorithm searches optimum posterior probability distribution of snow properties using a cost function between physically based emission simulations and Tb observations. A two-layer snowpack based on local snow cover knowledge is assumed to simulate emission using the Dense Media Radiative Transfer-Multi Layered (DMRT-ML) model. Overall, the results of this thesis reinforce the applicability of a physics-based emission model in SWE retrievals. This research highlights the necessity to consider the two-part emission characteristics of a slab-dominated tundra snowpack and suggests performing inversion in a Bayesian framework

Modeling the Observed Microwave Emission from Shallow Multi-Layer Tundra Snow Using DMRT-ML

The observed brightness temperatures (Tb) at 37 GHz from typical moderate density dry snow in mid-latitudes decreases with increasing snow water equivalent (SWE) due to volume scattering of the ground emissions by the overlying snow. At a certain point, however, as SWE increases, the emission from the snowpack offsets the scattering of the sub-nivean emission. In tundra snow, the Tb slope reversal occurs at shallower snow thicknesses. While it has been postulated that the inflection point in the seasonal time series of observed Tb V 37 GHz of tundra snow is controlled by the formation of a thick wind slab layer, the simulation of this effect has yet to be confirmed. Therefore, the Dense Media Radiative Transfer Theory for Multi Layered (DMRT-ML) snowpack is used to predict the passive microwave response from airborne observations over shallow, dense, slab-layered tundra snow. Airborne radiometer observations coordinated with ground-based in situ snow measurements were acquired in the Canadian high Arctic near Eureka, NT, in April 2011. The DMRT-ML was parameterized with the in situ snow measurements using a two-layer snowpack and run in two configurations: a depth hoar and a wind slab dominated pack. With these two configurations, the calibrated DMRT-ML successfully predicted the Tb V 37 GHz response (R correlation of 0.83) when compared with the observed airborne Tb footprints containing snow pits measurements. Using this calibrated model, the DMRT-ML was applied to the whole study region. At the satellite observation scale, observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) over the study area reflected seasonal differences between Tb V 37 GHz and Tb V 19 GHz that supports the hypothesis of the development of an early season volume scattering depth hoar layer, followed by the growth of the late season emission-dominated wind slab layer. This research highlights the necessity to consider the two-part emission characteristics of a slab-dominated tundra snowpack at 37 GHz Tb

Monitoring Seasonal Snow Density from Satellite Based Passive Microwave Remote Sensing and Automatic Weather Stations

Seasonal snow plays an important role in Earth’s systems and for hydrological applications one of the most important properties is the quantity of liquid water stored in the snowpack, referred to as snow water equivalent (SWE). SWE is related to the depth and density of a snowpack, so accurate estimates of both those properties are necessary to estimate SWE. However, the current understanding of snow density is limited to sparsely distributed in situ samples, which is especially limiting in an environment with restricted access like the Canadian tundra. Models can be used to estimate snow density in lieu of in situ sampling and there are a variety of such models available. However, it was determined that none of the available snow density models were entirely suitable for an environment like the Canadian tundra, each for their own reasons. A new remote sensing algorithm was proposed to estimate snow density from satellite based passive microwave observations and operational automatic weather station (AWS) networks. In this research, an experiment was designed to evaluate the potential for the remote sensing algorithm to monitor snow density in the Canadian Tundra. AWS data were used parametrize a two-layer snowpack model (representing a depth hoar layer underlying a wind slab) and 3D gradient descent machine learning was used to isolate the volume scattering contributions of each layer density independently. New components were added to the machine learning cost function to incorporate prior knowledge and constrain the model’s behaviour. The model was trained at the AWS site in Eureka, Nunavut and was then applied to AWS sites distributed across the Canadian tundra. Model performance was quite consistent at high arctic sites but began to degrade across the subarctic with increased distance from the training site, suggesting the need for more robust model training and forcing in the future. Estimation skill consistently improved over the course of algorithm runs and snow density estimates were often close to the ±10% uncertainty range of the in situ samples by the end of the season – showing good promise for estimating snow density at peak SWE accumulation, which could be useful for applications where total water storage in the snowpack is of concern

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

Amélioration de la caractérisation de la neige et du sol arctique afin d’améliorer la prédiction de l’équivalent en eau de la neige en télédétection micro-ondes

Le phénomène de l’amplification arctique consiste en une augmentation plus prononcée des températures de surface dans cette région que sur le reste du globe. Ce phénomène est notamment dû à la diminution marquée du couvert nival provoquant un déséquilibre dans le bilan d’énergie de surface via une réduction généralisée de l’albédo (rétroaction positive). L’accélération du réchauffement est jusqu’à trois fois plus élevée dans ces régions. Il est donc primordial, dans un contexte de changement climatique arctique, de poursuivre et d’améliorer le suivi à grande échelle du couvert nival afin de mieux comprendre les processus gouvernant la variabilité spatio-temporelle du manteau neigeux. Plus spécifiquement, l’Équivalent en Eau de la Neige (EEN) est généralement utilisé pour quantifier deux propriétés (hauteur et densité) de la neige. Son estimation à grande échelle dans les régions éloignées tel que l’Arctique provient actuellement essentiellement de produits en micro-ondes passives satellitaires. Cependant, il existe encore beaucoup d’incertitudes sur les techniques d’assimilation de l’ÉEN par satellite et ce projet vise une réduction de l’erreur liée à l’estimation de l’ÉEN en explorant deux des principales sources de biais tels que : 1) la variabilité spatiale de l’épaisseur et des différentes couches du manteau neigeux arctique liées à la topographie et la végétation au sol influençant l’estimation de l’ÉEN; et 2) les modèles de transfert radiatif micro-ondes de la neige et du sol ne bénéficient pas actuellement d’une bonne paramétrisation en conditions arctiques, là où les erreurs liées à l’assimilation de l’ÉEN sont les plus importantes. L’objectif global est donc d’analyser les propriétés géophysiques du couvert nival en utilisant des outils de télédétection et de modélisation pour diminuer l’erreur liée à la variabilité spatiale locale dans l’estimation du ÉEN à grande échelle, tout en améliorant la compréhension des processus locaux qui affectent cette variabilité. Premièrement, une analyse haute résolution à l’aide de l’algorithme Random Forest a permis de prédire la hauteur de neige à une résolution spatiale de 10 m avec une RMSE de 8 cm (23%) et d’en apprendre davantage sur les processus de distribution de la neige en Arctique. Deuxièmement, la variabilité du manteaux neigeux arctique (hauteur et microstructure) a été incorporée dans des simulations en transfert radiatif micro-ondes de la neige et comparée au capteur satellitaire SSMIS. L’ajout de variabilité améliore la RMSE des simulations de 8K par rapport à un manteau neigeux uniforme. Finalement, une paramétrisation du sol gelé est présentée à l’aide de mesures de rugosité provenant de photogrammétrie (Structure-from-Motion). Cela a permis d’investiguer trois modèles de réflectivité micro-ondes du sol ainsi que la permittivité effective du sol gelé avec une rugosité SfM d’une précision de 0.1 mm. Ces données de rugosité SfM avec une permittivité optimisée (ε'_19 = 3.3, ε'_37 = 3.6) réduisent significativement l’erreur des températures de brillance simulées par rapport à des mesures au sols (RMSE = 3.1K, R^2 = 0.71) pour toutes les fréquences et polarisations. Cette thèse offre une caractérisation des variables de surface (neige et sol) en Arctique en transfert radiatif micro-ondes qui bénéficie aux multiples modélisations (climatiques et hydrologiques) de la cryosphère