The source region of the Yangtze River (SRYR, 1.4 × 10⁵ km² above Zhimenda station) on the Tibetan Plateau (TP) has 78% permafrost coverage. The streamflow depth increased at a rate of 2.5 mm/a since 2000. Quantification of the water contribution brought by permafrost thawing is a difficult task. In this study, we used Sentinel‐1 data and the SBAS‐InSAR technique to monitor terrain deformation from September 2016 to December 2021, and then utilized the long‐term deformation rate to assess ground ice meltwater release and the seasonal deformation to evaluate water storage in the active layer. Results reveal that 55.3% of the terrain in the SRYR has subsidence >2.5 mm/a, indicating widespread ground ice melting. The release rate of ground ice meltwater is 4.3 mm/a in the entire SRYR, above 6 mm/a at the Dangqu and Tuotuo River subbasins. The water release rate is relatively small (∼3%) in comparison to the streamflow depth of 151 mm per year during the investigation period of 2017–2021. We did not detect a strong increasing or decreasing trend among the 5‐year seasonal deformation, which reflects that the total soil water content in the active layer did not change significantly during the short investigation period. The results provide a data basis for ground ice richness and loss information in the SRYR and help to understand the impact of permafrost thawing on the regional water cycle in the permafrost environment.

In this study, we applied small baseline subset-interferometric synthetic aperture radar (SBAS-InSAR) to monitor the ground surface deformation from 2017 to 2020 in the permafrost region within an ~400 km × 230 km area covering the northern and southern slopes of Mt. Geladandong, Tanggula Mountains on the Tibetan Plateau. During SBAS-InSAR processing, we inverted the network of interferograms into a deformation time series using a weighted least square estimator without a preset deformation model. The deformation curves of various permafrost states in the Tanggula Mountain region were revealed in detail for the first time. The study region undergoes significant subsidence. Over the subsiding terrain, the average subsidence rate was 9.1 mm/a; 68.1% of its area had a subsidence rate between 5 and 20 mm/a, while just 0.7% of its area had a subsidence rate larger than 30 mm/a. The average peak-to-peak seasonal deformation was 19.7 mm. There is a weak positive relationship (~0.3) between seasonal amplitude (water storage in the active layer) and long-term deformation velocity (ground ice melting). By examining the deformation time series of subsiding terrain with different subsidence levels, we also found that thaw subsidence was not restricted to the summer and autumn thawing times but could last until the following winter, and in this circumstance, the winter uplift was greatly weakened. Two import indices for indicating permafrost deformation properties, i.e., long-term deformation trend and seasonal deformation magnitude, were extracted by direct calculation and model approximations of deformation time series and compared with each other. The comparisons showed that the long-term velocity by different calculations was highly consistent, but the intra-annual deformation magnitudes by the model approximations were larger than those of the intra-annual highest-lowest elevation difference. The findings improve the understanding of deformation properties in the degrading permafrost environment.

Soil temperature plays an essential role in the permafrost thermal state and degradation process. Especially the soil temperatures at 10 cm and 50 cm depths in the active layer, which are much easier to be observed in situ, have great effects on the surface water cycles and vegetation, and could be used as the upper boundary for permafrost models to simulate the thermal state of the permafrost and active layer thicknesses. However, due to the limitations of the observation data, there are still large uncertainties in the soil temperature data, including at these two depths, in the permafrost region of Qinghai–Tibet Plateau (QTP). In this study, we evaluated and calibrated the applicability of four daily shallow soil temperature datasets (i.e., MERRA-2, GLDAS-Noah, ERA5-Land, and CFSR) by using the in situ soil temperature data from eight observation sites from 2004 to 2018 in the permafrost region along the Qinghai–Tibet Engineering Corridor. The results revealed that there were different uncertainties for all four sets of reanalysis data, which were the largest (Bias = −2.44 °C) in CFSR and smallest (Bias= −0.43 °C) in GLDAS-Noah at depths of 10 cm and 50 cm. Overall, the reanalysis datasets reflect the trends of soil temperature, and the applicability of reanalysis data at 50 cm depth is better than at 10 cm depth. Furthermore, the GLDAS-Noah soil temperatures were recalibrated based on our observations using multiple linear regression and random forest models. The accuracy of the corrected daily soil temperature was significantly improved, and the RMSE was reduced by 1.49 °C and 1.28 °C at the depth of 10 cm and 50 cm, respectively. The random forest model performed better in the calibration of soil temperature data from GLDAS-Noah. Finally, the warming rates of soil temperature were analyzed, which were 0.0994 °C/a and 0.1005 °C/a at 10 cm and 50 cm depth from 2004 to 2018, respectively.

Ground-penetrating radar (GPR) is a convenient geophysical technique for active-layer soil moisture detection in permafrost regions, which is theoretically based on the petrophysical relationship between soil moisture (θ) and the soil dielectric constant (ε). The θ–ε relationship varies with soil type and thus must be calibrated for a specific region or soil type. At present, there is lack of such a relationship for active-layer soil moisture estimation for the Qinghai–Tibet plateau permafrost regions. In this paper, we utilize the Complex Refractive Index Model to establish such a calibration equation that is suitable for active-layer soil moisture estimation with GPR velocity. Based on the relationship between liquid water, temperature, and salinity, the soil water dielectric constant was determined, which varied from 84 to 88, with an average value of 86 within the active layer for our research regions. Based on the calculated soil-water dielectric constant variation range, and the exponent value range within the Complex Refractive Index Model, the exponent value was determined as 0.26 with our field-investigated active-layer soil moisture and dielectric data set. By neglecting the influence of the soil matrix dielectric constant and soil porosity variations on soil moisture estimation at the regional scale, a simple active-layer soil moisture calibration curve, named CRIM, which is suitable for the Qinghai–Tibet plateau permafrost regions, was established. The main shortage of the CRIM calibration equation is that its calculated soil-moisture error will gradually increase with a decreasing GPR velocity and an increasing GPR velocity interpretation error. To avoid this shortage, a direct linear fitting calibration equation, named as υ-fitting, was acquired based on the statistical relationship between the active-layer soil moisture and GPR velocity with our field-investigated data set. When the GPR velocity interpretation error is within ±0.004 m/ns, the maximum moisture error calculated by CRIM is within 0.08 m³/m³. While when the GPR velocity interpretation error is larger than ±0.004 m/ns, a piecewise formula calculation method, combined with the υ-fitting equation when the GPR velocity is lower than 0.07 m/ns and the CRIM equation when the GPR velocity is larger than 0.07 m/ns, was recommended for the active-layer moisture estimation with GPR detection in the Qinghai–Tibet plateau permafrost regions.