在构造活动强烈的造山带，周期性复发的强震的同震垂直位移累积被认为是山脉持续隆升的一个重要驱动因素^[1].地震同时也导致了大量地表物质的迁移^[2-3]，改变了地形地貌^[4].当同震构造变形驱动的地表抬升量大于同震滑坡的地表剥蚀量，结果是山脉平均高程的增加，反之，是山脉平均高程的降低^[4].然而，这种关系的确立将取决于地震滑坡物质河流卸载时间与强震复发周期长短的比较^[5].因此定量化研究地震滑坡物质河流卸载时间对理解地震与造山带地貌演化的关系有着十分重要的意义. 2008年汶川地震最大地表同震垂直位移量达约9~10m^[6-8]，地震触发了大量的滑坡^{[7, 9-10]}.经过计算，Parker等^[4]认为汶川地震同震剥蚀的物质质量大于地震同震抬升导致的物质增加量，进而提出类似汶川地震的强震在龙门山并没有产生造山的作用.另外，Ouimet^[11]认为强震触发的滑坡增加了龙门山地区千年尺度的侵蚀速率.但是这些研究缺少对地震滑坡物质的河流卸载时间的估算，而估算地震滑坡物质的河流卸载时间将有助于对青藏高原东缘大地震、滑坡侵蚀和地貌演化之间关系的认识.

岷江为流经龙门山地区流量最大的河流，2008年汶川地震沿岷江干流两岸诱发了大量滑坡（图 1），这为研究地震滑坡物质河流卸载提供了良好机会.因此，岷江流域是研究地震活动与地形地貌演化关系的一个理想地区.本研究将定量地估算岷江对汶川地震滑坡物质卸载可能所需时间，并讨论大地震事件滑坡物质的卸载对龙门山晚新生代构造地貌演化的影响.

图 1

Fig. 1

图 1 沿岷江干流(图中蓝色粗线)发育了2008年汶川地震所诱发的大量滑坡 Fig. 1 Large-scale landslides induced by the 2008 Wenchuan earthquake along the Minjiang River (thick blue lines)

龙门山构造带是中国典型的推覆构造发育地区之一，形成于中生代和早新生代，主要由大致平行的三条逆冲、逆冲-走滑断裂带组成^[12]：自北西向南东发育汶川-茂县断裂（F1）、映秀-北川断裂（F2）和灌县-安县断裂（F3）.汶川地震是龙门山断裂带的映秀-北川断裂突发错动的结果^[13-14]，在地表上形成了300多公里长的地表破裂带，地表破裂带变形特征以逆冲为主兼具右旋走滑分量^[13-17].龙门山地貌上是青藏高原东缘与四川盆地的分界线，也是青藏高原周缘地形梯度变化最大的地区，在不到50km的范围内地形平均高程从4000m变化到500m.其中龙门山中段是龙门山地形地貌特征典型的地区^[18]，也是地震活动频发区^[19-20].

岷江发源于岷山南麓，向南流经岷山构造带和龙门山构造带，都江堰以上的上游段河道长约340km，在龙门山断裂带内表现为深切河谷、基岩型河道^[21].在龙门山构造带岷江沿汶川-茂县断裂带（F1）成北东-南西向流动，然后沿东南方向斜切映秀-北川断裂带（F2）和安县-灌县断裂带（F3），经都江堰流入四川盆地（图 2）.本文研究的岷江是指位于都江堰以上发育于龙门山构造带的岷江干流.

图 2

Fig. 2

图 2 研究区主要活动构造及位置图(蓝色虚线矩形为所研究的岷江河段) F1为汶川-茂县断裂, F2为映秀-北川断裂, F3为灌县-安县断裂, F4为岷江断裂. Fig. 2 Topography and active faults in the central Longmen Shan orogenic belt (The dashed rectangle shows location of study area) F1:Wenchuan-Maoxian fault, F2:Yingxiu-Beichuan fault, F3:Guanxian-Anxian fault, F4:Minjiang fault

