氣候變暖背景下青海三江源區季節凍土凍融特征研究

李萬志， 馬海玲， 龐昕瑋， 白文蓉， 李紅梅， 余 迪， 李漠雨

（1. 青海省氣候中心，青海 西寧 810001； 2. 高原科學與可持續發展研究院，青海 西寧 810001； 3. 常德市氣象局，湖南 常德 415000； 4. 中國科學院 大氣物理研究所，北京 100029）

0 引言

凍土是指溫度低于0 ℃且含冰的土壤或巖石，是一種對氣候變化十分敏感且性質極不穩定的土體［1-2］，其凍融過程是氣候變化的重要指示器［3-4］。凍土層作為一種特殊的區域性隔水層或弱透水層，它的存在及其對全球變化的響應改變了寒區水文徑流和生態環境［5］。影響凍土最大季節凍結深度變化的成因十分復雜，氣候因素、地形因素和人為活動等均可對季節凍結深度產生影響，其中氣候因素影響較大。受氣候變暖影響，季節凍土厚度不斷變淺，多年凍土層變薄甚至消失［6-7］。近50 年來，中國地表氣溫升高了1.1 ℃，增溫速率明顯高于全球或北半球同期的增溫速率［8］。IPCC 第六次報告（AR6）指出，氣候變暖的速度正在加快，尤其是全球地表平均溫度的上升速率十分明顯［9］。

在受氣候變化影響的各圈層中，冰凍圈首當其沖，是對氣候系統的反饋作用較為顯著和快速的圈層［10-13］。在極端天氣狀態下，青藏高原冰凍圈和相關聯的生態系統變化可能產生級聯放大的影響［14-15］。青藏高原是地球上面積最大的高海拔多年凍土區，占我國多年凍土總面積的70%左右［16-17］。Gao 等［18］、Guo 等［19］分別采用不同的物理模型模擬了流域尺度、青藏高原地區以及全球尺度的多年凍土變化，揭示出凍土的退化事實。尤其針對青藏高原的凍土變化，諸多研究表明，氣候變暖對該地區季節凍土影響主要表現為年最大凍結深度減小、凍結初始日推遲、融化終止日提前、凍結持續期縮短等變化趨勢［20-22］，導致該地區出現植被退化、冰川退縮等一系列生態環境問題［23-24］。因此，開展凍土變化特征研究，揭示凍土退化過程和機制，已成為氣候和冰凍圈科學研究迫在眉睫的任務。

青海三江源（本文簡稱三江源）地處青藏高原腹地，青海省南部，以山地地貌為主，山脈綿延，地形復雜，境內河流密布，湖泊、沼澤眾多，雪山、冰川廣布，是長江、黃河、瀾滄江三大河流的發源地，也是中國乃至東亞地區重要的水資源保護區，具有高寒生物資源豐富、生態環境脆弱以及對氣候變化敏感的特征［25-26］。作為青藏高原生態屏障極為重要的組成部分，高亢的地勢和嚴寒的氣候使其成為我國典型的冰凍圈核心作用區之一［27］，大面積高原多年凍土的存在孕育了廣泛的高寒草原、高寒草甸與高寒濕地［26］，其生態系統服務功能對中國的生態狀況及國民經濟發展起著重要的作用［29-30］。受氣候變暖和人為活動的共同影響，三江源正面臨著以“變暖變濕”為主的氣候變化，是全球氣候變化的敏感區和生態環境的脆弱區［31-33］?；诖?，本文在前人研究成果的基礎上，通過補充最近較長序列氣象臺站凍土深度和地溫等氣象觀測資料，對比分析氣候變暖前后三江源區的季節凍土凍融特征，研究成果可為三江源區高寒生態治理及寒區工程建設提供理論依據。

1 資料與方法

1.1 研究區概況

