高澤民, 劉興旺,2, 邵延秀,2, 鄒小波

(1. 蘭州地球物理國家野外科學(xué)觀測研究站, 甘肅 蘭州 730000;2. 中國地震局蘭州巖土地震研究所, 甘肅 蘭州 730000)

0 引言

新生代以來印度板塊與歐亞板塊的匯聚碰撞致使青藏塊體不斷隆升,與此同時(shí)在塊體北緣伴生發(fā)育了阿爾金斷裂帶、北祁連斷裂帶等大型走滑—逆沖型構(gòu)造帶,共同構(gòu)成了高原北部邊界[1-2]。其中三危山斷裂位于青藏高原以北前緣隆升地帶,是阿爾金斷裂帶的重要伴生斷裂[3]。活動(dòng)斷裂運(yùn)動(dòng)學(xué)和古地震研究表明,三危山斷裂運(yùn)動(dòng)方式以逆沖為主兼具走滑,局部伴有正斷特征[4],斷裂東北段存在多次古地震事件[5]。對(duì)于其最新活動(dòng)性,早期研究認(rèn)為該斷裂活動(dòng)時(shí)間為早、中更新世,晚更新世以來不活動(dòng)[6]。近些年研究指出,斷裂為一條晚更新世活動(dòng)斷裂,在晚更世晚期或全新世早期也有活動(dòng)[4-5]。然而,從地貌學(xué)角度將構(gòu)造地貌分析方法應(yīng)用到三危山地區(qū)活動(dòng)構(gòu)造研究中的實(shí)例并不多,對(duì)三危山斷裂不同段構(gòu)造活動(dòng)強(qiáng)度量化及其對(duì)地貌演化的影響問題仍缺乏研究。

采用高分辨率地形數(shù)據(jù)從地貌形態(tài)中提取定量化地形參數(shù)信息分析構(gòu)造活動(dòng)特征已成為地貌和構(gòu)造研究領(lǐng)域的重要發(fā)展方向,地貌參數(shù)能夠反映構(gòu)造活動(dòng)對(duì)區(qū)域河流發(fā)育及演化過程的影響[7]。在構(gòu)造運(yùn)動(dòng)活躍的隆升區(qū),長時(shí)間跨度下平均抬升速率的差異可由地貌參數(shù)有效的指示[8]。通過分析河流水系、地形特征與地貌參數(shù)值間的相關(guān)性,可揭示造山帶水系格局變化對(duì)構(gòu)造活動(dòng)的響應(yīng)[9]。由地貌參數(shù)指示的構(gòu)造活動(dòng)強(qiáng)度與斷裂實(shí)地調(diào)查的活動(dòng)性強(qiáng)度間存在很好的一致性[10],因此,本文嘗試從地貌學(xué)角度出發(fā),利用高程、坡度、河道陡峭指數(shù)與河道坡降指數(shù)方法對(duì)三危山地區(qū)地貌特征進(jìn)行定量化分析,以期揭示三危山斷裂不同段地貌發(fā)育特征與斷裂活動(dòng)性差異的關(guān)系。

1 區(qū)域概況

三危山所在的瓜州—敦煌地區(qū)處于新構(gòu)造時(shí)期以來以強(qiáng)烈坳陷為特征的地區(qū),而三危山斷裂在整體下沉區(qū)卻表現(xiàn)出明顯不同于整體的新構(gòu)造運(yùn)動(dòng)特征,其所控制的北東向三危山隆起區(qū)分割了敦煌盆地和踏實(shí)盆地(圖1)。三危山斷塊隆起區(qū)地層主要由前震旦紀(jì)敦煌群變質(zhì)巖和華力西期巖漿巖組成(圖1),第四系地層主要構(gòu)成了山前略有起伏的傾斜地帶及戈壁平原,以上更系統(tǒng)和全新統(tǒng)出露最為廣泛,中、下更新統(tǒng)地層局部保留(圖1)。三危山斷裂對(duì)不同微地貌單元具有控制作用,多數(shù)地段構(gòu)成前第四紀(jì)地層與第四紀(jì)地層的分界線,亦構(gòu)成了三危山斷塊隆起區(qū)與山前緩傾斜平原的分界線。

