Flight Dynamics Operations of the TanDEM-X Formation

Since end of 2010 the German TerraSAR-X and TanDEM-X satellites are routinely operated as the first configurable single-pass Synthetic Aperture Radar interferometer in space. The two 1340 kg satellites fly in a 514 km sun-synchronous orbit. In order to collect sufficient measurements for the generation of a global digital elevation model and to demonstrate new interferometric SAR techniques and applications, more than three years of formation flying are foreseen with flexible baselines ranging from 150 m to few kilometers. As a prerequisite for the close formation flight an extensive flight dynamics system was established at DLR/GSOC, which comprises of GPS-based absolute and relative navigation and impulsive orbit and formation control. Daily formation maintenance maneuvers are performed by TanDEM-X to counterbalance natural and artificial disturbances. The paper elaborates on the routine flight dynamics operations and its interactions with mission planning and ground-station network. The navigation and formation control concepts and the achieved control accuracy are briefly outlined. Furthermore, the paper addresses non-routine operations experienced during formation acquisition, frequent formation reconfiguration, formation maintenance problems and space debris collision avoidance, which is even more challenging than for single-satellite operations. In particular two close approaches of debris are presented, which were experienced in March 2011 and April 2012. Finally, a formation break-up procedure is discussed which could be executed in case of severe onboard failures

Utilization of bistatic TanDEM-X data to derive land cover information

Forests have significance as carbon sink in climate change. Therefore, it is of high importance to track land use changes as well as to estimate the state as carbon sink. This is useful for sustainable forest management, land use planning, carbon modelling, and support to implement international initiatives like REDD+ (Reducing Emissions from Deforestation and Degradation). A combination of field measurements and remote sensing seems most suitable to monitor forests. Radar sensors are considered as high potential due to the weather and daytime independence. TanDEM-X is a interferometric SAR (synthetic aperture radar) mission in space and can be used for land use monitoring as well as estimation of biophysical parameters. TanDEM-X is a X-band system resulting in low penetration depth into the forest canopy. Interferometric information can be useful, whereas the low penetration can be considered as an advantage. The interferometric height is assumable as canopy height, which is correlated with forest biomass. Furthermore, the interferometric coherence is mainly governed by volume decorrelation, whereas temporal decorrelation is minimized. This information can be valuable for quantitative estimations and land use monitoring. The interferometric coherence improved results in comparison to land use classifications without coherence of about 10% (75% vs. 85%). Especially the differentiation between forest classes profited from coherence. The coherence correlated with aboveground biomass in a R² of about 0.5 and resulted in a root mean square error (RSME) of 14%. The interferometric height achieved an even higher correlation with the biomass (R²=0.68) resulting in cross-validated RMSE of 7.5%. These results indicated that TanDEM-X can be considered as valuable and consistent data source for forest monitoring. Especially interferometric information seemed suitable for biomass estimation

In-depth verification of Sentinel-1 and TerraSAR-X geolocation accuracy using the Australian Corner Reflector Array

This article shows how the array of corner reflectors (CRs) in Queensland, Australia, together with highly accurate geodetic synthetic aperture radar (SAR) techniques—also called imaging geodesy—can be used to measure the absolute and relative geometric fidelity of SAR missions. We describe, in detail, the end-to-end methodology and apply it to TerraSAR-X Stripmap (SM) and ScanSAR (SC) data and to Sentinel-1interferometric wide swath (IW) data. Geometric distortions within images that are caused by commonly used SAR processor approximations are explained, and we show how to correct them during postprocessing. Our results, supported by the analysis of 140 images across the different SAR modes and using the 40 reflectors of the array, confirm our methodology and achieve the limits predicted by theory for both Sentinel-1 and TerraSAR-X. After our corrections, the Sentinel-1 residual errors are 6 cm in range and 26 cm in azimuth, including all error sources. The findings are confirmed by the mutual independent processing carried out at University of Zurich (UZH) and German Aerospace Center (DLR). This represents an improve�ment of the geolocation accuracy by approximately a factor of four in range and a factor of two in azimuth compared with the standard Sentinel-1 products. The TerraSAR-X results are even better. The achieved geolocation accuracy now approaches that of the global navigation satellite system (GNSS)-based survey of the CRs positions, which highlights the potential of the end-to-end SAR methodology for imaging geodesy

