張 威，趙 賀

（遼寧師范大學(xué)地理科學(xué)學(xué)院，遼寧大連 116029）

0 引言

冰川槽谷是冰川侵蝕作用下形成的一種典型地貌。研究冰川槽谷橫剖面的形態(tài)特征有助于區(qū)分不同谷地的成因以及分析冰川的侵蝕過程和侵蝕能力［1-5］。近幾十年來，針對青藏高原及其周邊山地冰川槽谷橫剖面形態(tài)的研究已取得了階段性進(jìn)展，并且相關(guān)研究已進(jìn)入定量化階段［5-10］。崔之久［11］在研究青藏高原主要山地冰川侵蝕地貌特征時指出，僅有唐古拉山和祁連山的冰川槽谷橫剖面屬于寬淺型（槽谷的深寬比超過1∶6）。李吉均［12］在研究橫斷山脈冰川侵蝕地貌時，指出此山地共保存4種槽谷剖面形態(tài)，分別是普通“U”形、深窄“U”形、懸鏈“U”形和“箱形”，并將前3種剖面形態(tài)用函數(shù)形式進(jìn)行了表達(dá)。后續(xù)學(xué)者［5，8，13］應(yīng)用冪函數(shù)模型（y=axb）［14］分別對天山、白馬雪山和螺髻山等地的部分冰川槽谷形態(tài)進(jìn)行了詳細(xì)的分析。近期，Zimmer等［15］在研究美國內(nèi)華達(dá)山脈第四紀(jì)冰川侵蝕地貌特征時，提出了定量化表達(dá)冰川作用谷地形態(tài)的V指數(shù)模型，并證明此方法在描述冰川作用谷地形態(tài)時較冪函數(shù)法、二次多項式法準(zhǔn)確率更穩(wěn)定、適用的槽谷形態(tài)更多。此前，筆者已將此方法應(yīng)用于他念他翁山脈的冰川槽谷橫剖面形態(tài)特征分析，取得了較好的結(jié)果［16］。

冰川侵蝕作用谷地具有多種形態(tài)特征，除了最常見的冰川作用谷地為“U”形谷地外，還發(fā)育箱型、倒梯形等谷地形態(tài)，在以往的研究中，研究者大多集中于對偏“U”形谷地形態(tài)特征進(jìn)行定性與定量化分析，而對其他形態(tài)的谷地形態(tài)特征的定量化研究很少；同時針對唐古拉山寬淺型槽谷的具體表現(xiàn)也可以進(jìn)一步深化。因此，運(yùn)用最新的V指數(shù)模型，對于冰川槽谷形態(tài)進(jìn)行定量分析，并探討其影響因素仍需進(jìn)一步深入。目前，探討影響冰川槽谷剖面形態(tài)特征的主要因素時，大多采用較為傳統(tǒng)的皮爾遜相關(guān)性分析法等［16-18］。這些方法的解釋能力有限，且僅能分析各個因素的單獨(dú)作用，而忽略了這些因素之間的相互作用。地理探測器可以用來探測地理屬性的空間差異性，并揭示其背后的驅(qū)動因素［19］，它的優(yōu)點(diǎn)是：（1）探測數(shù)據(jù)類型可以是定性的，也可以是定量的；（2）可以探測兩因素交互作用于因變量。此方法在多方面地理研究中都已有較好的應(yīng)用［20-21］。

作為青藏高原中部重要的冰川發(fā)育區(qū)，唐古拉山是研究冰川槽谷形態(tài)特征的理想山地。另外，由于區(qū)域內(nèi)地形、冰川發(fā)育及巖性等特征分異顯著，加之地理位置的特殊性，使得本研究區(qū)冰川槽谷橫剖面特征存在一定的空間異質(zhì)性。基于上述研究背景，本文選取唐古拉山中西段形態(tài)典型的冰川槽谷作為研究對象，基于V指數(shù)模型和地理探測器方法，分析研究區(qū)冰川槽谷的形態(tài)特征，以及影響冰川槽谷發(fā)育的主控因素。分析結(jié)果有助于了解大陸性冰川的發(fā)育特點(diǎn)、侵蝕程度以及寬淺型槽谷的具體特征。

1 研究區(qū)概況

