Spatio-temporal influence of tundra snow properties on Ku-band (17.2 GHz) backscatter

During the 2010/11 boreal winter, a distributed set of backscatter measurements was collected using a ground-based Ku-band (17.2 GHz) scatterometer system at 26 open tundra sites. A standard snow-sampling procedure was completed after each scan to evaluate local variability in snow layering, depth, density and water equivalent (SWE) within the scatterometer field of view. The shallow depths and large basal depth hoar encountered presented an opportunity to evaluate backscatter under a set of previously untested conditions. Strong Ku-band response was found with increasing snow depth and snow water equivalent (SWE). In particular, co-polarized vertical backscatter increased by 0.82 dB for every 1 cm increase in SWE (R2 = 0.62). While the result indicated strong potential for Ku-band retrieval of shallow snow properties, it did not characterize the influence of sub-scan variability. An enhanced snow-sampling procedure was introduced to generate detailed characterizations of stratigraphy within the scatterometer field of view using near-infrared photography along the length of a 5m trench. Changes in snow properties along the trench were used to discuss variations in the collocated backscatter response. A pair of contrasting observation sites was used to highlight uncertainties in backscatter response related to short length scale spatial variability in the observed tundra environment

Exploiting the ANN Potential in Estimating Snow Depth and Snow Water Equivalent From the Airborne SnowSAR Data at X- and Ku-Bands

Within the framework of European Space Agency (ESA) activities, several campaigns were carried out in the last decade with the purpose of exploiting the capabilities of multifrequency synthetic aperture radar (SAR) data to retrieve snow information. This article presents the results obtained from the ESA SnowSAR airborne campaigns, carried out between 2011 and 2013 on boreal forest, tundra and alpine environments, selected as representative of different snow regimes. The aim of this study was to assess the capability of X- and Ku-bands SAR in retrieving the snow parameters, namely snow depth (SD) and snow water equivalent (SWE). The retrieval was based on machine learning (ML) techniques and, in particular, of artificial neural networks (ANNs). ANNs have been selected among other ML approaches since they are capable to offer a good compromise between retrieval accuracy and computational cost. Two approaches were evaluated, the first based on the experimental data (data driven) and the second based on data simulated by the dense medium radiative transfer (DMRT). The data driven algorithm was trained on half of the SnowSAR dataset and validated on the remaining half. The validation resulted in a correlation coefficient R ≃ 0.77 between estimated and target SD, a root-mean-square error (RMSE) ≃ 13 cm, and bias = 0.03 cm. ANN algorithms specific for each test site were also implemented, obtaining more accurate results, and the robustness of the data driven approach was evaluated over time and space. The algorithm trained with DMRT simulations and tested on the experimental dataset was able to estimate the target parameter (SWE in this case) with R = 0.74, RMSE = 34.8 mm, and bias = 1.8 mm. The model driven approach had the twofold advantage of reducing the amount of in situ data required for training the algorithm and of extending the algorithm exportability to other test sites

Airborne investigation of quasi-specular Ku-band radar scattering for satellite altimetry over snow-covered Arctic sea ice

Surface-based Ku-band radar altimetry investigations indicate the radar signal is typically backscattered from well above the snow-sea ice interface. However, this would induce a bias in satellite altimeter sea ice thickness retrievals not reflected by buoy validation. Our study presents a mechanism to potentially explain this paradox: probabilistic quasi-specular radar scattering from the snow-ice interface. We introduce the theory for this mechanism before identifying it in airborne Ku-band radar observations collected over landfast first year Arctic sea ice near Eureka, Canada, in spring 2016. Based on SAR data, this study area likely represents level first year sea ice across the Arctic. Radar backscatter from the snow and ice interfaces were estimated by co-aligning laser scanner and radar observations with in situ measurements. On average, 4-5 times more radar power was scattered from the snow-ice than the air-snow interface over first-year ice. However, return power varied by up to 20 dB between consecutive radar echoes, particularly from the snow-ice interface, depending on local slope and roughness. Measured laser-radar snow depths were more accurate when radar returns were specular, but there was no systematic bias between airborne and in situ snow depths. The probability and strength of quasi-specular returns depend on the measuring height above and slope distribution of sea ice, so these findings have implications for satellite altimetry snow depth and freeboard estimates. This mechanism could explain the apparent differences in Ku-band radar penetration into snow on sea ice when observed from the range of a surface-, airborne- or satellite-based sensor

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

Developing Parameter Constraints for Radar-based SWE Retrievals

