基于DEM的丹霞地貌演化階段劃分*

章桂芳 ，陳凱倫，張浩然，張慧

(1. 中山大學地球科學與工程學院，廣東 廣州 510275;2. 廣東省地球動力作用與地質(zhì)災(zāi)害重點實驗室，廣東 廣州 510275;3. 東北農(nóng)業(yè)大學資源與環(huán)境學院，黑龍江 哈爾濱 150030)

丹霞地貌是指以陡崖坡為特征的紅層地貌[1]，目前中國已發(fā)現(xiàn)丹霞地貌共有1 000余處[2]。自20世紀30年代末，陳國達提出“丹霞地形”以來[3]，丹霞地貌的研究已經(jīng)經(jīng)歷了初創(chuàng)、成型和發(fā)展3個階段之后，2009年開始了國際化研究階段[1]。雖然丹霞地貌研究已經(jīng)取得不少成果[4-5]，并且隨著現(xiàn)代技術(shù)手段的發(fā)展，出現(xiàn)了一些涉及定量研究、微觀實驗和現(xiàn)代測試手段的精細研究實例[6-16]。但是，從總體上來說，丹霞地貌基礎(chǔ)研究仍然十分薄弱，多數(shù)研究偏重于宏觀、定性的描述和推論，缺乏系統(tǒng)、定量和對比研究[1]。

數(shù)字高程模型(Digital Elevation Model, DEM)是地表形態(tài)的數(shù)字表達和模擬，蘊含著豐富的地形信息[17]。高分辨率DEM數(shù)據(jù)的獲取以及分析技術(shù)的不斷突破使得DEM成為地形地貌研究領(lǐng)域不可或缺的重要技術(shù)手段[18-26]，通過對DEM地貌參數(shù)的提取與分析來解決構(gòu)造運動、階地提取以及地貌發(fā)育階段的確定等問題已被多數(shù)學者所接受[20,27-39]。在DEM數(shù)據(jù)集中，航天飛機雷達地形測繪項目(Shuttle Radar Topography Mission — SRTM)提供的DEM數(shù)據(jù)覆蓋了北緯60°至南緯56°范圍內(nèi)80%面積的區(qū)域[40]，原始數(shù)據(jù)經(jīng)處理后平面精度能夠達到±20 m，高程精度達到±16 m[41]，是目前使用廣泛的全球地形數(shù)據(jù)[42-43]。相對其在黃土地貌和喀斯特地貌等方面的大量應(yīng)用[44-49]，使用DEM進行丹霞地貌的相關(guān)研究則非常缺乏。

從本質(zhì)上說，地表形態(tài)來源于地球內(nèi)外力對地殼的綜合作用[50-51]，作為紅層侵蝕地貌，丹霞地貌的發(fā)育和演化也受控于內(nèi)動力和外動力兩大動力系統(tǒng)[1]，它始于紅層盆地的抬升，在流水、風化、重力等外動力作用下，溝谷不斷展寬，崖壁崩塌后退，山頂面積不斷縮小。黃進將其演化過程定性劃分為幼年期、壯年期、老年期3個階段及“回春期”[52]；彭華在“中國丹霞”申報世界自然遺產(chǎn)過程中，基于侵蝕量和保留量的對比關(guān)系，將“中國丹霞”6個世界自然遺產(chǎn)提名地組成青年早期、青年晚期、壯年早期、壯年晚期、老年早期和老年晚期六個階段的演化系列[1]。

本次研究擬采用DEM數(shù)據(jù)，對貴州赤水、福建泰寧、湖南崀山、廣東丹霞山、江西龍虎山和浙江江郎山等6個丹霞地貌遺產(chǎn)地計算面積高程積分曲線及積分值，以期對丹霞地貌的不同發(fā)展階段提供定量劃分依據(jù)，促進丹霞地貌演化模式的定量研究進程。

1 數(shù)據(jù)與方法

1.1 基礎(chǔ)數(shù)據(jù)

本次研究選取的是中國丹霞6個遺產(chǎn)地提名序列，分別為貴州赤水、福建泰寧、湖南崀山、廣東丹霞山、江西龍虎山和浙江江郎山，他們代表了丹霞地貌從青年到老年的演化序列，均分布在中國南方(圖 1a)。需要說明的是：赤水的東西部和泰寧的南北部的地貌發(fā)育存在較大的區(qū)域差異，因此對這兩個地區(qū)分別進行分區(qū)研究。