河流搬运物质通常以溶解质、悬移质和底移质三种形态存在.溶解质所占比例非常小^[22]，在特殊的气候和岩性条件下才会考虑溶解质对地貌演化的贡献.地震滑坡物质大量进入河道，以河流悬移质形式存在的地震滑坡物质随河水而流动，其河流卸载时间可能很快.Hovius等^[23]报道了台湾集集地震震后10年河流悬移质含量就恢复到了震前水平.河流底移质的移动速度相对悬移质来说较慢，以河流底移质形式存在的地震滑坡物质的河流卸载时间可能相对较长.在长期的地貌演化过程中，河流物质中底移质对河流过程和形态的影响最大^[24].此外，河流底移质一般占河流中总物质量的30%^[25]，甚至可以达到50%^[22].因此，我们选择定量计算以底移质形式存在的地震滑坡物质的河流卸载时间.同时假设地震滑坡物质进入河道后，30%的物质以河流底移质的形式存在.本文采用Yanites等^[26]的方法估算地震滑坡物质的河流卸载时间（T）.T可表示为区域滑坡物质的质量（M_L）与河流底移质的搬运能力（Q_T）的比值，公式为

(1)

区域滑坡物质的质量由公式（2）和（4）得到

(2)

式中，ρ_L为滑坡物质的密度，V_TL为地震滑坡总方量.龙门山中部第四纪地层坡度较缓，不易发生滑坡，汶川地震中大量的滑坡发生在较老的岩层中，92%的滑坡发生在第三纪以前的地层中^[27]，可认为汶川地震滑坡主要为岩质滑坡，而龙门山地表岩石的密度为2450kg/m^3[28]，所以我们假设地震滑坡物质的密度为2450kg/m³.

区域滑坡物质的单个滑坡体方量由公式（3）^[29]得到

(3)

其中V_L为单个滑坡方量，A_L为单个滑坡面积，R²为决定系数.Larsen等^[30]认为滑坡方量（V）和面积（A）具有幂律的关系V=aA^b（a、b是常数），对于岩质滑坡，b=1.3-1.6.因此我们选择此公式计算单个滑坡的方量.以此经验公式为基础，建立岷江流域地震滑坡总方量公式为

(4)

式中N是研究区内滑坡个数.

在基岩型河道剪切应力侵蚀模型中，用河床剪切应力表征河流对河床侵蚀能力大小^[31-32].河床剪切应力（τ_b）由公式（5）^[33]表示为

(5)

式中，W为河宽，河宽数据采集自1:50000地形图，部分地区来源于空间分辨率为28.5m的Landsat ETM遥感影像（http：//glovis.usgs.gov/）.S为河流的坡度，提取自空间分辨率为90m的SRTMDEM数据（http：//srtm.csi.cgiar.org/）.由于部分高程数据在90m距离内并无变化，使得相邻两个像元的高程差为0，因此计算的坡度也为零.经过反复试验，以水平距离1000m提取河流坡度可消除此效应.n是曼宁系数，对于山区河通常取值为0.04^[34]. Q_W为河流多年平均流量（m³/s）.

岷江流域是人类社会经济活动较频繁的地域，在岷江及其支流上建造了大量的水电站，势必影响河流的流量（图 3）.例如，岷江上游映秀湾水电站建于1965-1972年，是一座比较大的水电站.岷江出山口紫坪铺水文站记录的河流流量数据显示：映秀湾水电站建站前10年河流的平均流量为463m³/s，建站后10年河流的平均流量为442m³/s.由于河流的流量影响基岩型河道的剪切应力计算，河流的实际搬运能力可能会被低估.岷江流域大规模的水电建设开始于1980年以后^[35].为减少人类经济活动对计算带来的误差，因此，河流多年平均流量选用1980以前岷江各站水文数据，数据来源于长江流域水文资料^[36].

图 3

Fig. 3

