更全的杂志信息网

Surface-deformation monitoring in the permafrost regions over the Tibetan Plateau, using Sentinel-1 data

更新时间:2016-07-05

1 Introduction

A variety of traditional methods has been used to figure out the seasonal or long term changes of the active layer, for instance, mechanical probing (Brown et al., 2000; Cheng and Wu, 2007; Zorigt et al., 2016)and frost or thaw tubes (Harris, 1970; Warren and Brown, 1972; Mackay, 1973a; Mark Nixon and Taylor, 1998). Data acquired from temperature sensors can also be used to assess differences in thaw depth during freezing and thawing periods (Zhang et al., 1997; Brown et al., 2000; Wu et al., 2010). With the development of the technology in geophysical survey in the past decades, time-domain reflectometry(TDR) and GPS were also applied to determine the variation of the active layer (Baker et al., 1982; Little et al., 2003). However, some of these traditional measurements, such as drilling boreholes in wintertime, are not easy to obtain due to the tough fieldwork conditions on the Tibetan Plateau, with high altitude and cold weather. Furthermore, the greatest disadvantage of all these in-situ methods is that only sparse, point-based observations can be supplied,though these observational results have high accuracy.Therefore, mapping permafrost displacement over a broad area by D-InSAR should be an adaptive way for determining seasonal and long-term surface changes of permafrost.

Advanced InSAR methods such as Permanent Scatters InSAR (PS-InSAR) and SBAS-InSAR have been proposed for detecting surface change using D-InSAR technology (Ferretti et al., 2000; Usai,2001; Berardino et al., 2002). SBAS-InSAR has already been applied to retrieve long-term surface-deformation data over cold regions or mine areas, where permanent scatters are difficult to find (Castañeda et al., 2009; Liu et al., 2010; Liu et al., 2012; Chen et al., 2013; Short et al., 2014). Furthermore, SBAS-In-SAR can significantly reduce the de-correlation induced by a long spatial and temporal baseline because the least-square method is incorporated into this approach for improving temporal resolution (Berardino et al., 2002). These specific advantages show a potential application of the SBAS-InSAR method for deformation measurement over permafrost regions. Previous studies of deformation measurements with D-InSAR technology have been conducted in some cold regions (Li et al., 2004; Liu et al., 2010,2012; Daout et al., 2017), and C-band and L-band data also have been utilized for detecting the surface changes of permafrost over the Tibetan Plateau (Xie et al., 2010; Chen et al., 2013; Daout et al., 2017).However, baseline or temporal de-correlations still have great adverse impacts on D-InSAR results, due to large baselines between acquisitions and vegetation growth or snowfall; and radar results also are significantly affected by phase delay, which is in turn caused by atmospheric factors such as humidity(Zebker et al., 1997). Although some previous studies have inverted seasonal and long-term surface variation of the active layer on the Tibetan Plateau, some of these studies still lack the ground truth. In addition,few studies use Sentinel-1 data for inverting time-series displacements of ground surfaces in permafrost regions of the Tibetan Plateau and validating its applicability.

In this study, using the SBAS-InSAR method, we acquired and analyzed the seasonal and long-term deformations of the TGL and LDH regions during two freeze thaw cycles, from October 2014 to June 2016,with Sentinel-1 images. The objective of this study was to estimate the applicability of Sentinel-1 C-band SAR data over specific permafrost regions under different vegetation coverage and soil-moisture regimes by using in-situ leveling observations.

2 Study area

TGL (91.8°E–92.0°E, 33.0°N–33.2°N) is located in the interior of a continuous permafrost zone on the Tibetan Plateau (Figure 1), with a high altitude of around 5,100 m a.s.l.. This region is covered by degraded alpine meadow, with plant stature less than 10 cm and coverage of 20%–30%. The mean annual air temperature is approximately −4.9 °C, and the extreme maximum and the minimum values of daily average temperature are 17.6 °C and −29.6 °C, respectively (Jiao et al., 2014). Precipitation primarily occurs in summer, between May and September, with an annual precipitation of about 400 mm; and the maximum precipitation occurs in July and August. The active-layer thickness of the in-situ leveling observational field in the TGL region is about 3.15 m. The soils are composed of gravelly sandy loam, including around 5%–30% gravel, 65%–95% sand, and 3%–22% silt, as well as 2%–13% clay (Hu et al.,2015).