研究采用ASTER DEM數(shù)據(jù)(V2版)，來源于中國科學院計算機網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云平臺(http:∥www.gscloud.cn)，空間分辨率為30 m。根據(jù)世界自然遺產(chǎn)保護區(qū)劃定范圍裁剪研究區(qū)域DEM(圖 1b-i)，并對其面積和高程信息統(tǒng)計如下(表 1)：赤水東部、丹霞山和龍虎山的面積最大，而江郎山的面積最小，為13.05 km2；赤水地區(qū)因地處云貴高原邊緣，高程最低值、最高值和平均值均為所有地區(qū)的最大值，而泰寧地區(qū)的高程差最小(北部為384 m，南部為403 m)，說明其地勢起伏是所有地區(qū)中最小的。

表1 各個丹霞地貌區(qū)域面積和高程信息統(tǒng)計Table 1 Area and elevation statistics of the Danxia sites

圖1 各丹霞地貌研究區(qū)域的位置及DEMFig.1 Location of Danxia sites and DEM in this study

1.2 研究方法

本次研究主要采用Strahler提出的面積高程積分(Hypsometric Integral，簡稱HI)[53]，該方法通過面積高程積分將Davies的地貌發(fā)育模式定量化，分為3個階段：幼年階段(HI>0.6)，壯年階段(0.4

面積高程積分值的計算方法有積分曲線法、體積比例法和起伏比法等[60]，本次研究采用積分曲線法，計算過程如下：

① 將流域的高差無限細分，分別求取大于每一細分單元高程值的流域面積，并將其均一化后，作為x值，即xi=ai/A(xi為大于第i分級高程的面積百分比，ai為大于第i分級高程的面積，A為流域總面積)；

② 高差細分的高程均一化值作為y值，即yi=hi/H(yi為第i分級相對高差百分比，hi為第i分級高程與流域最低點的相對高差，H為流域內(nèi)高差)；

③ 將各點展布到直角坐標系中，形成一條直線，即y=f(x)；

2 結(jié)果與討論

2.1 面積高程積分結(jié)果

采用積分曲線法對各個丹霞區(qū)域的DEM數(shù)據(jù)計算獲得的面積高程積分曲線如圖2所示，面積高程積分則如表2所示，考慮到圖件的清晰表達，圖 2中僅顯示赤水西部和泰寧北部的面積高程積分曲線，赤水東部和泰寧南部的曲線將在后面的討論中呈現(xiàn)。由圖 2可知，赤水西部地區(qū)的曲線呈現(xiàn)凸形特征，結(jié)合其面積高程積分值為0.65，表明該地區(qū)處于較年輕的演化階段；泰寧北部、崀山和丹霞山的曲線呈現(xiàn)S狀，表明該區(qū)域地貌已經(jīng)演化至較為成熟的階段，且從泰寧北部到崀山到丹霞山的曲線的凹形特征逐步加強，則表明這3個區(qū)域的地貌演化越來越成熟；龍虎山的曲線在相對高程較高的部分(相對高程h/H在0.3～1.0左右)斜率很大，在相對高程較低的部分(相對高程h/H在0～0.3左右)則趨于平緩，這說明龍虎山地區(qū)地勢較高的區(qū)域面積占比較小，而地勢較低的區(qū)域占了很大比例。

圖2 各丹霞地貌區(qū)面積高程積分曲線匯總圖Fig.2 Hypsometric curves of the Danxia sites

從總體特征來看，從赤水西部到龍虎山，隨著丹霞地貌的發(fā)育程度逐漸成熟，面積高程積分曲線形態(tài)從凸形轉(zhuǎn)為凹形，且凹形特征逐步加強，說明隨著丹霞地貌的發(fā)育，地層侵蝕量在不斷增加，符合發(fā)育過程中侵蝕量的變化特征。分析表 2中各個丹霞區(qū)域的面積高程積分值可以發(fā)現(xiàn)赤水東西部和泰寧南北部的面積高程積分均具有較大的差異；與Strahler劃分的地貌發(fā)展階段標準進行對比，丹霞地貌的HI值整體偏小，即向老年期偏移。需要指出的是，龍虎山(老年早期)和丹霞山(壯年晚期)在丹霞地貌世界遺產(chǎn)提名地序列中的發(fā)育階段比江郎山(老年晚期)年輕，但是它們的曲線形態(tài)與此并不相符(圖 2)，而且龍虎山和丹霞山的面積高程積分值均低于江郎山(表 2)，似乎違背規(guī)律，其原因與龍虎山和丹霞山地貌發(fā)育的區(qū)域差異關(guān)系密切，即與研究區(qū)的范圍密切相關(guān)，接下來作詳細分析。

2.2 區(qū)域差異性對比

丹霞地貌的發(fā)育受控于內(nèi)動力和外動力2大因素[1]，即使在同一地區(qū)，影響丹霞地貌發(fā)育的因素及其強度也不盡相同，因此在同一地區(qū)保留了不同演化階段的地貌形態(tài)，呈現(xiàn)出區(qū)域的差異性。下面選取赤水、泰寧、丹霞山和龍虎山進行同一區(qū)域的差異性對比研究。

表2 各個丹霞地貌區(qū)域面積高程積分值Table 2 Hypsometric indexes of the Danxia sites

2.2.1赤水西部和東部 “赤水丹霞”是發(fā)育于亞熱帶巨型紅層盆地南緣的高原——峽谷型和山嶺型丹霞景觀，是侏羅系、白堊系的紅層經(jīng)過第三紀長期剝蝕夷平后在第四紀強烈抬升過程中因流水的快速切割而形成的[61-62]。貴州赤水的東部和西部地區(qū)由赤水河自然分割(圖 3a)，東部以流水深度侵蝕切割形成山原峽谷地貌，西部河流切割程度小，呈高原峽谷地貌，顯示其發(fā)育中的區(qū)域差異，因此將赤水分為東西兩片區(qū)域進行差異性對比研究。計算獲得東西兩區(qū)面積高程積分曲線及HI值如圖 3b所示：赤水東區(qū)和西區(qū)的曲線均呈現(xiàn)凸形特征，HI值均大于0.5，說明兩區(qū)發(fā)育程度都較為年輕，但是西區(qū)的凸形特征比東區(qū)更為顯著，HI值也較東區(qū)大得較多，表明西區(qū)地貌侵蝕程度比東區(qū)低， 符合東西兩區(qū)的現(xiàn)有的地貌特征。

2.2.2 泰寧北部和南部 泰寧盆地是在華夏古板塊武夷山隆起帶背景上發(fā)育的白堊紀紅色斷陷盆地，由朱口和梅口兩個北東向紅色盆地構(gòu)成，構(gòu)成了泰寧丹霞的北部和南部(圖 4a)。泰寧盆地的丹霞地貌是由白堊紀紅層經(jīng)過第三紀以來多次抬升和長期剝蝕夷平，以流水切割、侵蝕為主因而形成的。與赤水類似，泰寧的兩個盆地中的丹霞地貌也具有區(qū)域的差異性。

圖 4b為計算獲得的泰寧北部和南部的面積高程積分曲線及HI值，兩者的曲線均呈現(xiàn)凹形特征且HI值小于0.5(分別為北部0.43，南部0.36)，顯示兩者均處于侵蝕量大于保存量的階段，北部地區(qū)的曲線在南部之上且相距較遠，顯示其發(fā)育程度比南部地區(qū)年輕且差距較大。據(jù)已有研究，北部的朱口盆地紅層地殼抬升相對較晚，地貌發(fā)育時間較短，所以呈現(xiàn)以深切峽谷曲流等為特征的典型的青年期丹霞特征，而其南部的梅口盆地則以發(fā)展到壯年期的丹霞地貌為主，呈現(xiàn)丹霞峰叢景觀。所以本研究獲得的面積高程積分曲線及HI值印證了泰寧盆地的丹霞發(fā)育由于地殼抬升時間、抬升幅度和抬升速度的差異而存在較大的區(qū)域差異。

2.2.3 丹霞山東部和西部 丹霞山發(fā)育在南嶺褶皺帶的白堊紀構(gòu)造盆地丹霞盆地中，錦江由北到南穿越而過(圖 5a)，將丹霞山分為東西兩個部分。丹霞山主要由近水平巖層構(gòu)成，總體上處于地貌發(fā)育的壯年期階段，表現(xiàn)出疏密相間，組合有序的簇群式丹霞峰叢——峰林景觀。

圖 3 赤水東西兩區(qū)面積高程積分對比Fig.3 Hypsometric curves of the eastern and western parts of Chishui

圖 4 泰寧南北兩區(qū)面積高程積分對比Fig.4 Hypsometric curves of the northern and southern parts of Taining

圖5 丹霞山東西兩區(qū)面積高程積分對比Fig.5 Hypsometric curves of the eastern and western parts of Danxiashan

以丹霞山為整體計算獲得的面積高程積分曲線呈現(xiàn)強烈凹形特征(圖 2)，HI值為0.27(表 2)，比發(fā)育階段為老年晚期的江郎山的HI值更小，似乎并不合理。因此，本研究嘗試以錦江為界，分別計算東西兩區(qū)的面積積分曲線及HI值(圖 5b)，結(jié)果表明丹霞山東西兩區(qū)的HI值非常相近，分別為0.35和0.36，說明廣東丹霞山東西兩側(cè)丹霞地貌的總體發(fā)育程度較為相近。而兩者的積分值與包含錦江的整個區(qū)域計算獲得的0.27有明顯差距，說明錦江兩岸區(qū)域的地貌形態(tài)極大地影響了整個丹霞山區(qū)域的地貌發(fā)育程度。分析原因可能與區(qū)域內(nèi)的的錦江密切相關(guān)，因為錦江流經(jīng)丹霞山地區(qū)的主河道多為寬谷，造成兩岸區(qū)域地勢平緩，在面積高程積分中則表現(xiàn)為強烈的侵蝕。丹霞山的分析結(jié)果進一步表明不同的區(qū)域尺度內(nèi)獲得的HI值具有很大的差距，同時表明面積高程積分可以敏銳地指示丹霞地貌發(fā)育的多期性。

2.2.4 龍虎山全區(qū)和核心區(qū) 龍虎山發(fā)育在中生代陸相盆地——信江盆地邊緣，瀘溪河穿區(qū)而過，是由晚白堊紀紅色碎屑巖組成的疏散峰林寬谷型丹霞地貌區(qū)，而龍虎山地質(zhì)公園則為龍虎山丹霞地貌區(qū)域的核心(圖 6a)。

圖6 龍虎山全區(qū)和核心區(qū)域面積高程積分對比Fig.6 Hypsometric curves of the whole area and center part of Longhushan

與丹霞山類似，以龍虎山全區(qū)為整體計算獲得的面積高程積分值為0.24(表 2)，比發(fā)育階段為老年晚期的江郎山的HI值(0.35)更小，因此本研究嘗試只計算龍虎山的核心區(qū)域(龍虎山地質(zhì)公園)的面積高程積分，結(jié)果顯示其HI值達到了0.47，幾乎達到是龍虎山全區(qū)的2倍。究其原因，應(yīng)該與瀘溪河谷密切相關(guān)，瀘溪河谷兩側(cè)主要為準平原化的低丘崗地，龍虎山和龜峰等地保存的峰叢、峰林、孤峰、殘丘等為該區(qū)僅存的高地勢區(qū)域。龍虎山準平原化代表了地貌發(fā)育的老年階段，因此對于龍虎山地區(qū)的面積高程積分計算受到區(qū)域尺度的巨大影響，這也印證了面積高程積分在丹霞地貌演化定量研究中具有很好的指示作用。

3 結(jié) 論

本次研究嘗試將面積高程積分方法用于丹霞地貌演化階段的定量劃分，主要取得以下幾點認識：

1)不同丹霞區(qū)域的面積高程積分研究結(jié)果表明面積高程積分可以較好地指示丹霞地貌的不同演化階段：隨著丹霞地貌從年輕到成熟的發(fā)育，面積高程積分曲線形態(tài)從凸形轉(zhuǎn)為凹形，且凹形特征逐步加強，HI值也越來越小，符合不同發(fā)育階段所代表的侵蝕量的變化。