在三危山斷裂以南更大區(qū)域范圍內(nèi)是以左旋走滑著名的阿爾金斷裂帶,阿爾金斷裂帶的東段控制了祁連山斷裂系向北西方向的延伸,該新構(gòu)造運(yùn)動(dòng)活躍區(qū)以北西向的活動(dòng)斷裂最為發(fā)育,為區(qū)域地殼運(yùn)動(dòng)差異性強(qiáng)烈地區(qū),活動(dòng)斷裂控制了現(xiàn)今地貌的形成,中強(qiáng)震及小震較多。區(qū)域上該區(qū)曾發(fā)生過1932年昌馬7.6級(jí)地震,1785年玉門惠回堡6級(jí)地震,2002年玉門5.9級(jí)地震[11-13],這些地震的等震線長軸均為北西向,與發(fā)育在祁連山隆起區(qū)的北西向活動(dòng)斷裂密切相關(guān)。三危山斷裂以北為新構(gòu)造時(shí)期以來以整體性抬升弱的差異性運(yùn)動(dòng)為特征,該地區(qū)長期處于剝蝕夷平狀態(tài),形成現(xiàn)今荒涼的戈壁殘山,這一地區(qū)少見中強(qiáng)震,顯示出以整體緩慢抬升的新構(gòu)造運(yùn)動(dòng)為主,差異性運(yùn)動(dòng)非常微弱。

圖1 三危山斷裂地質(zhì)構(gòu)造簡圖Fig.1 Simple geological map of Sanweishan fault

2 三危山斷裂活動(dòng)特征

三危山斷裂位于三危山西北麓,斷裂西起黨河水庫附近,在鳴沙山一帶由于風(fēng)積砂的堆積呈隱伏狀,向北東穿過西水溝、東水溝,過樹溝子、漢峽、蘆草溝、瓜州口等地延伸至雙塔水庫一帶,斷裂全長約175 km,傾角50°～70°,傾向SE[4]。根據(jù)斷裂幾何特征的差異,三危山斷裂可分為三段:西水溝—樹溝子段(西段),長約60 km;樹溝子—十工口子段(中段),長約65 km;十工口子—雙塔段(東段)長約50 km(圖1)。三危山隆起區(qū)以三危山斷裂為北西邊界,其形態(tài)明顯受三危山斷裂控制,山體寬度介于4～10 km,海拔高度介于約1 450～1 950 m,總體形態(tài)西高東低、南緩北陡。

晚更新世以來,三危山地區(qū)運(yùn)動(dòng)特征表現(xiàn)出間歇性抬升的特點(diǎn)[6],三危山斷裂表現(xiàn)出分段活動(dòng)的特征[4]。根據(jù)洪積階地拔河高度及部分階地的熱釋光測年,張?jiān)Ｃ鞯萚6]得到斷裂西、中、東段抬升速率分別約為2.06 mm/a、0.47 mm/a和0.84 mm/a,具有明顯的分段活動(dòng)特征。需要指出的是,該抬升速率實(shí)際是河流下切速率,階地下切是個(gè)復(fù)雜的過程,受氣候和構(gòu)造雙重作用控制[14-16],在完全不考慮氣候作用的影響下,可以認(rèn)為該下切速率與抬升速率相當(dāng),也等同于斷裂垂直滑動(dòng)速率。但全新世以來,祁連山西段氣候在控制階地下切中起了決定性的作用[17],本區(qū)下切速率并不能等同于抬升速率,下切速率要遠(yuǎn)遠(yuǎn)大于抬升速率。

