Mapping Frozen Ground in the Qilian Mountains in 2004–2019 Using Google Earth Engine Cloud Computing

• Published in: Remote Sensing
• Main Authors: Yuan Qi, Shiwei Li, Youhua Ran, Hongwei Wang, Jichun Wu, Xihong Lian, Dongliang Luo
• Format: Text
• Language: English
• Published: Multidisciplinary Digital Publishing Institute 2021
• Subjects: Qilian mountains; permafrost; Stefan Equation; TTOP model; climate change; Google Earth Engine; Active layer thickness
• Online Access: https://doi.org/10.3390/rs13010149

The permafrost in the Qilian Mountains (QLMs), the northeastern margin of the Qinghai–Tibet Plateau, changed dramatically in the context of climate warming and increasing anthropogenic activities, which poses significant influences on the stability of the ecosystem, water resources, and greenhouse gas cycles. Yet, the characteristics of the frozen ground in the QLMs are largely unclear regarding the spatial distribution of active layer thickness (ALT), the maximum frozen soil depth (MFSD), and the temperature at the top of the permafrost or the bottom of the MFSD (TTOP). In this study, we simulated the dynamics of the ALT, TTOP, and MFSD in the QLMs in 2004–2019 in the Google Earth Engine (GEE) platform. The widely-adopted Stefan Equation and TTOP model were modified to integrate with the moderate-resolution imaging spectroradiometer (MODIS) land surface temperature (LST) in GEE. The N-factors, the ratio of near-surface air to ground surface freezing and thawing indices, were assigned to the freezing and thawing indices derived with MODIS LST in considerations of the fractional vegetation cover derived from MODIS normalized difference vegetation index (NDVI). The results showed that the GEE platform and remote sensing imagery stored in Google cloud could be quickly and effectively applied to obtain the spatial and temporal variation of permafrost distribution. The area with TTOP < 0 °C is 8.4 × 104 km2 (excluding glaciers and lakes) and accounts for 46.6% of the whole QLMs, the regional mean ALT is 2.43 ± 0.44 m, while the regional mean MFSD is 2.54 ± 0.45 m. The TTOP and ALT increase with the decrease of elevation from the sources of the sub-watersheds to middle and lower reaches. There is a strong correlation between TTOP and elevation (slope = −1.76 °C km−1, p < 0.001). During 2004–2019, the area of permafrost decreased by 20% at an average rate of 0.074 × 104 km2·yr−1. The regional mean MFSD decreased by 0.1 m at a rate of 0.63 cm·yr−1, while the regional mean ALT showed an exception of a decreasing ...