烏伊嶺凍土區(qū)活動(dòng)層厚度影響因素分析

呂姝琴

（哈爾濱師范大學(xué)）

0 引言

隨著全球氣候變暖，極端氣候事件頻發(fā)［1］.IPCC第六次評(píng)估報(bào)告指出，將在2100 年前全球平均溫升控制在2℃以內(nèi)［2］.寒區(qū)生態(tài)系統(tǒng)對(duì)氣候變化的響應(yīng)研究已成為熱點(diǎn)［3］.近地表面的多年凍土是寒區(qū)陸地生態(tài)系統(tǒng)重要的組成部分，對(duì)多年凍土的研究是生態(tài)、水文、地質(zhì)和工程建設(shè)研究者關(guān)注的重要議題［4］.

全球氣候變暖對(duì)冰凍圈的影響巨大，凍土作為冰凍圈的組成部分，也受其影響［5］.季節(jié)凍結(jié)是指地層年平均地溫高于土壤凍結(jié)溫度的地層的凍結(jié)，是地層在冬季放出熱量的結(jié)果［6］.季節(jié)融化則是地層的年平均地溫低于凍結(jié)溫度的凍土的融化，是由于熱量暖季進(jìn)入地層所引起.通常把每年冬季凍結(jié)、夏季融化的接近地表的土層稱為活動(dòng)層，活動(dòng)層是凍土地層內(nèi)水熱交換最主要的區(qū)域［7］，受經(jīng)緯度、高程、植被覆蓋度、地表溫度、土壤性質(zhì)及氣候環(huán)境等諸多因素的影響.因此，對(duì)于凍土活動(dòng)層厚度的監(jiān)測(cè)是了解凍土對(duì)環(huán)境響應(yīng)的重要基礎(chǔ).

小興安嶺森林是溫帶、寒溫帶的過渡地帶，屬于針闊葉混交林與北方針葉林的過渡區(qū)域，處于中國受氣候變暖影響最顯著的地區(qū)對(duì)環(huán)境變化具有很高的敏感性［8］.而烏伊嶺區(qū)作為最典型的、稀有的林間濕地區(qū)，其凍土活動(dòng)層的變化對(duì)濕地植被的演替以及植被的保護(hù)有著重要的意義.許多專家和學(xué)者的研究表明，有關(guān)凍土活動(dòng)層厚度變化及其影響因素的研究多集中于青藏高原，鮮少有涉及到小興安嶺地區(qū).因此，該文以小興安嶺地區(qū)烏伊嶺區(qū)為例探究?jī)鐾粱顒?dòng)層厚度與影響因子，土壤濕度、土壤溫度以及高程和經(jīng)緯度之間的關(guān)系.

1 研究區(qū)與研究方法

1.1 研究區(qū)

烏伊嶺區(qū)位于黑龍江省東北部小興安嶺頂峰（128°57′～129°44′E、48°33′～49°08′N）（如圖1所示）.氣候上屬溫帶大陸性季風(fēng)氣候，由于受海洋環(huán)流和西伯利亞冷空氣影響，四季氣候呈明顯特點(diǎn).冬季漫長(zhǎng)而寒冷，多西北風(fēng)；夏季短促而燥熱；多東南風(fēng)和大暴雨，降水量大而且集中.年平均氣溫－1.1 ℃，歷年平均最低氣溫－8.5 ℃，最高氣溫6.1 ℃.1 月份氣溫最低，平均氣溫－24.6 ℃.7月份氣溫最高，平均氣溫19.1 ℃.年平均降水量585.7 mm.地形上屬低山丘陵，地勢(shì)多東南，西北坡向，山峰沿著支脈兩側(cè)呈枝狀分布，頂峰起伏不大，坡勢(shì)平緩.最高點(diǎn)在翠峰林場(chǎng)與東克林林場(chǎng)分界處，海拔606 m；最低點(diǎn)在扎克大旗河匯入烏云河的河口處，海拔157 m.有中國北方特有的森林沼澤、灌叢沼澤和泥炭蘚沼澤［9］.此外，還有各種鳥類、魚類及國家重點(diǎn)保護(hù)動(dòng)物，他們共同形成了北方高緯度、多種類、復(fù)合型濕地自然生態(tài)系統(tǒng).

圖1 研究區(qū)域

1.2 研究方法