2)將各個丹霞區(qū)域的面積高程積分值與Strahler劃分的地貌發(fā)展階段標準進行對比，發(fā)現(xiàn)丹霞地貌的HI值整體往老年期偏移，這與中國丹霞演化系列是建立在幾個代表小區(qū)的基礎(chǔ)上有關(guān)，這提示我們丹霞地貌發(fā)育階段的定量化研究與研究區(qū)的劃定密切相關(guān)。

3)同一地區(qū)不同區(qū)域的面積高程積分對比研究表明其可以敏銳地反映丹霞地貌發(fā)育的多期性，對于定量研究丹霞地貌具有重要作用。丹霞地貌發(fā)育的區(qū)域差異是進行丹霞演化階段確定中極其重要的影響因素，在進行面積高程積分的計算中需要充分考慮。

致謝：感謝中山大學地理科學與規(guī)劃學院彭華教授在本次研究和論文撰寫中給予的指導和幫助。

參考文獻：

[1] 彭華, 潘志新, 閆羅彬, 等. 國內(nèi)外紅層與丹霞地貌研究述評 [J]. 地理學報, 2013, 64(9): 1170-1181.

PENG H, PAN Z X, YAN L B, et al. A review of the research on red beds and Danxia landform [J]. Acta Geographica Sinica, 2013, 64(9): 1170-1181.