唐古拉山脈位于青海省和西藏自治區(qū)的邊界，本區(qū)新生代發(fā)育受控于印度-歐亞板塊的大型逆沖推覆構(gòu)造帶，走向與山脈延伸方向一致，呈北西西向展布［22］，整個山脈的地勢特點(diǎn)是由西向東逐漸降低［23］，海拔由6 621 m過渡到4 000 m以上（圖1）。唐古拉山還是印度夏季風(fēng)北上的第二大屏障，也是我國海洋性冰川向大陸性冰川過渡的重要地帶［23］。唐古拉山中西段（31°91′～33°40′N，90°83′～93°43′E）地勢平緩，少數(shù)山頭海拔達(dá)6 100 m以上，現(xiàn)代冰川發(fā)育面積可達(dá)678 km2［23］。根據(jù)冰川發(fā)育類型可將本研究區(qū)分為東部（A區(qū)）、中部（B區(qū)）和西部（C區(qū)）（圖1）［23-24］。其中，東部是山谷冰川發(fā)育區(qū)；中部是山谷冰川向冰帽冰川過渡區(qū)，此區(qū)的冰川發(fā)育兼?zhèn)淞吮北ê蜕焦缺ǖ奶攸c(diǎn)；西部是冰帽冰川發(fā)育區(qū)。槽谷主要由冰帽冰川和山谷冰川侵蝕形成，其中一些主谷的兩側(cè)發(fā)育冰川懸谷，形成了大小各異的復(fù)合型山谷冰川。冰川性質(zhì)包括極大陸性冰川和亞大陸性冰川［25］。本區(qū)共發(fā)生六次冰川作用，包括倒數(shù)第二次冰期、末次冰期早冰階、中冰階、晚冰階以及全新世的新冰期和小冰期［26］。研究區(qū)屬于典型的季風(fēng)山地氣候，主要受太平洋東南季風(fēng)和印度洋西南季風(fēng)的影響，年降水量約為400 mm［26］。區(qū)內(nèi)出露石炭紀(jì)、二疊紀(jì)、侏羅紀(jì)、古近紀(jì)和新近紀(jì)地層，巖石類型主要為花崗巖、碳酸鹽巖、礫巖、頁巖、泥巖等［27］。

圖1 唐古拉山中西段研究區(qū)范圍與創(chuàng)建的冰川槽谷橫剖面示例圖（A：冰帽作用區(qū)；B：冰帽冰川至山谷冰川過渡作用區(qū)；C：山谷冰川作用區(qū)）Fig.1 Example map of study area extent and created glacial valley cross section in the middle and west Tanggula Mountain（A：ice cap action area；B：ice cap glacier to valley glacier transition action area；C：valley glacier action area）

2 材料與方法

2.1 數(shù)據(jù)來源與提取

以ArcGIS10.6作為技術(shù)平臺，用于生成河網(wǎng)、提取谷地橫截面、獲取谷底高程數(shù)據(jù)集，在工作空間中使用WGS-84坐標(biāo)系處理數(shù)據(jù)。本文使用的GDEMV2 30 m分辨率的數(shù)字高程數(shù)據(jù)可以從地理空間數(shù)據(jù)云網(wǎng)站（https：//www.gscloud.cn/）下載。根據(jù)現(xiàn)場實際考察和前人的研究成果并結(jié)合谷歌地圖，選取研究區(qū)受后續(xù)河流影響較小、谷地填充物較少且發(fā)育典型的冰川槽谷91條。由于末次冰期以來的冰川槽谷形態(tài)特征清晰，受后期其他營力影響較小，故本文以末次冰期以來的冰川作用末端作為邊界。部分槽谷的末端可以通過具有年代學(xué)的冰川遺跡確定；另一些槽谷可以首先根據(jù)CF法和TSAM法重建末次冰期冰川物質(zhì)平衡線高度，然后再應(yīng)用THAR法即可確定槽谷的末端。首先根據(jù)冰川地貌基本特征確定每條冰川槽谷的大致范圍，然后基于每條槽谷內(nèi)的DEM生成河網(wǎng)矢量數(shù)據(jù)集，為了準(zhǔn)確地提取流域水系，在谷歌地圖上對水系進(jìn)行了反復(fù)對照，最終選擇1 000為閾值。之后使用ArcMap中的ET Geowizards站線工具，即可自動生成指定寬度且垂直于河網(wǎng)的橫剖面線。本文選擇的寬度是2 km，因為這個寬度可以最大程度地覆蓋研究區(qū)內(nèi)所有山谷類型?？紤]到研究區(qū)冰川作用谷地特征的空間異質(zhì)性，每條槽谷生成一條位于槽谷中間位置的剖面線（圖1）。最后基于DEM數(shù)據(jù)提取每條橫截面上的高程數(shù)據(jù)，然后將數(shù)據(jù)導(dǎo)入MATLAB腳本中計算V指數(shù)。其中，腳本運(yùn)行前要求輸入的指定高度根據(jù)每條谷地實際谷肩距谷底的高程設(shè)定。本文的谷肩位置通過谷歌地圖判讀結(jié)合野外考察獲取。