在斷裂西段臭水溝南,斷層陡坎發(fā)育于晚更新世洪積扇上[圖2(a)],形成的斷層陡坎高度在(3.9±0.5) m左右,其下更新世洪積扇上的陡坎高度約(2±0.2) m[圖2(b)],但較低陡坎存在一定的不確定性,可能是不同期洪積扇形成的侵蝕坎,而(3.9±0.5) m的數(shù)據(jù)較為可靠。在千佛洞附近,洪積扇上發(fā)育斷層陡坎[圖2(c)]兩條剖面線測量得到的陡坎高度分別為(4.6±0.2) m和(5.4±0.5) m[圖2(d)]。從區(qū)域上看,兩處陡坎可能發(fā)育于同一級(jí)地貌面上,平均高度在(4.6±0.7) m。該地貌面可能形成于晚更新世中段,以60 ka計(jì)算,則得到斷裂垂直滑動(dòng)速率為(0.08±0.01) mm/a;根據(jù)張?jiān)Ｃ鞯萚6]在西段榆樹溝得到的年代40 ka計(jì)算,則斷裂垂直滑動(dòng)速率為(0.12±0.01) mm/a,綜合考慮三危山斷裂西段晚更新世以來的垂直滑動(dòng)速率在0.1 mm/a 左右。

在三危山斷裂的中、東段,斷層陡坎主要表現(xiàn)為基巖陡坎,晚更新世洪積扇上基本無斷層陡坎發(fā)育,很難確定斷裂的滑動(dòng)速率。在斷裂東段,僅在十工口子附近發(fā)現(xiàn)在相當(dāng)于晚更新世上段(約20 ka)的階地上發(fā)育反向陡坎,其高度僅有0.3 m[4],其滑動(dòng)速率小于西段。從滑動(dòng)速率上判斷,三危山斷裂活動(dòng)性差異明顯,其活動(dòng)性西段強(qiáng)于中東段。雖然前人對(duì)三危山斷裂垂直滑動(dòng)速率的研究存在高估的可能性[6],但其速率值差異也反映了相同的特征。

圖2 三危山斷裂西段斷層陡坎測量Fig.2 Measurements on the fault scarps in the western segment of the Sanweishan fault

3 三危山地區(qū)地形因子特征

數(shù)字高程模型(Digital Elevation Model,DEM)是利用一定范圍內(nèi)的地形高程數(shù)據(jù)將地表形態(tài)數(shù)字化并以模型化方式呈現(xiàn)的實(shí)體地面模型,其蘊(yùn)含的地形地貌數(shù)據(jù)豐富,被廣泛應(yīng)用于地學(xué)分析領(lǐng)域[18]。本文采用1∶5萬數(shù)字地形圖等高線生成的DEM數(shù)據(jù)(格網(wǎng)分辨率25 m),利用GIS技術(shù)及Matlab程序提取了三危山研究區(qū)各流域盆地內(nèi)的河流或沖溝67條(圖3),通過對(duì)高程、坡度、河道陡峭指數(shù)及縱剖面等地形因子的定量分析,綜合研究三危山地區(qū)地貌特征因子對(duì)區(qū)域隆升差異的響應(yīng)。

3.1 高程與坡度

高程是對(duì)地形點(diǎn)高低大小值的度量,而坡度則反映地表起伏的陡緩,二者在對(duì)地形發(fā)育特征的研究中均得到了廣泛應(yīng)用[2,19]。研究利用ArcGIS空間分析功能,設(shè)置廊帶剖面測線兩側(cè)的數(shù)據(jù)統(tǒng)計(jì)緩沖寬度為100 m,長度為140～160 km,分別提取和計(jì)算三危山北麓沿山前斷裂走向的平均高程剖面[圖4(a)]與坡度剖面[圖4(b)],進(jìn)而揭示出三危山地區(qū)山體隆升及地表起伏特征。

圖4 三危山高程與坡度因子提取結(jié)果Fig.4 Extraction of elevation and slope factor of Sanwei Shan

