Area-based measurements of snow water equivalent (SWE) are important for understanding earth system processes such as glacier mass balance, winter hydrological storage in drainage basins, and ground thermal regimes. Remote sensing techniques are ideally suited for wide-scale area-based mapping with the most commonly used technique to measure SWE being passive microwave, which is limited to coarse spatial resolutions of 25 km or greater and to areas without significant topographic variation. Passive microwave also has a negative bias for large SWE. Another method is repeat-pass synthetic aperture radar interferometry (InSAR) that allows measurement of SWE change at much higher spatial resolution. However, it has not been widely adopted because (1) the phase unwrapping problem has not been robustly addressed, especially for interferograms with poor coherence, and (2) SWE change maps scaled directly from repeat-pass interferograms are not an absolute measurement but contain unknown offsets for each contiguous coherent area. We develop and test a novel method for repeat-pass InSAR-based dry-snow SWE estimation that exploits the sensitivity of the dry-snow refraction-induced InSAR phase to topographic variations. The method robustly estimates absolute SWE change at spatial resolutions of < 1 km without the need for phase unwrapping. We derive a quantitative signal model for this new SWE change estimator and identify the relevant sources of bias. The method is demonstrated using both simulated SWE distributions and a 9-year RADARSAT-2 (C-band, 5.405 GHz) spotlight-mode dataset near Inuvik, Northwest Territories (NWT), Canada. SWE results are compared to in situ snow survey measurements and estimates from ERA5 reanalysis. Our method performs well in high-relief areas, thus providing complementary coverage to passive-microwave-based SWE estimation. Further, our method has the advantage of requiring only a single wavelength band and thus can utilize existing spaceborne synthetic aperture radar systems.

Interferometric synthetic aperture radar (InSAR) time series analysis of natural terrain allows for characterization of long-term geophysical trends over extended areas and, in the case of distributed scatterers (DSs), is significantly enhanced by methods that exploit the full complex-valued scattering statistics. Phase-linking (PL) estimators impose a phase-closure constraint in order to estimate the temporal wrapped-phase history of a DS directly from its complex backscatter sample coherence matrix. Some PL methods, such as the SqueeSAR and maximum-likelihood-estimator of Interferometric phase (EMI) estimators, rely on knowledge of the coherence magnitude matrix. The true coherence magnitude is a priori unknown and must therefore be estimated from the data. Bias in these estimated coherence magnitudes reduces PL performance when the true coherence magnitude is low. Many areas of the Earth are seasonally snow-covered and, for natural terrain, this leads to severe cross-season decorrelation. This poses a significant challenge for PL estimators due to bias of the near-zero cross-season coherence magnitude estimates. We introduce a clustering approach to mitigate the PL estimator bias problem that exploits the fact that in natural terrain, many DSs decorrelate similarly. This allows for averaging over large numbers of same-behaving DS, which provides robust debiasing of the coherence magnitudes used during PL. We apply our method to a RADARSAT-2 spotlight-mode InSAR dataset over a site in the western Canadian Arctic and demonstrate significant reductions in a posteriori phase variance when compared to existing PL methods.

We discuss enhanced processing methods for high resolution Synthetic Aperture Radar (SAR) interferometry (InSAR) to monitor small landslides with difficult spatial characteristics, such as very steep and rugged terrain, strong spatially heterogeneous surface motion, and coherence-compromising factors, including vegetation and seasonal snow cover. The enhanced methods mitigate phase bias induced by atmospheric effects, as well as topographic phase errors in coherent regions of layover, and due to inaccurate blending of high resolution discontinuous with lower resolution background Digital Surface Models (DSM). We demonstrate the proposed methods using TerraSAR-X (TSX) Staring Spotlight InSAR data for three test sites reflecting diverse challenging landslide-prone mountain terrains in British Columbia, Canada. Comparisons with corresponding standard processing methods show significant improvements with resulting displacement residuals that reveal additional movement hotspots and unprecedented spatial detail for active landslides/rockfalls at the investigated sites.