LDH (91.5°E–91.85°E, 31.6°N–32.1°N) is located on the south margin of a continuous permafrost zone of the Tibetan Plateau (Figure 1). Vegetation coverage, soil moisture, and topography in this region significantly differ from those in the TGL region.Vegetation type at LDH is alpine paludal meadow,with a vegetation coverage of more than 80%. The average altitude of this area is above 4,500 m a.s.l., the mean annual air temperature is around −1.3 °C, and the mean annual precipitation (from 1971 to 2004) is about 411 mm (Yang et al., 2012). The active-layer thickness is around 1 m. Soil in this area is much wetter than that at TGL, and a pure-ice layer exists at the bottom of the active layer.

Figure 1 Position of the study area

3 Data sources

3.1 Sentinel-1 C-band data and DEM

The Sentinel-1 program was established by the European Space Agency to replace the older Earth-observation mission of the European remote-sensing satellite (ERS) and ensure continuous data acquisition.The Sentinel-1 mission includes two polar-orbiting satellites with an altitude of 693 km, which means they can comprise a constellation and improve the observation efficiency. These two satellites both carry a C-band synthetic aperture radar (SAR) device that can operate day and night, and it also can acquire all-weather images (European Space Agency, 2013).Two sets of 27 and 42 Sentine-1A/B data acquired from October 2014 to June 2017 and covering TGL and LDH were utilized in this study; see the list in Tables 1 and 2, respectively. All acquisition images used in this study were collected in Interferometric Wide Swath (IW) mode, with a swath width of around 250 km and spatial resolution of 5m×20m.

Table 1 Sentinel-1 acquisitions in TGL

?

Table 2 Sentinel-1 acquisitions in LDH

?

Precise orbit files of remotely sensed data, which are necessary in interferometry processing were obtained from the Payload Data Ground Segment(PDGS). The Shuttle Radar Topography Mission DEM (SRTM DEM), with a spatial resolution of 30 m and absolute height error of less than 16 m (Farr et al.,2007), was also used in this study because the indispensable DEM is used for topography subtraction during differential interferometry processes.

3.2 In-situ leveling data

A leveling observational field was arranged at TGL in 2010, consisting of one leveling datum point and 25 benchmarks evenly distributed over an area of 16 m2. The leveling-measurement datum point is an iron tube, inserted into a 60-m-deep borehole, with about 15 cm protruding aboveground. The leveling observing site was also established at the LDH region in October 2015, including 16 benchmarks with an area of around 500 m2 and one datum point (same as the datum point in TGL). We hypothesized that the leveling datum points of the two observation sites kept stable during each freeze–thaw cycle because the iron tubes were deeply inserted into the ground(Mackay, 1973b, 1977). These leveling marks were measured every month or two with a SOKKIA SDL30 digital level, which performs well in adverse environments, with a height-observing accuracy of 1.0 mm (standard deviation for 1-km double-run leveling) using random bidirectional (RAB) code staves.

4 Methods

The SBAS method assumes an unwrapping phase of each interferogram as observation and then acquires the cumulative phase by solving a linear leastsquares problem.

整个这节课的安排实在是太妙了!不仅在潜移默化中渗透了小说的三要素,而且还扎扎实实培养了学生该如何阅读小说这一能力,对于学生今后的自主阅读小说提供了很好的思路。

Considering that N+1 images at the same orbit are collected in chronological sequence of t0, t1, …, tn,and hypothesizing that each image could serve as a master image and produce at least one interferogram with other acquisitions, total M interferograms are produced during the interfering process; and they satisfy the following equation:

将刚制得的高盐稀态酱醪进行淋油并充分搅拌,以日式高盐稀态酱油的发酵工艺进行发酵,先将其混匀并分装至9个容器中,置于15 ℃的生化培养箱中进行酱油前期发酵,每天搅拌1次至第15 d,按表1添加 2×106个/g(酱醪)酵母,随后每天升高发酵温度1 ℃,直至30 ℃并保持,持续发酵到第45 d;再按表1添加酵母2×106个/g(酱醪),以9号样品为空白对照;发酵约90 d时,控制酱醪发酵温度于35 ℃到发酵结束,整个发酵周期为180 d[23]。

Assuming that the j-th interferogram is calculated by acquisition at the time of tA and tB, so the azimuth and range coordinates of x and r can be described as δØj(x,r)=ØB(x,r)−ØA(x,r)where λ is the transmitted-signal central wavelength and j∈(1, …, M); d(tB, x, r) and d(tA, x, r) mean the cumulative deformation in the direction of line of sight (LOS) at times tA and tB, relating to the value of d(t0, x, r)=0.