[2] 黃進, 陳致均, 齊德利.中國丹霞地貌分布(下)[J]. 山地學報, 2015, 33(4): 649-673.

HUANG J, CHEN Z J, QI D L. The distribution of Danxia landform in China [J]. Mountain Research, 2015, 33(4): 649-673.

[3] 陳國達, 劉輝泗. 江西貢水流域地質(zhì)[J]. 江西地質(zhì)匯刊, 1939 (2): 164.

CHEN G D, LIU H S. Geology of Gongshui valley， Jiangxi[J]. Geological Transactions of Jiangxi, 1939 (2): 1-64.

[4] 彭華. 中國丹霞地貌研究進展[J]. 地理科學, 2000, 20(3): 203-211.

PENG H. A survey of the Danxia Landform research in China [J]. Scientia Geographica Sinica, 2000, 20(3): 203-211.

[5] 朱誠, 馬春梅, 張廣勝, 等. 中國典型丹霞地貌成因研究[M]. 北京: 科學出版社, 2016.

ZHU C, MA C M, ZHANG G S, et al. Development mechanisms of the typical Danxia landforms in China [M]. Beijing: Science Press, 2016.

[6] 郭國林, 郭福生, 劉曉東, 等. 丹霞地貌砂巖的微觀化學風化作用電子探針研究[J]. 中國巖溶, 2006, 25(2): 172-176.