Keywords: InSAR, Permafrost, Infrastructure Monitoring, Deformation Monitoring The combined effect of climate change and accelerated economic development in Northern regions increases the threat of permafrost related surface deformation to buildings and transportation infrastructure. Satellite based InSAR provides a means for monitoring infrastructure that may be both remote and spatially extensive. However, permafrost poses challenges for InSAR monitoring due to the complex temporal deformation patterns caused by both seasonal active layer fluctuations and long-term changes in permafrost thickness. These dynamics suggest a need for increasing the temporal resolution of multi-temporal InSAR methods. To address this issue we have developed a method that combines and jointly processes two or more same side geometry InSAR stacks to provide a high-temporal resolution estimate of surface deformation. The method allows for combining stacks from more than a single SAR sensor and for a combination of frequency bands. Data for this work have been collected and analysed for an area near the community of Umiujaq, Quebec in Northern Canada and include scenes from RADARSAT-2, TerraSAR-X and COSMO-SkyMed. Multiple stack based surface deformation estimates are compared for several cases including results from the three sensors individually and for all sensors combined. The test cases show substantially similar surface deformation results which correlate well with surficial geology. The best spatial coverage of coherent targets was achieved when data from all sensors were combined. The proposed multiple stack method is demonstrated to improve the estimation of surface deformation in permafrost affected areas and shows potential for deriving InSAR based permafrost classification maps to aid in the monitoring of Northern infrastructure.

Snow water equivalent (SWE) is an important surface parameter for understanding a number of Earth system processes. Synthetic aperture radar (SAR) has considerable potential for measuring SWE of dry-snow because SAR can penetrate through the snow to the ground surface and is both amplitude- and phase-sensitive to refraction from the snow. Previous work on refraction-based SWE measurement by SAR has utilized the repeat-pass InSAR phase signal to estimate changes in SWE that occur between SAR acquisitions. These are subject to temporal decorrelation effects and consider only the refraction that occurs along the SAR beam center rather than the entire synthetic aperture. This study examines the refractive effect of dry-snow along the synthetic aperture and its impact on SAR image formation including defocusing and phase bias of the system impulse response. Snow phase compensation during time domain processing to recover the snow-free impulse response function (IRF) is described and demonstrated. The feasibility of using the mapdrift and image sharpness autofocus methods to estimate SWE is examined, and the effect of key system parameters on the estimation performance is derived. Experimental validation of the method was conducted by acquiring L-band data with the Simon Fraser University (SFU) Airborne SAR System over a pair of corner reflectors installed on the Kluane icefield in northwestern Canada. Results from both simulations and the icefield experiment are presented and compared including an analysis of errors affecting the estimation.

Area-based measurements of snow water equivalent (SWE) are important for understanding earth system processes such as glacier mass balance, winter hydrological storage in drainage basins and ground thermal regimes. Remote sensing techniques are ideally suited for wide-scale area-based mapping with the most commonly used technique to measure SWE being passive-microwave, which is limited to coarse spatial resolutions of 25 km or greater, and to areas without significant topographic variation. Passive-microwave also has a negative bias for large SWE. Repeat-pass synthetic aperture radar interferometry (InSAR) as an alternate technique allows measurement of SWE change at much higher spatial resolution. However, it has not been widely adopted because: (1) the phase unwrapping problem has not been robustly addressed, especially for interferograms with poor coherence and; (2) SWE change maps scaled directly from repeat-pass interferograms are not an absolute measurement but contain unknown offsets for each contiguous coherent area. We develop and test a novel method for repeat-pass InSAR based dry-snow SWE estimation that exploits the sensitivity of the dry-snow refraction-induced InSAR phase to topographic variations. The method robustly estimates absolute SWE change at spatial resolutions of < 1 km, without the need for phase unwrapping. We derive a quantitative signal model for this new SWE change estimator and identify the relevant sources of bias. The method is demonstrated using both simulated SWE distributions and a 9-year RADARSAT-2 spotlight-mode dataset near Inuvik, NWT, Canada. SWE results are compared to in situ snow survey measurements and estimates from ERA5 reanalysis. Our method performs well in high-relief areas and in areas with high SWE (> 150 mm), thus providing complementary coverage to other passive- and active-microwave based SWE estimation methods. Further, our method has the advantage of requiring only a single wavelength band and thus can utilize existing spaceborne synthetic aperture radar systems. In application, a first order analysis of SWE trends within three drainage basins suggests that differences between basin-level accumulations are a function of major landcover types, and that re-vegetation following a forest-tundra fire that occurred over 50 years ago continues to affect the spatial distribution of SWE accumulation in the study area.