1.2.1 问卷调查 问卷调查名称为《甲状腺疾病和糖尿病全国调查-2014(TIDE)》,问卷调查内容包括研究对象的一般资料(姓名、年龄、民族、文化程度、职业、收入、生育史、吸烟、饮酒史)、碘营养状况(食盐来源、进食盐习惯、进食海带、紫菜情况、是否服用含碘药物、3个月内接受造影剂检查情况)、疾病史及家族史(甲状腺病史、甲状腺疾病家族史、糖尿病病史、糖尿病家族史、糖尿病并发症、高血压病史、高脂血症病史、高尿酸血症病史)。采用整体分层抽样调查法,发放问卷1 464份,收集问卷1 430份,有效回收率为97.68%。

“Objective”表明该测试是一种客观评量,体现在各站测试目标明确,得分标准清楚;“Structured”意为测试的各站内容经过详尽的规划、非随机的测验,具有特定的目的性(即考查临床技能的某一特定部分,设计考试一定要有蓝图,见表1);“Clinical”意为测验的是临床技能,而且使用真实的临床案例;“Examination”包含3层涵义,该测试同时可作为回馈性测验、总结性测验和执照考试。据此,OSCE可运用于教学、毕业生或新进人员评鉴及执照考试。

represents the topography phase, most of which can be removed by an available external DEM; therefore, the residual topography phase can be neglected during the interfering process because of itsindicates de-correlation noise, such as system noise(caused by system thermal and quantisation). Equation (2) can be simplified as below without taking

For the purpose of obtaining the cumulative deformation with physical meaning, the phase in Equation (3) can be expressed as the mean phase velocity multiplied by the acquisition-time difference:

Therefore, the phase value of the j-th interferogram can be given by

鄯善联合站25 MPa注水系统(温西六区域)的运行注水量为9.17 m3/h,注水站出口压力27.2 MPa,额定流量20.2 m3/h,额定压力是25 MPa,工频运行。通过运行数据与额定参数对比,注水泵效率行低的原因是实际注水量远低于额定排量,处于“大马拉小车”现象,同时没有采取任何节能措施,工频运行。考虑实际注水情况,建议对鄯善联合站25 MPa注水系统(温西六区域)B12#注水泵安装变频控制柜,降低运行频率,调节电动机转速,达到合理调配流量作用,提高注水泵机组效率;或者更换注水泵运行,使用小的额定排量注水泵,如额定排量为13.9 m3/h注水泵。

Figure 6 demonstrates the cumulative deformation of TGL during the two-year measurement period,i.e., the total ground displacement on November 3,2016 (Figure 6b) relative to the surface height on October 27, 2014 (Figure 6a). The study area of TGL is included in the black rectangle, and deformation of these two maps is in the satellite's LOS direction, as measured in mm. A positive value means upwards from the land surface moving towards the satellite;negative represents land surface moving away from the satellite (subsidence). Some pixels appear as white in the figures because the coherence of these areas in the interferograms is lower than the coherence threshold. Therefore, these pixels were not included in calculating the deformation during the SBAS process.Regions that surround the study area appear as either subsidence or uplift, which means ground-deformation change is heterogeneous on the Tibetan Plateau.As evident from Figure 6, the region of interest underwent a small subsidence of around 5 mm during the monitoring period from late October 2014 to early November 2016.

where B is always a rank-deficiency matrix with a size of M×N because many different images are chosen as the master image during the computing process. By contrast, a generalized matrix of B can be obtained after applying the singular-value decomposition (SVD) method, and then the velocity (v) can be calculated by employing the minimum-norm solution. Finally, deformation of various periods can be achieved according to the integration of phase velocity at the same time.

It is difficult to analyze the temporal variations in ground-surface deformation in detail due to the sparse in-situ measurements during the monitoring period;and thus the displacement observed by D-InSAR was analyzed in the following sections. Generally, the freeze–thaw processes over permafrost regions on the Tibetan Plateau can be divided into four stages: thawing settlement in summer, frost heave in autumn,cooling process in winter, and warming process in spring (Zhao et al., 2000). However, we found some anomalous deformations inconsistent with those stages were presented during the monitoring period(Figure 6). For example, the uplift of land surface occurred from May to mid-July 2015 and between June and October 2016, which can be observed from In-SAR and the terrestrial leveling method, respectively.These anomalies can be explained by the rainfall or meltwater infiltration in summer. Seasonal ground thaw commonly occurs from early May to late October each year (lasting around 110 days), and the greatest concentration of precipitation falls at this region during the same period, accounting for about 83% of the annual precipitation (Jiao et al., 2014).The soil type of the TGL region is gravelly sandy loam with large porosity, and thus the percolating rainwater migrates downwards to the freezing activelayer soils beneath the thawing surfaces during the thawing period, subsequently increasing the ice content of the frozen active layer and hence the frost heaving of the ground surface. These phenomena have been verified in laboratory and field experiments at various permafrost regions, such as on the Tibetan Plateau and along the western Arctic coast(Cheng, 1982; Mackay, 1983). Water infiltration during the thawing period may have great impact on hydrological processes and on latent and sensible changes over permafrost regions. Therefore, these implications should be considered in the studies on hydrological and energy exchanges, by relating these processes with land-surface displacement.