從研究所得結(jié)果來看,無論高程還是坡度在空間分布上均呈現(xiàn)西高東低的趨勢(shì),三危山平均高程最大值域在斷裂西段,大部分海拔高度大于1 400 m,斷裂中段高程最低,多在約1 200 m,西段略高于東段。坡度的分布與高程分布存在相一致的趨勢(shì),在山體的整個(gè)西段,平均坡度集中在20°以上,表明該段坡度變化陡峭,而中東段坡度值變化總體趨于平緩,平均坡度在10°左右,除樹溝子以西山體北麓部分表現(xiàn)出較高值外,段落其余大部分山體區(qū)均以低值分布,但東段山體北麓坡度高值分布比中段更為集中連續(xù),在距斷裂較近的山前沿走向延伸。

3.2 河道陡峭指數(shù)

河道陡峭指數(shù)(ks)是反映基巖河道巖石抗侵蝕能力與河流動(dòng)力侵蝕能力相對(duì)大小的系數(shù),與河道比降(S)和流域面積(A)間存在冪律關(guān)系[20-22]:

S=ksA-θ

(1)

式中:ks即為河道陡峭指數(shù);θ為河道凹度。通常將不同河道的陡峭指數(shù)進(jìn)行標(biāo)準(zhǔn)化計(jì)算(ksn)后再進(jìn)行比較。大量實(shí)踐研究指出,河道陡峭指數(shù)與基巖隆升速率呈正相關(guān)關(guān)系[21,23],高ks值指示基巖區(qū)基巖隆升速率較快,而低ks值對(duì)應(yīng)的基巖隆升速率則較慢。因此,可根據(jù)ks值的變化判斷其對(duì)構(gòu)造隆升差異響應(yīng)的敏感程度。

研究通過計(jì)算各沖溝的標(biāo)準(zhǔn)化河道陡峭指數(shù)ksn,對(duì)三危山67條支流的河道陡峭指數(shù)空間分布規(guī)律進(jìn)行繪制。如圖5所示,三危山北麓各支流ksn值在東西走向上變化趨勢(shì)明顯,總體來看,絕大部分西段支流ksn值大于中、東段支流,在斷裂分段點(diǎn)處有ksn值激增的特點(diǎn)。其中,ksn最大值位于斷裂西段,達(dá)110,最小值位于斷裂東段,ksn值僅14;斷裂西段流域內(nèi)平均ksn值為56,東段次之,為31,最小值為中段流域,為29。

圖5 67條河流的河道陡峭指數(shù)Fig.5 The ksn value of 67 rivers in the study area

3.3 河道坡降指數(shù)

河道坡降指數(shù)(Stream gradient index,SL)是單位距離的河段坡度與河段中點(diǎn)距源頭距離的乘積[24],可表示為關(guān)系:

(2)

式中:▽H為單位河段高程變化值;▽L為單位河段長度;L為河段中心距河源的距離。此外,河道剖面上河源與出水口兩點(diǎn)間相接的直線代表全河段處于動(dòng)力平衡時(shí)均衡狀態(tài),其斜率為均衡坡降指標(biāo),用K表示。利用均衡坡降指標(biāo)值(K)對(duì)各河流坡降指標(biāo)進(jìn)行標(biāo)準(zhǔn)化(SL/K)計(jì)算,可進(jìn)一步比較不同河流間的河段坡降指標(biāo)差異。河道坡降指數(shù)對(duì)河流各河段坡度變化響應(yīng)敏感,對(duì)巖性差別、構(gòu)造活動(dòng)的強(qiáng)弱具有很好的指示作用[19]。

