State-of-the-art in studies of glacial isostatic adjustment for the British Isles: a literature review

Understanding the effects of glacial isostatic adjustment (GIA) of the British Isles is essential for the assessment of past and future sea-level trends. GIA has been extensively examined in the literature, employing different research methods and observational data types. Geological evidence from palaeo-shorelines and undisturbed sedimentary deposits has been used to reconstruct long-term relative sea-level change since the Last Glacial Maximum. This information derived from sea-level index points has been employed to inform empirical isobase models of the uplift in Scotland using trend surface and Gaussian trend surface analysis, as well as to calibrate more theory-driven GIA models that rely on Earth mantle rheology and ice sheet history. Furthermore, current short-term rates of GIA-induced crustal motion during the past few decades have been measured using different geodetic techniques, mainly continuous GPS (CGPS) and absolute gravimetry (AG). AG-measurements are generally employed to increase the accuracy of the CGPS estimates. Synthetic aperture radar interferometry (InSAR) looks promising as a relatively new technique to measure crustal uplift in the northern parts of Great Britain, where the GIA-induced vertical land deformation has its highest rate. This literature review provides an in-depth comparison and discussion of the development of these different research approaches

Monitoring permafrost environments with Synthetic Aperture Radar (SAR) sensors

Permafrost occupies approximately 24% of the exposed land area in the Northern Hemisphere. It is an important element of the cryosphere and has strong impacts on hydrology, biological processes, land surface energy budget, and infrastructure. For several decades, surface air temperatures in the high northern latitudes have warmed at approximately twice the global rate. Permafrost temperatures have increased in most regions since the early 1980s, the averaged warming north of 60°N has been 1-2°C. In-situ measurements are essential to understanding physical processes in permafrost terrain, but they have several limitations, ranging from difficulties in drilling to the representativeness of limited single point measurements. Remote sensing is urgently needed to supplement ground-based measurements and extend the point observations to a broader spatial domain. This thesis concentrates on the sub-arctic permafrost environment monitoring with SAR datasets. The study site is selected in a typical discontinuous permafrost region in the eastern Canadian sub-Arctic. Inuit communities in Nunavik and Nunatsiavut in the Canadian eastern sub-arctic are amongst the groups most affected by the impacts of climate change and permafrost degradation. Synthetic Aperture Radar (SAR) datasets have advantages for permafrost monitoring in the Arctic and sub-arctic regions because of its high resolution and independence of cloud cover and solar illumination. To date, permafrost environment monitoring methods and strategies with SAR datasets are still under development. The variability of active layer thickness is a direct indication of permafrost thermal state changes. The Differential SAR Interferometry (D-InSAR) technique is applied in the study site to derive ground deformation, which is introduced by the thawing/freezing depth of active layer and underlying permafrost. The D-InSAR technique has been used for the mapping of ground surface deformation over large areas by interpreting the phase difference between two signals acquired at different times as ground motion information. It shows the ability to detect freeze/thaw-related ground motion over permafrost regions. However, to date, accuracy and value assessments of D-InSAR applications have focused mostly on the continuous permafrost region where the vegetation is less developed and causes fewer complicating factors for the D-InSAR application, less attention is laid on the discontinuous permafrost terrain. In this thesis, the influencing factors and application conditions for D-InSAR in the discontinuous permafrost environment are evaluated by using X- band and L-band data. Then, benefit from by the high-temporal resolution of C-band Sentinel-1 time series, the seasonal displacement is derived from small baseline subsets (SBAS)-InSAR. Landforms are indicative of permafrost presence, with their changes inferring modifications to permafrost conditions. A permafrost landscape mapping method was developed which uses multi-temporal TerraSAR-X backscatter intensity and interferometric coherence information. The land cover map is generated through the combined use of object-based image analysis (OBIA) and classification and regression tree analysis (CART). An overall accuracy of 98% is achieved when classifying rock and water bodies, and an accuracy of 79% is achieved when discriminating between different vegetation types with one year of single-polarized acquisitions. This classification strategy can be transferred to other time-series SAR datasets, e.g., Sentinel-1, and other heterogeneous environments. One predominant change in the landscape tied to the thaw of permafrost is the dynamics of thermokarst lakes. Dynamics of thermokarst lakes are developed through their lateral extent and vertical depth changes. Due to different water depth, ice cover over shallow thermokarst ponds/lakes can freeze completely to the lake bed in winter, resulting in grounded ice; while ice cover over deep thermokarst ponds/lakes cannot, which have liquid water persisting under the ice cover all winter, resulting in floating ice. Winter ice cover regimes are related to water depths and ice thickness. In the lakes having floating ice, the liquid water induces additional heat in the remaining permafrost underneath and surroundings, which contributes to further intensified permafrost thawing. SAR datasets are utilized to detect winter ice cover regimes based on the character that liquid water has a remarkably high dielectric constant, whereas pure ice has a low value. Patterns in the spatial distribution of ice-cover regimes of thermokarst ponds in a typical discontinuous permafrost region are first revealed. Then, the correlations of these ice-cover regimes with the permafrost degradation states and thermokarst pond development in two historical phases (Sheldrake catchment in the year 1957 and 2009, Tasiapik Valley 1994 and 2010) were explored. The results indicate that the ice-cover regimes of thermokarst ponds are affected by soil texture, permafrost degradation stage and permafrost depth. Permafrost degradation is difficult to directly assess from the coverage area of floating-ice ponds and the percentage of all thermokarst ponds consisting of such floating-ice ponds in a single year. Continuous monitoring of ice-cover regimes and surface areas is recommended to elucidate the hydrological trajectory of the thermokarst process. Several operational monitoring methods have been developed in this thesis work. In the meanwhile, the spatial distribution of seasonal ground thaw subsidence, permafrost landscape, thermokarst ponds and their winter ice cover regimes are first revealed in the study area. The outcomes help understand the state and dynamics of permafrost environment.Der Permafrostboden bedeckt etwa 24% der exponierten Landfläche in der nördlichen Hemisphäre. Es ist ein wichtiges Element der Kryosphäre und hat starke Auswirkungen auf die Hydrologie, die biologischen Prozesse, das Energie-Budget der Landoberfläche und die Infrastruktur. Seit mehreren Jahrzehnten erhöhen sich die Oberflächenlufttemperaturen in den nördlichen hohen Breitengraden etwa doppelt so stark wie die globale Rate. Die Temperaturen der Permafrostböden sind in den meisten Regionen seit den frühen 1980er Jahren gestiegen. Die durchschnittliche Erwärmung nördlich von 60° N beträgt 1-2°C. In-situ-Messungen sind essentiell für das Verständnis der physischen Prozesse im Permafrostgelände. Es gibt jedoch mehrere Einschränkungen, die von Schwierigkeiten beim Bohren bis hin zur Repräsentativität begrenzter Einzelpunktmessungen reichen. Fernerkundung ist dringend benötigt, um bodenbasierte Messungen zu ergänzen und punktuelle Beobachtungen auf einen breiteren räumlichen Bereich auszudehnen. Diese Dissertation konzentriert sich auf die Umweltbeobachtung der subarktischen Permafrostböden mit SAR-Datensätzen. Das Untersuchungsgebiet wurde in einer typischen diskontinuierlichen Permafrostzone in der kanadischen östlichen Sub-Arktis ausgewählt. Die Inuit-Gemeinschaften in den Regionen Nunavik und Nunatsiavut in der kanadischen östlichen Sub-Arktis gehören zu den Gruppen, die am stärksten von den Auswirkungen des Klimawandels und Permafrostdegradation betroffen sind. Synthetische Apertur Radar (SAR) Datensätze haben Vorteile für das Permafrostmonitoring in den arktischen und subarktischen Regionen aufgrund der hohen Auflösung und der Unabhängigkeit von Wolkendeckung und Sonnenstrahlung. Bis heute sind die Methoden und Strategien mit SAR-Datensätzen für Umweltbeobachtung der Permafrostböden noch in der Entwicklung. Die Variabilität der Auftautiefe der aktiven Schicht ist eine direkte Indikation der Veränderung des thermischen Zustands der Permafrostböden. Die Differential-SAR-Interferometrie(D-Insar)-Technik wird im Untersuchungsgebiet zur Ableitung der Bodendeformation, die durch Auftau- / und Gefriertiefe der aktiven Schicht und des unterliegenden Permafrostbodens eingeführt wird, eingesetzt. Die D-InSAR-Technik wurde für Kartierung der Landoberflächendeformation über große Flächen verwendet, indem der Phasenunterschied zwischen zwei zu verschiedenen Zeitpunkten als Bodenbewegungsinformation erfassten Signalen interpretiert wurde. Es zeigt die Fähigkeit, tau- und gefrierprozessbedingte Bodenbewegungen über Permafrostregionen zu detektieren. Jedoch fokussiert sich die Genauigkeit und Wertschätzung der D-InSAR-Anwendung bis heute hauptsächlich auf kontinuierliche Permafrostregion, wo die Vegetation wenig entwickelt ist und weniger komplizierte Faktoren für D-InSAR-Anwendung verursacht. Das diskontinuierliche Permafrostgelände wurde nur weniger berücksichtigt. In dieser Dissertation wurden die Einflussfaktoren und Anwendungsbedingungen für D-InSAR im diskontinuierlichen Permafrostgebiet mittels X-Band und L-Band Daten ausgewertet. Dann wurde die saisonale Verschiebung dank der hohen Auflösung der C-Band Sentinel-1 Zeitreihe von „Small Baseline Subsets (SBAS)-InSAR“ abgeleitet. Landformen weisen auf die Präsenz des Permafrosts hin, wobei deren Veränderungen auf die Modifikation der Permafrostbedingungen schließen. Eine Kartierungsmethode der Permafrostlandschaft wurde entwickelt, dabei wurde Multi-temporal TerraSAR-X Rückstreuungsintensität und interferometrische Kohärenzinformationen verwendet. Die Landbedeckungskarte wurde durch kombinierte Anwendung objektbasierter Bildanalyse (OBIA) und Klassifikations- und Regressionsbaum Analyse (CART) generiert. Eine Gesamtgenauigkeit in Höhe von 98% wurde bei Klassifikation der Gesteine und Wasserkörper erreicht. Bei Unterscheidung zwischen verschiedenen Vegetationstypen mit einem Jahr einzelpolarisierte Akquisitionen wurde eine Genauigkeit von 79% erreicht. Diese Klassifikationsstrategie kann auf andere Zeitreihen der SAR-Datensätzen, z.B. Sentinel-1, und auch anderen heterogenen Umwelten übertragen werden. Eine vorherrschende Veränderung in der Landschaft, die mit dem Auftauen des Permafrosts verbunden ist, ist die Dynamik der Thermokarstseen. Die Dynamik der Thermokarstseen ist durch Veränderungen der seitlichen Ausdehnung und der vertikalen Tiefe entwickelt. Aufgrund der unterschiedlichen Wassertiefen kann die Eisdecke über den flachen Thermokarstteichen/-seen im Winter bis auf den Wasserboden vollständig gefroren sein, was zum geerdeten Eis führt, während die Eisdecke über den tiefen Thermokarstteichen/-seen es nicht kann. In den tiefen Thermokarstteichen/-seen bleibt den ganzen Winter flüssiges Wasser unter der Eisdecke bestehen, was zum Treibeis führt. Das Wintereisdeckenregime bezieht sich auf die Wassertiefe und die Eisdicke. In den Seen mit Treibeis leitet das flüssige Wasser zusätzliche Wärme in den restlichen Permafrost darunter oder in der Umgebung, was zur weiteren Verstärkung des Permafrostauftauen beiträgt. Basiert auf den Charakter, dass das flüssige Wasser eine bemerkenswert hohe Dielektrizitätskonstante besitzt, während reines Eis einen niedrigen Wert hat, wurden die SAR Datensätzen zur Erkennung des Wintereisdeckenregimes verwendet. Zunächst wurden Schemen in der räumlichen Verteilung der Eisdeckenregimes der Thermokarstteiche in einer typischen diskontinuierlichen Permafrostregion abgeleitet. Dann wurden die Zusammenhänge dieser Eisdeckenregimes mit dem Degradationszustand des Permafrosts und der Entwicklung der Thermokarstteiche in zwei historischen Phasen (Sheldrake Einzugsgebiet in 1957 und 2009, Tasiapik Tal in 1994 und 2010) erforscht. Die Ergebnisse deuten darauf, dass die Eisdeckenregimes der Thermokarstteiche von der Bodenart, dem Degradationszustand des Permafrosts und der Permafrosttiefe beeinflusst werden. Es ist schwer, die Permafrostdegradation in einem einzelnen Jahr direkt durch den Abdeckungsbereich der Treibeis-Teiche und die Prozentzahl aller aus solchen Treibeis-Teichen bestehenden Thermokarstteiche abzuschätzen. Ein kontinuierliches Monitoring der Eisdeckenregimes und -oberflächen ist empfehlenswert, um den hydrologischen Verlauf des Thermokarstprozesses zu erläutern. In dieser Dissertation wurden mehrere operativen Monitoringsmethoden entwickelt. In der Zwischenzeit wurden die räumliche Verteilung der saisonalen Bodentauabsenkung, die Permafrostlandschaft, die Thermokarstteiche und ihre Wintereisdeckenregimes erstmals in diesem Untersuchungsgebiet aufgedeckt. Die Ergebnisse tragen dazu bei, den Zustand und die Dynamik der Permafrostumwelt zu verstehen