凍土活動(dòng)層厚度主要通過釬探法、坑探法、鉆探法和土壤溫度線性插值法［10］進(jìn)行測(cè)量，該文主要通過金屬釬探法采集活動(dòng)層厚度數(shù)據(jù)，并輔以探地雷達(dá)數(shù)據(jù)進(jìn)行驗(yàn)證.釬探法運(yùn)用“T”型金屬探測(cè)棒垂直于地面向下插入泥土之中，直至金屬探測(cè)棒不再繼續(xù)往下深入，記錄當(dāng)前金屬探測(cè)棒深度，得到該點(diǎn)的活動(dòng)層厚度數(shù)據(jù).探地雷達(dá)數(shù)據(jù)采集則使用500 MHz和250 MHz屏蔽天線探地雷達(dá)，采用中心頻率分別為500 MHz 和250 MHz的兩組設(shè)備.具體數(shù)據(jù)采集與處理如圖2所示.

土壤溫度和土壤濕度由便攜式土壤參數(shù)速測(cè)儀進(jìn)行測(cè)定，經(jīng)緯度及高程由GPS 工具箱進(jìn)行定位和確定.

1.3 統(tǒng)計(jì)方法

多元線性回歸方程是對(duì)于給定的實(shí)際問題，找出某一個(gè)因素及對(duì)其有顯著影響的其他若干幾個(gè)因素，并通過從一組已知觀測(cè)值出發(fā)，運(yùn)用相關(guān)理論及統(tǒng)計(jì)軟件對(duì)其進(jìn)行方差分析后建立一個(gè)數(shù)學(xué)模型表達(dá)關(guān)系方式，進(jìn)而對(duì)該模型的顯著性、可信度和擬合優(yōu)度及標(biāo)準(zhǔn)誤差等進(jìn)行統(tǒng)計(jì)檢驗(yàn)［11］.

在該文中將凍土活動(dòng)層厚度作為因變量，土壤溫度、土壤濕度、高程和植被類型作為自變量，對(duì)其進(jìn)行相關(guān)性分析.原始方程如下：

式中：b0，b1，b2，…，bn為回歸系數(shù)，ω 為隨機(jī)誤差.

2 結(jié)果與分析

2.1 烏伊嶺區(qū)凍土活動(dòng)層厚度的分布特征

烏伊嶺區(qū)凍土活動(dòng)層厚度的整體變化范圍在0.4 ～6.21 m，最小凍土活動(dòng)層厚度出現(xiàn)在灌叢草甸中，最大凍土活動(dòng)層厚度出現(xiàn)在溪流邊.在野外調(diào)查采樣的10 組探地雷達(dá)樣線，120 個(gè)插釬點(diǎn)數(shù)據(jù)中，灌叢草甸的平均活動(dòng)層厚度最小為0.59 m，溪邊的平均活動(dòng)層厚度最深為2.69 m.且不同植被類型活動(dòng)層厚度差異明顯，以不同植被類型區(qū)平均活動(dòng)層厚度作為標(biāo)準(zhǔn)進(jìn)行對(duì)比，得出：興安杜鵑（0.97 m）＞小葉章（0.854 m）＞白樺（0.684 m）＞落葉松（0.645 m）＞薹草（0.634 m）＞篤斯越桔（0.569 m）＞杜香（0.453 m）.

從整體來看，凍土活動(dòng)層厚度與土壤濕度以及不同植被類型有著緊密的聯(lián)系，土壤濕度最大的區(qū)域活動(dòng)層厚度也是最大，即溪邊凍土活動(dòng)層厚度最深.而從植被類型來看，白樺林與落葉松林土壤濕度差異不大，因此他們的活動(dòng)層厚度差異也不太明顯.于此相反的是落葉松和興安杜鵑，兩者生長(zhǎng)區(qū)域土壤濕度差異不明顯，但他們活動(dòng)層厚度差異明顯.白樺和薹草生長(zhǎng)區(qū)域活動(dòng)層厚度差異不明顯，土壤濕度數(shù)據(jù)差異明顯，土壤溫度差異不明顯（見表1）.

典型植被 凍土活動(dòng)層厚度／m 土壤溫度／℃ 土壤濕度／% 高程／m白樺 0.684 20.08 13.87 368.13篤斯越桔 0.569 21.17 10.97 354.53薹草 0.634 19.0 20.79 366.33杜香 0.453 20.49 12.26 350.89落葉松 0.645 17.93 14.5 367.36興安杜鵑 0.97 17.95 15.44 393.57小葉章 0.854 20.29 14.17 376.13