研究依次對(duì)西水溝—雙塔段的67條沖溝或河流的坡降指標(biāo)進(jìn)行了提取。結(jié)果表明[圖6(a)],研究區(qū)SL值大部分處于0～100區(qū)間,均值81,最大值位于西水溝—樹溝子段,達(dá)587;SL指標(biāo)值整體西段與中、東段差異較大,西段SL平均值達(dá)185,而中、東段SL平均值分別僅為53、54。為進(jìn)一步比較不同河流間的河段坡降指標(biāo)差異,研究利用均衡坡降指標(biāo)值(K)對(duì)各河流坡降指標(biāo)進(jìn)行了標(biāo)準(zhǔn)化(SL/K),結(jié)果顯示[圖6(b)],各河流標(biāo)準(zhǔn)化河道坡降指標(biāo)(SL/K)與坡降指標(biāo)分布趨勢(shì)一致,斷裂西段河道最陡,SL/K平均值3.3,東段次之,SL/K均值2.3,中段河道坡降最緩,SL/K均值僅1.6。

圖6 67條河流的坡降因子提取結(jié)果Fig.6 Extraction of SL factor of 67 rivers in the study area

4 討論

河道陡峭指數(shù)與河道坡降指數(shù)值的變化受巖性抗侵蝕能力影響,在抗侵蝕能力強(qiáng)的基巖區(qū),河道陡峭指數(shù)與河道坡降指數(shù)均在高值區(qū)間占有較大比例[25-26]。研究區(qū)前第四系地層主要為前震旦系敦煌群和華力西期巖漿巖,敦煌群變質(zhì)巖與侵入巖體構(gòu)成了三危山隆升山地(圖1)。巖體類型主要為花崗質(zhì)片麻巖、片巖、石英閃長巖及白云質(zhì)大理巖等。第四系地層以更新世和全新世的風(fēng)積、沖積、沖洪積及殘坡的砂巖、礫石等松散堆積為主。綜合研究區(qū)巖性分布特征認(rèn)為,三危山斷裂西、中、東段一帶所在基巖區(qū)巖石抗侵蝕能力均較強(qiáng),差異性并不顯著,ksn、SL/K值在樹溝子—十工口子、十工口子—雙塔堡一帶應(yīng)偏大,而研究結(jié)果并不符合這一特征,ksn及SL/K高值均分布于西水溝—樹溝子段,樹溝子—十工口子段主要為低值,分布特征與巖性抗侵蝕能力所在區(qū)不一致,研究認(rèn)為流域地貌參數(shù)受巖性影響有限。

降水可以增大河道徑流量進(jìn)而提升對(duì)河道侵蝕能力,對(duì)流域高程、河道縱剖面坡度變化有一定影響[23]。據(jù)1981—2005年甘肅西北高原地區(qū)的降雨量數(shù)據(jù)統(tǒng)計(jì)結(jié)果顯示[27],三危山地區(qū)降水量正處于水平年穩(wěn)態(tài)階段,降雨量少(38.28～51.1 mm/a)且穩(wěn)定。因此,在降雨量少且?guī)r性接近的條件下,降雨對(duì)地表的侵蝕作用程度相似,地貌參數(shù)應(yīng)相同或接近。實(shí)際上,河道陡峭指數(shù)(圖5)及河流的坡降指數(shù)與標(biāo)準(zhǔn)化坡降指數(shù)變化起伏較大(圖6),其值未因侵蝕條件相同而表現(xiàn)出一致性,表明降雨對(duì)地表剝蝕改造能力不強(qiáng)。綜合認(rèn)為降雨不是影響三危山地區(qū)地貌參數(shù)的主因,ksn、SL/K指標(biāo)差異分布可能更多的與構(gòu)造活動(dòng)相關(guān)。

