Exploring the potential of high temporal resolution X-band SAR time series for various permafrost applications with ground truth observations in the Lena River Delta, Siberia.

Permafrost is a subsurface phenomenon that cannot be directly monitored with satellite remote sensing. A variety of indirect approaches are currently being developed which aim to measure permafrost-related processes and environmental variables. Results of these studies aid the planning of future satellite missions which will allow large-scale permafrost monitoring. This thesis contributes to this ongoing effort by assessing the potential of repeat-pass TerraSAR-X (TSX) time series for permafrost-related applications. For the first time, multi-year Synthetic Aperture Radar (SAR) data with high temporal (11 days) and spatial (3 m) resolution was analysed for a region characterized by continuous permafrost in the Siberian Arctic. Extensive in situ data was collected during three summer and winter expeditions to validate and interpret remote sensing results. Three case studies were carried out: (i) the detection of land surface changes (e.g. ground freezing and thawing, surface wetness variations, snow cover onset and melt); (ii) monitoring bedfast lake ice and ice phenology (freeze-up, melt onset, break-up); and (iii) differential SAR interferometry (DInSAR) for thaw subsidence monitoring. For the first two case studies, time series of both backscatter intensity and 11-day interferometric coherence (i.e. a measure of phase stability between two SAR images) were investigated. Backscatter intensity was generally shown to be insensitive to the land surface changes but responded to events that occurred at the time of TSX acquisition (rain, snow shower, melt/freeze crust on snow). Interferometric coherence decreased dramatically across the entire image upon snow cover onset and melt, permitting the possible use of coherence for the monitoring of these events. Backscatter intensity was found to be an excellent tool for the detection and monitoring of bedfast lake ice due in part to improved temporal resolution compared to previously used SAR systems. Ice phenology was mostly well tracked with backscatter intensity. Interferometric coherence was found to be sensitive to the lake ice grounding and to the onset of surface melt on the lakes with bedfast ice. The investigation of coherence was a useful preparative step for the following DInSAR analysis. For the third case study, coherent 11-day and 22-day interferograms were available only for one summer of the two-year TSX time series. The cumulative DInSAR displacement strongly underestimated the subsidence observed on the ground. In situ observations revealed high variability of subsidence, which likely caused errors in phase unwrapping. Conventional DInSAR processing might therefore not be suitable for the accurate representation of permafrost thaw subsidence. This study highlights the importance of field measurements for the quantification of thaw subsidence with DInSAR, which were mostly omitted in the previous studies. All in all, this thesis shows the limitations and potential of TSX time series to spatially and temporally monitor permafrost. It thus provides an important contribution to the methodological development of a long-term permafrost monitoring scheme

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

High resolution spatial variability in spring snowmelt for an Arctic shrub-tundra watershed

Arctic tundra environments are characterized by spatially heterogeneous end-of-winter snow cover because of high winds that erode, transport and deposit snow over the winter. This spatially variable end-of-winter snow cover subsequently influences the spatial and temporal variability of snowmelt and results in a patchy snowcover over the melt period. Documenting changes in both snow cover area (SCA) and snow water equivalent (SWE) during the spring melt is essential for understanding hydrological systems, but the lack of high-resolution SCA and SWE datasets that accurately capture micro-scale changes are not commonly available, and do not exist for the Canadian Arctic. This study applies high-resolution remote sensing measurements of SCA and SWE using a fixed-wing Unmanned Aerial System (UAS) to document snowcover changes over the snowmelt period for an Arctic tundra headwater catchment. Repeat measurements of SWE and SCA were obtained for four dominant land cover types (tundra, short shrub, tall shrub, and topographic drift) to provide observations of spatially distributed snowmelt patterns and basin-wide declines in SWE. High-resolution analysis of snowcover conditions over the melt reveal a strong relationship between land cover type, snow distribution, and snow ablation rates whereby shallow snowpacks found in tundra and short shrub regions feature rapid declines in SWE and SCA and became snow-free approximately 10 days earlier than deeper snowpacks. In contrast, tall shrub patches and topographic drift regions were characterized by large initial SWE values and featured a slow decline in SCA. Analysis of basin-wide declines in SCA and SWE reveal three distinct melt phases characterized by 1) low melt rates across a large area resulting in a minor change in SCA, but a very large decline in SWE with, 2) high melt rates resulting in drastic declines in both SCA and SWE, and 3) low melt rates over a small portion of the basin, resulting in little change to either SCA or SWE. The ability to capture high-resolution spatio-temporal changes to tundra snow cover furthers our understanding of the relative importance of various land cover types on the snowmelt timing and amount of runoff available to the hydrological system during the spring freshet