2.2 凍土活動(dòng)層變化影響因素分析

2.2.1 烏伊嶺地區(qū)植被覆蓋類型與凍土活動(dòng)層變化影響分析

在采樣區(qū)域的選擇上根據(jù)凍土區(qū)優(yōu)先的原則，根據(jù)不同植被類型，選擇凍土分布的地區(qū)的典型植被，白樺（A1）、篤斯越桔（A2）、薹草（A3）、杜香（A4）、落葉松（A5）、興安杜鵑（A6）、小葉章作為不同的植被覆蓋類型.由于不同的植被類型在變量的屬性上屬于分類變量，因此需要先將變量轉(zhuǎn)換為虛擬變量再進(jìn)行多元線性回歸分析.又因?yàn)檫@7個(gè)變量的取值存在著互斥的關(guān)系，同時(shí)將所有虛擬變量引入模型會(huì)產(chǎn)生完全共線性問題.因此，只引入6個(gè)虛擬變量（A1～A6），剩下的一個(gè)分類變量作為參照組.由于小葉章作為優(yōu)勢(shì)種的區(qū)域活動(dòng)層厚度平均值為0.854m，僅次于興安杜鵑為優(yōu)勢(shì)種區(qū)域的活動(dòng)層厚度，且小葉章的分布范圍廣于興安杜鵑，故將小葉章作為此次分析的參照組.

由表2 可知，調(diào)整過的R2說明不同的植被類型能解釋凍土活動(dòng)層厚度變化的26.8%，表明不同的植被覆蓋類型與凍土活動(dòng)層厚度變化存在著線性關(guān)系，但不是影響凍土活動(dòng)層厚度變化的主要因素.其中德賓－沃森（DW）檢驗(yàn)是推斷殘差是否存在自相關(guān)的統(tǒng)計(jì)檢驗(yàn)方法，根據(jù)DW檢驗(yàn)表可得本分析dL、dU 值分別為1.404 和1.805，即dL＜DW＜dU，說明本次分析殘差獨(dú)立性檢驗(yàn)基本符合，樣本間獨(dú)立性檢驗(yàn)基本符合，但回歸關(guān)系不確定.

以小葉章作為標(biāo)準(zhǔn)變量進(jìn)行分析，杜香和篤斯越桔對(duì)凍土活動(dòng)層厚度的影響程度要高與小葉章，且與凍土活動(dòng)層厚度變化呈顯著相關(guān)（P ＝0.001／0.012＜0.05）.白樺、薹草、落葉松和興安杜鵑與凍土活動(dòng)層厚度變化無顯著相關(guān)性（顯著性均大于0.05），且以上4 種植被覆蓋類型與小葉章對(duì)凍土活動(dòng)層厚度的影響差距較大，但落葉松和薹草對(duì)凍土活動(dòng)層厚度影響效果較為相似.由此可知，杜香和篤斯越桔對(duì)凍土活動(dòng)層厚度變化產(chǎn)生的影響較大，呈顯著負(fù)相關(guān)，杜香對(duì)凍土活動(dòng)層厚度影響最大，而其余植被覆蓋類型與凍土活動(dòng)層厚度變化無顯著性，見表3.

模型 R R2 調(diào)整后R2 標(biāo)準(zhǔn)估算的錯(cuò)誤 德賓－沃森1 0.583a 0.340 0.268 23.83160 1.632 a.預(yù)測(cè)變量：（常量），A6，A4，A5，A3，A2，A1 b.因變量：活動(dòng)層

2.2.2 烏伊嶺地區(qū)土壤溫濕度和高程對(duì)凍土活動(dòng)層變化影響分析