GUO G L, GUO F S, LIU X D, et al. Study by Epma on microcosmic chemical weathering to sandstone in Danxia landform [J]. Carsologica Sinica, 2006, 25(2): 172-176.

[7] 姜勇彪, 郭福生, 劉林清, 等. 龍虎山丹霞地貌區(qū)河流階地地貌面的熱釋光測年研究[J]. 東華理工大學學報, 2006, 29(3): 225-228.

JIANG Y B, GUO F S, LIU L Q, et al. TL dating research of geomorphic surface on lower terraces of rives in Longhushang Danxia landform areas [J]. Journal of East China Institute of Technology, 2006, 29(3): 225-228.

[8] 黃進. 丹霞地貌區(qū)地殼上升速率公式改進及其意義[J]. 經(jīng)濟地理, 2006, 26(增刊): 8-13.

HUANG J. Improvement on the formula used to calculate tectonic uplift rates in Danxia landforms areas and its significance [J]. Economic Geography, 2006, 26(Suppl):8-13.

《林黛玉進賈府》的人物語言特點是“含而不露，意在言外”，簡而言之就是話中有話。而《老人與海》人物語言簡潔凝練，在平實地敘述中彰顯主題。現(xiàn)以《林黛玉進賈府》中王熙鳳和《老人與?！飞５貋喐绲娜宋镎Z言描寫為例，分析比較兩篇課文語言風格。

[9] 陳姝, 朱誠, 彭華, 等. 廣東丹霞山洞穴景觀巖體穩(wěn)定性的抗壓試驗研究 [J]. 安徽師范大學學報(自然科學版), 2010, 33(2): 170-174.

CHEN S, ZHU C, PENG H, et al. Anti-Pressure experimental study on rocks in level cave of Jinshiyan of Danxiashan in Guangdong province [J]. Journal of Anhui Normal University (Natural Science), 2010, 33(2): 170-174.

[10] 朱誠, 彭華, 李中軒, 等. 浙江江郎山丹霞地貌發(fā)育的年代與成因[J]. 地理學報, 2009, 64(1): 21-32.

ZHU C, PENG H, LI Z X, et al. Age and genesis of the Danxia Landfom on Jianglang mountains, Zhejiang Province [J]. Journal of Geographical Science, 2009, 64(1): 21-32.

[11] 朱誠, 彭華, 歐陽杰, 等.浙江方巖丹霞地貌發(fā)育的年代、成因與特色研究[J]. 地理科學, 2009b, 29(2): 229-237.

ZHU C, PENG H, OUYANG J, et al. Age genesis and characteristics of Danxia landform of Fangyan in Zhejiang province, China[J]. Scientia Geographica Sinica, 2009b, 29(2): 229-237.

[12] ZHU C, PENG H, OUYANG J, et al. Rock resistance and the development of horizontal grooves on Danxia slopes [J]. Geomorphology, 2010, 123(1): 84-96.

[13] 歐陽杰, 朱誠, 彭華, 等.湖南崀山丹霞地貌巖體抗酸侵蝕脆弱性的實驗研究[J]. 地球科學進展, 2011, 26(9): 965-970.

OUYANG J, ZHU C, PENG H, et al. Experimental research on vulnerability of Danxia rocks to resistance against acid erosion in Langshan, Hunan Province [J]. Advances in Earth Science, 2011, 26(9): 965-970.

[14] 彭華, 邱卓煒, 潘志新. 丹霞山順層洞穴風化特征的試驗研究[J]. 地理科學, 2014, 34(4): 454-463.

PENG H, QIU Z W, PAN Z X. Experimental study on the weathering features of bedding caves at Mt Danxiashan [J]. Scientia Geographica Sinica, 2014, 34(4): 454-463.

[15] 齊德利, 陳致均, 王隨繼, 等.崆峒山丹霞地貌地層歸屬演化及地貌年齡[J]. 山地學報, 2015, 33(4): 408-415.

QI D L，CHEN Z J，WANG S J, et al. Stratigraphic classification, evolution stage and geomorphologic age of Kongtongshan Danxia landform in Pingliang, Gansu, China [J]. Mountain Research, 2015, 33(4): 408-415.

[16] 陳智, 彭華, GREIF V等.紅層軟巖夾層的物質(zhì)組成與結(jié)構(gòu)特征對其力學性質(zhì)影響的定量研究——以崀山世界自然遺產(chǎn)地雷劈石為例[J]. 中山大學學報(自然科學版), 2015, 54(4): 139-150.