2.2 V指數(shù)模型

山區(qū)河流和冰川侵蝕機(jī)制的差異會造成兩者谷地橫剖面形態(tài)的不同，通常河谷呈“V”形，而冰川谷呈“U”形。該模型據(jù)此研究了谷底上方特定冰川作用深度的谷面面積（Ax）與具有相同深度和寬度的理想“V”形截面面積（Av）間的關(guān)系。V指數(shù)（V）計算公式為：

V指數(shù)描述谷地形態(tài)與“V”形谷的偏差程度，理想的“V”形谷的V指數(shù)為0，“U”形谷的V指數(shù)>0；而凸谷的V指數(shù)<0（圖2）。

圖2 V指數(shù)和相應(yīng)的谷形的對應(yīng)關(guān)系Fig.2 Correspondence between V and corresponding valley shape

2.3 V指數(shù)與b值之間的關(guān)系

假設(shè)冰川槽谷某一側(cè)橫剖面的形態(tài)為半“U”形，根據(jù)冪函數(shù)模型［14］和V指數(shù)模型的計算公式，這個半“U”形的V指數(shù)可以表示為：

進(jìn)而得到V指數(shù)與b值的關(guān)系式：

如果槽谷剖面的兩側(cè)是對稱的，那么這個完整槽谷橫剖面的V指數(shù)與冪函數(shù)b值之間的關(guān)系就可用公式（3）表示，且兩者呈正相關(guān)關(guān)系。另外，典型的冰川槽谷橫剖面b值主要集中在1.5～2.5［6-7，28-29］，那么對應(yīng)的典型冰川槽谷橫剖面V指數(shù)應(yīng)介于0.20～0.43。即V指數(shù)小于0.20的槽谷呈現(xiàn)由“V”形谷向典型“U”形谷過渡的形態(tài)，侵蝕程度?。唤橛?.20～0.43之間屬于具有典型“U”形形態(tài)的冰川槽谷；而大于0.43的槽谷較典型“U”形谷形態(tài)更復(fù)雜。在實地考察中發(fā)現(xiàn)V指數(shù)大于0.43的槽谷一般谷底寬闊，谷壁與谷底有夾角，有些甚至接近90°，且谷壁弧度不明顯，總體呈現(xiàn)出“倒梯形”或近似“箱形”的形態(tài)。雖然部分槽谷形態(tài)存在不完全對稱的情況，但一般相差不大，所以在實際計算中即可用上述關(guān)系將V指數(shù)與槽谷剖面特征進(jìn)行對應(yīng)。

2.4 地理探測器

分異探測：當(dāng)某個影響因素對某一地理特征影響較大時，那么這個因素X和地理特征Y的空間分布應(yīng)該呈現(xiàn)出相似性，具體結(jié)果q通過GeoDetector軟件可以獲?。?9］。此模型的表達(dá)式是：

式中：L是X和Y的分層，也就是分類或分區(qū)；Ni和N分別為層h和全區(qū)的單元數(shù)；σ2i和σ2分別是X、Y的方差。q的取值范圍為［0，1］，q值越大表示影響因素X對地理屬性Y的解釋力越強(qiáng)。極端情況下，q值為1表明因子X完全控制了Y的空間分布，q值為0則表明因子X與Y沒有任何關(guān)系。

交互探測：用于識別各個影響因素Xs間的交互作用，即評估在因素X1和X2的共同作用下對地理特征Y的解釋能力是增大或減小，或這些因素對Y的影響是相互獨(dú)立的。具體方法是先分別計算X1和X2因素對Y的q值：q（X1）、q（X2），其次計算兩個因素交互下的q值：q（X1∩X2），最后把q（X1）、q（X2）與q（X1∩X2）進(jìn)行比較。交互探測關(guān)系包含以下幾種：（1）當(dāng)存在關(guān)系q（X1∩X2）>q（X1）+q（X2）時表示非線性增強(qiáng)；（2）當(dāng)存在關(guān)系q（X1∩X2）=q（X1）+q（X2）時表示相互獨(dú)立；（3）當(dāng)存在關(guān)系q（X1∩X2）>max［q（X1），q（X2）］時表示雙因子增強(qiáng)；（4）當(dāng)存在關(guān)系min［q（X1），q（X2）］

