Description: In permafrost areas, seasonal freeze-thaw cycles of active layer result in upward and downward movements of the ground. Additionally, relatively uniform thawing of the ice-rich layer at the permafrost table can contribute to net long-term surface lowering. We use a simple method to quantify surface lowering (subsidence) and uplift in a yedoma area of the Lena River Delta, Siberian Arctic (Kurungnakh Island), using reference rods (metal pipes and fiberglass rods) installed deeply in permafrost. The metal pipes were 2 m long and 3 cm in diameter and were anchored at least 1 m below the typical active layer. The fiberglass rods were 2 m long and 1 cm in diameter and were anchored at least 70 m below the typical active layer. We assume, therefore, that the rods were motionless relative to the permafrost. The plexiglass plate with a size of 10 by 10 cm was fixed in its horizontal position by the rod but could move freely with the surface vertically along the rod. We repeatedly measured distance between the top of a rod and a plexiglass plate resting on the ground surface. Several distance measurements around each rod were taken at each visit and averaged. Altogether 12 metal pipes were installed at the study site in April 2013 and 19 fiberglass rods were installed in April 2014. Measurements were conducted during field campaigns from spring 2013 to summer 2017 with some gaps. We provide here the measured distances between the top of a rod and a plexiglass plate. To obtain the ground displacement, the user have to define the period of interest and calculate the displacement.

Description: In permafrost areas, seasonal freeze-thaw cycles of active layer result in upward and downward movements of the ground. Additionally, relatively uniform thawing of the ice-rich layer at the permafrost table can contribute to net long-term surface lowering. We use a simple method to quantify surface lowering (subsidence) and uplift in a yedoma area of the Lena River Delta, Siberian Arctic (Kurungnakh Island), using reference rods (metal pipes and fiberglass rods) installed deeply in permafrost. The metal pipes were 2 m long and 3 cm in diameter and were anchored at least 1 m below the typical active layer. The fiberglass rods were 2 m long and 1 cm in diameter and were anchored at least 70 m below the typical active layer. We assume, therefore, that the rods were motionless relative to the permafrost. The plexiglass plate with a size of 10 by 10 cm was fixed in its horizontal position by the rod but could move freely with the surface vertically along the rod. We repeatedly measured distance between the top of a rod and a plexiglass plate resting on the ground surface. Several distance measurements around each rod were taken at each visit and averaged. Altogether 12 metal pipes were installed at the study site in April 2013 and 19 fiberglass rods were installed in April 2014. Measurements were conducted during field campaigns from spring 2013 to summer 2017 with some gaps. We provide here the measured distances between the top of a rod and a plexiglass plate. To obtain the ground displacement, the user have to define the period of interest and calculate the displacement.

Description: Differential SAR interferometry (DInSAR) uses the phase difference between two SAR signals acquired on two dates over the same area to measure small-scale ground motion. During the last decade the method has been adapted for monitoring permafrost-related ground motion. Here we perform DInSAR on TerraSAR-X data to assess its viability for seasonal thaw subsidence detection in a yedoma landscape of the Lena River Delta. TerraSAR-X is a right-looking SAR satellite launched in 2007, operating in the X-band (wavelength 3.1 cm, frequency 9.6 GHz), with a revisit time of eleven days. All data that we used were acquired in StripMap mode with HH polarization from a descending orbit at 08:34 local acquisition time (22:34 UTC). The incidence angle of the track we use is approximately 31 degrees. The scene size covered an area of approximately 18 km x 56 km. The slant range and azimuth pixel spacing were approximately 0.9 m and 2.4 m, respectively. Based on the ground temperature data we roughly estimated the beginning and the end of thaw season in 2013. The corresponding TerraSAR-X time series used for this study includes nine Single-Look Slant Range Complex (SSC) images taken between 7 June and 14 September 2013. The time span between the acquisitions that we used for interferometry was 11 days, with one exception when the time span was 22 days due to a missing acquisition. The data were processed using the Gamma radar software. The SSC data were converted to Gamma Single Look Complex (SLC) format and the SLC data were then consecutively co-registered with subpixel accuracy (typically better than 0.2 pixels) in such a way that the co-registered slave image became the master for the next image. This way of co-registering also ensures subpixel co-registration accuracy for all interferometric combinations of the nine images. Multilooking was performed with the factor 4 in the range and factor 3 in the azimuth directions to reduce the noise and obtain roughly square ground range pixels. The ground size of the multilooked pixel is approximately 7 m. We removed the topographic phase term using ArcticDEM that is a freely available high-resolution (5 m) circum-Arctic DEM produced from optical stereographic WorldView imagery acquired from 2012–2016. Obtained differential interferograms were then filtered with an adaptive filter based on the local fringe spectrum with the filtering window size of 128 pixels and an alpha exponent of 0.4. Interferograms, featuring especially low coherence, were additionally filtered with a window size of 64 pixels. For the phase unwrapping we used a branch-cut algorithm with the seeding point located approximately in the middle of the study area with relatively high coherence. We did not attempt to unwrap the areas, separated from the main study area by the river channels. The influence of atmospheric phase delays was evident in the unwarpped interferograms. In order to enhance the displacement signal and reduce atmospheric noise, all eight unwrapped interferograms were summed up in a time-continuous stack. Phase rate per day was calculated from the stack. A strong linear ramp was present across the phase rate map. To remove the trend, a 2D linear function was fit to the data and then subtracted from the phase rate map. The phase rate was then converted to vertical displacement rate in meters, under the assumption that the ground movement is purely vertical. The resulting displacement rate map was geocoded using ArcticDEM to the Universal Transverse Mercator (UTM) projection, zone 52N WGS84 with a pixel size of 5 m. The map was finally converted to the displacement magnitude by multiplying the rate by 99 days (from 7 June to 14 September 2013) and converted to centimeters. As opposed to the results, published in the related paper, here we did not start the unwrapping from the known bedrock position, as it was partly affected by low coherence as well as rather remote from the main area of interest and only weakly connected to the rest of the map over a small and noisy area of valid pixels. It means that the displacement map published here, features only displacement values relative to each other, without a fixed reference point. The spatial pattern of the signal, however, did not change with this alteration in processing. The DInSAR map showed a distinct subsidence in most of the thermokarst basins relative to the upland. Moreover, the spatial pattern of DInSAR signal was in high agreement with the surface wetness in the basins, identified with the near infra-red band of a high-resolution optical image. Drier parts of the basins were clearly separated from wetter parts that showed a prominent subsidence. In general, low coherence in combination with atmospheric effects as well as remoteness of a reference ground point were severe obstacles for the retrieval of a wide-area seasonal thaw subsidence map with TerraSAR-X data.

