Soil Moisture Calibration Equations for Active Layer GPR Detection—a Case Study Specially for the Qinghai–Tibet Plateau Permafrost Regions

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 m3/m3. 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

Vegetation Mapping in the Permafrost Region: A Case Study on the Central Qinghai-Tibet Plateau

An accurate and detailed vegetation map is of crucial significance for understanding the spatial heterogeneity of subsurfaces, which can help to characterize the thermal state of permafrost. The absence of an alpine swamp meadow (ASM) type, or an insufficient resolution (usually km-level) to capture the spatial distribution of the ASM, greatly limits the availability of existing vegetation maps in permafrost modeling of the Qinghai-Tibet Plateau (QTP). This study generated a map of the vegetation type at a spatial resolution of 30 m on the central QTP. The random forest (RF) classification approach was employed to map the vegetation based on 319 ground-truth samples, combined with a set of input variables derived from the visible, infrared, and thermal Landsat-8 images. Validation using a train-test split (i.e., 70% of the samples were randomly selected to train the RF model, while the remaining 30% were used for validation and a total of 1000 runs) showed that the average overall accuracy and Kappa coefficient of the RF approach were 0.78 (0.68–0.85) and 0.69 (0.64–0.74), respectively. The confusion matrix showed that the overall accuracy and Kappa coefficient of the predicted vegetation map reached 0.848 (0.844–0.852) and 0.790 (0.785–0.796), respectively. The user accuracies for the ASM, alpine meadow, alpine steppe, and alpine desert were 95.0%, 83.3%, 82.4%, and 86.7%, respectively. The most important variables for vegetation type prediction were two vegetation indices, i.e., NDVI and EVI. The surface reflectance of visible and shortwave infrared bands showed a secondary contribution, and the brightness temperature and the surface temperature of the thermal infrared bands showed little contribution. The dominant vegetation in the study area is alpine steppe and alpine desert. The results of this study can provide an accurate and detailed vegetation map, especially for the distribution of the ASM, which can help to improve further permafrost studies

Vegetation Mapping in the Permafrost Region: A Case Study on the Central Qinghai-Tibet Plateau

An accurate and detailed vegetation map is of crucial significance for understanding the spatial heterogeneity of subsurfaces, which can help to characterize the thermal state of permafrost. The absence of an alpine swamp meadow (ASM) type, or an insufficient resolution (usually km-level) to capture the spatial distribution of the ASM, greatly limits the availability of existing vegetation maps in permafrost modeling of the Qinghai-Tibet Plateau (QTP). This study generated a map of the vegetation type at a spatial resolution of 30 m on the central QTP. The random forest (RF) classification approach was employed to map the vegetation based on 319 ground-truth samples, combined with a set of input variables derived from the visible, infrared, and thermal Landsat-8 images. Validation using a train-test split (i.e., 70% of the samples were randomly selected to train the RF model, while the remaining 30% were used for validation and a total of 1000 runs) showed that the average overall accuracy and Kappa coefficient of the RF approach were 0.78 (0.68–0.85) and 0.69 (0.64–0.74), respectively. The confusion matrix showed that the overall accuracy and Kappa coefficient of the predicted vegetation map reached 0.848 (0.844–0.852) and 0.790 (0.785–0.796), respectively. The user accuracies for the ASM, alpine meadow, alpine steppe, and alpine desert were 95.0%, 83.3%, 82.4%, and 86.7%, respectively. The most important variables for vegetation type prediction were two vegetation indices, i.e., NDVI and EVI. The surface reflectance of visible and shortwave infrared bands showed a secondary contribution, and the brightness temperature and the surface temperature of the thermal infrared bands showed little contribution. The dominant vegetation in the study area is alpine steppe and alpine desert. The results of this study can provide an accurate and detailed vegetation map, especially for the distribution of the ASM, which can help to improve further permafrost studies

Ground Deformation and Permafrost Degradation in the Source Region of the Yellow River, in the Northeast of the Qinghai-Tibet Plateau

The source region of the Yellow River (SRYR) is situated on the permafrost boundary in the northeast of the Qinghai-Tibet Plateau (QTP), which is an area highly sensitive to climate change. As a result of increasing global temperatures, the permafrost in this region has undergone significant degradation. In this study, we utilized Sentinel-1 to obtain ground surface deformation data in the SRYR from June 2017 to January 2022. We then analyzed the differences in terrain deformation under various environmental conditions. Our findings indicated an overall subsidence trend in the SRYR, with a long-term deformation velocity of −4.2 mm/a and seasonal deformation of 8.85 mm. Furthermore, the results showed that terrain deformation varied considerably from region to region, and that the Huanghe’ yan sub-basin with the highest permafrost coverage among all sub-basins significantly higher subsidence rates than other regions. Topography strongly influenced ground surface deformation, with flat slopes exhibiting much higher subsidence rates and seasonal deformation. Moreover, the ground temperature and ground ice richness played a certain role in the deformation pattern. This study also analyzed regional deformation details from eight boreholes and one profile line covering different surface conditions, revealing the potential for refining the permafrost boundary. Overall, the results of this study provide valuable insights into the evolution of permafrost in the SRYR region