2.5 冰川槽谷形態(tài)空間分異的影響因素

根據(jù)冰川的發(fā)育和運(yùn)動特征，影響冰川槽谷形態(tài)特征的因素主要有地形因素、冰川動力因素、冰川性質(zhì)因素、巖性因素和構(gòu)造因素［4］。其中，地形因素主要是指坡度、地形起伏度和槽谷朝向；冰川動力因素是指冰川作用區(qū)面積和冰川作用正差。具體而言，坡度是指選取的每個橫剖面位置處的坡度值；地形起伏度是指選取橫剖面位置以上山地的最高點(diǎn)海拔高度與最低點(diǎn)海拔高度的差值；本文的槽谷朝向是指每條槽谷選取的中間位置橫剖面的方向，采用0°及360°為正北方向，-22.5°<β≤22.5°用N表示，為北向；22.5°<β≤67.5°用NE表示，為東北向；以此類推順時針依次用E、SE、S、SW、W、NW表示其他槽谷朝向。研究區(qū)槽谷朝向以南、北方向為主，東、西方向最少（圖3）。由于冰川槽谷的相關(guān)ELA數(shù)據(jù)獲取較易，且計算更便捷，所以冰川作用區(qū)面積可由冰川物質(zhì)平衡線（ELA）以上的冰川發(fā)育面積（積累區(qū)面積）替代［17］，其中，末次冰期平衡線高度通過CF法和TSAM法重建，均值約為5 200 m，與前人的研究結(jié)果相似［23］；冰川作用正差是指山地最大海拔高度與ELA海拔高度之差［17］。如圖4所示，研究區(qū)內(nèi)冰川性質(zhì)分界線以西是極大陸型冰川，以東是亞大陸型冰川［25］。區(qū)域內(nèi)選取的槽谷剖面覆蓋的巖性主要包括七大類（圖4）。而構(gòu)造特征無法轉(zhuǎn)換成數(shù)據(jù)，所以本文未利用地理探測器探測此因素，僅用文字討論。

圖3 冰川槽谷朝向玫瑰圖Fig.3 Glacial valley heading rose map

圖4 唐古拉山中西部V指數(shù)與巖性分布圖［27］Fig.4 V and lithological distribution map in the middle and west of Tanggula Mountain［27］

3 結(jié)果分析

3.1 冰川槽谷橫剖面形態(tài)分布規(guī)律

研究區(qū)冰川槽谷長2.5～12.5 km，平均寬度1.8 km，深100～200 m，海拔介于5 690～4 710 m。計算結(jié)果表明，研究區(qū)冰川槽谷橫剖面V指數(shù)的最小值為0.04，最大值為0.95，均值為0.42。各谷地剖面的V指數(shù)數(shù)值都>0，其中，V指數(shù)<0.20的有17條，占19%；V指數(shù)介于0.20～0.43之間的有44條，占48%；V指數(shù)>0.43的有30條，占33%（圖5）。

圖5 唐古拉山中西段V指數(shù)箱型圖Fig.5 V box diagram of the middle and west section of Tanggula Mountain

分別統(tǒng)計北坡和南坡的冰川槽谷剖面V指數(shù)以分析其發(fā)育特征的差異，選取的冰川槽谷位于北坡的有52條，位于南坡的有39條。比較每條冰川槽谷V指數(shù)發(fā)現(xiàn)，無論北坡還是南坡槽谷發(fā)育特征差異都較大，北坡冰川槽谷V指數(shù)<0.20的有6條，占12%；V指數(shù)介于0.20～0.43之間的有31條，占60%；V指數(shù)>0.43的有15條，占28%。南坡冰川槽谷V指數(shù)<0.20的有11條，約占29%；V指數(shù)介于0.20～0.43之間的有13條，占33%；V指數(shù)>0.43的有15條，占38%（圖4、圖6）。由此可見，北坡的槽谷主要呈現(xiàn)出典型的“U”形形態(tài)，其次是近似“箱形”的形態(tài)，少部分呈現(xiàn)由“V”形向典型“U”形過渡的形態(tài)；而南坡的槽谷三種形態(tài)發(fā)育數(shù)量相差不大。其中，北坡發(fā)育的冰川槽谷V指數(shù)最小值為0.10，最大值為0.95，均值為0.42；南坡發(fā)育的冰川槽谷V指數(shù)最小值為0.04，最大值為0.95，均值為0.47。