該文選擇土壤濕度、土壤溫度和高程［4，12，13］作為影響因子，分析其對(duì)凍土活動(dòng)層厚度變化的影響.由表4 可知，調(diào)整后的R2為53.6%，說明所選環(huán)境因子土壤溫濕度和高程可以解釋凍土活動(dòng)層厚度變化的53.6%，是影響活動(dòng)層厚度的變化的主要因素，且呈線性關(guān)系.德賓－沃森（DW）檢驗(yàn)值從DW檢驗(yàn)表可知dU ＝1.689＜DW＜4－dU，即殘差序列不存在自相關(guān)，該回歸方程模型可以被充分解釋.由表5 可知，海拔會(huì)影響凍土活動(dòng)層厚度的變化，且與凍土活動(dòng)層厚度變化呈正比；而土壤溫度和土壤濕度則不會(huì)影響凍土活動(dòng)層厚度變化，土壤溫度和土壤濕度的顯著性均大于0.05.共線性統(tǒng)計(jì)包括方差膨脹因子（VIF）和容差兩個(gè)指標(biāo)，事實(shí)上，VIF ＝容差的倒數(shù)（1／容差），只需要判斷其中一個(gè)指標(biāo)即可.由圖3 可知，在本次分析中，所有容忍度值都大于0.1（最小值為0.632），說明該研究自變量多重共線不嚴(yán)重.該次分析殘差直方圖服從正態(tài)分布，標(biāo)準(zhǔn)差接近于1（標(biāo)準(zhǔn)正態(tài)分布），說明該次線性回歸符合正態(tài)性條件.

圖3 土壤溫濕度和高程與凍土活動(dòng)層厚度線性關(guān)系標(biāo)準(zhǔn)化殘差直方圖

通過上述分析，可以得到海拔與凍土活動(dòng)層厚度的多元線性回歸方程為：活動(dòng)層＝0.987×海拔－316.061，其關(guān)系圖如圖4所示.

模型 R R2 調(diào)整后R2 標(biāo)準(zhǔn)估算的錯(cuò)誤 德賓－沃森1 0.747a 0.559 0.536 18.98063 1.811 a.預(yù)測(cè)變量：（常量），濕度，海拔，溫度b.因變量：活動(dòng)層

圖4 活動(dòng)層厚度與海拔線性關(guān)系圖

2.2.3 烏伊嶺地區(qū)不同植被類型和土壤溫濕度與活動(dòng)層厚度的關(guān)系

由2.2.2的結(jié)論可知，土壤溫濕度與活動(dòng)層厚度變化無顯著相關(guān)性，但根據(jù)研究土壤溫濕度皆是影響凍土活動(dòng)層變化的重要因素［3，14－18］，不同植被類型的分布也與土壤溫濕度的變化息息相關(guān)［19－20］，因此該文試圖將植被類型作為控制變量，探究土壤溫濕度與凍土活動(dòng)層厚度的關(guān)系，運(yùn)用分層線性回歸的方式，對(duì)三者進(jìn)行分析，由于植被類型屬于分類變量，因此將其轉(zhuǎn)換為虛擬變量后代入方程進(jìn)行分析，分析結(jié)果見表6.

變量 模型1 模型2（常量）85.4 74.522 A1－17－15.601 A2－28.511*－24.682*A3－22.025－29.793**A4－40.114***－37.913***A5－20.9－21.862 A6 11.6 9.587溫度 —－0.25濕度 —1.125*R2 0.34 0.422 ΔR2 0.34 0.082 F 4.724*** 4.84***

由表6可知，模型1 是控制變量即植被類型對(duì)因變量?jī)鐾粱顒?dòng)層厚度的影響；模型2 是在控制變量（植被類型）的基礎(chǔ)上，加入因變量土壤溫濕度對(duì)凍土活動(dòng)層厚度的影響.其中模型1 的擬合度為0.34，說明植被類型作為控制變量對(duì)作為因變量的凍土活動(dòng)層厚度的解釋度為34%，F(xiàn) 值是4.724，在0.001水平上顯著，說明控制變量中有一個(gè)變量與凍土活動(dòng)層厚度之間存在影響關(guān)系.模型2是在模型1 的基礎(chǔ)上加入了土壤溫濕度作為自變量，以反應(yīng)土壤溫濕度對(duì)凍土活動(dòng)層厚度的影響.模型2的擬合度為0.442，說明在植被類型作為控制變量的條件下，土壤溫濕度對(duì)凍土活動(dòng)層厚度的解釋度為44.2%，F(xiàn) 值為4.84，在0.001 水平上顯著，說明該模型有意義，R2的變化量為0.082則說明該模型有效；由變量的回歸系數(shù)的顯著性可看出控制變量中的篤斯越桔、薹草和杜香對(duì)模型起到了干擾作用，因此在排除了這三種植被的干擾之后，自變量土壤濕度與凍土活動(dòng)層厚度變化呈正相關(guān)（B ＝1.125；P ＜0.05），而土壤溫度與凍土活動(dòng)層厚度變化無顯著相關(guān)性（P＞0.05）.