Atmospheric artifacts correction for InSAR using empirical model and numerical weather prediction models

lnSAR has been proved its unprecedented ability and merits of monitoring ground deformation on large scale with centimeter to millimeter scale accuracy. However, several factors affect the reliability and accuracy of its applications. Among them, atmospheric artifacts due to spatial and temporal variations of atmosphere state often pose noise to interferograms. Therefore, atmospheric artifacts m itigalion remains one of the biggest challenges to be addressed in the In SAR community. State-of-the-art research works have revealed atmospheric artifacts can be partially compensated with empirical models, temporal-spatial filtering approach in lnSAR time series, pointwise GPS zenith path delay and numerical weather prediction models. In this thesis, firstly, we further develop a covariance weighted linear empirical model correction method. Secondly, a realistic LOS direction integration approach based on global reanalysis data is employed and comprehensively compared with the conventional method that integrates along zenith direction. Finally, the realistic integration method is applied to local WRF numerical forecast model data. l'vbreover, detailed comparisons between different global reanalysis data and local WRF model are assessed. In terms of empirical models correcting methods, many publications have studied correcting stratified tropospheric phase delay by assuming a linear model between them and topography. However, most of these studies ha\19 not considered the effect of turbulent atmospheric artefacts when adjusting the linear model to data. In this thesis, an improved technique that minimizes the influence of turbulent atmosphere in the model adjustment has been presented. In the proposed algorithm, the model is adjusted to the phase differences of pixels instead of using the unwrapped phase of each pixel. In addition, the different phase differences are weighted as a function of its APS covariance estimated from an empirical variogram to reduce in the model adjustment the impact of pixel pairs with significant turbulent atmosphere. The performance of the proposed method has been validated with both simulated and real Sentinel-1 SAR data in Tenerife island, Spain. Considering methods using meteorological observations to mitigate APS, an accurate realistic com puling strategy utilizing global atmospheric reanalysis data has been implemented. With the approach, the realistic LOS path along satellite and the monitored points is considered, rather than converting from zenith path delay. Com pared with zenith delay based method, the biggest advantage is that it can avoid errors caused by anisotropic atmospheric behaviour. The accurate integration method is validated with Sentinel-1 data in three test sites: Tenerife island, Spain, Almeria, Spain and Crete island, Greece. Compared to conventional zenith method, the realistic integration method shows great improvement. A variety of global reanalysis data are available from different weather forecasting organizations, such as ERA-Interim, ERAS, MERRA2. In this study, the realistic integration mitigation method is assessed on these different reanalysis data. The results show that these data are feasible to mitigate APS to some extent in most cases. The assessment also demonstrates that the ERAS performs the best statistically, compared to other global reanalysis data. l'vbreover, as local numerical weather forecast models have the ability to predict high spatial resolution atmospheric parameters, by using which, it has the potential to achieve APS mitigation. In this thesis, the realistic integration method is also employed on the local WRF model data in Tenerife and Almeria test s ites. However, it turns out that the WRF model performs worse than the original global reanalysis data.Las técnicas lnSAR han demostrado su capacidad sin precedentes y méritos para el monitoreo de la deformaci6n del suelo a gran escala con una precisión centimétrica o incluso milimétrica. Sin embargo, varios factores afectan la fiabilidad y precisión de sus aplicaciones. Entre ellos, los artefactos atmosféricos debidos a variaciones espaciales y temporales del estado de la atm6sfera a menudo añaden ruido a los interferogramas. Por lo tanto, la mitigación de los artefactos atmosféricos sigue siendo uno de los mayores desafíos a abordar en la comunidad lnSAR. Los trabajos de investigaci6n de vanguardia han revelado que los artefactos atmosféricos se pueden compensar parcialmente con modelos empíricos, enfoque de filtrado temporal-espacial en series temporales lnSAR, retardo puntual del camino cenital con GPS y modelos numéricos de predicción meteorológica. En esta tesis, en primer lugar, desarrollamos un método de corrección de modelo empírico lineal ponderado por covarianza. En segundo lugar, se emplea un enfoque realista de integracion de dirección LOS basado en datos de reanálisis global y se compara exhaustivamente con el método convencional que se integra a lo largo de la dirección cenital. Finalmente, el método de integraci6n realista se aplica a los datos del modelo de pronóstico numérico WRF local. Ademas, se evalúan las comparaciones detalladas entre diferentes datos de reanálisis global y el modelo WRF local. En términos de métodos de corrección con modelos empíricos, muchas publicaciones han estudiado la corrección del retraso estratificado de la fase troposférica asumiendo un modelo lineal entre ellos y la topografía. Sin embargo, la mayoría de estos estudios no han considerado el efecto de los artefactos atmosféricos turbulentos al ajustar el modelo lineal a los datos. En esta tesis, se ha presentado una técnica mejorada que minimiza la influencia de la atm6sfera turbulenta en el ajuste del modelo. En el algoritmo propuesto, el modelo se ajusta a las diferencias de fase de los pixeles en lugar de utilizar la fase sin desenrollar de cada pixel. Además, las diferentes diferencias de fase se ponderan en función de su covarianza APS estimada a partir de un variograma empírico para reducir en el ajuste del modelo el impacto de los pares de pixeles con una atm6sfera turbulenta significativa. El rendimiento del método propuesto ha sido validado con datos SAR Sentinel-1 simulados y reales en la isla de Tenerife, España. Teniendo en cuenta los métodos que utilizan observaciones meteorológicas para mitigar APS, se ha implementado una estrategia de computación realista y precisa que utiliza datos de reanálisis atmosférico global. Con el enfoque, se considera el camino realista de LOS a lo largo del satélite y los puntos monitoreados, en lugar de convertirlos desde el retardo de la ruta cenital. En comparación con el método basado en la demora cenital, la mayor ventaja es que puede evitar errores causados por el comportamiento atmosférico anisotrópico. El método de integración preciso se valida con los datos de Sentinel-1 en tres sitios de prueba: la isla de Tenerife, España, Almería, España y la isla de Creta, Grecia. En comparación con el método cenital convencional, el método de integración realista muestra una gran mejora.Postprint (published version

Atmospheric artifacts correction for InSAR using empirical model and numerical weather prediction models

lnSAR has been proved its unprecedented ability and merits of monitoring ground deformation on large scale with centimeter to millimeter scale accuracy. However, several factors affect the reliability and accuracy of its applications. Among them, atmospheric artifacts due to spatial and temporal variations of atmosphere state often pose noise to interferograms. Therefore, atmospheric artifacts m itigalion remains one of the biggest challenges to be addressed in the In SAR community. State-of-the-art research works have revealed atmospheric artifacts can be partially compensated with empirical models, temporal-spatial filtering approach in lnSAR time series, pointwise GPS zenith path delay and numerical weather prediction models. In this thesis, firstly, we further develop a covariance weighted linear empirical model correction method. Secondly, a realistic LOS direction integration approach based on global reanalysis data is employed and comprehensively compared with the conventional method that integrates along zenith direction. Finally, the realistic integration method is applied to local WRF numerical forecast model data. l'vbreover, detailed comparisons between different global reanalysis data and local WRF model are assessed. In terms of empirical models correcting methods, many publications have studied correcting stratified tropospheric phase delay by assuming a linear model between them and topography. However, most of these studies ha\19 not considered the effect of turbulent atmospheric artefacts when adjusting the linear model to data. In this thesis, an improved technique that minimizes the influence of turbulent atmosphere in the model adjustment has been presented. In the proposed algorithm, the model is adjusted to the phase differences of pixels instead of using the unwrapped phase of each pixel. In addition, the different phase differences are weighted as a function of its APS covariance estimated from an empirical variogram to reduce in the model adjustment the impact of pixel pairs with significant turbulent atmosphere. The performance of the proposed method has been validated with both simulated and real Sentinel-1 SAR data in Tenerife island, Spain. Considering methods using meteorological observations to mitigate APS, an accurate realistic com puling strategy utilizing global atmospheric reanalysis data has been implemented. With the approach, the realistic LOS path along satellite and the monitored points is considered, rather than converting from zenith path delay. Com pared with zenith delay based method, the biggest advantage is that it can avoid errors caused by anisotropic atmospheric behaviour. The accurate integration method is validated with Sentinel-1 data in three test sites: Tenerife island, Spain, Almeria, Spain and Crete island, Greece. Compared to conventional zenith method, the realistic integration method shows great improvement. A variety of global reanalysis data are available from different weather forecasting organizations, such as ERA-Interim, ERAS, MERRA2. In this study, the realistic integration mitigation method is assessed on these different reanalysis data. The results show that these data are feasible to mitigate APS to some extent in most cases. The assessment also demonstrates that the ERAS performs the best statistically, compared to other global reanalysis data. l'vbreover, as local numerical weather forecast models have the ability to predict high spatial resolution atmospheric parameters, by using which, it has the potential to achieve APS mitigation. In this thesis, the realistic integration method is also employed on the local WRF model data in Tenerife and Almeria test s ites. However, it turns out that the WRF model performs worse than the original global reanalysis data.Las técnicas lnSAR han demostrado su capacidad sin precedentes y méritos para el monitoreo de la deformaci6n del suelo a gran escala con una precisión centimétrica o incluso milimétrica. Sin embargo, varios factores afectan la fiabilidad y precisión de sus aplicaciones. Entre ellos, los artefactos atmosféricos debidos a variaciones espaciales y temporales del estado de la atm6sfera a menudo añaden ruido a los interferogramas. Por lo tanto, la mitigación de los artefactos atmosféricos sigue siendo uno de los mayores desafíos a abordar en la comunidad lnSAR. Los trabajos de investigaci6n de vanguardia han revelado que los artefactos atmosféricos se pueden compensar parcialmente con modelos empíricos, enfoque de filtrado temporal-espacial en series temporales lnSAR, retardo puntual del camino cenital con GPS y modelos numéricos de predicción meteorológica. En esta tesis, en primer lugar, desarrollamos un método de corrección de modelo empírico lineal ponderado por covarianza. En segundo lugar, se emplea un enfoque realista de integracion de dirección LOS basado en datos de reanálisis global y se compara exhaustivamente con el método convencional que se integra a lo largo de la dirección cenital. Finalmente, el método de integraci6n realista se aplica a los datos del modelo de pronóstico numérico WRF local. Ademas, se evalúan las comparaciones detalladas entre diferentes datos de reanálisis global y el modelo WRF local. En términos de métodos de corrección con modelos empíricos, muchas publicaciones han estudiado la corrección del retraso estratificado de la fase troposférica asumiendo un modelo lineal entre ellos y la topografía. Sin embargo, la mayoría de estos estudios no han considerado el efecto de los artefactos atmosféricos turbulentos al ajustar el modelo lineal a los datos. En esta tesis, se ha presentado una técnica mejorada que minimiza la influencia de la atm6sfera turbulenta en el ajuste del modelo. En el algoritmo propuesto, el modelo se ajusta a las diferencias de fase de los pixeles en lugar de utilizar la fase sin desenrollar de cada pixel. Además, las diferentes diferencias de fase se ponderan en función de su covarianza APS estimada a partir de un variograma empírico para reducir en el ajuste del modelo el impacto de los pares de pixeles con una atm6sfera turbulenta significativa. El rendimiento del método propuesto ha sido validado con datos SAR Sentinel-1 simulados y reales en la isla de Tenerife, España. Teniendo en cuenta los métodos que utilizan observaciones meteorológicas para mitigar APS, se ha implementado una estrategia de computación realista y precisa que utiliza datos de reanálisis atmosférico global. Con el enfoque, se considera el camino realista de LOS a lo largo del satélite y los puntos monitoreados, en lugar de convertirlos desde el retardo de la ruta cenital. En comparación con el método basado en la demora cenital, la mayor ventaja es que puede evitar errores causados por el comportamiento atmosférico anisotrópico. El método de integración preciso se valida con los datos de Sentinel-1 en tres sitios de prueba: la isla de Tenerife, España, Almería, España y la isla de Creta, Grecia. En comparación con el método cenital convencional, el método de integración realista muestra una gran mejora

Monitoring permafrost environments with Synthetic Aperture Radar (SAR) sensors

Permafrost occupies approximately 24% of the exposed land area in the Northern Hemisphere. It is an important element of the cryosphere and has strong impacts on hydrology, biological processes, land surface energy budget, and infrastructure. For several decades, surface air temperatures in the high northern latitudes have warmed at approximately twice the global rate. Permafrost temperatures have increased in most regions since the early 1980s, the averaged warming north of 60°N has been 1-2°C. In-situ measurements are essential to understanding physical processes in permafrost terrain, but they have several limitations, ranging from difficulties in drilling to the representativeness of limited single point measurements. Remote sensing is urgently needed to supplement ground-based measurements and extend the point observations to a broader spatial domain. This thesis concentrates on the sub-arctic permafrost environment monitoring with SAR datasets. The study site is selected in a typical discontinuous permafrost region in the eastern Canadian sub-Arctic. Inuit communities in Nunavik and Nunatsiavut in the Canadian eastern sub-arctic are amongst the groups most affected by the impacts of climate change and permafrost degradation. Synthetic Aperture Radar (SAR) datasets have advantages for permafrost monitoring in the Arctic and sub-arctic regions because of its high resolution and independence of cloud cover and solar illumination. To date, permafrost environment monitoring methods and strategies with SAR datasets are still under development. The variability of active layer thickness is a direct indication of permafrost thermal state changes. The Differential SAR Interferometry (D-InSAR) technique is applied in the study site to derive ground deformation, which is introduced by the thawing/freezing depth of active layer and underlying permafrost. The D-InSAR technique has been used for the mapping of ground surface deformation over large areas by interpreting the phase difference between two signals acquired at different times as ground motion information. It shows the ability to detect freeze/thaw-related ground motion over permafrost regions. However, to date, accuracy and value assessments of D-InSAR applications have focused mostly on the continuous permafrost region where the vegetation is less developed and causes fewer complicating factors for the D-InSAR application, less attention is laid on the discontinuous permafrost terrain. In this thesis, the influencing factors and application conditions for D-InSAR in the discontinuous permafrost environment are evaluated by using X- band and L-band data. Then, benefit from by the high-temporal resolution of C-band Sentinel-1 time series, the seasonal displacement is derived from small baseline subsets (SBAS)-InSAR. Landforms are indicative of permafrost presence, with their changes inferring modifications to permafrost conditions. A permafrost landscape mapping method was developed which uses multi-temporal TerraSAR-X backscatter intensity and interferometric coherence information. The land cover map is generated through the combined use of object-based image analysis (OBIA) and classification and regression tree analysis (CART). An overall accuracy of 98% is achieved when classifying rock and water bodies, and an accuracy of 79% is achieved when discriminating between different vegetation types with one year of single-polarized acquisitions. This classification strategy can be transferred to other time-series SAR datasets, e.g., Sentinel-1, and other heterogeneous environments. One predominant change in the landscape tied to the thaw of permafrost is the dynamics of thermokarst lakes. Dynamics of thermokarst lakes are developed through their lateral extent and vertical depth changes. Due to different water depth, ice cover over shallow thermokarst ponds/lakes can freeze completely to the lake bed in winter, resulting in grounded ice; while ice cover over deep thermokarst ponds/lakes cannot, which have liquid water persisting under the ice cover all winter, resulting in floating ice. Winter ice cover regimes are related to water depths and ice thickness. In the lakes having floating ice, the liquid water induces additional heat in the remaining permafrost underneath and surroundings, which contributes to further intensified permafrost thawing. SAR datasets are utilized to detect winter ice cover regimes based on the character that liquid water has a remarkably high dielectric constant, whereas pure ice has a low value. Patterns in the spatial distribution of ice-cover regimes of thermokarst ponds in a typical discontinuous permafrost region are first revealed. Then, the correlations of these ice-cover regimes with the permafrost degradation states and thermokarst pond development in two historical phases (Sheldrake catchment in the year 1957 and 2009, Tasiapik Valley 1994 and 2010) were explored. The results indicate that the ice-cover regimes of thermokarst ponds are affected by soil texture, permafrost degradation stage and permafrost depth. Permafrost degradation is difficult to directly assess from the coverage area of floating-ice ponds and the percentage of all thermokarst ponds consisting of such floating-ice ponds in a single year. Continuous monitoring of ice-cover regimes and surface areas is recommended to elucidate the hydrological trajectory of the thermokarst process. Several operational monitoring methods have been developed in this thesis work. In the meanwhile, the spatial distribution of seasonal ground thaw subsidence, permafrost landscape, thermokarst ponds and their winter ice cover regimes are first revealed in the study area. The outcomes help understand the state and dynamics of permafrost environment.Der Permafrostboden bedeckt etwa 24% der exponierten Landfläche in der nördlichen Hemisphäre. Es ist ein wichtiges Element der Kryosphäre und hat starke Auswirkungen auf die Hydrologie, die biologischen Prozesse, das Energie-Budget der Landoberfläche und die Infrastruktur. Seit mehreren Jahrzehnten erhöhen sich die Oberflächenlufttemperaturen in den nördlichen hohen Breitengraden etwa doppelt so stark wie die globale Rate. Die Temperaturen der Permafrostböden sind in den meisten Regionen seit den frühen 1980er Jahren gestiegen. Die durchschnittliche Erwärmung nördlich von 60° N beträgt 1-2°C. In-situ-Messungen sind essentiell für das Verständnis der physischen Prozesse im Permafrostgelände. Es gibt jedoch mehrere Einschränkungen, die von Schwierigkeiten beim Bohren bis hin zur Repräsentativität begrenzter Einzelpunktmessungen reichen. Fernerkundung ist dringend benötigt, um bodenbasierte Messungen zu ergänzen und punktuelle Beobachtungen auf einen breiteren räumlichen Bereich auszudehnen. Diese Dissertation konzentriert sich auf die Umweltbeobachtung der subarktischen Permafrostböden mit SAR-Datensätzen. Das Untersuchungsgebiet wurde in einer typischen diskontinuierlichen Permafrostzone in der kanadischen östlichen Sub-Arktis ausgewählt. Die Inuit-Gemeinschaften in den Regionen Nunavik und Nunatsiavut in der kanadischen östlichen Sub-Arktis gehören zu den Gruppen, die am stärksten von den Auswirkungen des Klimawandels und Permafrostdegradation betroffen sind. Synthetische Apertur Radar (SAR) Datensätze haben Vorteile für das Permafrostmonitoring in den arktischen und subarktischen Regionen aufgrund der hohen Auflösung und der Unabhängigkeit von Wolkendeckung und Sonnenstrahlung. Bis heute sind die Methoden und Strategien mit SAR-Datensätzen für Umweltbeobachtung der Permafrostböden noch in der Entwicklung. Die Variabilität der Auftautiefe der aktiven Schicht ist eine direkte Indikation der Veränderung des thermischen Zustands der Permafrostböden. Die Differential-SAR-Interferometrie(D-Insar)-Technik wird im Untersuchungsgebiet zur Ableitung der Bodendeformation, die durch Auftau- / und Gefriertiefe der aktiven Schicht und des unterliegenden Permafrostbodens eingeführt wird, eingesetzt. Die D-InSAR-Technik wurde für Kartierung der Landoberflächendeformation über große Flächen verwendet, indem der Phasenunterschied zwischen zwei zu verschiedenen Zeitpunkten als Bodenbewegungsinformation erfassten Signalen interpretiert wurde. Es zeigt die Fähigkeit, tau- und gefrierprozessbedingte Bodenbewegungen über Permafrostregionen zu detektieren. Jedoch fokussiert sich die Genauigkeit und Wertschätzung der D-InSAR-Anwendung bis heute hauptsächlich auf kontinuierliche Permafrostregion, wo die Vegetation wenig entwickelt ist und weniger komplizierte Faktoren für D-InSAR-Anwendung verursacht. Das diskontinuierliche Permafrostgelände wurde nur weniger berücksichtigt. In dieser Dissertation wurden die Einflussfaktoren und Anwendungsbedingungen für D-InSAR im diskontinuierlichen Permafrostgebiet mittels X-Band und L-Band Daten ausgewertet. Dann wurde die saisonale Verschiebung dank der hohen Auflösung der C-Band Sentinel-1 Zeitreihe von „Small Baseline Subsets (SBAS)-InSAR“ abgeleitet. Landformen weisen auf die Präsenz des Permafrosts hin, wobei deren Veränderungen auf die Modifikation der Permafrostbedingungen schließen. Eine Kartierungsmethode der Permafrostlandschaft wurde entwickelt, dabei wurde Multi-temporal TerraSAR-X Rückstreuungsintensität und interferometrische Kohärenzinformationen verwendet. Die Landbedeckungskarte wurde durch kombinierte Anwendung objektbasierter Bildanalyse (OBIA) und Klassifikations- und Regressionsbaum Analyse (CART) generiert. Eine Gesamtgenauigkeit in Höhe von 98% wurde bei Klassifikation der Gesteine und Wasserkörper erreicht. Bei Unterscheidung zwischen verschiedenen Vegetationstypen mit einem Jahr einzelpolarisierte Akquisitionen wurde eine Genauigkeit von 79% erreicht. Diese Klassifikationsstrategie kann auf andere Zeitreihen der SAR-Datensätzen, z.B. Sentinel-1, und auch anderen heterogenen Umwelten übertragen werden. Eine vorherrschende Veränderung in der Landschaft, die mit dem Auftauen des Permafrosts verbunden ist, ist die Dynamik der Thermokarstseen. Die Dynamik der Thermokarstseen ist durch Veränderungen der seitlichen Ausdehnung und der vertikalen Tiefe entwickelt. Aufgrund der unterschiedlichen Wassertiefen kann die Eisdecke über den flachen Thermokarstteichen/-seen im Winter bis auf den Wasserboden vollständig gefroren sein, was zum geerdeten Eis führt, während die Eisdecke über den tiefen Thermokarstteichen/-seen es nicht kann. In den tiefen Thermokarstteichen/-seen bleibt den ganzen Winter flüssiges Wasser unter der Eisdecke bestehen, was zum Treibeis führt. Das Wintereisdeckenregime bezieht sich auf die Wassertiefe und die Eisdicke. In den Seen mit Treibeis leitet das flüssige Wasser zusätzliche Wärme in den restlichen Permafrost darunter oder in der Umgebung, was zur weiteren Verstärkung des Permafrostauftauen beiträgt. Basiert auf den Charakter, dass das flüssige Wasser eine bemerkenswert hohe Dielektrizitätskonstante besitzt, während reines Eis einen niedrigen Wert hat, wurden die SAR Datensätzen zur Erkennung des Wintereisdeckenregimes verwendet. Zunächst wurden Schemen in der räumlichen Verteilung der Eisdeckenregimes der Thermokarstteiche in einer typischen diskontinuierlichen Permafrostregion abgeleitet. Dann wurden die Zusammenhänge dieser Eisdeckenregimes mit dem Degradationszustand des Permafrosts und der Entwicklung der Thermokarstteiche in zwei historischen Phasen (Sheldrake Einzugsgebiet in 1957 und 2009, Tasiapik Tal in 1994 und 2010) erforscht. Die Ergebnisse deuten darauf, dass die Eisdeckenregimes der Thermokarstteiche von der Bodenart, dem Degradationszustand des Permafrosts und der Permafrosttiefe beeinflusst werden. Es ist schwer, die Permafrostdegradation in einem einzelnen Jahr direkt durch den Abdeckungsbereich der Treibeis-Teiche und die Prozentzahl aller aus solchen Treibeis-Teichen bestehenden Thermokarstteiche abzuschätzen. Ein kontinuierliches Monitoring der Eisdeckenregimes und -oberflächen ist empfehlenswert, um den hydrologischen Verlauf des Thermokarstprozesses zu erläutern. In dieser Dissertation wurden mehrere operativen Monitoringsmethoden entwickelt. In der Zwischenzeit wurden die räumliche Verteilung der saisonalen Bodentauabsenkung, die Permafrostlandschaft, die Thermokarstteiche und ihre Wintereisdeckenregimes erstmals in diesem Untersuchungsgebiet aufgedeckt. Die Ergebnisse tragen dazu bei, den Zustand und die Dynamik der Permafrostumwelt zu verstehen

Land cover and forest mapping in boreal zone using polarimetric and interferometric SAR data

Remote sensing offers a wide range of instruments suitable to meet the growing need for consistent, timely and cost-effective monitoring of land cover and forested areas. One of the most important instruments is synthetic aperture radar (SAR) technology, where transfer of advanced SAR imaging techniques from mostly experimental small test-area studies to satellites enables improvements in remote assessment of land cover on a global scale. Globally, forests are very suitable for remote sensing applications due to their large dimensions and relatively poor accessibility in distant areas. In this thesis, several methods were developed utilizing Earth observation data collected using such advanced SAR techniques, as well as their application potential was assessed. The focus was on use of SAR polarimetry and SAR interferometry to improve performance and robustness in assessment of land cover and forest properties in the boreal zone. Particular advances were achieved in land cover classification and estimating several key forest variables, such as forest stem volume and forest tree height. Important results reported in this thesis include: improved polarimetric SAR model-based decomposition approach suitable for use in boreal forest at L-band; development and demonstration of normalization method for fully polarimetric SAR mosaics, resulting in improved classification performance and suitable for wide-area mapping purposes; establishing new inversion procedure for robust forest stem volume retrieval from SAR data; developing semi-empirical method and demonstrating potential for soil type separation (mineral soil, peatland) under forested areas with L-band polarimetric SAR; developing and demonstrating methodology for simultaneous retrieval of forest tree height and radiowave attenuation in forest layer from inter-ferometric SAR data, resulting in improved accuracy and more stable estimation of forest tree height

The SAR Handbook: Comprehensive Methodologies for Forest Monitoring and Biomass Estimation

This Synthetic Aperture Radar (SAR) handbook of applied methods for forest monitoring and biomass estimation has been developed by SERVIR in collaboration with SilvaCarbon to address pressing needs in the development of operational forest monitoring services. Despite the existence of SAR technology with all-weather capability for over 30 years, the applied use of this technology for operational purposes has proven difficult. This handbook seeks to provide understandable, easy-to-assimilate technical material to remote sensing specialists that may not have expertise on SAR but are interested in leveraging SAR technology in the forestry sector

ALOS-2/PALSAR-2 Calibration, Validation, Science and Applications

Twelve edited original papers on the latest and state-of-art results of topics ranging from calibration, validation, and science to a wide range of applications using ALOS-2/PALSAR-2. We hope you will find them useful for your future research

Biomass Representation in Synthetic Aperture Radar Interferometry Data Sets

This work makes an attempt to explain the origin, features and potential applications of the elevation bias of the synthetic aperture radar interferometry (InSAR) datasets over areas covered by vegetation. The rapid development of radar-based remote sensing methods, such as synthetic aperture radar (SAR) and InSAR, has provided an alternative to the photogrammetry and LiDAR for determining the third dimension of topographic surfaces. The InSAR method has proved to be so effective and productive that it allowed, within eleven days of the space shuttle mission, for acquisition of data to develop a three-dimensional model of almost the entire land surface of our planet. This mission is known as the Shuttle Radar Topography Mission (SRTM). Scientists across the geosciences were able to access the great benefits of uniformity, high resolution and the most precise digital elevation model (DEM) of the Earth like never before for their a wide variety of scientific and practical inquiries. Unfortunately, InSAR elevations misrepresent the surface of the Earth in places where there is substantial vegetation cover. This is a systematic error of unknown, yet limited (by the vertical extension of vegetation) magnitude. Up to now, only a limited number of attempts to model this error source have been made. However, none offer a robust remedy, but rather partial or case-based solutions. More work in this area of research is needed as the number of airborne and space-based InSAR elevation models has been steadily increasing over the last few years, despite strong competition from LiDAR and optical methods. From another perspective, however, this elevation bias, termed here as the “biomass impenetrability”, creates a great opportunity to learn about the biomass. This may be achieved due to the fact that the impenetrability can be considered a collective response to a few factors originating in 3D space that encompass the outermost boundaries of vegetation. The biomass, presence in InSAR datasets or simply the biomass impenetrability, is the focus of this research. The report, presented in a sequence of sections, gradually introduces terminology, physical and mathematical fundamentals commonly used in describing the propagation of electromagnetic waves, including the Maxwell equations. The synthetic aperture radar (SAR) and InSAR as active remote sensing methods are summarised. In subsequent steps, the major InSAR data sources and data acquisition systems, past and present, are outlined. Various examples of the InSAR datasets, including the SRTM C- and X-band elevation products and INTERMAP Inc. IFSAR digital terrain/surface models (DTM/DSM), representing diverse test sites in the world are used to demonstrate the presence and/or magnitude of the biomass impenetrability in the context of different types of vegetation – usually forest. Also, results of investigations carried out by selected researchers on the elevation bias in InSAR datasets and their attempts at mathematical modelling are reviewed. In recent years, a few researchers have suggested that the magnitude of the biomass impenetrability is linked to gaps in the vegetation cover. Based on these hints, a mathematical model of the tree and the forest has been developed. Three types of gaps were identified; gaps in the landscape-scale forest areas (Type 1), e.g. forest fire scares and logging areas; a gap between three trees forming a triangle (Type 2), e.g. depending on the shape of tree crowns; and gaps within a tree itself (Type 3). Experiments have demonstrated that Type 1 gaps follow the power-law density distribution function. One of the most useful features of the power-law distributed phenomena is their scale-independent property. This property was also used to model Type 3 gaps (within the tree crown) by assuming that these gaps follow the same distribution as the Type 1 gaps. A hypothesis was formulated regarding the penetration depth of the radar waves within the canopy. It claims that the depth of penetration is simply related to the quantisation level of the radar backscattered signal. A higher level of bits per pixels allows for capturing weaker signals arriving from the lower levels of the tree crown. Assuming certain generic and simplified shapes of tree crowns including cone, paraboloid, sphere and spherical cap, it was possible to model analytically Type 2 gaps. The Monte Carlo simulation method was used to investigate relationships between the impenetrability and various configurations of a modelled forest. One of the most important findings is that impenetrability is largely explainable by the gaps between trees. A much less important role is played by the penetrability into the crown cover. Another important finding is that the impenetrability strongly correlates with the vegetation density. Using this feature, a method for vegetation density mapping called the mean maximum impenetrability (MMI) method is proposed. Unlike the traditional methods of forest inventories, the MMI method allows for a much more realistic inventory of vegetation cover, because it is able to capture an in situ or current situation on the ground, but not for areas that are nominally classified as a “forest-to-be”. The MMI method also allows for the mapping of landscape variation in the forest or vegetation density, which is a novel and exciting feature of the new 3D remote sensing (3DRS) technique. Besides the inventory-type applications, the MMI method can be used as a forest change detection method. For maximum effectiveness of the MMI method, an object-based change detection approach is preferred. A minimum requirement for the MMI method is a time-lapsed reference dataset in the form, for example, of an existing forest map of the area of interest, or a vegetation density map prepared using InSAR datasets. Preliminary tests aimed at finding a degree of correlation between the impenetrability and other types of passive and active remote sensing data sources, including TerraSAR-X, NDVI and PALSAR, proved that the method most sensitive to vegetation density was the Japanese PALSAR - L-band SAR system. Unfortunately, PALSAR backscattered signals become very noisy for impenetrability below 15 m. This means that PALSAR has severe limitations for low loadings of the biomass per unit area. The proposed applications of the InSAR data will remain indispensable wherever cloud cover obscures the sky in a persistent manner, which makes suitable optical data acquisition extremely time-consuming or nearly impossible. A limitation of the MMI method is due to the fact that the impenetrability is calculated using a reference DTM, which must be available beforehand. In many countries around the world, appropriate quality DTMs are still unavailable. A possible solution to this obstacle is to use a DEM that was derived using P-band InSAR elevations or LiDAR. It must be noted, however, that in many cases, two InSAR datasets separated by time of the same area are sufficient for forest change detection or similar applications

Monitoring Snow Cover and Snowmelt Dynamics and Assessing their Influences on Inland Water Resources

Snow is one of the most vital cryospheric components owing to its wide coverage as well as its unique physical characteristics. It not only affects the balance of numerous natural systems but also influences various socio-economic activities of human beings. Notably, the importance of snowmelt water to global water resources is outstanding, as millions of populations rely on snowmelt water for daily consumption and agricultural use. Nevertheless, due to the unprecedented temperature rise resulting from the deterioration of climate change, global snow cover extent (SCE) has been shrinking significantly, which endangers the sustainability and availability of inland water resources. Therefore, in order to understand cryo-hydrosphere interactions under a warming climate, (1) monitoring SCE dynamics and snowmelt conditions, (2) tracking the dynamics of snowmelt-influenced waterbodies, and (3) assessing the causal effect of snowmelt conditions on inland water resources are indispensable. However, for each point, there exist many research questions that need to be answered. Consequently, in this thesis, five objectives are proposed accordingly. Objective 1: Reviewing the characteristics of SAR and its interactions with snow, and exploring the trends, difficulties, and opportunities of existing SAR-based SCE mapping studies; Objective 2: Proposing a novel total and wet SCE mapping strategy based on freely accessible SAR imagery with all land cover classes applicability and global transferability; Objective 3: Enhancing total SCE mapping accuracy by fusing SAR- and multi-spectral sensor-based information, and providing total SCE mapping reliability map information; Objective 4: Proposing a cloud-free and illumination-independent inland waterbody dynamics tracking strategy using freely accessible datasets and services; Objective 5: Assessing the influence of snowmelt conditions on inland water resources