图 3 映秀湾水电站建设前后河流流量变化图(灰色阴影区为水电站建设期间, 数据资料来源自http://www.mwr.gov.cn/) Fig. 3 The inter-annual change of flow discharge before and after the construction of the Yingxiuwan dam.The grey shadow region is construction period of the dam (http://www.mwr.gov.cn/)

河流底移质搬运能力（Q_T）由公式（6）^[37]表示：

(6)

式中，ρ_s为底移质的密度，值为2.65g/cm³，ρ为水的密度，值为1g/cm³，g为重力加速度，值为9.8m/s². D为底移质的中值粒径（cm），岷江河床底移质值中值粒径为10cm^[38].τ_b为河床剪切应力，τ_c*为河床临界启动剪切应力，对于基岩型河道，河床临界启动剪切应力数值通常变化范围为0.03~0.08^[39].我们取常数0.03作为岷江的河床临界启动剪切应力.

图 4显示岷江搬运能力整体变化平稳，往出山口方向具有变大的趋势.在河流斜穿断层F2和F3的位置，河流有着较高的搬运能力.因为断层所处位置也是龙门山地形梯度变化较大的地方，河流的比降变化较快（图 5a）.河流在比降大的地方有着较大的重力势能，因而对沉积物具有较强的搬运能力.另外，河床剪切力空间分布也说明了位于龙门山构造带内岷江河段的水动力条件具有类似的空间分布特征（图 5b）.当上游的河流搬运能力小于下游，地震滑坡物质则由河流下游相对上游的“剩余”的搬运能力来搬运，此时根据公式（1）计算的结果为河流搬运时间的下限^[26].根据公式（2）计算研究区内滑坡物质质量约为17885Mt，根据公式（6）计算岷江上下游搬运能力的差值为1.73×10⁹kg/yr.如果岷江保持现有的搬运能力并且有能力完全搬运地震滑坡物质，那么岷江流域的所有地震滑坡物质被河流搬运出龙门山至少需要3100yr.

图 4

Fig. 4

图 4 龙门山构造带中段岷江河段搬运能力的空间分布特征 F1为汶川-茂县断裂, F2为映秀-北川断裂, F3为灌县-安县断裂. Fig. 4 The spatial distribution of the sediment transport capacity along the Minjiang River in the central Longmen Shan orogenic belt F1:Wenchuan-Maoxian fault, F2:Yingxiu-Beichuan fault, F3:Guanxian-Anxian fault

图 5

Fig. 5

图 5 龙门山构造带中段岷江河段的河流比降(a)和河床剪切应力(b)空间分布特征 F1为汶川-茂县断裂, F2为映秀-北川断裂, F3为灌县-安县断裂. Fig. 5 Spatial distribution of the river slopes (a) and the shear stress (b) along the Minjiang River in the central Longmen Shan orogenic belt F1:Wenchuan-Maoxian fault, F2:Yingxiu-Beichuan fault, F3:Guanxian-Anxian fault

前人研究结果表明龙门山断裂带中段强震的复发周期可能约为2000~3000yr^{[19, 40-41]}，与地震滑坡物质河流卸载时间相比，暗示岷江对汶川地震滑坡物质搬运时间可能略长于一个地震周期.换句话说，至少在龙门山中段在一个地震周期内绝大部分地震滑坡物质可以被卸载出龙门山.需要说明的是我们定量估算的结果仍受到一些不确定因素的影响，比如地质历史时期气候变化对河流流量的可能影响，地震滑坡堵江形成堰塞湖对河流底移质搬运能力的影响，滑坡物质从坡面进入河道的时间，以及现代大坝建设对底移物质搬运的影响.尽管如此，在龙门山晚新生代长期的地貌演化过程中，地震滑坡物质被卸载的速度可能是很快的.因为在龙门山山谷和河谷中缺少大面积裸露的基岩和滞留过厚的（厚度大于100m）沉积物^{[21, 42]}.另外，龙门山山前快速堆积的一套厚度大于500m的晚新生代沉积物可能与长期重复发生的大地震活动有关^[43]，从侧面也说明了岷江的河流卸载作用很强.龙门山中段类似2008年汶川大地震复发周期是根据经验公式计算出来，存在一定误差，地震物质河流卸载时间结果也存在一些不确定因素的影响，但是从我们的估算结果来看二者的量级几乎一致.

在龙门山造山带晚新生代隆升历史过程中，龙门山地区每次逆冲型大地震都会形成不同规模的地表垂直位移和同震滑坡，如果大地震产生的滑坡物质被侵蚀和剥蚀所平衡，那么为什么目前龙门山仍然会保持高大而陡峻的地形？我们认为与河流对龙门山剥蚀物质的长期快速的卸载作用有关.由于河流对龙门山山脉体系内剥蚀物质的快速卸载，引起壳幔边界Moho面均衡反弹^[44-45]，促使地表抬升^[7]，使得龙门山主峰高出青藏高原内部，从而形成龙门山现今的高陡山峰.但是剥蚀卸载引起的均衡反弹只是使山峰的高度比剥蚀前更高，山脉的平均高程应该是在逐步降低的^[46].值得注意的是，龙门山地区古地形重建结果显示，在3.6Ma以前，龙门山与青藏高原内部高程相近，但小于现代平均高程^[47].说明可能还有别的地球动力学作用机制在维持着龙门山整体持续的隆升.我们推断可能是由于地表剥蚀物质的快速卸载诱发了龙门山下部深部物质的上涌导致了龙门山持续的隆升（图 6）.横跨龙门山断裂带的被动源地震探测数据和大地电磁数据揭示^{[44-45, 48-49]}：龙门山及其以西部分地区下部为低速异常，指示可能存在由于岩石圈深部物质上涌；四川盆地西部与松潘-甘孜为高速异常，指示可能有地幔物质的下潜.这些地球物理探测结果暗示青藏高原东部可能存在局部的地幔物质垂向对流.因此，地幔岩浆物质上涌可能是高原东部抬升的主要因素之一.另外，下地壳的流变特性在地壳内物质均衡重置过程中也起到了重要作用^[50].层析成像结果揭示青藏高原东部中下地壳刚性较小^[51-52]，具备流动的条件^[53-57]，所以下地壳的部分物质也可能受到了上地幔垂向运动的作用而在龙门山造山带下方蓄积.因此，地表剥蚀（或侵蚀）物质的快速卸载驱动地壳均衡反弹，形成了龙门山构造带现今的高陡山峰；另一方面，诱发龙门山构造带下方的深部物质上涌，形成顶托作用，有利于龙门山保持高陡的地形特征.

图 6

Fig. 6

图 6 龙门山晚新生代隆升的地球动力学机制示意图 (a)伴随着周期性大地震的复发, 龙门山下方地壳缩短增厚; 与此同时, 河流长期对地震滑坡及其它地表剥蚀物质的快速卸载引起龙门山地壳均衡反弹, 形成高陡山峰; (b)均衡反弹诱发深部物质的上涌, 从而维持了龙门山高陡的地形特征. Fig. 6 Geodynamic model for the uplift processes of Longmen Shan orogenic belt during late Cenozoic (a) with the repeated large seismic events, the crust under the Longmen Shan is shortened.At the same time, the long-term fluvial unloading of erosion material produced by earthquake-related landslides may induce isostatic crutsal-mantle boundary rebound, which lead to formation of the high and steep orogenic summit; (b)the isostatic rebound further triggers the deep mantle upwelling and maintains the present-day steep topography

山脉如果失去持续的内动力作用，将很快会被地表作用夷平^[58].然而龙门山在地壳水平变形速率很小的情况下^[59-60]，依然保持了高陡的地形特征，并没有被水系和其它地表作用消减，说明内动力在持续的作用于龙门山的抬升.而这很可能与大地震的频繁发生有关，大地震的周期性复发标志了龙门山构造变形强烈，所以周期性大地震造成的地表抬升的累积可能是龙门山晚新生代以来高陡地貌形成的主要因素^[7].尽管如此，地表剥蚀快速卸载诱发的深部物质活动在龙门山高陡地貌特征的形成中也起到了重要作用.诚然，对这一过程仍需作进一步的研究.需要说明的是，大地震引发的大范围、大规模的滑坡是龙门山地表过程中的重要力量，但是除了地震滑坡侵蚀，其它成因的地表剥蚀物质的快速卸载对龙门山的均衡反弹也应有贡献.均衡反弹是形成和维持造山带的高地形特征的重要深部过程^{[58, 61]}，所以对龙门山地貌上的高陡地形梯度带成因的研究不应忽视均衡反弹的作用^[62].因此，地表作用和深部地球物理过程的共同作用塑造了龙门山现今独特的地貌特征.

根据上述分析，得出如下初步结论和认识：

（1） 根据经验公式估算出汶川地震滑坡物质的总质量约为17885Mt.汶川地震滑坡物质进入岷江河道后，如果30%的地震滑坡物质以底移质的形式存在，在现今龙门山中段岷江的河流水动力条件下，汶川地震滑坡物质被岷江卸载出龙门山至少需要3100yr.

（2） 龙门山中段类似2008年汶川大地震的复发周期约为3000yr，暗示大地震产生的同震滑坡物质可能在相同时间段内被河流卸载作用平衡，并被卸载出龙门山构造带.

（3） 我们推断在龙门山长期的地貌演化过程中，除了周期性大地震造成的地表构造抬升累积外，周期性大地震相关的地表侵蚀物质快速卸载所驱动的地壳均衡反弹及诱发的深部地幔对流也造成龙门山的隆升作用，这可能是龙门山晚新生代隆升的重要地球动力学机制之一.

感谢评审人提出的建设性修改意见和建议，感谢美国科罗拉多大学的BrianJ.Yanites博士的有益讨论.