三江源區地處青藏高原腹地，地理位置為31°39′～37°10′ N、89°24′～102°27′ E，海拔為2 800～6 564 m，面積為13.0×104km2，包括長江、黃河和瀾滄江源區等三大源區，其面積分別為13.0×104km2、11.8×104km2和5.4×104km2。境內河流密布，湖泊、沼澤眾多，雪山冰川廣布，是中國面積最大的天然濕地分布區，被譽為“中華水塔”［34-35］。這里氣候屬典型高原大陸性氣候，年溫差小、日溫差大，年平均氣溫為-5.1～8.3 ℃；年降雨量為253.1～732.6 mm。

圖1 三江源區地形及氣象站點分布圖Fig. 1 Topography and distribution of meteorological stations in the Three-River Source Region on northeastern Qinghai-Tibet Plateau

1.2 數據來源

該研究數據主要采用三江源區18 個國家氣象站1961—2021 年的凍土、地溫以及氣溫觀測資料。資料均來自青海省氣象信息中心，并經過嚴格的質量控制（表1）。由于各氣象站季節凍結深度開始觀測年份不同，所以在做最大季節凍結深度特征研究時，剔除五道梁、沱沱河站點。其次，為保證時間年限的一致性，剩余各站最大季節凍結深度研究時間均從1981 年開始，清水河站1981—1989 年凍土觀測數據采用5 年滑動平均值進行補充。所以，該研究主要基于16 個國家氣象站1981—2021 年凍土最大季節凍結深度觀測資料分析其變化特征。平均地表溫度、氣溫等觀測數據較為完整，在進行季節凍結時間研究時采用1961—2021 年觀測數據。文中涉及的地圖是基于自然資源部標準地圖服務系統（http：//bzdt. ch. mnr. gov. cn/）下載的審圖號為GS（2019）3266的青海省地圖制作，底圖無修改。

表1 三江源區氣象站點基本信息Table 1 Basic information of meteorological stations in the Three-River Source Region

1.3 研究方法

文中主要根據日平均地溫（0 cm）<0 ℃的開始和結束時間代表土的凍結和融化時間，地表土的凍結初日是上一年9 月1 日后淺表層土（0～3.2 m）第一次凍結日期，土的凍結終日是當年8 月31 日前最后一次出現淺表層土（0～3.2 m）凍結的日期，年內日平均地表溫度（0 cm）<0 ℃天數反映了凍土凍結時間，1 年中土的凍結深度的最大值為年最大凍結深度［36］。在進行統計分析時采用了線性趨勢、累計距平［37］和顯著性t檢驗［38］等方法。累積距平曲線是確定氣候突變的一種方法，計算公式為：

式中：Xi為氣象要素歷年值為氣象要素多年平均值，若氣象要素累積距平指標絕對值達到最大時，所對應的t為突變年份。同時，通過計算信噪比［39］檢驗轉折是否達到氣候突變的標準。計算公式如下：

2 結果與分析

2.1 平均氣溫變化

1961—2021 年，三江源區年平均氣溫為-0.34 ℃，總體表現為顯著的升高趨勢，升溫速率為0.38 ℃·（10a）-1［圖2（a）］。氣溫的變化過程并不是持續、均勻的單調增暖，而是呈現出冷暖相間、波動上升的特點。20 世紀70 年代至90 年代后期平均氣溫相對較低，且變化較為平穩，在1983 年出現最低值（-1.7 ℃），進入21 世紀后氣溫開始呈現顯著增溫趨勢，并在2020 年達到最高值（2.0 ℃）。從三江源區年平均氣溫累積距平變化曲線來看，年平均氣溫累積距平呈明顯的“單谷型”分布，在1997年出現轉折點，1997 年之前呈顯著的下降趨勢，并在1997年達到“谷底”后開始回升［圖2（b）］，計算其S/N值為1.2，大于1.0，說明1997 年為三江源區平均氣溫的突變年份，這與整個青藏高原年平均氣溫突變年份一致［40-41］。所以，在研究三江源區變暖前后凍土變化特征時，以1997 年為時間節點，1997 年以前為變暖前期，之后為變暖后期。