The Zonation of Mountain Frozen Ground under Aspect Adjustment Revealed by Ground-Penetrating Radar Survey—A Case Study of a Small Catchment in the Upper Reaches of the Yellow River, Northeastern Qinghai–Tibet Plateau

Permafrost distribution is of great significance for the study of climate, ecology, hydrology, and infrastructure construction in high-cold mountain regions with complex topography. Therefore, updated high-resolution permafrost distribution mapping is necessary and highly demanded in related fields. This case study conducted in a small catchment in the northeast of the Qinghai Tibet Plateau proposes a new method of using ground-penetrating radar (GPR) to detect the stratigraphic structure, interpret the characteristics of frozen ground, and extract the boundaries of permafrost patches in mountain areas. Thus, an empirical–statistical model of mountain frozen ground zonation, along with aspect (ASP) adjustment, is established based on the results of the GPR data interpretation. The spatial mapping of the frozen ground based on this model is compared with a field survey dataset and two existing permafrost distribution maps, and their consistencies are all higher than 80. In addition, the new map provides more details on the distribution of frozen ground. In this case, the influence of ASP on the distribution of permafrost in mountain areas is revealed: the adjustment of ASP on the lower limit of continuous and discontinuous permafrost is 180–200 m, the difference in the annual mean ground temperature between sunny and shady slopes is up to 1.4–1.6 °C, and the altitude-related temperature variation and uneven distribution of solar radiation in different ASPs comprehensively affect the zonation of mountain frozen ground. This work supplements the traditional theory of mountain permafrost zonation, the results of which are of value to relevant scientific studies and instructive to engineering construction in this region

The Zonation of Mountain Frozen Ground under Aspect Adjustment Revealed by Ground-Penetrating Radar Survey—A Case Study of a Small Catchment in the Upper Reaches of the Yellow River, Northeastern Qinghai–Tibet Plateau

Permafrost distribution is of great significance for the study of climate, ecology, hydrology, and infrastructure construction in high-cold mountain regions with complex topography. Therefore, updated high-resolution permafrost distribution mapping is necessary and highly demanded in related fields. This case study conducted in a small catchment in the northeast of the Qinghai Tibet Plateau proposes a new method of using ground-penetrating radar (GPR) to detect the stratigraphic structure, interpret the characteristics of frozen ground, and extract the boundaries of permafrost patches in mountain areas. Thus, an empirical–statistical model of mountain frozen ground zonation, along with aspect (ASP) adjustment, is established based on the results of the GPR data interpretation. The spatial mapping of the frozen ground based on this model is compared with a field survey dataset and two existing permafrost distribution maps, and their consistencies are all higher than 80. In addition, the new map provides more details on the distribution of frozen ground. In this case, the influence of ASP on the distribution of permafrost in mountain areas is revealed: the adjustment of ASP on the lower limit of continuous and discontinuous permafrost is 180–200 m, the difference in the annual mean ground temperature between sunny and shady slopes is up to 1.4–1.6 °C, and the altitude-related temperature variation and uneven distribution of solar radiation in different ASPs comprehensively affect the zonation of mountain frozen ground. This work supplements the traditional theory of mountain permafrost zonation, the results of which are of value to relevant scientific studies and instructive to engineering construction in this region

Spatiotemporal Variations of Soil Temperature at 10 and 50 cm Depths in Permafrost Regions along the Qinghai-Tibet Engineering Corridor

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

Mapping Surficial Soil Particle Size Fractions in Alpine Permafrost Regions of the Qinghai–Tibet Plateau

Spatial information of particle size fractions (PSFs) is primary for understanding the thermal state of permafrost in the Qinghai-Tibet Plateau (QTP) in response to climate change. However, the limitation of field observations and the tremendous spatial heterogeneity hamper the digital mapping of PSF. This study integrated log-ratio transformation approaches, variable searching methods, and machine learning techniques to map the surficial soil PSF distribution of two typical permafrost regions. Results showed that the Boruta technique identified different covariates but retained those covariates of vegetation and land surface temperature in both regions. Variable selection techniques effectively decreased the data redundancy and improved model performance. In addition, the spatial distribution of soil PSFs generated by four log-ratio models presented similar patterns. Isometric log-ratio random forest (ILR-RF) outperformed the other models in both regions (i.e., R2 ranged between 0.36 to 0.56, RMSE ranged between 0.02 and 0.10). Compared with three legacy datasets, our prediction better captured the spatial pattern of PSFs with higher accuracy. Although this study largely improved the accuracy of spatial distribution of soil PSFs, further endeavors should also be made to improve model accuracy and interpretability for a better understanding of the interaction and processes between environmental predictors and soil PSFs at permafrost regions