青藏高原多年凍土區(qū)地溫分布特征

青藏高原多年凍土區(qū)地溫分布特征

0局域因素的影響青藏高原是世界上最富有的多年冷氣候地區(qū)。隨著世界氣候的持續(xù)變化和人類(lèi)活動(dòng)的增多，以及青藏鐵路等工程的建設(shè)和運(yùn)營(yíng)，青藏高原的多年凍土逐漸惡化。沿著青藏公路由楚瑪爾河高平原至風(fēng)火山,是青藏高原工程走廊重要的一段,在公路兩側(cè)寬約10km的區(qū)域,分布有公路、鐵路、通信光纜、輸油管道、天然氣管道、青藏±500kV直流電網(wǎng)等多條工程管線(xiàn)。該區(qū)域近40年來(lái)氣候存在變暖的趨勢(shì),升溫率為0.25℃/10a,20世紀(jì)90年代比60年代氣溫升高1.05℃,持續(xù)增溫將引起凍土區(qū)環(huán)境和凍土工程特性的顯著變化,只是受局地因素的影響存在一定的地區(qū)差異。局地因素直接影響地面接收太陽(yáng)輻射的強(qiáng)度,而青藏高原強(qiáng)烈的太陽(yáng)輻射使局地因素影響作用增強(qiáng)。因此,局地因素對(duì)小區(qū)域多年凍土的影響不容忽視。目前,國(guó)外的一些凍土空間模型中,有直接將凍土與氣候、地形和坡向進(jìn)行統(tǒng)計(jì)分析建模的,也有將坡度、坡向?qū)μ?yáng)輻射收入的影響折算成等效緯度進(jìn)行計(jì)算的,但國(guó)內(nèi)有關(guān)凍土空間模型的文獻(xiàn)中,尚未考慮局地因素的影響。然而氣候變暖引起多年凍土的退化,引發(fā)各種熱融災(zāi)害,如熱融滑塌引起的斜坡失穩(wěn),數(shù)量逐漸增多和面積逐漸變大的熱融湖塘、公路路基的沉陷變形與翻漿、鐵路路橋過(guò)渡段的變形沉陷與開(kāi)裂等給工程的正常運(yùn)營(yíng)帶來(lái)一定的危害。為了保護(hù)區(qū)域內(nèi)的凍土工程設(shè)施和進(jìn)行災(zāi)害評(píng)價(jià)、防治及預(yù)測(cè),建立考慮局地因素的凍土分布模型和繪制較大比例尺凍土分布圖具有重要意義。本文利用GIS良好的空間分析和制圖功能,提出了考慮局地因素坡向的凍土分布模型,輸出研究區(qū)內(nèi)凍土工程地溫區(qū)劃圖,并結(jié)合氣候轉(zhuǎn)暖對(duì)地溫的影響關(guān)系,預(yù)測(cè)了研究區(qū)50年后凍土的變化趨勢(shì)。1氣候變化影響多年凍土分布的研究點(diǎn)樣本點(diǎn)研究區(qū)以青藏公路為基準(zhǔn)線(xiàn),向兩側(cè)外延約5km,北起楚瑪爾河,南至風(fēng)火山,縱向長(zhǎng)約110km,海拔高程4420～5120m,對(duì)應(yīng)公路基樁號(hào)K2964～3073,介于北緯34°41′～35°26′、東經(jīng)92°50′～93°30′,面積約1100km2。29個(gè)樣本點(diǎn)主要分布在楚瑪爾河、五道梁、北麓河、風(fēng)火山等地的青藏公路沿線(xiàn),從北向南縱穿該區(qū)域內(nèi)的多年凍土區(qū),基本體現(xiàn)了高原凍土的高海拔特征,樣本點(diǎn)具備較好的代表性。青藏高原多年凍土是地氣系統(tǒng)長(zhǎng)期相互作用的產(chǎn)物,大的氣候背景決定了大范圍多年凍土分布的宏觀(guān)格局,但局地因素(如坡向、坡度、植被、巖性、雪蓋厚度及持續(xù)期、有機(jī)層和土壤特性等)的影響在一定條件下會(huì)超過(guò)大氣候背景,造成同樣氣候環(huán)境下局地多年凍土的異常分布。局地因素通過(guò)影響太陽(yáng)輻射、大氣對(duì)流和熱量傳導(dǎo)等過(guò)程進(jìn)而影響凍土區(qū)的地溫,從而影響局地區(qū)域內(nèi)多年凍土的存在與否。研究區(qū)內(nèi)地形地貌多樣,分布有高平原、山丘、河谷、盆地等,考慮局地因素的影響,方能得到更精準(zhǔn)的凍土分布結(jié)果。本文選取了研究區(qū)內(nèi)較易量化的植被覆蓋度、坡度、坡向3個(gè)局地因素進(jìn)行分析對(duì)比,篩選出影響最大的因子參與建模。1.1研究區(qū)植被覆蓋度植被能削減降雨侵蝕能量,調(diào)節(jié)地表徑流,阻滯泥沙遷移,還能保持水分,改良土壤,增強(qiáng)土壤滲透性能,提高其抗沖蝕性。高原上的植被對(duì)下伏的多年凍土具有較好的保護(hù)作用,并促進(jìn)有機(jī)質(zhì)的積累和土壤的發(fā)育,而多年凍土的隔水作用,又使植被得以在高寒干旱地區(qū)生存。例如,青藏高原的高寒草甸在冬季可形成隔熱層,阻礙土壤熱量散失,夏季可減少地面受熱,能使地表年溫差降低4～5℃。研究區(qū)內(nèi)植被主要分布于楚瑪爾河高平原、五道梁盆地和北麓河盆地,總體生長(zhǎng)稀疏,植株低矮,生存周期短,覆蓋度較低。通過(guò)對(duì)2006年8月覆蓋研究區(qū)的Quickbird全色影像運(yùn)用歸一化植被指數(shù)NDVI進(jìn)行計(jì)算,結(jié)合2011年9月的野外考察數(shù)據(jù)進(jìn)行校核修正,獲取研究區(qū)內(nèi)的植被覆蓋度圖層(圖1)。NDVI值的范圍是-1～1,一般綠色植被區(qū)的范圍是0.2～0.8,而研究區(qū)NDVI值為-0.378～0.394,植被覆蓋度較低。1.2坡位與地面亮說(shuō)坡度對(duì)植被的生長(zhǎng)影響較大,坡度越大越不利于植被生長(zhǎng)。不同的坡度之間,物種多樣性與豐富度也有較大差異。坡地上輻射分布有隨坡度而變化的規(guī)律,從而影響著坡地上的溫度分布。地面亮溫和地面濕度與坡度呈負(fù)相關(guān)。圖2是運(yùn)用ArcGIS的空間分析得出的坡度圖層,由圖2可知,研究區(qū)地勢(shì)整體較平緩,驟然抬升的區(qū)域較少。1.3地下冰發(fā)育情況不同坡向的地域接受的太陽(yáng)輻射差異極大,影響了水蒸發(fā)量、植被覆蓋、坡面侵蝕等諸多因素,從而造成相同海拔高度上不同坡向地域多年凍土的分布明顯不同。祁連山南北坡凍土下界分布高程差210m左右,214國(guó)道K359+100陰、陽(yáng)坡面下0.3m處年均地溫差達(dá)4.22℃。在風(fēng)火山上,陰、陽(yáng)坡年均地溫差1.7～2.4℃,凍土厚度差50～70m,地下冰發(fā)育情況也明顯不同。青藏鐵路DK1139+670處實(shí)測(cè)數(shù)據(jù)表明,陰、陽(yáng)坡面下0.5m處年均溫差為3℃以上。天山坡向?qū)Χ嗄陜鐾恋挠绊懢薮?在海拔高度相近的情況下,陰坡和陽(yáng)坡處年均地溫差達(dá)4.6℃。如果氣溫遞減率按0.5℃/100m計(jì)算,則坡向的影響相當(dāng)于920m高程差產(chǎn)生的結(jié)果。King提出了以坡向?yàn)橹饕罁?jù)的多年凍土分帶方案:將南坡無(wú)多年凍土的最高海拔定義為連續(xù)多年凍土下界,而將北坡出現(xiàn)多年凍土的最低海拔定義為不連續(xù)多年凍土下界。因此,在青藏高原進(jìn)行較小區(qū)域凍土空間建模時(shí),坡向是局地因素中不容忽視的一個(gè)影響因素。坡向通過(guò)影響太陽(yáng)輻射、大氣對(duì)流和熱量傳導(dǎo)等過(guò)程而影響凍土區(qū)的地溫,從而較大程度地影響著局地區(qū)域內(nèi)多年凍土的存在與否。本文主要考慮研究區(qū)內(nèi)坡向所造成的南北方向的地溫差異,故假定:南北方向坡向影響大小相等,方向相反;東西方向影響大小相等,方向相同。坡向數(shù)據(jù)的獲取如下:在ArcGIS平臺(tái)上,對(duì)研究區(qū)的數(shù)字高程數(shù)據(jù)DEM進(jìn)行空間分析獲取坡向圖層(圖3);參考天山的觀(guān)測(cè)數(shù)據(jù),運(yùn)用算式2.3cosθ(θ為坡向,0°～360°)將坡向量度轉(zhuǎn)化為線(xiàn)性量度,其值域?yàn)?2.3～2.3,結(jié)果如圖4。1.4年際變化的主要影響因素對(duì)研究區(qū)內(nèi)29個(gè)鉆孔點(diǎn)的經(jīng)緯度、高程、坡度、歸一化坡向及其植被覆蓋度與對(duì)應(yīng)的年均地溫值進(jìn)行Pearson相關(guān)性分析(表1),發(fā)現(xiàn)大尺度因素緯度、高程都與年均地溫在0.01顯著水平下呈負(fù)顯著相關(guān),經(jīng)度與年均地溫在0.05顯著水平下為正相關(guān),而且相關(guān)系數(shù)|高程|>|緯度|>|經(jīng)度|。局地因素坡向(歸一化數(shù)據(jù))與年均地溫在0.05顯著水平下的相關(guān)系數(shù)為0.268,在局地因素中最大,體現(xiàn)了小區(qū)域中陰陽(yáng)坡效應(yīng)的重要影響?；诖?選取緯度、高程及坡向作為影響因子,構(gòu)建研究區(qū)多年凍土分布模型。2結(jié)合品計(jì)算模型的年凍土分布模型青藏高原凍土模型或凍土對(duì)氣候的響應(yīng)模型總體可分為建立在傳熱學(xué)基礎(chǔ)上的物理模型和適合于凍土制圖的經(jīng)驗(yàn)?zāi)Ｐ?。?jīng)驗(yàn)?zāi)Ｐ痛蠖贾皇褂糜邢薜淖兞?與GIS結(jié)合緊密,具有空間性的優(yōu)點(diǎn),青藏高原應(yīng)用較廣泛的是高程模型和基于年均地溫的回歸模型。本文建立考慮局地因子坡向影響的年均地溫回歸模型,進(jìn)行多年凍土的分布與制圖研究。利用研究區(qū)內(nèi)29個(gè)鉆孔點(diǎn)2000-2010年的監(jiān)測(cè)數(shù)據(jù),對(duì)鉆孔點(diǎn)的年均地溫與高程、緯度和歸一化坡向進(jìn)行線(xiàn)性多元回歸統(tǒng)計(jì),得到考慮局地因子坡向的關(guān)系式:T=65.461-1.222φ-0.005H-0.299cosθ(1)式中:T為年均地溫(℃),φ為十進(jìn)制表示的緯度(°),H為高程(m),θ為坡向(°)。年均地溫與高程、緯度及坡向之間的復(fù)相關(guān)系數(shù)為0.936,具有較好的相關(guān)性。3研究區(qū)地溫區(qū)劃分凍土環(huán)境影響建筑物的穩(wěn)定性,線(xiàn)性工程的設(shè)計(jì)施工與維護(hù)保養(yǎng)尤其要考慮凍土因素。研究區(qū)是高原上管線(xiàn)工程較為密集的地區(qū),故進(jìn)行多年凍土分布模擬與制圖具有重要的意義。本文利用“回歸方程模擬值+修正值”的方法計(jì)算研究區(qū)的年均地溫,依據(jù)高原工程地溫分區(qū)標(biāo)準(zhǔn)進(jìn)行劃分,獲取相應(yīng)的工程地溫分區(qū)圖。具體操作如下:1)運(yùn)用ArcGIS平臺(tái)的柵格計(jì)算模塊,將高程、緯度和歸一化坡向因子代入公式(1)得到年均地溫模擬值圖層;2)對(duì)鉆孔樣本點(diǎn)的地溫實(shí)測(cè)值與模擬值求殘差,用反距離權(quán)重法將殘差插值到整個(gè)研究區(qū),得到相應(yīng)的殘差值圖層;3)將年均地溫模擬值圖層與殘差值圖層疊置求和,得到研究區(qū)修正后的年均地溫圖層。反距離權(quán)重插值法是以樣本點(diǎn)與插值點(diǎn)間的距離為權(quán)重的一種加權(quán)平均法,離插值點(diǎn)越近的樣本點(diǎn)賦予的估值權(quán)重越大,其表達(dá)式如下:z=[∑i=1nzid2i]/[∑i=1n1d2i](2)式中:z為所要估計(jì)的網(wǎng)格點(diǎn)的年平均地溫殘差值,zi為第i個(gè)鉆孔點(diǎn)的殘差值,di為插值點(diǎn)到第i個(gè)鉆孔點(diǎn)的距離,n為用于插值的鉆孔點(diǎn)數(shù)目。結(jié)合ArcGIS良好的數(shù)據(jù)分析和制圖功能,依據(jù)表2進(jìn)行重分類(lèi)和矢量化,產(chǎn)生的圖5即為研究區(qū)多年凍土工程地溫區(qū)劃圖。計(jì)算輸出的結(jié)果表明,研究區(qū)內(nèi)多年凍土地溫帶劃分與實(shí)際分布基本相符。根據(jù)分帶指標(biāo)進(jìn)行面積統(tǒng)計(jì)可知,多年凍土區(qū)占總面積的94.06%,其中,低溫穩(wěn)定帶占多年凍土區(qū)面積的15.94%,主要分布在風(fēng)火山和可可西里的高山基巖區(qū);低溫基本穩(wěn)定帶占16.97%,主要分布在風(fēng)火山及可可西里丘陵地帶;高溫不穩(wěn)定帶占48%,主要分布于可可西里和北麓河盆地東緣;高溫極不穩(wěn)地帶占19.09%,主要分布于北麓河盆地和楚瑪爾河高平原;季節(jié)凍土區(qū)或融區(qū)因較難用模型模擬,本文結(jié)合歷史資料和專(zhuān)家經(jīng)驗(yàn)進(jìn)行數(shù)字化修正,疊加到多年凍土工程地溫分區(qū)圖上。4凍土分布的立地條件通過(guò)Pearson相關(guān)性分析,篩選出局地因素——坡向參與研究區(qū)凍土分布模型的構(gòu)建。該模型充分反映了研究區(qū)內(nèi)多年凍土的分布特征,輸出符合實(shí)際的工程地溫區(qū)劃圖。多年凍土區(qū)占整個(gè)區(qū)域的94.06%,其中,低溫穩(wěn)定帶占多年凍土區(qū)面積的15.94%,主要分布在風(fēng)火山和可可西里的高山基巖區(qū);低溫基本穩(wěn)定帶占16.97%,主