圖2 1961—2021年三江源區年平均氣溫、年平均氣溫累計距平變化曲線（單位：℃）Fig. 2 Changing curve of annual average temperature （a） and cumulative annual average temperature anomaly （b）in the Three-River Source Region on northeastern Qinghai-Tibet Plateau from 1961 to 2021（unit： ℃）

通過三江源區國家基準站年平均氣溫空間分布特征來看，年平均氣溫為-5.0～4.6 ℃，總體呈東高西低的趨勢，其中囊謙站年平均氣溫最高，五道梁站最低［圖3（a）］。從年平均氣溫變率空間分布來看，三江源區年平均氣溫變率為-0.03～0.83 ℃·（10a）-1，除河南站呈下降趨勢外，其他站均呈升溫趨勢，通過0.001顯著性水平檢驗，其中同德站升溫速率最高［圖3（b）］。對比變暖前后各站平均氣溫差值得出，變暖后年平均氣溫較變暖前升高0.3～3.0 ℃，其中同德站增溫幅度最大，河南站增溫幅度最?。蹐D3（c）］。

圖3 1961—2021年三江源區年均氣溫（Tavg）、年均氣溫變率（Tslope）增暖前后氣溫差值（ΔT）空間分布圖Fig. 3 Spatial distribution of annual average temperature（Tavg） （a）， annual average temperature variability （Tslope）（b）and temperature difference before and after warming（ΔT） （c） in the Three-River Source Region on northeastern Qinghai-Tibet Plateau from 1961 to 2021

2.2 年最大季節凍結深度變化

1961—2021 年三江源區各站平均年最大季節凍結深度為142.5 cm，受氣溫變暖影響呈退化趨勢，退化速率為2.4 cm·（10a）-1。在年代變化上，20世紀80 年代后期至90 年代為季節凍結深度最大的時段，最大值出現在1986 年（172.7 cm）。進入21世紀后，季節凍結深度開始出現明顯退化，在2019年出現最低值（108.8 cm）［圖4（a）］，且變暖后年最大季節凍結深度較變暖前減少了11 cm。

圖4 1961—2021年三江源區最大季節凍結深度（單位：cm）、最大季節凍結深度累積距平變化曲線及最大季節凍結深度（單位：cm）、最大季節凍結深度變率［單位：cm·（10a）-1］空間分布圖Fig. 4 Changing curve of annual maximum depth of frost penetration （unit： mm） （a）， cumulative anomaly of the mean annual maximun seasonal freezing depth （b）， spatial distribution map of the mean annual maximun seasonal freezing depth （unit： mm） （c），and change rate of the mean annual maximun seasonal freezing depth change rate ［unit： cm·（10a）-1］ （d） in the Three-River Source Region on northeastern Qinghai-Tibet Plateau from 1961 to 2021

從年最大季節凍結深度累積距平值來看，三江源區年最大季節凍結深度有兩個轉折點，1979 年之前季節凍結深度變化相對平穩，之后呈顯著增加趨勢，在2004 年達到峰值后開始呈減小趨勢［圖4（b）］。計算兩個轉折點S/N值，其中1979 年的S/N值為0.6，小于1.0，2004 年的S/N值為1.1，大于1.0，說明年最大季節凍結深度在2004年發生突變，突變后年最大季節凍結深度較之前減小了16.1 cm。年最大季節凍結深度雖然受氣溫變化影響，但其突變年份與氣溫突變年份相比有一定推后，說明季節凍結深度由于受其他因素影響，對氣溫變化有一定的滯后性，這與李林等［32］研究結論較為一致。

空間分布上，三江源區國家基準站年最大季節凍結深度自西北向東南減小，最大值位于瑪多站（211.9 cm），最小值位于囊謙站（66.7 cm）［圖4（c）］?？臻g變率上，玉樹站、瑪多站、河南站年最大季節凍結深度呈增加趨勢，增加速率為0.07～4.0 cm·（10a）-1，其中玉樹站增加速率最大，其余13 站均呈減少趨勢，減少速率為1.1～16.8 cm·（10a）-1，曲麻萊站減少趨勢最顯著［圖4（d）］。對比變暖前后各站年最大季節凍結深度得出，玉樹站、瑪多站變暖后較變暖前季節凍結深度分別增加了6.1 cm 和8.6 cm，其余各站均出現退化，退化幅度為4.4～48.2 cm，其中清水河站退化幅度最大，治多站最小。

