A Monte-Carlo-Based Method for Probabilistic Permafrost Mapping Across Northeast China During 2003 to 2022

Permafrost degradation under climate warming has profound implications for ecological processes, hydrology, and human activities. Northeast China, characterized by sporadic and isolated patch permafrost near the southern limit of latitudinal permafrost (SLLP), represents one of the most sensitive and complex permafrost regions. This study aims to improve the reliability of permafrost mapping by incorporating parameter uncertainty into simulations. We developed a Monte Carlo–Temperature at the Top of Permafrost (TTOP) (MC–TTOP) framework that integrates an equilibrium model with Monte Carlo sampling to quantify parameter sensitivity and model uncertainty. Using all-sky daily air temperature data and land use and land cover information, we generated probabilistic estimates of mean annual ground temperature (MAGT), permafrost occurrence probability (PZI), and associated uncertainties. Model validation against borehole observations demonstrated improved accuracy compared with global-scale simulations, with a reduced bias and RMSE. Results reveal that permafrost in Northeast China was relatively stable during 2003–2010 but experienced pronounced degradation after 2011, with the total area decreasing to ~2.79 × 105 km2 by 2022. Spatial uncertainty was greatest in transitional zones near the southern boundary, where PZI-based delineations tended to overestimate permafrost extent. Regional comparisons further showed that permafrost in Northeast China is more fragmented and uncertain than that on the Tibetan Plateau, owing to complex snow–vegetation–topography interactions and intensive human disturbances. Overall, the MC-TTOP simulations indicate an accelerated permafrost degradation after 2011, with the highest uncertainty concentrated in southern transitional zones near the SLLP, demonstrating that the MC-TTOP framework provides a robust tool for probabilistic permafrost mapping, offering improved reliability for regional-scale assessments and important insights for future risk evaluation under climate change.