IMPROVING AND EXPANDING PRECISION ORBIT DERIVED ATMOSPHERIC DENSITIES

Atmospheric drag is the most uncertain non-conservative force acting on a low Earth orbiting satellite. The existing atmospheric density models are not accurate enough to model the variations in density, which significantly affect the drag on satellites since drag is directly proportional to atmospheric density. In this research, precision orbit ephemerides (POE) are used as measurements in an optimal orbit determination scheme to estimate corrections to baseline atmospheric density models. These corrections improve the drag estimates, which in turn improve orbit determination and prediction and also provide a better understanding of the upper atmosphere. The POE are used as measurements in a sequential measurement and filtering scheme using the Orbit Determination Tool Kit (ODTK) software, which provides the orbit determination. Five atmospheric density models are available in ODTK, which are used as baseline atmospheric density models to which corrections are made in the orbit determination. These density models are Jacchia 1971, Jacchia-Roberts, CIRA 1972, MSISE 1990, and NRLMSISE 2000. The user has the option to specify the ballistic coefficient (BC) correlated half-life and density correlated half-life. These half-lives are usually given values of 1.8, 18, or 180 minutes. If all five baseline density models are used along with three different combinations of ballistic coefficient and density correlated half-lives, then this would result in forty-five different cases. All the forty-five cases are examined in some studies and only a selected few are examined in others, the details of which are given in the appropriate sections. The POE derived densities are validated by comparing them with accelerometer derived densities for satellites which have accelerometers onboard, such as the Challenging Minisatellite Payload (CHAMP) and the Gravity Recovery and Climate Experiment (GRACE). The trend in the variation is compared quantitatively by calculating the cross correlation between the POE and accelerometer derived densities, and the magnitude is compared by calculating the root mean square between the two. The accelerometer derived densities for both CHAMP and GRACE are available from Sean Bruinsma of CNES and also from Eric Sutton of the United States Air Force Research Laboratory, and are used in this research. The effect of different functions of geomagnetic planetary amplitude (ap) as an input in orbit determination to estimate atmospheric density was investigated. The three different functions of input are 3-hourly ap step functions, linear interpolated ap functions, and ap osculating spline functions. These three different types of functions were used as inputs for all the forty-five different combinations obtained by using the five different baseline atmospheric density models and three different combinations of ballistic coefficient and density correlated half-lives as stated earlier, and POE derived density was estimated for both CHAMP and GRACE. The POE derived densities were compared with the accelerometer derived densities by calculating the CC and RMS. To create continuous data sets of POE derived densities that span a period of one week, the linear weighted blending technique was used to blend the 14 hour POE derived densities in their overlap periods. CIRA 1972 was used as the baseline atmospheric density model and a BC correlated half-life of 1.8 minutes and density correlated half-life of 180 minutes were used as inputs in ODTK to generate these POE derived density estimates. These one week continuous POE derived densities showed better correlation with accelerometer derived densities than HASDM densities for both CHAMP and GRACE. The average cross-sectional area of the satellite that is normal to the velocity vector, the area facing the Sun, and the area facing the Earth, were determined so that these areas could be used to estimate the atmospheric drag, the force due to solar radiation pressure, and the force due to Earth radiation pressure (infrared and Earth albedo). This was done for both TerraSAR-X and ICESat. For TerraSAR-X, the area normal to the velocity vector was assumed be a constant and equal to the frontal area, and the area facing the Earth was also assumed to be constant. However, the area facing the Sun varied with time. The average area facing the Sun for a period of 14 hours and also the annual average area were calculated and used to calculate the POE derived densities. The POE derived densities calculated using these two different average areas facing the Sun were found to be very similar. Since TerraSAR-X does not have an accelerometer onboard, the POE derived densities could not be compared with accelerometer derived densities, but instead were compared with Jacchia-71 densities since this was also one of the outputs from ODTK. The POE derived densities were also compared with NRLMSISE 2000 densities. The attitude of ICESat as a function of beta angle was given in the literature and so was the average area of each side of the satellite when it was modeled as a rectangular box with two solar panels. This information was used to estimate the 30-hour average area normal to the velocity vector, area facing the Earth, and area facing the Sun, for ICESat. The POE derived densities using these areas were estimated by ODTK and compared with the Jacchai-71 density model