CHEN Z, PENG H, GREIF V, et al. Analysis of mineralogical composition and texture to quantify the mechanical properties of red bed soft rock: a case study of Leipi Rock at world heritage site, Mt Langshan in Hunan province [J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2015, 54(4): 139-150.

[17] SUMMERFIELD M A. Geomorphology and global tectonics [M]. London: John Wiley & Sons Ltd Press, 2000.

[18] PHILIP T G. Geomorphologic signatures: Classification of aggregated slope unit objects from digital elevation and remote sensing data [J]. Earth Surface Processes and Landforms, 1998, 23(7): 581-594.

[19] 張會平, 楊農(nóng), 張岳橋, 等. 岷江水系流域地貌特征及其構(gòu)造指示意義[J]. 第四紀研究, 2006, 26(1): 126-135.

ZHANG H P, YANG N, ZHANG Y Q, et al. Geomorphology of the Minjiang drainage system (Sichuan, China) and its structural implications [J]. Quaternary Sciences, 2006, 26 (1): 126-135.

[20] DEMOULIN A. An automated method to extract fluvial terraces from digital elevation models: The Vesdre valley, a case study in eastern Belgium [J].Geomorphology, 2007, 91(1/2): 51-64.

[21] 李利波, 徐剛, 胡健民, 等. 基于DEM渭河上游流域的活動構(gòu)造量化分析 [J]. 第四紀研究,2012, 32 (5): 866-879.

LI L B, XU G, HU J M, et al. Quantitative analysis of relative active tectonics of the upstream region of Weihe River based on DEM [J]. Quaternary Sciences, 2012, 32 (5): 866-879.

[22] TAROLLI P. High-resolution topography for understanding Earth surface processes: Opportunities and challenges[J]. Geomorphology, 2014, 216(1): 295-312.

[23] 湯國安. 我國數(shù)字高程模型與數(shù)字地形分析研究進展[J]. 地理學報, 2014, 69(09) : 1305-1325.

TANG G A. Progress of DEM and digital terrain analysis in China[J].Acta Geographica Sinica,2014,69 (09): 1305-1325.

[24] TAUFIK M, PUTRA Y S, HAYATI N. The utilization of global digital elevation model for watershed management: A case study: Bungbuntu sub watershed, pamekasan [J]. Procedia Environmental Sciences, 2015, 24: 297-302.

[25] TARQUINI S, NANNIPIERI L. The 10 m-resolution TINITALY DEM as a trans-disciplinary basis for the analysis of the Italian territory: Current trends and new perspectives [J]. Geomorphology, 2017, 281: 108-115.

[26] 謝超, 周本剛, 李正芳. 東喜馬拉雅構(gòu)造結(jié)地貌形態(tài)及其構(gòu)造指示意義[J]. 地震地質(zhì), 2017, 39 (2): 276-286.

XIE C, ZHOU B G, LI Z F. Geomorphic features of Eastern Himalayan syntaxis and its tectonic implications [J]. Seismology and Geology, 2017, 39 (2): 276-286.

[27] 施煒. 黃河中游晉陜峽谷的DEM 流域特征分析及其新構(gòu)造意義[J]. 第四紀研究, 2008, 28(2): 288-298.

SHI W. DEM drainage analysis of the Shanxi-Shaanxi Gorge in the middle reaches of the Huanghe River and its neotectonic implications [J]. Quaternary Sciences, 2008, 28 (2): 288-298.

[28] FONT M, AMORESE D, LAGARDE J L. DEM and GIS analysis of the stream gradient index to evaluate effects of tectonics: the Normandy intraplate area (NW France) [J]. Geomorphology, 2010, 119(3): 172-180.

[29] 張會平, 張培震, 鄭德文, 等. 祁連山構(gòu)造地貌特征: 青藏高原東北緣晚新生代構(gòu)造變形和地貌演化過程的啟示[J]. 第四紀研究, 2012, 32(5): 907-920.

ZHANG H P, ZHANG P Z, ZHENG D W, et al. Tectonic geomorphology of the Qilian Shan: insights into the Late Cenozoic landscape evolution and deformation in the north eastern Tibetan Plateau[J]. Quaternary Sciences, 2012, 32 (5): 907-917.

[30] 宋效東, 劉學軍, 湯國安, 等. DEM 與地形分析的并行計算[J]. 地理與地理信息科學, 2012, 28(4): 1-7.

SONG X D, LIU X J, TANG G A, et al. Parallel computing of the digital elevation model and digital terrain analysis [J]. Geography and Geo-Information Science, 2012, 28(4): 1-7.

[31] 王巖, 劉少峰, 高明星, 等. 洮河水系流域地貌特征及其構(gòu)造指示意義[J]. 地學前緣, 2010, 17(4): 43-49.

WANG Y, LIU S F, GAO M X, et al. Geomorphology of the Taohe River drainage system and its structural implications [J]. Earth Science Frontiers, 2010, 17 (4): 43-49.

[32] 胡小飛, 潘保田, KIRBY E, 等. 河道陡峭指數(shù)所反映的祁連山北翼抬升速率的東西差異[J]. 科學通報, 2010, 55(23): 2329-2338.

HU X F, PAN B T, KIRBY E, et al. Spatial differences in rock uplift rates inferred from channel steepness indices along the northern flank of the Qilian Mountain, northeast Tibetan Plateau[J]. Science China Press, 2010, 55(23): 2329- 2338 .

[33] 宋卓沁. 青藏高原東緣典型河流地貌及其活動構(gòu)造指示[D]. 北京：中國地震局地震預(yù)測研究所, 2014.

[34] MAZZOLI S, ASCIONE A, BUSCHER J T, et al. Low-angle normal faulting and focused exhumation associated with late Pliocene change in tectonic style in the southern Apennines (Italy) [J]. Tectonics, 2014, 33(9): 1802-1818.

[35] ROBERTS G P, MESCHIS M, HOUGHTON S, et al. The implications of revised Quaternary palaeoshoreline chronologies for the rates of active extension and uplift in the upper plate of subduction zones [J]. Quaternary Science Reviews, 2014, 78 (11): 169-187.

[36] SCHELLEKENS J, BROLSMA R J, DAHM R J, et al. Rapid setup of hydrological and hydraulic models using OpenStreetMap and the SRTM derived digital elevation model [J]. Environmental Modelling & Software, 2014, 61(C): 98-105.

[37] ZHANG H P, ZHANG P Z, ZHENG D W, et al. Transforming the Miocene Altyn Tagh fault slip into shortening of the northeastern Qilian Shan: Insights from the drainage basin geometry [J]. Terra Nova, 2014, 26(3): 216-221.

[38] DOMENEGHETTI A, CARISI F, CASTELLARIN A, et al. Evolution of flood risk over large areas: quantitative assessment for the Po river [J]. Journal of Hydrology, 2015, 527: 809-823.

[39] 蘇琦, 袁道陽, 謝虹. 祁連山東段石羊河流域及鄰區(qū)地貌特征及其構(gòu)造意義[J]. 地質(zhì)論評, 2017, 63(1): 7-20.

SU Q, YUAN D Y, XIE H. Geomorphic features of the Shiyang river drainage basin and adjacent aea in Eastern Qilian Mountains and its insight into tectonic implications [J]. Geological Review, 2017, 63(1): 7-20.

[40] FARR T G, KOBRICK M. Shuttle radar topography mission produces a wealth of data [J]. EOS Transactions of the American Geophysical Union, 2000, 81(48): 583-585.

[41] BERNHARD R, MICHAEL E, ACHIM R, RICHARD B. The Shuttle Radar Topography Mission: a new class of digital elevation models acquired by spaceborne radar [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2003, 57(4): 241-262.

[42] FARR T G, ROSEN P A, CARO E, et al. The shuttle radar topography mission [J]. Reviews of Geophysics, 2007, 45(2): 1-34.

[43] O’LOUGHLIN F E, PAIVA R C D, DURAND M, et al. A multi-sensor approach towards a global vegetation corrected SRTM DEM product [J]. Remote Sensing of Environment, 2016, 182: 49-59.

[44] 周毅, 湯國安, 王春, 等. 基于高分辨率DEM的黃土地貌正負地形自動分割技術(shù)研究[J]. 地理科學, 2010, 30(2): 261-266.

ZHOU Y, TANG G AN, WANG C, et al. Automatic segmentation of Loess positive and negative terrains based on high resolution grid DEMs [J]. Scientia Geographica Sinica, 2010, 30(2): 261-266.

[45] 羅偉, 鄢志武, 薛重生. 基于DEM的巖溶地區(qū)地形地貌特征提取與分析[J]. 華中師范大學學報(自然科學版), 2013, 47(3): 436-439.

LUO W, YAN Z W, XUE C S. Extraction and analysis of the topography in Karst area based on DEM [J]. Journal of Huazhong Normal University (Natural Sciences), 2013, 47(3): 436-439.

[46] 祝世杰, 湯國安, 李發(fā)源, 等. 基于DEM的黃土高原面積高程積分研究[J]. 地理學報, 2013, 68 (71):921-932.

ZHU S J, TANG G AN, LI F Y, XIONG L Y, et al. Spatial variation of hypsometric integral in the Loess Plateau based on DEM [J]. Acta Geographica Sinica, 2013, 68 (71): 921-932.

[47] 熊禮陽, 湯國安, 袁寶印, 等. 基于DEM的黃土高原(重點流失區(qū))地貌演化的繼承性研究[J]. 中國科學(地球科學), 2014, 57(2): 313-321.

XIONG L Y, TANG G A, YUAN B Y, et al. Geomorphological inheritance for loess landform evolution in a severe soil erosion region of Loess Plateau of China based on digital elevation models [J]. Science China(Earth Sciences), 2014, 57(2): 313-321.

[48] 李飛, 楊小平, 郝紅科. 基于DEM對喀斯特地貌山地及水系信息的提取——以桂林陽朔縣為例[J]. 測繪與空間地理信息, 2015, 38(8): 73-76.

LI F, YANG X P, HAO H K. Based on the DEM of Karst landform mountain and water system of information extraction: with Yangshuo County as an example [J]. Geomatics & Spatial Information Technology, 2015, 38(8): 73-76.

[49] 楊先武, 錢葉青, 鄭春霞. 喀斯特峰林峰叢地貌形態(tài)表達研究綜述[J].地理與地理信息科學, 2017, 33(4): 22-27.

YANG X W, QIAN Y Q, ZHENG C X. Research review on the morphology of Karst fenglin and fengcong landforms[J]. Geography and Geo-Informatjion Science, 2017, 33(4): 22-27.

[50] HENGL T, REUTER H I. Geomorphometry: concepts, software, applications-developments in soil science,Vol 33 [M]. Hungary: Elsevier, 2009.

[51] HENGL T, REUTER H I. Geomorphometry: concepts, software, applications-developments in soil science, Vol 33[M]. Hungary, Elsevier, 2009.

[52] 黃進. 丹霞地貌坡面發(fā)育的一種基本方式[J]. 熱帶地理, 1982, 3(2): 107-134.

HUANG J. A basic model of the slope development of Danxia landform[J]. Tropical Geography, 1982, 3(2): 107-134.

[53] STRAHLER A N. Hypsometric (Area-Altitude) analysis of erosional topography [J]. Bulletin of the Geological Society of America, 1952, 63: 1117-1142.

[55] BROCKLEHUST S H, WHIPPLE K X. Hypsometry of glaciated landscapes [J]. Earth Surface Processes and Landforms,2004, 29(7): 907-926.

[56] ANSAN V, MANGOLD N. New observations of Warrego Valleys, Mars: evidence for precipitation and surface runoff [J]. Planetary and Space Science, 2006 (54): 219-242.

[57] 張敬春, 李川川, 張梅, 等. 格爾木河流域面積-高程積分值的地貌學分析[J]. 山地學報, 2011, 29(3): 257-268.

ZHANG J C, LI C C, ZHANG M, et al. Geomorphologic analysis of the Golmud River drainage basin based on hypsometric integral value [J]. Journal of Mountain Science, 2011, 29(3): 257-268.

[58] 祝世杰, 湯國安, 李發(fā)源, 等. 基于DEM的黃土高原面積高程積分研究[J]. 地理學報, 2013, 68 (71): 921-932.

ZHU S J, TANG G A, LI F Y, et al. Spatial variation of hypsometric integral in the Loess Plateau based on DEM [J]. Acta Geographica Sinica, 2013, 68 (71): 921-932.

[59] 高明星, 陳桂華, 徐錫偉. 地貌參數(shù)指示的臨潭-宕昌斷裂帶最新構(gòu)造隆升差異與地震活動[J]. 地震地質(zhì), 2015, 37(3): 709-718.

GAO M X, CHEN G H, XU X W. Geomorphic indices indicated recent differential tectonic uplift of the Lintan-Dangchang Fault and the Minxian-Zhangxian earthquake [J]. Seismology and Geology, 2015, 37(3): 709-718.

[60] 蘇琦, 袁道陽, 謝虹. 祁連山-河西走廊黑河流域地貌特征及其構(gòu)造意義[J]. 地震地質(zhì), 2016, 38(3): 560-581.

SU Q, YUAN D Y, XIE H. Geomorphic features of the Heihe river drainage basin in western Qilian shan-hexi corridor and its tectonic implications [J]. Seismology and Geology, 2016, 38(3): 560-581.

[61] 常直楊, 王建, 白世彪, 等. 面積高程積分值計算方法的比較[J]. 干旱區(qū)資源與環(huán)境, 2015, 29(3):171-175.

CHANG Z Y, WANG J, BAI S B, et al. Comparison of hypsometric integral methods[J]. Journal of Arid Land Resources and Environment, 2015, 29 (3):171-175.

[62] 熊康寧, 肖時珍. 世界遺產(chǎn)與赤水丹霞奇觀[M]. 北京: 高等教育出版社, 2012.

XIONG K N, XIAO S Z. World heritage and Chishui Danxia landscape [M]. Beijing: Higher Education Press, 2012.