Figure 2 Flow chart of data processing

In the TGL region, simple linear regressions were applied on InSAR-derived and ground-based measurement results for detecting long-term variation in the active-layer thickness. Both the InSAR results and leveling measurements showed a decreasing trend during the study period, with a larger slope of D-In-SAR (Figure 5). The main reason for this trend should be the low sampling frequency of the terrestrial leveling observations, which may cause the loss of some larger values of deformation during the observing period. For instance, surface height was not recorded between 18 April 2015 and 30 October 2015 or between 31 October 2016 or 8 March 2017 by leveling measurement. Another possible explanation is that D-InSAR results included some residual phase during the interfering process due to DEM errors or atmospheric delay, and this question should be solved in future studies.

新的核磁共振技术,无一例外地面临着诸多挑战,而要成功地将这些技术应用于临床,仍有大量工作要做。DTI的研究已经产生了最实质性的成果,但在这项技术被广泛应用于临床之前,仍然需要进一步的改进,然后进行标准化。目前,国内外各学者对 DTI相关参数在病理基础表现上持有不同的观点,在预后判定方面还未做系统的研究。有关ADC值、FA值变化的病理机制、临床意义及与临床预后的研判方面的研究将是在未来的几年里研究的重点与热点。脊髓成像界无疑会取得许多伟大的成就,未来的研究转化将使临床从中受益。

Figure 3 The spatial temporal baselines of the acquisitions in TGL (black-filled circles and black lines mean the radar acquisitions and interferograms, respectively)

Figure 4 The spatial temporal baselines of the acquisitions in LDH (black-filled circles and black lines mean the radar acquisitions and interferograms, respectively)

During SBAS processing, the threshold of coherence and minimum number of interferograms were set as 0.23 and four separately in TGL, which means that pixels would be included in the time-series processing with a coherence higher than 0.23 in more than four interferograms. In the LDH region, pixels'coherence were to be no lower than 0.2 in no fewer than 50 interferograms during the time-series calculating process. The lower coherence threshold was employed for avoiding seriously incoherent results due to LDH's high vegetation coverage and abundant soilwater content; more interfering results restraint was used because a large number of interferograms were produced during D-InSAR process.

式中,p(θ)为单个粒子散射相位函数,r为粒子半径,Qsca为散射效率因子,n(r)为粒子谱分布函数.

5 Results

In Figure 5, we compared the results observed from the SBAS-InSAR method with terrestrial leveling measurements at TGL. The deformation obtained from the two methods showed a consistent trend and significant seasonal variations during the study period. As observed by D-InSAR and terrestrial leveling,subsidence of around 5 mm occurred from late October to December 2014, though the soils should have been uplifted due to the effects of frost heave during that period. Similar phenomenon can also be identified during the same period of 2015 and 2016, which means that land subsidence did occur during the autumn freezing stage in the study area.

Figure 5 Surface deformation obtained from SBAS-InSAR method and terrestrial leveling at TGL.D-InSAR observations have been converted into vertical deformation

As InSAR results have shown, the land surface in TGL rose by about 13 mm from mid-December 2014 to mid-July 2015; and it increased by 4 mm, with remarkable fluctuations from late-January to August 2016 (Figure 5). Meanwhile, the active-layer thickness increased 6 mm, according to the leveling measurement. The terrestrial leveling observation indicated no subsidence during summer 2016. On the contrary, the ground rose around 15 mm by October 2016, as compared with June 2016. The in-situ observations are sparser than the satellite acquisitions in temporal distribution. Displacement results of InSAR measurements were larger than the leveling observations before July 2016, but these values were less than the ground truth after that time.

That is the integration of phase velocity in each period between master and slave images. With this equation organized in matrix form, it can eventually be expressed as:

As shown in Figure 7, the trend and magnitude of deformation obtained from D-InSAR in the LDH region were inconsistent with in-situ observations. The leveling results were much higher than the measurements of radar technology, indicating an insensitivity of the deformation measured by Sentinel-1 C-band data to the freeze–thaw cycles of active-layer soils, indicating the Sentinel-1 data set performs poorly in the LDH region.

6 Discussion