2.3 地表土的凍結初日、凍結終日變化

1961—2021 年，三江源區平均地表凍結初日為10 月24 日，呈顯著推遲趨勢，推遲速率為1.0 d·（10a）-1，通過0.001 顯著性水平檢驗。在年代變化上，地表凍結初日在21 世紀前變化并不明顯，21世紀后明顯開始推遲，尤其近十年推遲趨勢明顯［圖5（a）］。

圖5 1961—2021年三江源區地表土的凍結初日、凍結終日變化曲線及氣候變暖前后各觀測站地表土的凍結初日和凍結終日變化Fig.5 Changes in of freezing onset date （a） and freezing end date （b） of ground surface soil in the Three-River Source Region on northeastern Qinghai-Tibet Plateau from 1961 to 2021 and changes in freezing onset date and freezing end date of ground surface soil at each station before and after climate warming （c）

平均地表凍結終日為5 月18 日，呈顯著提前趨勢，提前速率為3.3 d·（10a）-1，通過0.001 顯著性水平檢驗。其中，地表凍結終日最晚的日期為1990年6月6日，最早日期為2021年4月21日［圖5（b）］。地表凍結初日和凍結終日在氣候變暖前后差異較大，變暖前凍結終日為5月22日，變暖后為5月10日，提前了12 d；而凍結初日在變暖前為9月25日，變暖后為10月9日，推遲了14 d，與氣溫變化有著很好的響應。

從三江源區各站季節性凍土的凍結初日和凍結終日變化來看，澤庫站、治多站、瑪多站、清水河站、甘德站、達日站、河南站7 站在9 月就開始凍結，其中甘德站最早（9 月8 日），其余11 站在10 月才開始凍結，囊謙站最晚（10月24日）。同德站、玉樹站、囊謙站、班瑪站4 站在4 月就已經開始消融，其中囊謙站最早（4 月17 日）；澤庫站、甘德站6 月開始消融，清水河站7 月開始消融，消融日期最晚，其余11站均是5月開始消融。從變暖前后各站地表的凍結初日和凍結終日變化來看，凍結初日除瑪多站提前外，其他各站均推遲。其中，瑪多站增暖后較增暖前凍結初日提前14 d，其余各站推遲2.1～28.6 d，曲麻萊站推遲時間最長。凍結終日各站提前0.3～45.1 d，其中瑪多站提前時間最長，班瑪站最短［圖5（c）］。

2.4 季節凍土凍結時間變化

三江源區季節性凍土的平均凍結時間為133.9 d，呈顯著減少趨勢，減少速率為1.9 d·（10a）-1。在年代變化上，20 世紀60 至70 年代凍結時間變化不明顯，80年代持續時間最長，90年代開始減少；進入21世紀后減少趨勢尤為顯著，其中凍結時間最長的年份為1983 年（150.8 d），最短為2015 年（117.2 d）［圖6（a）］，且變暖前凍結時間為137.5 d，變暖后為128.7 d，相差8.8 d。從凍結時間累積距平值來看，三江源區季節性凍土的凍結時間累計距平呈“單峰型”分布，在2002 年之前呈增多趨勢，到2002 年達到頂峰后開始減少［圖6（b）］。計算其S/N值為1.2，大于1.0，說明凍結時間在2002年發生了突變，突變年份與年最大季節凍結深度突變年份一致，且受氣溫變化影響均具有一定的滯后性。

圖6 1961—2021年三江源區季節性凍土的凍結時間、凍結時間累計距平變化曲線Fig.6 Changing curve of seasonally frozen soil duration （a） and accumulated anomaly of seasonally frozen soil duration （b）in the Three-River Source Region on northeastern Qinghai-Tibet Plateau from 1961 to 2021