One Decade of Glacier Mass Changes on the Tibetan Plateau Derived from Multisensoral Remote Sensing Data

The Tibetan Plateau (TP) with an average altitude of 4,500 meters above sea level is characterized by many glaciers and ice caps. Glaciers are a natural indicator for climate variability in this high mountain environment where meteorological stations are rare or non-existent. In addition, the melt water released from the Tibetan glaciers is feeding the headwaters of the major Asian river systems and contributes to the rising levels of endorheic lakes on the plateau. As many people directly rely on the glacier melt water a continuous glacier monitoring program is necessary in this region. In situ measurements of glaciers are important, but are spatial limited due to large logistical efforts, physical constrains and high costs. Remote sensing techniques can overcome this gap and are suitable to complement in situ measurements on a larger scale. In the last decade several remote sensing studies dealt with areal changes of glaciers on the TP. However, glacier area changes only provide a delayed signal to a changing climate and the amount of melt water released from the glaciers cannot be quantified. Therefore it is important to measure the glacier mass balance. In order to estimate glacier mass balances and their spatial differences on the TP, several remote sensing techniques and sensors were synthesized in this thesis. In a first study data from the Ice Cloud and Elevation Satellite (ICESat) mission were employed. ICESat was in orbit between 2003 and 2009 and carried a laser altimeter which recorded highly accurate surface elevation measurements. As in mid-latitudes these measurements are rather sparse glaciers on the TP were grouped into eight climatological homogeneous sub-regions in order to perform a statistical sound analysis of glacier elevation changes. To assess surface elevation changes of a single mountain glacier from ICESat data, an adequate spatial sampling of ICESat measurements need to be present. This is the case for the Grosser Aletschgletscher, located in the Swiss Alps which served as a test site in this thesis. In another study data from the current TanDEM-X satellite mission and from the Shuttle Radar Topography Mission (SRTM) conducted in February 2000 were employed to calculate glacier elevation changes. In a co-authored study, these estimates could be compared with glacier elevation changes obtained from the current French Pléiades satellite mission. In order to calculate glacier mass balances, the derived elevation changes were combined with assumptions about glacier area and ice density in all studies. In this thesis contrasting patterns of glacier mass changes were found on the TP. With an ICESat derived estimate of -15.6±10.1 Gt/a between 2003 and 2009 the average glacier mass balance on the TP was clearly negative. However, some glaciers in the central and north-western part of the TP showed a neutral mass balance or a slightly positive anomaly which was also confirmed by data from the current TanDEM-X satellite mission. A possible explanation of this anomaly in mass balance could be a compensation of the temperature driven glacier melt due to an increase in precipitation

Monitoring active open-pit mine stability in the Rhenish coalfields of Germany using a coherence-based SBAS method