The deformation process over cold regions is relatively intricate. On a local scale, permafrost displacement is determined by many factors, such as water/ice content of the active layer, its thickness, vegetation coverage, ground-surface roughness, and geomorphological processes, among many others (Liu et al., 2010; Chen et al., 2013). In permafrost regions,seasonal deformation is mainly affected by freeze thaw processes of the active layer, due to water–ice phase change, along with temporal variation, i.e.,ground-surface heaving in cold seasons and subsiding in warm seasons. In addition, long-term surface variation is caused by ice–water phase changes in soils,which might be a useful proxy for permafrost agradation or degradation.

6.1 Seasonal and long-term changes of ground-surface deformations in the TGL region

In this study, InSAR Scientific Computing Environment (ISCE) software was implemented in D-In-SAR processing because it is freely attainable and widely used in the InSAR scientific community (Rosen et al., 2012). An adaptive filter (Goldstein-Werner filter) was applied in this software to smooth interferograms (Goldstein and Werner, 1998), and statisticalcost, network-flow algorithm for phase unwrapping(SNAPHU) was employed in phase unwrapping in the final stage of interfering. Genetic InSAR Analysis Toolbox (GIAnT) software was applied in the timeseries analyzing and deformation computing process(Agram et al., 2012, 2013). A high-coherence area in most interfering results was assumed to be stale during the monitoring period, and this area was selected as the phase reference area during the time-series processing, which means that all inverted deformations referred to this area. Deformations of a specific pixel(from SBAS-InSAR time-series results) that corresponded with an in-situ measurement area was extracted and compared with ground truth, to validate the D-InSAR results. A flow chart of data processing is shown in Figure 2.

Figure 6 Cumulative deformation map of the TGL region (a) Reference surface height on October 27, 2014; (b) Cumulative deformation map on November 3, 2016 at TGL. A positive value represents surface uplift, and negative means surface settlement; white indicates areas of seriously incoherent signal (Different colors means different precipitation)

Figure 7 Surface deformation obtained from the SBAS-InSAR method and terrestrial leveling at LDH

Before time-series InSAR analysis, interferograms should first be produced. In the TGL region, all total 103 interferograms (Figure 3) were produced when setting the spatial perpendicular baseline less than 100 m and the interval less than 200 days between master and slave acquisitions. All the imagery data was cropped as 0.2°×0.2° grids to save computing time. A total of 324 interferograms(Figure 4) were produced in the LDH region. The largest perpendicular baseline and time separation were set as 100 m and 300 days, respectively. Similarly with the TGL region, only a small region in LDH, around 0.50°×0.35° in terms of latitude and longitude, was calculated in the interfering process.

Land-surface displacement over cold regions might be used to indicate permafrost changes, but it is difficult to monitor slight changes of active-layer thickness. For example, it is impossible to detect a surface settlement of 5 mm using ground temperature or drilling methods. In the TGL region, the decreased displacement trend might signify permafrost degradation because the melting of ice in the top of permafrost can induce surface settlement. Some studies have shown that significant variation of active-layer thickness occurred in the past three decades in the eastern and inland Tibetan Plateau (including the TGL region) because of remarkable climate warming (Zhao et al., 2004). A similar trend of permafrost degradation was also obtained on the Alaska North Slope between 1992 and 2000, with a magnitude around 10–40 mm per decade (Liu et al., 2012), while the subsidence of about 1–3 mm per year occurred, according to terrestrial leveling measurement and D-In-SAR results in the TGL region, separately. This phenomenon demonstrates a permafrost degradation in these regions under climate warming.

6.2 Applicability estimation of D-InSAR method over the TGL and LDH regions

Although it is difficult to estimate the accuracy of D-InSAR observation with statistical methods due to limited ground measurements, the differences between D-InSAR-derived displacements and leveling observations at TGL were less than 10 mm.However, irregular and incorrect deformations were obtained by the D-InSAR method over LDH. Such a large difference between TGL and LDH might be attributed to the differences in vegetation coverage,soil-water content, or a combination of the two.

Table 3 shows the porosity, vegetation coverage,and water content of all layers sampled at the two study sites. Details from this table illustrate that water content and porosity of soil at LDH is larger than at TGL in various sampled layers. We have discussed the frost heave in TGL during the groundthawing period due to large soil porosity and concentrated precipitation in summer. However, the same phenomenon did not display in LDH according to in-situ measurements, though porosity in LDH is larger than that in TGL; and precipitation there concentrated at the same period. Maybe normal deformation was related to exorbitant soil-water content in the LDH region during the thawing period, and longer leveling measuring periods would be necessary to validate this hypothesis.