選取的槽谷中位于東部、中部和西部的數(shù)量分別為37條、31條和23條。東部冰川槽谷V指數(shù)<0.20的有10條，占27%；V指數(shù)介于0.20～0.43之間的有20條，占54%；V指數(shù)>0.43的有7條，占19%。中部冰川槽谷V指數(shù)<0.20的有5條，占16%；V指 數(shù) 介 于0.20～0.43之 間 的 有19條，占61%；V指數(shù)>0.43的有7條，占23%。西部冰川槽谷V指數(shù)<0.20的有2條，占9%；V指數(shù)介于0.20～0.43之間的有5條，占21%；V指數(shù)>0.43的有16條，占70%（圖4、圖6）。其中，東部冰川槽谷的V指數(shù)最小值為0.04，最大值為0.95，均值為0.36；中部的冰川槽谷V指數(shù)最小值為0.04，最大值為0.95，均值為0.40；西部的冰川槽谷V指數(shù)最小值為0.04，最大值為0.95，均值為0.64。上述結(jié)果表明，冰帽冰川發(fā)育區(qū)的槽谷V指數(shù)均值在三個區(qū)域中最大，且槽谷以趨近“箱形”形態(tài)為主，侵蝕程度較強(qiáng)；而山谷冰川作用區(qū)和過渡區(qū)的槽谷V指數(shù)均值相似，且槽谷形態(tài)都以典型“U”形為主，侵蝕程度較弱。呈現(xiàn)上述現(xiàn)象的原因可能與地形因素、冰川動力因素密切相關(guān)。具體而言，西部夷平面廣布，冰帽冰川發(fā)育，在這種地形影響下冰川物質(zhì)平衡線的少量下降也會造成冰川的強(qiáng)烈擴(kuò)張。因此，此區(qū)冰川規(guī)模較大，冰川動力較強(qiáng)，冰川侵蝕能力更強(qiáng)，槽谷形態(tài)更接近于“箱形”。而東部和中部地形切割程度較西部強(qiáng)，所以此區(qū)冰川的擴(kuò)張不如西部強(qiáng)烈，冰川侵蝕能力也不如西部，槽谷形態(tài)更接近于典型“U”形。除此之外，箱形槽谷的形成可能也與冰前地形的非“V”形形態(tài)有關(guān)。

圖6 南北坡V指數(shù)分布特征（a）和東、中、西三部V指數(shù)分布特征（b）Fig.6 V distribution characteristics of north-south slope（a）and V distribution characteristics of east，middle and west slope（b）

3.2 槽谷橫剖面形態(tài)的影響因素分析

地理探測器可以分析各因素的空間分異程度。由表1可知，各個因素在不同程度上解釋了槽谷形態(tài)的空間差異。地理探測器q值顯示的各因素與槽谷形態(tài)間的相關(guān)性程度如下：冰川作用區(qū)面積（0.3695）>冰川作用正差（0.2143）>巖性（0.1816）>坡度（0.1784）>地形起伏度（0.1221）>冰川性質(zhì)（0.0973）>槽谷朝向（0.0680）。因此，上述定量的因素中冰川作用區(qū)面積和冰川作用正差因素是本研究區(qū)冰川槽谷形態(tài)空間分異最主要的控制因素，其次是巖性因素，再次是坡度和地形起伏度因素，最后是冰川性質(zhì)和槽谷朝向因素。另外，q值對應(yīng)的p值代表了每個因子的顯著性。一般而言，如果p值越小，說明因素X對因變量Y的影響越顯著。但有些數(shù)據(jù)總體差異較?。ㄈ缙露群偷匦纹鸱龋?，分區(qū)后平均值差異不大也會造成p值較大的情況。

表1 地理探測器各因素探測結(jié)果Table 1 Detection results of each factor of geographical detector