With the recent progress in synthetic aperture radar (SAR) technology, especially the new generation of SAR satellites (Sentinel-1 and TerraSAR-X), our ability to assess slope stability in open-pit mines has significantly improved. The main objective of this work is to map ground displacement and slope instability over three open-pit mines, namely, Hambach, Garzweiler and Inden, in the Rhenish coalfields of Germany to provide long-term monitoring solutions for open-pit mining operations and their surroundings. Three SAR datasets, including Sentinel-1A data in ascending and descending orbits and TerraSAR-X data in a descending orbit, were processed by a modified small baseline subset (SBAS) algorithm, called coherence-based SBAS, to retrieve ground displacement related to the three open-pit mines and their surroundings. Despite the continuously changing topography over these active open-pit mines, the small perpendicular baselines of both Sentinel-1A and TerraSAR-X data were not affected by DEM errors and hence could yield accurate estimates of surface displacement. Significant land subsidence was observed over reclaimed areas, with rates exceeding 500 mm/yr, 380 mm/yr, and 310 mm/yr for the Hambach, Garzweiler and Inden mine, respectively. The compaction process of waste materials is the main contributor to land subsidence. Land uplift was found over the areas near the active working parts of the mines, which was probably due to excavation activities. Horizontal displacement retrieved from the combination of ascending and descending data was analysed, revealing an eastward movement with a maximum rate of ∼120 mm/yr on the western flank and a westward movement with a maximum rate of ∼ 60 mm/yr on the eastern flank of the pit. Former open-pit mines Fortuna-Garsdorf and Berghein in the eastern part of Rhenish coalfields, already reclaimed for agriculture, also show subsidence, at locations reaching 150 mm/yr. The interferometric results were compared, whenever possible, with groundwater information to analyse the possible reasons for ground deformation over the mines and their surroundings

Satellite Formation-Flying and Rendezvous

GNSS has come to play an increasingly important role in satellite formation-flying and rendezvous applications. In the last decades, the use of GNSS measurements has provided the primary technique for determining the relative position of cooperative co-orbiting satellites in low Earth orbit

Assessment of high resolution SAR imagery for mapping floodplain water bodies: a comparison between Radarsat-2 and TerraSAR-X

Flooding is a world-wide problem that is considered as one of the most devastating natural hazards. New commercially available high spatial resolution Synthetic Aperture RADAR satellite imagery provides new potential for flood mapping. This research provides a quantitative assessment of high spatial resolution RADASAT-2 and TerraSAR-X products for mapping water bodies in order to help validate products that can be used to assist flood disaster management. An area near Dhaka in Bangladesh is used as a test site because of the large number of water bodies of different sizes and its history of frequent flooding associated with annual monsoon rainfall. Sample water bodies were delineated in the field using kinematic differential GPS to train and test automatic methods for water body mapping. SAR sensors products were acquired concurrently with the field visits; imagery were acquired with similar polarization, look direction and incidence angle in an experimental design to evaluate which has best accuracy for mapping flood water extent. A methodology for mapping water areas from non-water areas was developed based on radar backscatter texture analysis. Texture filters, based on Haralick occurrence and co-occurrence measures, were compared and images classified using supervised, unsupervised and contextual classifiers. The evaluation of image products is based on an accuracy assessment of error matrix method using randomly selected ground truth data. An accuracy comparison was performed between classified images of both TerraSAR-X and Radarsat-2 sensors in order to identify any differences in mapping floods. Results were validated using information from field inspections conducted in good conditions in February 2009, and applying a model-assisted difference estimator for estimating flood area to derive Confidence Interval (CI) statistics at the 95% Confidence Level (CL) for the area mapped as water. For Radarsat-2 Ultrafine, TerraSAR-X Stripmap and Spotlight imagery, overall classification accuracy was greater than 93%. Results demonstrate that small water bodies down to areas as small as 150m² can be identified routinely from 3 metre resolution SAR imagery. The results further showed that TerraSAR-X stripmap and spotlight images have better overall accuracy than RADARSAT-2 ultrafine beam modes images. The expected benefits of the research will be to improve the provision of data to assess flood risk and vulnerability, thus assisting in disaster management and post-flood recovery

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

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