Soil-water content and vegetation coverage may be the control factors for deformation-estimating in LDH. Soil with higher moisture content may result in specular scattering and coherence loss, which leads to incorrect deformation records because the phase values might be smoothed by the adaptive phase filter,and the unwarping phase may be produced incorrectly in these areas (Short et al., 2014). By contrast,soil-water content changed substantially from April to July (Table 3) in the LDH region. This change could cause complex phase changes, which has been demonstrated in laboratory and field experiments(Castañeda et al., 2009; Morrison et al., 2011). Another reason for the difference is high vegetation coverage in the LDH region, with some research indicating that vegetational variations show a significant consistency with the phase diversity; and it could adversely influence the final observational results of the ground-surface deformations (Zwieback and Hajnsek,2016).

AP病程中,外泌体在胰腺组织间液pH值的变化中也发挥了重要的调控作用。部分外泌体依靠本身的V-ATPase(一种细胞膜上的质子泵)主动提高囊泡内H+的浓度,导致组织液酸化,进而对周围组织造成一定的损伤[16]。

Table 3 Comparison of vegetation coverage and soil-water content over two study sites

?

7 Conclusions

The D-InSAR method can supply estimates of seasonal and long-term displacements of land surface over cold regions with broader coverage, and this approach has even higher temporal resolution than our ground-based leveling method. Compared with in-situ leveling measurements, Sentinel-1 C-band data performed well in the sparsely vegetated and drier TGL region. Some detailed variation, such as frost heave in summer time, can be identified by the D-InSAR method; and a long-term surface-subsidence tendency also can be detected by this approach. Seasonal change of land surface at the TGL site was less than 20 mm,which was mainly caused by the freeze–thaw cycle of active-layer soils; and long-term subsidence was around 1 to 3 mm per year, due to the thawing of ground ice. However, Sentinel-1 data showed poor results in the LDH region; and deformation was greatly undervalued by the D-InSAR method, mainly due to higher vegetation coverage and greater soil-water content.

Ground-based results are obtained by measuring limited points on the ground, but InSAR-observed deformations displayed changes over a broad area.Therefore, it is impossible to analyze accuracy between InSAR and leveling measurements based on statistical methods. In the future, some corner reflectors should be installed to enable acquiring more detailed comparative information between InSAR and ground-based measurements. Radar data with different bands also should be tested in the LDH region for the collection of surface deformations. For example,L-band data may produce better results in this region.

This work was supported by the Innovation Groups of the National Natural Science Foundation of China(41421061), the Chinese Academy of Sciences(KJZD-EW-G03-02), the project of the State Key Laboratory of Cryosphere Science (SKLCS-ZZ-2017), and CUHK Direct Grant (4053206). The Sentinel-1 SAR data were provided by the European Space Agency (ESA) through Sentinels Scientific Data Hub. The authors would like to thank the two anonymous reviewers, as well as the editor for their helpful comments and suggestions.

路线起点位于兰州市西固区新城镇黄河新桥南桥头,与已建的西固至新城一级公路终点顺接,沿黄河布线,经兰州维尼纶厂、河口、八盘峡、扶河、盐锅峡、恐龙湾、朱家台、孔家寺、大川,终点位于永靖县古城村,与临夏折桥至兰州达川二级公路及永靖县城市道路顺接,总体走向由北向南,全长48.24 km。

References:

Agram P, Jolivet R, Simons M, et al., 2012. GIAnT-generic InSAR analysis toolbox. American Geophysical Union, G43A-0897.

Agram PS, Jolivet R, Riel B, et al., 2013. New radar interferometric time series analysis toolbox released. EOS Transactions American Geophysical Union, 94(7): 69–70. DOI: 10.1002/2013EO070001.

Baker THW, Davis JL, Hayhoe HN, et al., 1982. Locating the frozenunfrozen interface in soils using time-domain reflectometry. Canadian Geotechnical Journal, 19(4): 511–517. DOI: 10.1139/t82-056.

Berardino P, Fornaro G, Lanari R, et al., 2002. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms. IEEE Transactions on Geoscience and Remote Sensing, 40(11): 2375–2383. DOI: 10.1109/Tgrs.2002.803792.

Brown J, Hinkel KM, Nelson FE, 2000. The circumpolar active layer monitoring (CALM) program: Research designs and initial results.Polar Geography, 24(3): 166–258. DOI: 10.1080/10889370009377698.

Castañeda C, Gutiérrez F, Manunta M, et al., 2009. DInSAR measurements of ground deformation by sinkholes, mining subsidence, and landslides, Ebro River, Spain. Earth Surface Processes and Landforms, 34(11): 1562–1574. DOI: 10.1002/esp.1848.

Chen FL, Lin H, Zhou W, et al., 2013. Surface deformation detected by ALOS PALSAR small baseline SAR interferometry over permafrost environment of Beiluhe section, Tibet Plateau, China. Remote Sensing of Environment, 138: 10–18. DOI: 10.1016/j.rse.2013.07.006.