構(gòu)造抬升可以改變區(qū)內(nèi)山體、河道高程起伏變化塑造地貌形態(tài),從而能夠影響河道坡降指數(shù)、河道陡峭指數(shù)分布趨勢(shì)。為進(jìn)一步明確區(qū)域地貌參數(shù)值與構(gòu)造活動(dòng)性的相關(guān)性,研究對(duì)所得地貌特征因子與斷裂滑動(dòng)速率進(jìn)行了綜合對(duì)比和分析。由前所述,三危山斷裂各個(gè)段落滑動(dòng)速率具有明顯的差異,西段大于中東段。地貌高程和坡度參數(shù)均顯示與滑動(dòng)速率分布相同的特征(圖4)。對(duì)比不同段落ksn、SL、SL/K地貌參數(shù)值,西段較大,中、東段較為接近(圖5,6)?？紤]到前述巖性與氣候因素作用對(duì)地貌參數(shù)的影響較小,地貌參數(shù)的變化主要受斷裂活動(dòng)性的控制,地貌因子的分布特征與斷層滑動(dòng)速率差異的一致性反映三危山的隆升主要受斷裂活動(dòng)性控制。

三危山斷裂地處青藏高原北緣向北擴(kuò)展的前緣地帶,屬于阿爾金斷裂帶向北擴(kuò)展的派生斷裂,三危山斷裂在一定程度上遠(yuǎn)程吸收了阿爾金斷裂走滑速率向東的衰減,調(diào)節(jié)著該地區(qū)構(gòu)造應(yīng)力分布[28]。但三危山斷裂整體滑動(dòng)速率不高,左旋滑動(dòng)速率僅約0.3 mm/a[29],構(gòu)造活動(dòng)整體不強(qiáng),這也與本區(qū)無強(qiáng)震記錄相符合。青藏高原北緣與阿爾金斷裂走向相近或交角較小的斷裂往往具有左旋走滑性質(zhì)[30],本區(qū)三危山斷裂和陽關(guān)斷裂與阿爾金斷裂走向近平行,均具有左旋走滑特征。三危山斷裂和陽關(guān)斷裂以左旋右階的形式展布(圖1),階區(qū)以擠壓隆升為主要特征,形成了鳴沙山及控制其隆升的鳴沙山斷裂。鳴沙山斷裂的新活動(dòng)及鳴沙山地區(qū)的擠壓隆起加速了三危山西段的隆升,此外,巴兔山斷裂與三危山斷裂作為一個(gè)統(tǒng)一的構(gòu)造變形體系吸收青藏高原向外圍的推擠變形,其分擔(dān)了三危山中東段的變形量,使三危山西段與中東段地貌發(fā)育差異性進(jìn)一步擴(kuò)大,顯示出西段強(qiáng)而中東段弱的特點(diǎn)。

5 結(jié)語

研究基于1∶5萬DEM數(shù)據(jù),利用GIS技術(shù)及Matlab程序提取了三危山研究區(qū)各流域盆地內(nèi)的河流或沖溝67條,通過對(duì)高程、坡度、河道陡峭指數(shù)及縱剖面等地形因子參數(shù)值的系統(tǒng)提取,綜合區(qū)內(nèi)地層發(fā)育概況、降水條件及野外地質(zhì)地貌特征調(diào)查,定量分析三危山地區(qū)地貌發(fā)育特征對(duì)斷裂構(gòu)造活動(dòng)的響應(yīng)。研究表明:

(1) 三危山各段地貌發(fā)育程度存在明顯差異,綜合地貌參數(shù)值(高程、坡度、河道陡峭指數(shù)值和河道坡降指數(shù)值)空間分布不均,有斷裂西段遠(yuǎn)大于中東段特征,揭示三危山抬升速率西段大于中東段,差異性抬升主要受控于三危山斷裂活動(dòng)性差異。

(2) 鳴沙山的隆升及鳴沙山斷裂的新活動(dòng)加速了三危山西段的隆升,此外,巴兔山斷裂與三危山斷裂作為一個(gè)統(tǒng)一的構(gòu)造變形體系吸收青藏高原向外圍的推擠變形,其分擔(dān)了三危山中東段的變形量,對(duì)三危山西段與中東段地貌發(fā)育差異性擴(kuò)大具有促進(jìn)作用。