冰川作用區(qū)面積和冰川作用正差因素是本研究區(qū)冰川槽谷形態(tài)空間分異最主要的控制因素，具體而言：研究區(qū)內(nèi)冰川作用面積越大、冰川動力越強(qiáng)時，槽谷的形態(tài)越接近典型“U”形或“箱形”。而巖性因素也是研究區(qū)內(nèi)槽谷形態(tài)變化不可忽視的因素，具體表現(xiàn)在：區(qū)內(nèi)的巖石類型中，花崗巖與閃長巖的硬度最強(qiáng)，其次是碳酸鹽巖，再次是礫巖、砂巖和頁巖，而泥巖等最弱。從代表區(qū)內(nèi)冰川槽谷剖面形態(tài)的V指數(shù)結(jié)果來看，冰川槽谷在不同巖性處的形態(tài)表現(xiàn)基本符合上述侵蝕規(guī)律。針對地形因素（包括坡度和地形起伏度因素）對區(qū)內(nèi)冰川槽谷形態(tài)的影響可以解釋為：因本文選取地形數(shù)據(jù)是指每條槽谷中的地形，也就是槽谷沿程的坡度和地形起伏度，此時的地形會通過影響基底剪切應(yīng)力從而改變冰底滑動速率的方式改變冰川槽谷形態(tài)［30-31］。

一般來說，極大陸型冰川的侵蝕能力比亞大陸型冰川弱［25］，但研究區(qū)內(nèi)槽谷V指數(shù)普遍較大的冰帽冰川發(fā)育區(qū)（A區(qū)）屬于極大陸型冰川，而槽谷V指數(shù)普遍較小的山谷冰川發(fā)育區(qū)（C區(qū)）屬于亞大陸型冰川，這與侵蝕規(guī)律不符，所以冰川性質(zhì)不是槽谷形態(tài)的主控因素。此外，槽谷朝向也不是本研究區(qū)槽谷形態(tài)的主控因素，因為槽谷總共有八個朝向，每個朝向的V指數(shù)均值不存在規(guī)律，且每個朝向均保存了多種形態(tài)的谷地。構(gòu)造活動主要通過降低冰川物質(zhì)平衡線、改變山地坡度、影響地?zé)嵬康确矫嬗绊懕ú酃绕拭孀兓?5，32］。青藏高原雖然處于不斷的隆升過程［33］，但整體性的抬升并不會造成小范圍研究區(qū)的冰川發(fā)育差異。并且根據(jù)前人的調(diào)查結(jié)果，唐古拉山脈整體處于同一個地層分區(qū)中，受大斷裂的影響較?。?4］，山脈內(nèi)斷層也較少，所以構(gòu)造因素在本研究區(qū)冰川槽谷形態(tài)特征方面的影響可以忽略不計。綜上所述，影響本研究區(qū)冰川槽谷形態(tài)差異最主要的因素是冰川作用區(qū)面積和冰川作用正差因素，其次是巖性因素，再次是坡度和地形起伏度因素；而冰川性質(zhì)因素、槽谷朝向因素和構(gòu)造因素對本研究區(qū)槽谷形態(tài)的影響可以忽略不計。

實際上，冰川槽谷形態(tài)特征往往是由多個因素共同作用的結(jié)果，所以為了進(jìn)一步確定不同因素的組合作用對槽谷形態(tài)的解釋力是加強(qiáng)還是削弱，需要進(jìn)行交互探測?；诒?～2中q值結(jié)果，冰川作用區(qū)面積與冰川作用正差、冰川作用區(qū)面積與冰川性質(zhì)、冰川作用區(qū)面積與巖性、冰川作用正差與冰川性質(zhì)、冰川作用正差與巖性、冰川性質(zhì)與巖性的因素組合關(guān)系都屬于雙因子增強(qiáng)，而其他兩兩因素的組合交互作用都大于其疊加作用，屬于非線性增強(qiáng)關(guān)系。另外，根據(jù)探測結(jié)果，冰川作用區(qū)面積因素和坡度因素的交互作用對冰川槽谷形態(tài)特征的影響最大。這兩因素可以通過既增大冰川規(guī)模又增強(qiáng)冰川基底滑動速率的方式使得冰川槽谷形態(tài)更趨近于典型“U”形或“箱形”。