在空間分布上，三江源區各站季節性凍土的凍結時間為88.3～180.6 d，總體呈西高東低的趨勢，其中五道梁站凍結時間最長，囊謙站最短［圖7（a）］。在長期變化上，三江源區各站土的凍結時間均呈減少趨勢，速率為1.7～7.8 d·（10a）-1，均通過顯著性水平檢驗。其中，河南站減少速率最小，同德站最大［圖7（b）］。對比變暖前后各站凍結時間變化得出，變暖后凍結時間較變暖前減少了8.6～27.4 d，其中同德站減少幅度最大，河南站最?。蹐D7（c）］。

圖7 三江源區季節性凍土的凍結時間、凍結時間變率、增暖前后凍結時間差空間分布圖Fig. 7 Spatial distribution of seasonally frozen ground duration （a）， seasonal ground freezing time variability （b） and seasonally frozen ground freezing time difference （c） before and after warming in the Three-River Source Region on northeastern Qinghai-Tibet Plateau

3 結論

通過1961—2021 年三江源區18 個國家氣象站的地溫、氣溫、凍土等觀測資料統計結果，計算得出了三江源區氣溫、年最大季節凍結深度以及淺表層土的凍結時間的突變年份，并對比分析了氣候變暖前后季節凍土的變化特征，得到如下結論：

（1）三江源區年平均氣溫為-0.34 ℃，呈東高西低分布，總體以0.38 ℃·（10a）-1速率增溫。各站增溫速率在-0.03～0.83 ℃·（10a）-1之間，除河南站呈降低趨勢外，其他均呈升溫趨勢。三江源區年平均氣溫在1997年發生突變，且1997年為氣候變暖前后時間節點，與變暖前相比，各站增溫在0.3～3.0 ℃。

（2）三江源區平均年最大季節凍結深度為142.5 cm，呈西北高東南低分布，總體以2.4 cm·（10a）-1速率退化。除玉樹站、瑪多站、河南站呈增加趨勢外，其他各站均呈退化趨勢。較變暖前相比，年最大季節凍結深度退化了11 cm。除玉樹站、瑪多站季節凍結深度增加外，其他各站均出現退化。三江源區年最大季節凍結深度在2004 年發生了突變，突變年份與氣溫突變年份不一致，且具有一定的滯后性。

（3）三江源區平均地表凍結初日為10 月24 日，以1.0 d ·（10a）-1速率推遲，平均凍結終日為5 月18日，以3.3 d·（10a）-1速率提前?？臻g上，凍結初日甘德站最早、囊謙站最晚，凍結終日為囊謙站最早、清水河站最晚。與變暖前相比，凍結終日提前12 d，凍結初日推遲14 d，各站凍結初日除瑪多站提前外，其他均推遲，凍結終日則均呈出現提前。

（4）三江源區季節性凍土的平均凍結時間為133.9 d，呈西多東少分布，總體以1.9 d·（10a）-1速率減少，空間上各站也均呈減少趨勢。相比變暖前，各站凍結時間減少8.6～27.4 d，同德站減少幅度最大。三江源區土的凍結時間在2002 年發生了突變，突變年份與平均氣溫并不一致，說明土的凍結時間受氣溫變化影響具有一定的滯后性。

受氣候變暖影響，三江源區季節性凍土的凍結深度及凍結時間均呈減少趨勢，凍結終日呈提前趨勢，初日呈推遲趨勢。但通過計算季節凍結深度及凍結時間的突變年份，發現其與氣溫的突變年份并不一致，說明季節凍土變化在受氣候變暖影響的同時，還受人類活動、下墊面等其他因素影響，后期研究中還需進一步加強該方面的研究，厘清季節凍土退化事實。同時，三江源區地廣人稀，站點稀疏，現有觀測站點數據并不能很好地反映三江源區季節凍結作用弱化的現象，后期還需加密觀測，使分析結果更加精確。