Terrestrial snow is an important freshwater reservoir with significant influence on the climate and energy balance. It exhibits natural spatiotemporal variability which has been enhanced by climate change, thus it is important to monitor on a large scale. Active microwave, or radar remote sensing has shown frequency-dependent promise in this regard, however, interpretation remains a challenge. The aim of this thesis was to develop constraints for radar based SWE retrievals which characterize and limit uncertainty with a focus on the underlying physical processes, snowpack stratigraphy, the influence of vegetation, and effects of background scattering. The University of Waterloo Scatterometer (UWScat) was used to make measurements at 9.6 and 17.2 GHz of snow and bare ground in a series of field-based campaigns in Maryhill and Englehart, ON, Grand Mesa, CO (NASA SnowEx campaign, year 1), and Trail Valley Creek, NT. Additional measurements from Tobermory, ON, and Churchill, MB (Canadian Snow and Ice Experiment) were included. The Microwave Emission Model for Layered Snowpacks, Version 3, adapted for backscattering (MEMLS3&a) was used to explore snowpack parameterization and SWE retrieval and the Freeman-Durden three component decomposition (FD3c) was used to leverage the polarimetric response. Physical processes in the snow accumulation environment demonstrated influence on regional snowpack parameterization and constraints in a SWE retrieval context with a single-layer snowpack parameterization for Maryhill, ON and a two-layer snowpack parameterization for Englehart, ON resulting in a retrieval RMSE of 21.9 mm SWE and 24.6 mm SWE, respectively. Use of in situ snow depths improved RMSE to 12.0 mm SWE and 10.9 mm SWE, while accounting for soil scattering effects further improved RMSE by up to 6.3 mm SWE. At sites with vegetation and ice lenses, RMSE improved from 60.4 mm SWE to 21.1 mm SWE when in situ snow depths were used. These results compare favorably with the common accuracy requirement of RMSE ≤ 30 mm and underscore the importance of understanding the driving physical processes in a snow accumulation environment and the utility of their regional manifestation in a SWE retrieval context. A relationship between wind slab thickness and the double-bounce component of the FD3c in a tundra snowpack was introduced for incidence angles ≥ 46° and wind slab thickness ≥ 19 cm. Estimates of wind slab thickness and SWE resulted in an RMSE of 6.0 cm and 5.5 mm, respectively. The increased double-bounce scattering was associated with path length increase within a growing wind slab layer. Signal attenuation in a sub-canopy SWE retrieval was also explored. The volume scattering component of the FD3c yielded similar performance to forest fraction in the retrieval with several distinct advantages including a real-time description of forest condition, accounting for canopy geometry without ancillary information, and providing coincident information on forest canopy in remote locations. Overall, this work demonstrated how physical processes can manifest regional outcomes, it quantified effects of natural inclusions and background scattering on SWE retrievals, it provided a means to constrain wind slab thickness in a tundra environment, and it improved characterization of coniferous forest in a sub-canopy SWE retrieval context. Future work should focus on identifying ice and vegetation conditions prior to SWE retrieval, testing the spatiotemporal validity of the methods developed herein, and finally, improving the integration of snowpack attenuation within retrieval efforts

Multiple Volume Scattering in Random Media and Periodic Structures with Applications in Microwave Remote Sensing and Wave Functional Materials

The objective of my research is two-fold: to study wave scattering phenomena in dense volumetric random media and in periodic wave functional materials. For the first part, the goal is to use the microwave remote sensing technique to monitor water resources and global climate change. Towards this goal, I study the microwave scattering behavior of snow and ice sheet. For snowpack scattering, I have extended the traditional dense media radiative transfer (DMRT) approach to include cyclical corrections that give rise to backscattering enhancements, enabling the theory to model combined active and passive observations of snowpack using the same set of physical parameters. Besides DMRT, a fully coherent approach is also developed by solving Maxwell’s equations directly over the entire snowpack including a bottom half space. This revolutionary new approach produces consistent scattering and emission results, and demonstrates backscattering enhancements and coherent layer effects. The birefringence in anisotropic snow layers is also analyzed by numerically solving Maxwell’s equation directly. The effects of rapid density fluctuations in polar ice sheet emission in the 0.5~2.0 GHz spectrum are examined using both fully coherent and partially coherent layered media emission theories that agree with each other and distinct from incoherent approaches. For the second part, the goal is to develop integral equation based methods to solve wave scattering in periodic structures such as photonic crystals and metamaterials that can be used for broadband simulations. Set upon the concept of modal expansion of the periodic Green’s function, we have developed the method of broadband Green’s function with low wavenumber extraction (BBGFL), where a low wavenumber component is extracted and results a non-singular and fast-converging remaining part with simple wavenumber dependence. We’ve applied the technique to simulate band diagrams and modal solutions of periodic structures, and to construct broadband Green’s functions including periodic scatterers.PHDElectrical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/135885/1/srtan_1.pd

Multiple Volume Scattering in Random Media and Periodic Structures with Applications in Microwave Remote Sensing and Wave Functional Materials

The objective of my research is two-fold: to study wave scattering phenomena in dense volumetric random media and in periodic wave functional materials. For the first part, the goal is to use the microwave remote sensing technique to monitor water resources and global climate change. Towards this goal, I study the microwave scattering behavior of snow and ice sheet. For snowpack scattering, I have extended the traditional dense media radiative transfer (DMRT) approach to include cyclical corrections that give rise to backscattering enhancements, enabling the theory to model combined active and passive observations of snowpack using the same set of physical parameters. Besides DMRT, a fully coherent approach is also developed by solving Maxwell’s equations directly over the entire snowpack including a bottom half space. This revolutionary new approach produces consistent scattering and emission results, and demonstrates backscattering enhancements and coherent layer effects. The birefringence in anisotropic snow layers is also analyzed by numerically solving Maxwell’s equation directly. The effects of rapid density fluctuations in polar ice sheet emission in the 0.5~2.0 GHz spectrum are examined using both fully coherent and partially coherent layered media emission theories that agree with each other and distinct from incoherent approaches. For the second part, the goal is to develop integral equation based methods to solve wave scattering in periodic structures such as photonic crystals and metamaterials that can be used for broadband simulations. Set upon the concept of modal expansion of the periodic Green’s function, we have developed the method of broadband Green’s function with low wavenumber extraction (BBGFL), where a low wavenumber component is extracted and results a non-singular and fast-converging remaining part with simple wavenumber dependence. We’ve applied the technique to simulate band diagrams and modal solutions of periodic structures, and to construct broadband Green’s functions including periodic scatterers.PHDElectrical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/137141/1/srtan_1.pd