目前，已經(jīng)對青藏高原部分山地的冰川槽谷形態(tài)特征及其影響因素做出了詳細(xì)的研究［16-18］。這些槽谷發(fā)育的冰川性質(zhì)包括海洋性冰川、亞大陸性冰川和極大陸性冰川。結(jié)合本文研究結(jié)果來看，無論哪種冰川性質(zhì)的槽谷形態(tài)都可能受到坡度因素和地形起伏度的影響。而從緯向分布的角度，這些槽谷分布在青藏高原的中高緯和中低緯度。其中，分布于中低緯度的槽谷形態(tài)主要受冰川規(guī)模的影響，而冰川規(guī)模的增大通常通過冰川作用區(qū)面積和冰川作用正差等因素的改變而實現(xiàn)，這與本文得到的冰川作用區(qū)面積和冰川作用正差因素是影響唐古拉山中西段冰川槽谷形態(tài)變化的主要因素結(jié)論相符。

表2 地理探測器的交互探測Table 2 Interactive detection of geographic detectors

4 對比分析

本文V指數(shù)與Zimmer在內(nèi)華達(dá)山脈的研究［15］和筆者在他念他翁山的研究［16］有較大的差異。內(nèi)華達(dá)山脈存在V指數(shù)>0.43的槽谷但最大值是0.67，他念他翁山也存在V指數(shù)>0.43的槽谷但最大值是0.53，而本文首次計算得到有V指數(shù)結(jié)果接近于1。那么V指數(shù)近似于1時代表了怎樣的冰川槽谷形態(tài)呢？首先，根據(jù)上述討論可知V指數(shù)>0.43的槽谷一般呈現(xiàn)出“倒梯形”或“箱形”形態(tài)。其次，基于V指數(shù)模型的計算原理［式（1）］，當(dāng)計算得到V指數(shù)結(jié)果接近于1時，面積（Ax）的大小應(yīng)為面積（Av）大小的近2倍，所以此時的冰川槽谷剖面形態(tài)接近于矩形或者說“箱形”。李吉均［12］在研究橫斷山脈范圍內(nèi)的貢嘎山和稻城冰帽冰川侵蝕地貌時也指出，這兩個山地保存著“箱形”形態(tài)的冰川槽谷，這與本文得到的結(jié)果相符。而唐古拉山保存的大多數(shù)“寬淺型”槽谷可能就是“箱形”形態(tài)。從V指數(shù)的結(jié)果來看，箱形槽谷相對于冰前河流“V”形谷偏差程度較大，可能代表了侵蝕程度較強(qiáng)。另一方面箱形槽谷的形成可能也與冰前地形的非“V”形形態(tài)密切相關(guān)。另外，根據(jù)“箱形”槽谷的形態(tài)特征和模型計算原理，V指數(shù)模型可以更好地識別“箱形”槽谷。因為冪函數(shù)、二次多項式等模型在計算時的假設(shè)前提就是槽谷橫剖面要呈拋物線或“U”形形態(tài)，而V指數(shù)模型只需輸入相應(yīng)的數(shù)據(jù)直接計算，不用考慮此假設(shè)前提。

5 結(jié)論

（1）冰川槽谷橫剖面的V指數(shù)與冪函數(shù)b值之間可以相互替代，且典型冰川槽谷橫剖面的V指數(shù)介于0.20～0.43之間。

（2）唐古拉山中西段保存著“箱形”形態(tài)的冰川槽谷，“箱形”形態(tài)的槽谷可以通過V指數(shù)（>0.43）進(jìn)行識別和定量描述。

（3）根據(jù)V指數(shù)計算結(jié)果，本研究區(qū)選取的冰川槽谷呈典型“U”形的數(shù)量最多，其次是近似于“箱形”的，最少的是由“V”形向典型“U”形過渡的槽谷。北坡的槽谷形態(tài)主要呈現(xiàn)典型“U”形，其次是近似“箱形”的形態(tài)，少部分呈現(xiàn)由“V”形向典型“U”形過渡的形態(tài)；而南坡的槽谷三種形態(tài)發(fā)育的數(shù)量基本相同。冰帽冰川發(fā)育區(qū)的槽谷以近似“箱形”為主；而山谷冰川作用區(qū)和過渡區(qū)的槽谷都以典型“U”形為主。

（4）冰川作用區(qū)面積和冰川作用正差因素是造成冰川槽谷形態(tài)差異的主要原因，而從交互關(guān)系來看冰川作用區(qū)面積因素與坡度因素的交互作用對冰川槽谷形態(tài)特征的影響最大。其他各個因素（地形因素、冰川動力因素、冰川性質(zhì)因素、巖性因素和構(gòu)造因素）對唐古拉中西段冰川槽谷形態(tài)的空間分異有著不同程度的影響。