Improving flood forecasting using multi-source remote sensing data – Report of the Floodfore project

Current remote sensing satellites can provide valuable information relevant to hydrological monitoring. And by using available in situ measurements together with the satellite data the information can be even more valuable. The FloodFore project developed new methods to estimate hydrological parameters from multi source remote sensing and in situ data. These hydrological parameters are important input to the watershed simulation model in order to improve the accuracy of its forecasts. In the project several new methods were either developed or demonstrated: satellite based snow water equivalent (SWE) estimation, weather radar based accumulated precipitation estimation, satellite based soil freezing state determination, and SWE estimation with high spatial resolution using both microwave radiometer and SAR data. Also a visualisation system for multi source information was developed to demonstrate the new products to users. The effect of the snow remote sensing estimates to the hydrological forecasting accuracy was studied for the Kemijoki river basin. The commercialisation possibilities of the results of the project were also studied

Détection des cycles de gel/dégel de la couche active du sol en toundra arctique à partir d’imageries radar à synthèse d’ouverture (RSO) multicapteur en bande C

L’augmentation de la température de l’air moyenne annuelle, chiffrée à +2,3 °C pour les régions de l’arctique Canadien entre 1948 et 2016, a des impacts considérables sur le couvert nival arctique et sur la végétation en place. Ces deux paramètres influencent le régime thermique du sol et donc, les cycles de gel/dégel de sa couche active dans l’écosystème arctique. L’importance du suivi de ces cycles réside dans leur influence sur plusieurs paramètres de la cryosphère tels que le cycle hydrologique et du carbone, la saison de croissance de la végétation, l’état du pergélisol sous-jacent ainsi que l’épaisseur de sa couche active. L’utilisation de données ponctuelles ou provenant de capteurs micro-onde passive à basse résolution présente un enjeu pour le suivi spatial et temporel de ces cycles. Le projet vise à développer un algorithme de détection des cycles de gel/dégel du sol en toundra arctique à partir d’imageries RSO multicapteur (i.e., Sentinel-1 et RADARSAT-2) ayant une couverture temporelle quasi journalière en bande C, afin d’évaluer l’impact de la variabilité spatiale et temporelle des paramètres influençant le régime thermique du sol tel que, les écosystèmes terrestres (i.e., écotype) et la présence de neige. L’étude se concentre sur une zone à l’intérieur du bassin versant du lac Greiner à proximité de la ville de Cambridge Bay au Nunavut. La normalisation de l’angle d’incidence a permis de diminuer le bruit dans les séries temporelles ainsi que de rendre possible l’utilisation d'images acquises à l'intérieur de plusieurs orbites d’observation. Cela a aussi permis d’uniformiser les données des deux capteurs pour les combiner en une seule série temporelle. Deux algorithmes de détections ont été utilisés, soit un algorithme de seuil saisonnier (STA) ainsi qu’un algorithme de détection de changement (CPD). La validation s’est faite à partir des données spatialement distribuées de température du sol et de l’air indépendamment sous forme de précision (%) et de délai (#jours) de détection. Les deux algorithmes ont permis d’obtenir une précision de détection de plus de 97% sur les sites de référence. Une spatialisation, pixel par pixel, de la méthode STA a permis la création de cartes de jour de gel/dégel pour le site d’étude. La combinaison des cartes de jour de transition avec la carte d’écotype a permis de modéliser l’impact des caractéristiques des écotypes sur le jour de transition. Les résultats obtenus dans ce projet démontrent clairement le potentiel de l’utilisation des données RSO en bande C pour la détection des cycles de gel/dégel, ce qui constitue un résultat important en raison de la quantité grandissante de données à cette fréquence (e.g., RCM, Sentinel-1A-C-D). La méthode présentée dans ce projet pourrait permettre de créer des cartes de transition pour tout le bassin versant du lac Greiner à partir de données RSO en bande C.Abstract : The observed average annual surface temperature increase of 2.3°C in the Canadian Arctic regions between 1948 and 2016 has significant effects on the Arctic snow cover and on the vegetation in place. Those two parameters influence the thermal regime of the ground and therefore the freeze and thaw (F/T) cycles of the soil active layer in the Arctic tundra ecosystem. The importance of monitoring these cycles lies in their influence on several parameters of the cryosphere such as the hydrological and carbon cycle, the vegetation growing season, the state of the underlying permafrost and the thickness of its active layer. The use of punctual data or low-resolution passive microwave sensors presents a challenge for the spatial and temporal monitoring of these cycles. The project aims to develop an algorithm for soil freeze/thaw cycles detection in arctic tundra from multisensor C-band imagery (i.e., Sentinel-1 and RADARSAT-2) to assess the impact of the spatial and temporal variability of the parameters influencing the thermal regime of the ground, such as the terrestrial ecosystems (i.e., ecotype) and the snow cover. The study focused on a region of the Greiner lake watershed on Victoria Island in Nunavut. An incidence angle normalization was applied to the backscatter time series to remove influence of the acquisition angle on backscatter and to allow for the use of images acquired within several orbits of observation. This also standardized the data from the two sensors to combine them into a single time series. Two detection algorithms were used on the normalized backscatter coefficient data, namely a seasonal threshold algorithm (STA) and a change point detection algorithm (CPD). A spatially distributed network of soil and air temperature were used for validation in the form of accuracy (%) and delay (#days) of detection. Both algorithms achieved a detection accuracy of more than 97% for the entire analysis period on the reference sites. A pixel-by-pixel spatialization of the STA method allowed to create F/T transition maps for the extended study site. The combination of the transition maps with the ecotype data made it possible to model the impact of ecotype characteristics on the day of transition. The results obtained in this project clearly demonstrate the potential of using C-band for the detection of F/T cycles, which is an important aspect due to the increasing number of data at this frequency (e.g., RCM, Sentinel -1A-C-D). The method presented in this project could then make it possible to create transition maps for the entire Greiner Lake watershed from C-band SAR data and thus improve the integration of this parameter in climate models

An Investigation into the Effects of Variable Lake Ice Properties on Passive and Active Microwave Measurements Over Tundra Lakes Near Inuvik, N.W.T.

The accurate estimation of snow water equivalent (SWE) in the Canadian sub-arctic is integral to climate variability studies and water availability forecasts for economic considerations (drinking water, hydroelectric power generation). Common passive microwave (PM) snow water equivalent (SWE) algorithms that utilize the differences in brightness temperature (Tb) at 37 GHz – 19 GHz falter in lake-rich tundra environments because of the inclusion of lakes within PM pixels. The overarching goal of this research was to investigate the use of multiple platforms and methodologies to observe and quantify the effects of lake ice and sub-ice water on passive microwave emission for the purpose of improving snow water equivalent (SWE) retrieval algorithms. Using in situ snow and ice measurements as input, the Helsinki University of Technology (HUT) multi-layer snow emission model was modified to include an ice layer below the snow layer. Emission for 6.9, 19, 37 and 89 GHz were simulated at horizontal and vertical polarizations, and were validated by high resolution airborne passive microwave measurements coincident with in situ sampling sites over two lakes near Inuvik, Northwest Territories (NWT). Overall, the general magnitude of brightness temperatures were estimated by the HUT model for 6.9 and 19 GHz H/V, however the variability was not. Simulations produced at 37 GHz exhibited the best agreement relative to observed temperatures. However, emission at 37 GHz does not interact with the radiometrically cold water, indicating that ice properties controlling microwave emission are not fully captured by the HUT model. Alternatively, active microwave synthetic aperture radar (SAR) measurements can be used to identify ice properties that affect passive microwave emission. Dual polarized X-band SAR backscatter was utilized to identify ice types by the segmentation program MAGIC (MAp Guided Ice Classification). Airborne passive microwave transects were grouped by ice type classes and compared to backscatter measurements. In freshwater, where there were few areas of high bubble concentration at the ice/water interface Tbs exhibited positive correlations with cross-polarized backscatter, corresponding to ice types (from low to high emission/backscatter: clear ice, transition zone between clear and grey ice, grey ice and rafted ice). SWE algorithms were applied to emission within each ice type producing negative or near zero values in areas of low 19 GHz Tbs (clear ice, transition zone), but also produced positive values that were closer to the range of in situ measurements in areas of high 19 GHz Tbs (grey and rafted ice). Therefore, cross-polarized X-band SAR measurements can be used as a priori ice type information for spaceborne PM algorithms, providing information on ice types and ice characteristics (floating, frozen to bed), integral to future tundra-specific SWE retrieval algorithms

Forward Modelling of Multifrequency SAR Backscatter of Snow-Covered Lake Ice: Investigating Varying Snow and Ice Properties Within a Radiative Transfer Framework

Lakes are a key feature in the Northern Hemisphere landscape. The coverage of lakes by ice cover has important implications for local weather conditions and can influence energy balance. The presence of lake ice is also crucial for local economies, providing transportation routes, and acting as a source of recreation/tourism and local customs. Both lake ice cover, from which ice dates and duration can be derived (i.e., ice phenology), and ice thickness are considered as thematic variables of lakes as an essential climate variable by the Global Climate Observing System (GCOS) for understanding how climate is changing. However, the number of lake ice phenology ground observations has declined over the past three decades. Remote sensing provides a method of addressing this paucity in observations. Active microwave remote sensing, in particular synthetic aperture radar (SAR), is popular for monitoring ice cover as it does not rely on sunlight and the resolution allows for the monitoring of small and medium-sized lakes. In recent years, our understanding of the interaction between active microwave signals and lake ice has changed, shifting from a double bounce mechanism to single bounce at the ice-water interface. The single bounce, or surface scattering, at the ice-water interface is due to a rough surface and high dielectric contrast between ice and water. However, further work is needed to fully understand how changes in different lake ice properties impact active microwave signals. Radiative transfer modelling has been used to explore these interactions, but there are a variety of limitations associated with past experiments. This thesis aimed to faithfully represent lake ice using a radiative transfer framework and investigate how changes in lake ice properties impact active microwave backscatter. This knowledge was used to model backscatter throughout ice seasons under both dry and wet conditions. The radiative transfer framework used in this thesis was the Snow Microwave Radiative Transfer (SMRT) model. To investigate how broad changes in ice properties impact microwave backscatter, SMRT was used to conduct experiments on ice columns representing a shallow lake with tubular bubbles and a deep lake without tubular bubbles at L/C/X-band frequencies. The Canadian Lake Ice Model (CLIMo) was used to parameterize SMRT. Ice properties investigated included ice thickness, snow ice bubble radius and porosity, root mean square (RMS) height of the ice-water interface, correlation length of the ice-water interface, and tubular bubble radius and porosity. Modelled backscatter indicated that changes in ice thickness, snow ice porosity, and tubular bubble radius and porosity had little impact on microwave backscatter. The property that had the largest impact on backscatter was RMS height at the ice-water interface, confirming the results of other recent studies. L and C-band frequencies were found to be most sensitive to changes in RMS height. Bubble radius had a smaller impact on backscatter, but X-band was found to be most sensitive to changes in this property and would be a valuable frequency for studying surface ice conditions. From the results of these initial experiments, SMRT was then used to simulate the backscatter from lake ice for two lakes during different winter seasons. Malcolm Ramsay Lake near Churchill, Manitoba, represented a shallow lake with dense tubular bubbles and Noell lake near Inuvik, Northwest Territories, represented a deep lake with no tubular bubbles. Both field data and CLIMo simulations for the two lakes were used to parameterize SMRT. Because RMS height was determined to be the ice property that had the largest impact on backscatter, simulations focused on optimizing this value for both lakes. Modelled backscatter was validated using C-band satellite imagery for Noell Lake and L/C/X-band imagery for Malcolm Ramsay Lake. The root mean square error values for both lakes ranged from 0.38 to 2.33 dB and Spearman’s correlation coefficient (ρ) values >0.86. Modelled backscatter for Noell Lake was closer to observed values compared to Malcolm Ramsay Lake. Optimized values of RMS height provided a better fit compared to a stationary value and indicated that roughness likely increases rapidly at the start of the ice season but plateaus as ice growth slows. SMRT was found to model backscatter from ice cover well under dry conditions, however, modelling backscatter under wet conditions is equally important. Detailed field observations for Lake Oulujärvi in Finland were used to parameterize SMRT during three different conditions. The first was lake ice with a dry snow cover, the second with an overlying layer of wet snow, and the third was when a slush layer was present on the ice surface. Experiments conducted under dry conditions continued to support the dominance of scattering from the ice-water interface. However, when a layer of wet snow or slush layer was introduced the dominant scattering interface shifted to the new wet layer. Increased roughness at the boundary of these wet layers resulted in an increase in backscatter. The increase in backscatter is attributed to the higher dielectric constant value of these layers. The modelled backscatter was found to be representative of observed backscatter from Sentinel-1. The body of work of this thesis indicated that the SMRT framework can be used to faithfully represent lake ice and model backscatter from ice covers and improved understanding of the interaction between microwave backscatter and ice properties. With this improved understanding inversion models can be developed to retrieve roughness of the ice-water interface, this could be used to build other models to estimate ice thickness based on other remote sensing data. Additionally, insights into the impact of wet conditions on radar backscatter could prove useful in identifying unsafe ice locations

The seasonal dynamics of Arctic surface hydrology in permafrost environments

Climate-induced landscape evolution is resulting in changes to biogeochemical and hydrologi- cal cycling. In the Arctic and sub-Arctic permafrost zones, rising air temperatures are warming, and in some regions even thawing, the frozen ground. Permafrost is a carbon sink. The thermal state of the ground therefore has important implications on carbon exchange with the atmo- sphere. Permafrost thaw mobilises previously sequestered carbon stocks, potentially turning these high latitude regions into a net carbon source. Borehole temperature and active layer depth measurements are the traditional means for monitoring permafrost, however these point measurements cannot easily be extrapolated to the landscape-scale; Earth Observation (EO) data may be used for such purposes. It is widely recognised that changes in the thermal state of permafrost may be associated with longterm changes in surface hydrology. As the ground shifts from a frozen to a thawed state, Arctic lakes display changes in surface extent. Therefore, it has become common practice to explore lake dynamics, using these as indicators of permafrost change; dynamics being the keyword. Surface hydrology is a dynamic process. Discharge studies in the Arctic and sub-Arctic regions are associated with flashy hydrographs. Currently, however, remote sensing of permafrost lake change is done on the scale of decades without explicitly taking seasonality and rapid hydrolog- ical phenology into consideration. To examine the seasonal changes in Arctic surface hydrology on the landscape scale high temporal resolution data are necessary. Synthetic aperture radar instruments are exemplary for such a task. The PhD research focuses on establishing operational techniques for mapping open surface water using synthetic aperture radar data, investigating straightforward raster classification methods and exploring their feasibility by undertaking map accuracy and sensitivity studies (chapter 3). The results are then used to justify error propagation when developing an auto- mated procedure that creates temporal composites of water body extent. These temporal water body classifications are the main EO product used to identify and image seasonal surface water change in Arctic permafrost environments (chapter 4). Furthermore, a terrain-based hydrolog- ical study is undertaken to explore the context of the detected changes and possible links to relief and stream channel network (chapter 5). The aim of this PhD is to demonstrate a new method of dynamic monitoring using the Euro- pean Space Agency’s Envisat Advanced Synthetic Aperture Radar, recommending its incorpo- ration in longterm lake change studies. Technical feasibility is explored, the inherent trade-off vii between spatial and temporal resolution discussed. An automated surface water change de- tection algorithm is developed and its applicability to monitoring spring floods is assessed; noting possible modifications to the drainage system given present-day land-use and land- cover changes that are taking place in the study area, the hydrocarbon-rich Yamalo-Nenets Autonomous District in the North of West Siberia (chapter 6). The key significance of this research is to improve the current knowledge of Arctic lake change by including spring flood events and seasonality in the equation. Therefore, it is strongly believed that this research is of benefit to the entire permafrost community.AM Trofaier was a recipient of a DOC-fFORTE Fellowship [Women in Research and Technology] of the Austrian Academy of Sciences