3 討論

3.1 植被類型與凍土活動(dòng)層影響關(guān)系分析

通過對(duì)不同植被類型的分析表明，不同植被類型對(duì)凍土活動(dòng)層厚度變化有較大影響，尤其是杜香和篤斯越桔，且杜香作為優(yōu)勢(shì)種的區(qū)域，其凍土活動(dòng)層厚度明顯小于篤斯越桔.杜香和篤斯越桔分別作為優(yōu)勢(shì)種的區(qū)域，活動(dòng)層厚度均小于其他與凍土活動(dòng)層厚度變化無顯著相關(guān)性的植被區(qū)域.而現(xiàn)有研究多為植被NDVI 或植被群落類型與凍土活動(dòng)層厚度變化之間的關(guān)系，少有將具體植被種類與凍土活動(dòng)層厚度變化情況進(jìn)行分析［21－22］.因此，白樺、薹草以及興安杜鵑等凍土區(qū)指示植物與凍土活動(dòng)層厚度變化無顯著相關(guān)性的原因還需進(jìn)一步探究.具體原因有可能是未將植被的葉面積指數(shù)、積雪厚度、苔蘚層厚度、土壤導(dǎo)熱系數(shù)等因素納入分析.

3.2 土壤溫濕度和高程對(duì)凍土活動(dòng)層厚度影響分析

研究表明，非地形因子土壤水分和溫度［23］和地形因子高程共同作用于凍土活動(dòng)層，令活動(dòng)層厚度發(fā)生改變.該文分析表明，高程與凍土活動(dòng)層厚度呈顯著正相關(guān).這與學(xué)者基于整個(gè)東北的區(qū)域尺度上得出的小興安嶺高海拔有利于多年凍土的發(fā)育和保存［24］.兩種研究結(jié)果不同，可能是由于研究的尺度不同所造成的.相較于整個(gè)東北地區(qū)，烏伊嶺區(qū)受到全球變暖以及其他因素的影響更為劇烈.

研究證實(shí)，土壤濕度越高會(huì)導(dǎo)致活動(dòng)層厚度越淺，也有研究證實(shí)土壤濕度越高會(huì)導(dǎo)致活動(dòng)層厚度越深.兩種結(jié)果均表明土壤濕度是導(dǎo)致活動(dòng)層變化的重要因素.該文將土壤溫濕度與高程作為影響因子進(jìn)行相關(guān)性分析時(shí)，土壤溫濕度均與活動(dòng)層厚度變化無顯著相關(guān)性，而植被類型與活動(dòng)層厚度具有顯著相關(guān)性，故而將植被類型作為控制變量，對(duì)土壤溫濕度進(jìn)行分析.得出結(jié)論：土壤濕度與活動(dòng)層厚度呈顯著正相關(guān)，土壤溫度與活動(dòng)層厚度無顯著相關(guān)性.這與杜二計(jì)的研究結(jié)果不同，可能是由于研究地區(qū)不同導(dǎo)致的.青藏高原由于地處中低緯度且平均海拔在4000 m 以上，有強(qiáng)烈的地表輻射［25－26］.烏伊嶺區(qū)地處小興安嶺，緯度高且植被種類繁多，所以土壤濕度對(duì)活動(dòng)層厚度加深有促進(jìn)作用.

4 結(jié)論

該文基于伊春市烏伊嶺凍土區(qū)的凍土活動(dòng)層厚度、植被類型、土壤溫濕度以及高程的實(shí)測(cè)數(shù)據(jù)，對(duì)實(shí)測(cè)數(shù)據(jù)進(jìn)行多元線性回歸和分層線性回歸分析，探究了凍土活動(dòng)層厚度的影響因素，結(jié)論如下：

（1）與小葉章相比，杜香和篤斯越桔對(duì)凍土活動(dòng)層厚度的影響更大，且均為負(fù)相關(guān)，而其他植被類型對(duì)凍土活動(dòng)層厚度變化無顯著影響；高程與凍土活動(dòng)層厚度的變化呈極顯著正相關(guān)，土壤溫度和土壤濕度與凍土活動(dòng)層厚度變化無顯著關(guān)系.

（2）將植被類型作為控制變量，在植被類型一定的情況下土壤濕度與凍土活動(dòng)層厚度變化顯著相關(guān).