Description: In permafrost areas, seasonal freeze-thaw cycles result in upward and downward movements of the ground. For some permafrost areas, long-term downward movements were reported during the last decade. We measured seasonal and multi-year ground movements in a yedoma region of the Lena River Delta, Siberia, in 2013–2017, using reference rods installed deep in the permafrost. The seasonal subsidence was 1.7 ± 1.5 cm in the cold summer of 2013 and 4.8 ± 2 cm in the warm summer of 2014. Furthermore, we measured a pronounced multi-year net subsidence of 9.3 ± 5.7 cm from spring 2013 to the end of summer 2017. Importantly, we observed a high spatial variability of subsidence of up to 6 cm across a sub-meter horizontal scale. In summer 2013, we accompanied our field measurements with Differential Synthetic Aperture Radar Interferometry (DInSAR) on repeat-pass TerraSAR-X (TSX) data from the summer of 2013 to detect summer thaw subsidence over the same study area. Interferometry was strongly affected by a fast phase coherence loss, atmospheric artifacts, and possibly the choice of reference point. A cumulative ground movement map, built from a continuous interferogram stack, did not reveal a subsidence on the upland but showed a distinct subsidence of up to 2 cm in most of the thermokarst basins. There, the spatial pattern of DInSAR-measured subsidence corresponded well with relative surface wetness identified with the near infra-red band of a high-resolution optical image. Our study suggests that (i) although X-band SAR has serious limitations for ground movement monitoring in permafrost landscapes, it can provide valuable information for specific environments like thermokarst basins, and (ii) due to the high sub-pixel spatial variability of ground movements, a validation scheme needs to be developed and implemented for future DInSAR studies in permafrost environments.

Description: Knowledge about the coverage and characteristics of glaciers in High Mountain Asia is still incomplete and heterogeneous. However, several applications such as modelling of past or future glacier development, runoff or glacier volume, rely on the existence and accessibility of complete datasets. In particular, precise outlines of glacier extent are required to spatially constrain glacier-specific calculations such as length, area and volume changes or flow velocities. As a contribution to the Randolph Glacier Inventory (RGI) and the Global Land Ice Measurements from Space (GLIMS) glacier database, we have produced a homogeneous inventory of the Pamir and the Karakoram mountain ranges using 28 Landsat TM and ETM+ scenes acquired around the year 2000. We applied a standardized method of automated digital glacier mapping and manual correction using coherence images from ALOS-1 PALSAR-1 as an additional source of information; we then separated the glacier complexes into individual glaciers using drainage divides derived by watershed analysis from the ASTER GDEM2, and separately delineated all debris-covered areas. Assessment of uncertainties was performed for debris-covered and clean-ice glacier parts using the buffer method and independent multiple digitizing of three glaciers representing key challenges such as shadows and debris cover. Indeed, along with seasonal snow at high elevations, shadow and debris cover represent the largest uncertainties in our final dataset. In total, we mapped more than 27'800 glaciers 〉0.02 km² covering an area of 35'520 ±1948 km² and an elevation range from 2260 m to 8600 m. Regional median glacier elevations vary from 4150 m (Pamir Alai) to almost 5400 m (Karakoram), which is largely due to differences in temperature and precipitation. Supraglacial debris covers an area of 3587 ±662 km², i.e. 10% of the total glacierised area. Larger glaciers have a higher share in debris-covered area (up to 〉20%), making it an important factor to be considered in subsequent applications.