多年冻土作为高寒生态系统的重要组成，其变化对土壤水热状况、植被生长发育、生态系统碳氮循环、工程建设的稳定性产生重要影响。但目前关于多年冻土的变化研究大多集中在1960年以后，对更长时间尺度的冻土变化研究还存在一定的不足。本研究基于1901‒2020年CRU月平均气温数据，结合气象站点和土壤类型数据，通过TTOP模型和Stefan公式模拟了三江源区1901‒2020多年冻土顶板温度(temperature at the top of the perennially frozen/unfrozen ground, TTOP)、活动层厚度 (active layer thickness , ALT)、季节冻结深度(seasonal freezing depth, SFD)的时空变化特征（图1，2，3）。同时，通过气象站点月平均气温数据、已有冻土分布图和实际钻孔数据验证模拟结果。
研究表明，在过去120年里，多年冻土温度并非一直处于退化状态，而是随气候的冷暖交替呈现出波动变化。具体表现为：在1900s到1930s，多年冻土面积由263.9 × 10³ km²减少至233.3 × 10³ km² ；而在1940s到1970s，从232.3 × 10³ km² 增加至260.9 × 10³ km²；但到2010时，再次减少到228.9 × 10³ km²。TTOP也经历了上升‒下降‒再上升的过程，总体上表现出从1900s的–1.34 ± 2.74 °C上升至2010s的–0.48 ± 2.69 °C，其中以TTOP< −5 °C的变化最为显著，面积占比从8%减少至1%。在气候相对温暖的时期，多年冻土顶板温度升高、多年冻土面积减少、活动层厚度增加。但在多年冻土退化为季节冻土的区域，其季节冻结深度大多大于3.5m，因而，虽气候变暖会导致季节冻结深度变浅，但整个季节冻土区的季节冻结深度变化不大。

三江源区,多年冻土,TTOP模型,Stefan公式