Cheng KT, 1982. The forming process of thick layered ground ice. Science China Chemistry, 25(7): 777–788. DOI: 10.1360/yb1982-25-7-777.

Cheng GD, Wu TH, 2007. Responses of permafrost to climate change and their environmental significance, Qinghai Tibet Plateau. Journal of Geophysical Research: Earth Surface, 112(F2): F05S03. DOI:10.1029/2006jf000631.

Daout S, Doin MP, Peltzer G, et al., 2017. Large-scale InSAR monitoring of permafrost freeze-thaw cycles on the Tibetan Plateau.Geophysical Research Letters, 44(2): 901–909. DOI:10.1002/2016gl070781.

European Space Agency, 2013. Sentinel-1 User Handbook (2013-09-01)[2018-02-11]. https://earth.esa.int/documents/247904/685163/Sentinel-1_User_Handbook.

Farr TG, Rosen PA, Caro E, et al., 2007. The shuttle radar topography mission. Reviews of Geophysics, 45(2): RG2004. DOI: 10.1029/2005rg000183.

Ferretti A, Prati C, Rocca F, 2000. Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry.IEEE Transactions on Geoscience and Remote Sensing, 38(5):2202–2212. DOI: 10.1109/36.868878.

Goldstein RM, Werner CL, 1998. Radar interferogram filtering for geophysical applications. Geophysical Research Letters, 25(21):4035–4038. DOI: 10.1029/1998gl900033.

Harris AR, 1970. Direct Reading Frost Gage is Reliable, Inexpensive.Research Note NC-89 Paul, MN: Forest Service, North Central Forest Experiment Station.

Hinzman LD, Kane DL, Gieck RE, et al., 1991. Hydrologic and thermal properties of the active layer in the Alaskan Arctic. Cold Regions Science and Technology, 19(2): 95–110. DOI: 10.1016/0165-232x(91)90001-W.

Hu GJ, Zhao L, Wu XD, et al., 2015. Modeling permafrost properties in the Qinghai-Xizang (Tibet) Plateau. Science China Earth Sciences, 58(12): 2309–2326. DOI: 10.1007/s11430-015-5197-0.

Jiao YL, Li R, Zhao L, et al., 2014. Processes of soil thawing-freezing and features of soil moisture migration in the permafrost active layer. Journal of Glaciology and Geocryology, 36(2): 237–247. DOI:10.7522/j.issn.1000-0240.2014.0030.

Kane DL, Hinzman LD, Zarling JP, 1991. Thermal response of the active layer to climatic warming in a permafrost environment. Cold Regions Science and Technology, 19(2): 111–122. DOI: 10.1016/0165-232x(91)90002-X.

Li Z, Li XW, Liu YZ, et al., 2004. Detecting the displacement field of thaw settlement by means of SAR interferometry. Journal of Glaciology and Geocryology, 26(4): 389–396. DOI: 10.3969/j.issn.1000-0240.2004.04.003.

Little JD, Sandall H, Walegur MT, et al., 2003. Application of differential global positioning systems to monitor frost heave and thaw settlement in tundra environments. Permafrost and Periglacial Processes, 14(4): 349–357. DOI: 10.1002/ppp.466.

Liu L, Zhang TJ, Wahr J, 2010. InSAR measurements of surface deformation over permafrost on the North Slope of Alaska. Journal of Geophysical Research: Earth Surface, 115(F3): F03023. DOI:10.1029/2009jf001547.

Liu L, Schaefer K, Zhang TJ, et al., 2012. Estimating 1992–2000 average active layer thickness on the Alaskan North Slope from remotely sensed surface subsidence. Journal of Geophysical Research: Earth Surface, 117(F1): F01005. DOI: 10.1029/2011JF002041.

Liu XD, Chen BD, 2000. Climatic warming in the Tibetan Plateau during recent decades. International Journal of Climatology, 20(14):1729–1742. DOI: 10.1002/1097-0088(20001130)20:14<1729::AID-JOC556>3.0.CO;2-Y.

Mackay JR, 1973a. A frost tube for the determination of freezing in the active layer above permafrost. Canadian Geotechnical Journal,10(3): 392–396. DOI: 10.1139/t73-033.

Mackay JR, 1973b. The growth of Pingos, western arctic coast,Canada. Canadian Journal of Earth Sciences, 10(6): 979–1004.DOI: 10.1139/e73-086.

Mackay JR, 1977. Pulsating pingos, Tuktoyaktuk Peninsula, N. W. T.Canadian Journal of Earth Sciences, 14(2): 209–222. DOI:10.1139/e77-023.

Mackay JR, 1983. Downward water movement into frozen ground,western arctic coast, Canada. Canadian Journal of Earth Sciences,20(1): 120–134. DOI: 10.1139/e83-012.

Mark Nixon F, Taylor AE, 1998. Regional active layer monitoring across the sporadic, discontinuous and continuous permafrost zones, Mackenzie Valley, northwestern Canada. Proceedings of the 1998 Seventh International Conference on Permafrost. Canada,Yellowknife: Collection Nordicana, 815–820.

Morrison K, Bennett JC, Nolan M, et al., 2011. Laboratory measurement of the DInSAR response to spatiotemporal variations in soil moisture. IEEE Transactions on Geoscience and Remote Sensing,49(10): 3815–3823. DOI: 10.1109/Tgrs.2011.2132137.

Qin YH, Wu TH, Li R, et al., 2016. Using ERA-Interim reanalysis dataset to assess the changes of ground surface freezing and thawing condition on the Qinghai-Tibet Plateau. Environmental Earth Sciences, 75(9): 1–13. DOI: 10.1007/s12665-016-5633-2.

Rosen PA, Gurrola E, Sacco GF, et al., 2012. The InSAR scientific computing environment. Proceedings of the 9th European Conference on Synthetic Aperture Radar. Nuremberg, Germany: VDEAssociation for Electrical, Electronics and Information Technologies, 730–733.

Short N, LeBlanc AM, Sladen W, et al., 2014. RADARSAT-2 D-In-SAR for ground displacement in permafrost terrain, validation from Iqaluit Airport, Baffin Island, Canada. Remote Sensing of Environment, 141: 40–51. DOI: 10.1016/j.rse.2013.10.016.

Usai S, 2001. A New Approach for Longterm Monitoring of Deformations by Differential SAR Interferometry. Delft, Holland: Delft University Press.

Warren R, Brown J, 1972. The performance of a frost-tube for the determination of soil freezing and thawing depths. Soil Science,113(2): 149–154.

Wu QB, Zhang TJ, Liu YZ, 2010. Permafrost temperatures and thickness on the Qinghai-Tibet Plateau. Global and Planetary Change,72(1–2): 32–38. DOI: 10.1016/j.gloplacha.2010.03.001.

Xie C, Li Z, Xu J, et al., 2010. Analysis of deformation over permafrost regions of Qinghai-Tibet plateau based on permanent scatterers. International Journal of Remote Sensing, 31(8): 1995–2008.DOI: 10.1080/01431160902929255.

Yang YH, Zhu BZ, Jiang FQ, et al., 2012. Prevention and management of wind-blown sand damage along Qinghai-Tibet Railway in Cuonahu Lake area. Sciences in Cold and Arid Regions, 4(2):132–139. DOI: 10.3724/SP.J.1226.2012.00132.

Zebker HA, Rosen PA, Hensley S, 1997. Atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps. Journal of Geophysical Research: Solid Earth,102(B4): 7547–7563. DOI: 10.1029/96jb03804.

Zhang T, Osterkamp TE, Stamnes K, 1997. Effects of climate on the active layer and permafrost on the north slope of Alaska, USA. Permafrost and Periglacial Processes, 8(1): 45–67. DOI: 10.1002/(SICI)1099-1530(199701)8:1<45::AID-PPP240>3.0.CO;2-K.

Zhao L, Cheng GD, Li SX, et al., 2000. Thawing and freezing processes of active layer in Wudaoliang region of Tibetan Plateau.Chinese Science Bulletin, 45(23): 2181–2187. DOI: 10.1007/Bf02886326.

Zhao L, Ping CL, Yang DQ, et al., 2004. Changes of climate and seasonally frozen ground over the past 30 years in Qinghai-Xizang(Tibetan) Plateau, China. Global and Planetary Change,43(1–2): 19–31. DOI: 10.1016/j.gloplacha.2004.02.003.

Zorigt M, Kwadijk J, Van Beek E, et al., 2016. Estimating thawing depths and mean annual ground temperatures in the Khuvsgul region of Mongolia. Environmental Earth Sciences, 75: 897. DOI:10.1007/s12665-016-5687-1.

Zwieback S, Hajnsek I, 2016. Influence of vegetation growth on the polarimetric zero-baseline DinSAR phase diversity-implications for deformation studies. IEEE Transactions on Geoscience and Remote Sensing, 54(5): 3070–3082. DOI: 10.1109/Tgrs.2015.2511118.

ZhenMingWu,LinZhao,LinLiu,RuiZhu,ZeShenGao,YongPingQiao,LiMingTian,HuaYunZhou,MeiZhenXie
《Sciences in Cold and Arid Regions》2018年第2期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号