New Reconstruction of Glacier Debris Flows Based on Tree Ring Response
-
摘要: 冰川型泥石流在瓤打曲流域频发，对当地安全造成了严重威胁，重建泥石流的暴发历史、流动范围和成因，可以为工程建设提供防灾减灾参数.基于树木年轮的愈伤组织和生长抑制响应，给出了年轮生长抑制的评判阈值和年轮定灾的Wit指数的新计算方法，并重建了1890—2021年的泥石流的暴发时间和流动范围；用泥石流实际暴发时间验证了年轮生长抑制的评判阈值和Wit新计算方法的合理性. 结果表明：（1）冰川活跃区年轮生长抑制的阈值修订为25%；（2）Wit指数计算需要排除生长释放因子，利用年轮的创伤组织和生长抑制会增大泥石流灾害事件的定年准确率.Abstract: Glacial debris flows occur frequently in the Rangdaqu River Basin, posing a serious threat to local safety. Reconstructing the eruption history, flow range, and causes of debris flows can provide disaster prevention and mitigation parameters for engineering construction. Based on the cambial injury tissue and growth inhibition response of trees, this study provides a new calculation method for the evaluation threshold of growth inhibition in tree rings and the Wit index of debrisflow identification in tree rings, and reconstructs the outbreak time and flow range of debris flows from 1890 to 2021. The results show that: (1) the threshold for growth inhibition in the active zone of glaciers has been revised to 25%; (2) the calculation of the Wit index needs to exclude growth release factors, and the use of traumatic tissue and growth inhibition in tree rings can improve the accuracy of dating debris flow events.
-
Key words:
- debris flows /
- glaciers /
- damaged tissue /
- growth inhibition /
- disaster reconstruction /
- engineering geology
-
0. 引言
冰川型泥石流在帕隆藏布流域分布较广，受气候变化影响显著（吕立群，2017），往往在冰川退缩期和强烈地震之后出现泥石流活动的高潮（李尧等，2022）. 泥石流成灾过程中，对流动路径上的树木生长造成扰动，年轮出现伤疤，生长释放或者抑制，统称为生长扰动（GD）. 根据树木年轮宽度变化或愈伤组织可以重建泥石流灾害历史（吕立群等，2024），并反映泥石流灾害对区域气候变化的响应（黄涛，2019）. 树木年代法相对于光释光、14C、宇宙成因核素暴露、火山灰测年等方法，其分辨率高和连续性强，被广泛用于泥石流灾害历史重建（铁永波等，2014；赖忠平等，2021）.
目前利用年轮重建灾害历史的指标主要有3种：（1）生长扰动（GD）数量；（2）响应指数It（生长扰动的树木与存活的样本数量的比率）；（3）Wit指数（对树木所受不同程度的扰动进行加权计算，以剔除部分除灾害外的影响）. 响应指数It是最早的年轮定灾方法，在国内外早期的研究中应用较多（Vădean et al.，2015；Šilhán，2021；王志兰等，2022）. 近几十年，开始通过采用树木年轮偏心圆和树木的正常生长年轮进行对比重建灾害事件（铁永波等，2014），或者通过年轮生长抑制和释放百分比重建灾害事件（吴佳亮等，2021），但都无法精准定灾. 直到Kogelnig-Mayer et al.（2011）提出用加权指数Wit来重建灾害事件，发现单独采用加权指数Wit（孟哲等，2022）或与其他指标（Peitzsch et al.，2013）一起进行灾害事件定年，误差年限大大降低. 国内部分学者利用Wit不仅重建了泥石流灾害事件的时间，还重建了其流动范围，并对树龄进行了敏感性分析（吕立群等，2024）. 由于年轮重建灾害史存在较强区域性，并没有统一的定灾标准，这与地理位置、气候条件和灾害类型有关（表 1），因此利用加权指数Wit并不能照搬前人标准，需根据地域特征重新尝试确认新的参数阈值.
表 1 基于树木年轮的灾害定年方法Table Supplementary Table Standard for tree ring disaster dating作者 国家 地区 样品数量 时期(年份) 海拔(m) 年均降水(mm) GD总数 定灾标准 灾害类型 It GD数量 Wit Mayera et al.（2010） 奥地利 Gratzental 227 1800-2007 2 106~1 166 1 526 1 155 It≥4% 未提及 未提及 泥石流 Corona et al.（2012） 法国 Arve 209 1771-1998 1 100~3 650 1 262 645 It > 1.7 GD数 > 5% 未提及 雪崩 Peitzsch et al.(2013) 美国 Rocky 673 1636-2017 1 100~2 700 未提及 2 134 样品数=10~20：
GD数≥3个，It≥15%高：Wit≥0.3 雪崩 样品数=21~50：
GD数≥5个，It≥10%中：Wit≥0.2 样品数=51~100：
GD数≥7个，It≥7%样品数 > 100：
GD数≥9个，It≥4.5%低：Wit < 0.2 Vădean et al.（2015） 罗马尼亚 Apuseni 20 1940-2010 平均835 800~1 000 72 It≥15% GD数≥3个 未提及 泥石流 Šilhán et al.（2021） 捷克共和国 Hostynsko-vsetınskahornatina Mts. 153 1920-2017 平均550 700~1 200 179 可能发生：5%~10%
一定发生：It≥10%未提及 未提及 滑坡 Franco-Ramos et al.（2019） 墨西哥 La Malinche 65 1933-2017 3 000~4 400 754 205 未提及 确定：GD≥4和Wit≥2
“潜在”泥石流：Wit≥1.5但GD≤4
噪声：GD数≤2个，Wit≤2泥石流 Schneuwly-Bollschweilera et al.（2013） 瑞士 Zermatt 385 1600-2008 1 420~4 545 690 581 未提及 未提及 Wit > 1 泥石流 注：GD. 生长扰动；It. 响应指数；Wit. 加权指数 本文以西藏瓤打曲流域的一条冰川型泥石流支沟为例，通过指数Wit重建1890—2021年间泥石流灾害活动，并对指数Wit的准确性进行验证，以期为冰川型泥石流灾害历史重建提供新的参数阈值和计算方法.
1. 研究区域
瓤打曲是帕隆藏布二级支流，流域面积304.2 km2，主沟长28.4 km（张金山等，2015），沟道两岸分布有大量冲洪积和冰碛物等第四纪堆积体（曾宪阳等，2019）. 瓤打曲地形陡峭，降水集中在5—9月，占年降水量的70%以上（苗晓岐，2022），流域内分布着大面积冰川，冰川融水补给地表径流，为泥石流形成提供有利的水源条件. 瓤打曲右侧支沟（图 1），位于舍托曲与瓤打曲的入汇处，流域面积8.0 km2，平均海拔3 952 m，下伏地层为灰岩、大理岩、砂质板岩，夹有厚层的变质石英砂岩. 该区域曾经在2002、2007、2009、2012年暴发过泥石流，频率较高，对当地安全造成了严重威胁. 沟道内2000年发生过大面积森林火灾，沟道左岸曾发生过大规模崩塌，为泥石流发育提供了大量物源. 将研究区划分为物源区AB段、流通区BC段和堆积区CDEF段，其中堆积区按照灾害流动方向划分为CD段、CE段和CF段（图 1）.
图 1 研究区及树木采集样点左上图底图审图号：GS（2020）4619号Fig. 1. Study area and distribution of tree samples2. 研究方法
2.1 树木年轮的采集与处理
2.1.1 年轮采样方法
使用5.15 mm生长锥在受泥石流影响的伤疤木、掩埋木和根系暴露木的树干进行取样，将取出的树芯样本进行分类编号，并存放于纸管中. 取样时记录样本位置，胸径，邻近树木状态. 树木采样点均位于泥石流沟道内和扇体上，共采集树木标本72棵，并采集未受泥石流影响的树木20棵作为参考木（图 1）.
2.1.2 样本处理
将所采树芯样品带回实验室后，粘在特制木质凹槽中，用胶带固定好两端且做好标记；待其自然晾干后，使用树木年轮打磨机打磨至光滑平整且年轮清晰可见. 使用精度为0.001 mm的LinTabTSAP Win 6.0分析系统来测量年轮宽度；测量完毕后，利用COFECHA程序对交叉定年结果进行检验，过滤掉年轮序列中的低频，并以99%的置信区间进行检验（王珂，2016）；将这些交叉定年的样品与参考树木年轮进行比对，减小误差，最终制成生长扰动年代分布表.
2.2 泥石流发生时间
2.2.1 生长抑制和生长释放
损伤和愈伤组织（图 2b）可以准确确定泥石流发生时间（Stoffel et al.，2006；Lundström et al.，2008），年轮中愈伤组织的第一层细胞形成的时间定义为泥石流灾害发生的时间（丁苗等，2016）. 年轮宽度突然减小所造成树木的生长抑制也是判断泥石流事件发生的重要依据（Strunk et al.，1997；Schneuwly et al.，2009）. 生长抑制一方面是由于泥石流对树木的巨大重量引起的物理反应，另一方面是这些泥沙颗粒对树干和根系造成了破坏. 当泥石流将某个区域的树木摧毁时，邻近的树木会形成一个竞争更少，光线、养分和水分更足的新环境，年轮的宽度也随之增加，即生长释放. 生长释放和生长抑制的评判依据为：
I=ri−ri−1ri−1， (1) 式中：i年的年轮宽度与i-1年的年轮宽度相比，I增加50%及以上体现了生长释放，I减少40%及以下体现了生长抑制（图 2a）.
2.2.2 Wit加权指数
极端气候、人类活动、动物的破坏对树木年轮的扰动也造成影响，Kogelnig-Mayer et al.（2011）提出用加权指数Wit来剔除这些干扰：
Wit=[(∑ni=1Ti×7)+(∑ni=1Ts×5)+ (∑ni=1Tm×3)+(∑ni=1Tw)]×∑ni=1Rt∑ni=1At， (2) 式中：Ti为具有损伤和愈伤组织的树木；Ts为具有强烈生长扰动的树木；Tm为具有中等生长扰动的树木；Tw为具有弱生长扰动的树木；Rt为在t年具有生长扰动的树木数量；At为在t年获取样本的树木数量. 扰动强度划分阈值见表 2.
表 2 生长扰动强度的分类阈值Table Supplementary Table Classification of growth disturbances (GDs)Ts Tm Tw Ti 年轮变化(%) 持续时间 年轮变化(%) 持续时间 年轮变化(%) 持续时间 生长抑制 ≥60% ≥5年 ≥60% < 5年且≥2年 ≥40%且 < 60% ≥2年 生长释放 ≥60% ≥5年 ≥60% < 5年且≥2年 ≥50%且 < 60% ≥2年 创伤及愈伤组织 存在创伤及愈伤组织 2.3 泥石流活动范围解译
为了解泥石流活动范围，使用了28幅遥感影像对泥石流活动区域进行了反演（表 3）. 利用ENVI5.3软件解译遥感数据，通过辐射定标、大气校正后可以观察到泥石流灾害发生后新形成的沟道和扩宽的河道. 因计算机对于阴影区域的边界难以处理，易导致分类误差大，故本文采取目视解译法. 根据不同的波段组合（Landsat5：波段4、3、2；Landsat 8：波段5、4、3），保证泥石流流动范围提取准确性.
表 3 遥感数据集Table Supplementary Table Collected Satellite datasets in this work编号 数据来源 日期(日/月/年) 分辨率(m) 编号 数据来源 日期(日/月/年) 分辨率(m) 1 Landsat5 17/01/1987 30 15 Landsat5 17/02/2004 30 2 Landsat5 05/02/1988 30 16 Landsat5 23/12/2006 30 3 Landsat5 08/12/1989 30 17 Landsat5 30/04/2007 30 4 Landsat5 11/12/1990 30 18 Landsat5 26/11/2008 30 5 Landsat5 14/12/1991 30 19 Landsat5 29/11/2009 30 6 Landsat5 30/11/1992 30 20 Landsat5 18/12/2010 30 7 Landsat5 01/11/1993 30 21 Landsat5 04/02/2011 30 8 Landsat5 06/12/1994 30 22 Landsat8 24/11/2013 30 9 Landsat5 07/11/1995 30 23 Landsat8 29/12/2014 30 10 Landsat5 07/11/1996 30 24 Landsat8 30/11/2015 30 11 Landsat5 15/11/1998 30 25 Landsat8 21/12/2017 30 12 Landsat5 18/01/1999 30 26 Landsat8 22/11/2018 30 13 Landsat5 22/12/2000 30 27 Landsat8 11/11/2020 30 14 Landsat5 15/12/2003 30 28 Landsat8 06/03/2022 30 3. 原方法的缺陷和原因分析
利用表 2的分类阈值计算得到1890—2021年的Wit值，与已知的泥石流发育时间（2002、2007、2009、2012）作对比，2002年Wit值为3.71，2007年Wit值为2.28，2009年Wit值为2.10，2012年Wit值为5.49（图 3），Wit值并不集中. 以Wit的最小值2.10来确定泥石流事件，灾害频率过高，以Wit的最大值5.49来确定泥石流事件，灾害频率过低，说明表 2的分类阈值和公式2的计算方法并不合理. 树木年轮在受到扰动后，首先生长抑制，然后紧接着出现生长释放，而且在冰川活跃区树木年轮生长缓慢，生长释放具有明显延迟性，这是导致利用Wit值重建泥石流灾害时间出现偏差的可能原因（吕立群等，2022）.
图 3 基于原阈值计算的Wit值重建的泥石流暴发时间Fig. 3. Reconstruction of debris flow occurrence by weighted index Wit不同树龄的年轮生长抑制/释放对泥石流扰动的响应不同（图 4）：树龄越小，生长抑制/释放对泥石流事件的扰动越敏感. 1958、1994、2002、2004、2007、2009、2012年泥石流事件中，生长抑制/释放集中在50 a树龄以下的树木中. 年轻树木树皮薄并且根系不稳固，容易受到泥石流等地质灾害的伤害，树木年龄和树木记录灾害事件的能力二者呈负相关，这与西藏迫龙沟利用树木年代学重建泥石流灾害事件中的规律基本一致（吕立群等，2022）. 树龄较小的树木在研究区居多，因为掩埋或者根系暴露可以存活，对大石块造成的撞击，后期很难愈合和存活（图 4）. 因此研究区扰动树木年轮生长抑制/释放比创伤及愈伤组织对Wit值的影响要大，导致利用Wit值重建泥石流灾害时间与实际暴发时间并不吻合.
图 4 不同树龄树木对泥石流扰动的响应S. 强烈扰动；M. 中等扰动；W. 弱扰动；I. 创伤Fig. 4. Tree sensitivity analysis to debris flows4. 基于新的Wit计算方法的泥石流灾害历史重建
4.1 新Wit计算方法的提出
由于研究区扰动树木年龄偏小，生长释放现象普遍且具有明显延迟性，为提高灾害重建的准确度，在公式（2）Wit计算中将生长释放去除；冰川活跃区树木生长缓慢，公式（1）中生长抑制的评判阈值修订为年轮宽度减少25%. 经过重新计算后，共得到967个扰动值，最早的生长抑制出现在1896年（表 4）. 表 2中的生长扰动强度分类阈值修订为表 5.
表 4 生长扰动年代分布表(GDs)Table Supplementary Table Distribution of growth disturbance (GDs)
生长抑制 创伤及
愈伤组织
生长抑制 创伤及
愈伤组织
生长抑制 创伤及
愈伤组织数量 占比 数量 占比 数量 占比 数量 占比 数量 占比 数量 占比 1896 1 1.00 0 0 1938 2 0.25 0 0.00 1980 8 0.17 0 0.00 1897 0 0.00 0 0 1939 2 0.25 0 0.00 1981 11 0.23 1 0.02 1898 0 0.00 0 0 1940 0 0.00 0 0.00 1982 12 0.25 1 0.02 1899 0 0.00 0 0 1941 4 0.50 0 0.00 1983 13 0.27 2 0.04 1900 0 0.00 0 0 1942 1 0.13 0 0.00 1984 4 0.08 3 0.06 1901 1 1.00 0 0 1943 2 0.22 0 0.00 1985 13 0.25 2 0.04 1902 0 0.00 0 0 1944 2 0.22 0 0.00 1986 10 0.19 2 0.04 1903 1 0.33 0 0 1945 4 0.44 0 0.00 1987 9 0.17 2 0.04 1904 1 0.33 0 0 1946 5 0.56 0 0.00 1988 17 0.30 2 0.04 1905 2 0.67 0 0 1947 0 0.00 0 0.00 1989 10 0.18 1 0.02 1906 0 0.00 0 0 1948 3 0.27 0 0.00 1990 16 0.29 1 0.02 1907 1 0.33 0 0 1949 2 0.18 0 0.00 1991 13 0.23 2 0.04 1908 0 0.00 0 0 1950 4 0.36 0 0.00 1992 11 0.19 2 0.03 1909 1 0.33 0 0 1951 1 0.09 0 0.00 1993 14 0.24 3 0.05 1910 1 0.33 0 0 1952 3 0.25 0 0.00 1994 18 0.31 4 0.07 1911 1 0.33 0 0 1953 3 0.25 0 0.00 1995 8 0.13 3 0.05 1912 1 0.33 0 0 1954 3 0.25 1 0.08 1996 22 0.35 3 0.05 1913 1 0.33 0 0 1955 5 0.31 1 0.06 1997 11 0.17 3 0.05 1914 0 0.00 0 0 1956 2 0.13 0 0.00 1998 18 0.29 2 0.03 1915 1 0.33 0 0 1957 7 0.44 0 0.00 1999 18 0.28 3 0.05 1916 1 0.33 0 0 1958 4 0.24 1 0.06 2000 14 0.22 3 0.05 1917 1 0.33 0 0 1959 4 0.21 1 0.05 2001 20 0.30 4 0.06 1918 0 0.00 0 0 1960 4 0.20 0 0.00 2002 11 0.16 4 0.06 1919 1 0.33 0 0 1961 0 0.00 0 0.00 2003 17 0.25 2 0.03 1920 1 0.33 0 0 1962 7 0.33 1 0.05 2004 18 0.26 4 0.06 1921 0 0.00 0 0 1963 3 0.14 1 0.05 2005 12 0.17 2 0.03 1922 1 0.25 0 0 1964 3 0.13 2 0.09 2006 19 0.27 3 0.04 1923 1 0.25 0 0 1965 7 0.27 0 0.00 2007 20 0.29 4 0.06 1924 1 0.25 0 0 1966 6 0.23 1 0.04 2008 15 0.21 4 0.06 1925 2 0.50 0 0 1967 2 0.07 1 0.04 2009 18 0.25 4 0.06 1926 0 0.00 0 0 1968 8 0.28 2 0.07 2010 25 0.35 3 0.04 1927 1 0.25 0 0 1969 4 0.13 2 0.07 2011 10 0.14 5 0.07 1928 2 0.40 0 0 1970 7 0.22 0 0.00 2012 15 0.21 6 0.08 1929 2 0.40 0 0 1971 8 0.24 1 0.03 2013 17 0.24 4 0.06 1930 0 0.00 0 0 1972 8 0.24 0 0.00 2014 16 0.22 3 0.04 1931 2 0.40 0 0 1973 9 0.26 0 0.00 2015 9 0.13 1 0.01 1932 0 0.00 0 0 1974 10 0.29 0 0.00 2016 24 0.33 2 0.03 1933 1 0.17 0 0 1975 10 0.26 0 0.00 2017 12 0.17 1 0.01 1934 1 0.14 0 0 1976 5 0.12 0 0.00 2018 15 0.21 1 0.01 1935 1 0.14 0 0 1977 11 0.26 0 0.00 2019 23 0.32 1 0.01 1936 2 0.29 0 0 1978 7 0.17 0 0.00 2020 17 0.24 2 0.03 1937 1 0.14 0 0 1979 11 0.24 0 0.00 2021 14 0.19 2 0.03 表 5 生长扰动强度的新分类阈值Table Supplementary Table New classification of growth disturbances (GDs)Ts Tm Tw Ti 年轮变化(%) 持续时间 年轮变化(%) 持续时间 年轮变化(%) 持续时间 生长抑制 ≥60% ≥5年 ≥60% < 5年且≥2年 ≥25%且 < 40% ≤2年 创伤及愈伤组织 存在创伤及愈伤组织 4.2 基于新Wit计算方法的泥石流暴发时间重建
利用表 5新的分类阈值得到1890-2021年的Wit值，经过数次计算发现Wit的定灾阈值与年轮的样本数量有较大关系：年轮样本数量At < 30时，Wit > 2大概率确定泥石流灾害时间；30≤At < 60时，Wit > 3大概率确定泥石流灾害时间；年轮样本数量At≥60时，Wit > 4大概率确定泥石流灾害时间（图 5）. 样本数量越大，排除极端气候、人类活动、动物的破坏对树木年轮的扰动的可能性越大，因此Wit值的越高.
图 5 基于新的Wit值计算方法的泥石流暴发时间重建Fig. 5. Reconstruction of debris flow occurrence by weighted index Wit根据新的生长抑制阈值（表 5）和新Wit计算方法重建的泥石流灾害年份（图 5）：1958年、1994年、2002年、2004年、2007年、2009年、2010年、2012年，与实际泥石流时间（2002、2007、2009、2012）吻合度高. 2010年Wit值较高，造成泥石流发生的假象，原因是：2009年泥石流规模巨大，造成大量树木生长抑制，年轮恢复正常生长普遍滞后一年，计算所得2010年Wit值较高.
按照新方法重建堆积区3条流动路径的灾害时间，ABCD段灾害年份：1958年、1968年、2002年、2004年、2009年；ABCE段灾害年份：2002年、2004年、2007年、2009年、2010年、2012年；ABCF段灾害年份：1958年、1994年、2009年、2010年（图 6）. 将3条路径灾害年份与全区域灾害年份对比，只有1968年灾害未被记录到，证明了新方法的合理性.
图 6 基于新的Wit值计算方法的泥石流暴发时间重建Fig. 6. Reconstruction of debris flow burst time based on different paths based on index Wit4.3 基于新Wit计算方法的泥石流流动范围重建
根据生长扰动树木的位置重建灾害年份泥石流活动范围（图 7）：可以看出1958-1994年期间堆积区CE段地势比CD段和CF段高，泥石流未对CE段树木造成明显扰动；2002年和2004年泥石流造成CD段和CE段生长扰动明显，同时这两年灾害时间较近，可能导致CD段堆积物叠加并堵江；2009年的泥石流覆盖整个堆积区，是规模最大的一年，这与泥石流的实际淤积范围基本一致.
5. 结论
研究区冰川运动活跃，树木年轮生长缓慢，尤其在受到泥石流的扰动之后，生长释放具有明显的延迟性，导致基于原方法计算的Wit指数与实际泥石流暴发时间不一致. 各个年龄段树木年轮对冰川型泥石流的扰动响应不同，生长释放在年轻树木中更加常见，树龄越小，生长释放对泥石流事件的扰动越敏感. 研究区年轻的扰动树木占据多数，创伤及愈伤组织并不常见，进一步增加了生长释放在Wit指数计算过程中的权重. 论文基于大量的年轮生长变化分析，将新的生长抑制评判阈值由原来的40%降低到25%. 利用Wit重建冰川地区泥石流灾害过程中，如果年轻树木较多，生长释放会大大增加Wit的误差，本文在重构新的Wit指数计算方法中排除了生长释放这一因子. 基于树木年轮创伤及愈伤组织，生长抑制因子构建的新的Wit指数，反映出研究区在1958、1994、2002、2004、2007、2009、2012年暴发过泥石流，与历史记录基本吻合，反映出本文提出的新Wit指数计算方法重建冰川型泥石流的合理性.
-
图 1 研究区及树木采集样点
左上图底图审图号：GS（2020）4619号
Fig. 1. Study area and distribution of tree samples
图 3 基于原阈值计算的Wit值重建的泥石流暴发时间
Fig. 3. Reconstruction of debris flow occurrence by weighted index Wit
图 4 不同树龄树木对泥石流扰动的响应
S. 强烈扰动；M. 中等扰动；W. 弱扰动；I. 创伤
Fig. 4. Tree sensitivity analysis to debris flows
图 5 基于新的Wit值计算方法的泥石流暴发时间重建
Fig. 5. Reconstruction of debris flow occurrence by weighted index Wit
图 6 基于新的Wit值计算方法的泥石流暴发时间重建
Fig. 6. Reconstruction of debris flow burst time based on different paths based on index Wit
表 1 基于树木年轮的灾害定年方法
Table 1. Standard for tree ring disaster dating
作者 国家 地区 样品数量 时期(年份) 海拔(m) 年均降水(mm) GD总数 定灾标准 灾害类型 It GD数量 Wit Mayera et al.（2010） 奥地利 Gratzental 227 1800-2007 2 106~1 166 1 526 1 155 It≥4% 未提及 未提及 泥石流 Corona et al.（2012） 法国 Arve 209 1771-1998 1 100~3 650 1 262 645 It > 1.7 GD数 > 5% 未提及 雪崩 Peitzsch et al.(2013) 美国 Rocky 673 1636-2017 1 100~2 700 未提及 2 134 样品数=10~20：
GD数≥3个，It≥15%高：Wit≥0.3 雪崩 样品数=21~50：
GD数≥5个，It≥10%中：Wit≥0.2 样品数=51~100：
GD数≥7个，It≥7%样品数 > 100：
GD数≥9个，It≥4.5%低：Wit < 0.2 Vădean et al.（2015） 罗马尼亚 Apuseni 20 1940-2010 平均835 800~1 000 72 It≥15% GD数≥3个 未提及 泥石流 Šilhán et al.（2021） 捷克共和国 Hostynsko-vsetınskahornatina Mts. 153 1920-2017 平均550 700~1 200 179 可能发生：5%~10%
一定发生：It≥10%未提及 未提及 滑坡 Franco-Ramos et al.（2019） 墨西哥 La Malinche 65 1933-2017 3 000~4 400 754 205 未提及 确定：GD≥4和Wit≥2
“潜在”泥石流：Wit≥1.5但GD≤4
噪声：GD数≤2个，Wit≤2泥石流 Schneuwly-Bollschweilera et al.（2013） 瑞士 Zermatt 385 1600-2008 1 420~4 545 690 581 未提及 未提及 Wit > 1 泥石流 注：GD. 生长扰动；It. 响应指数；Wit. 加权指数 
下载: 导出CSV
表 2 生长扰动强度的分类阈值
Table 2. Classification of growth disturbances (GDs)
Ts Tm Tw Ti 年轮变化(%) 持续时间 年轮变化(%) 持续时间 年轮变化(%) 持续时间 生长抑制 ≥60% ≥5年 ≥60% < 5年且≥2年 ≥40%且 < 60% ≥2年 生长释放 ≥60% ≥5年 ≥60% < 5年且≥2年 ≥50%且 < 60% ≥2年 创伤及愈伤组织 存在创伤及愈伤组织
下载: 导出CSV
表 3 遥感数据集
Table 3. Collected Satellite datasets in this work
编号 数据来源 日期(日/月/年) 分辨率(m) 编号 数据来源 日期(日/月/年) 分辨率(m) 1 Landsat5 17/01/1987 30 15 Landsat5 17/02/2004 30 2 Landsat5 05/02/1988 30 16 Landsat5 23/12/2006 30 3 Landsat5 08/12/1989 30 17 Landsat5 30/04/2007 30 4 Landsat5 11/12/1990 30 18 Landsat5 26/11/2008 30 5 Landsat5 14/12/1991 30 19 Landsat5 29/11/2009 30 6 Landsat5 30/11/1992 30 20 Landsat5 18/12/2010 30 7 Landsat5 01/11/1993 30 21 Landsat5 04/02/2011 30 8 Landsat5 06/12/1994 30 22 Landsat8 24/11/2013 30 9 Landsat5 07/11/1995 30 23 Landsat8 29/12/2014 30 10 Landsat5 07/11/1996 30 24 Landsat8 30/11/2015 30 11 Landsat5 15/11/1998 30 25 Landsat8 21/12/2017 30 12 Landsat5 18/01/1999 30 26 Landsat8 22/11/2018 30 13 Landsat5 22/12/2000 30 27 Landsat8 11/11/2020 30 14 Landsat5 15/12/2003 30 28 Landsat8 06/03/2022 30
下载: 导出CSV
表 4 生长扰动年代分布表(GDs)
Table 4. Distribution of growth disturbance (GDs)
生长抑制 创伤及
愈伤组织
生长抑制 创伤及
愈伤组织
生长抑制 创伤及
愈伤组织数量 占比 数量 占比 数量 占比 数量 占比 数量 占比 数量 占比 1896 1 1.00 0 0 1938 2 0.25 0 0.00 1980 8 0.17 0 0.00 1897 0 0.00 0 0 1939 2 0.25 0 0.00 1981 11 0.23 1 0.02 1898 0 0.00 0 0 1940 0 0.00 0 0.00 1982 12 0.25 1 0.02 1899 0 0.00 0 0 1941 4 0.50 0 0.00 1983 13 0.27 2 0.04 1900 0 0.00 0 0 1942 1 0.13 0 0.00 1984 4 0.08 3 0.06 1901 1 1.00 0 0 1943 2 0.22 0 0.00 1985 13 0.25 2 0.04 1902 0 0.00 0 0 1944 2 0.22 0 0.00 1986 10 0.19 2 0.04 1903 1 0.33 0 0 1945 4 0.44 0 0.00 1987 9 0.17 2 0.04 1904 1 0.33 0 0 1946 5 0.56 0 0.00 1988 17 0.30 2 0.04 1905 2 0.67 0 0 1947 0 0.00 0 0.00 1989 10 0.18 1 0.02 1906 0 0.00 0 0 1948 3 0.27 0 0.00 1990 16 0.29 1 0.02 1907 1 0.33 0 0 1949 2 0.18 0 0.00 1991 13 0.23 2 0.04 1908 0 0.00 0 0 1950 4 0.36 0 0.00 1992 11 0.19 2 0.03 1909 1 0.33 0 0 1951 1 0.09 0 0.00 1993 14 0.24 3 0.05 1910 1 0.33 0 0 1952 3 0.25 0 0.00 1994 18 0.31 4 0.07 1911 1 0.33 0 0 1953 3 0.25 0 0.00 1995 8 0.13 3 0.05 1912 1 0.33 0 0 1954 3 0.25 1 0.08 1996 22 0.35 3 0.05 1913 1 0.33 0 0 1955 5 0.31 1 0.06 1997 11 0.17 3 0.05 1914 0 0.00 0 0 1956 2 0.13 0 0.00 1998 18 0.29 2 0.03 1915 1 0.33 0 0 1957 7 0.44 0 0.00 1999 18 0.28 3 0.05 1916 1 0.33 0 0 1958 4 0.24 1 0.06 2000 14 0.22 3 0.05 1917 1 0.33 0 0 1959 4 0.21 1 0.05 2001 20 0.30 4 0.06 1918 0 0.00 0 0 1960 4 0.20 0 0.00 2002 11 0.16 4 0.06 1919 1 0.33 0 0 1961 0 0.00 0 0.00 2003 17 0.25 2 0.03 1920 1 0.33 0 0 1962 7 0.33 1 0.05 2004 18 0.26 4 0.06 1921 0 0.00 0 0 1963 3 0.14 1 0.05 2005 12 0.17 2 0.03 1922 1 0.25 0 0 1964 3 0.13 2 0.09 2006 19 0.27 3 0.04 1923 1 0.25 0 0 1965 7 0.27 0 0.00 2007 20 0.29 4 0.06 1924 1 0.25 0 0 1966 6 0.23 1 0.04 2008 15 0.21 4 0.06 1925 2 0.50 0 0 1967 2 0.07 1 0.04 2009 18 0.25 4 0.06 1926 0 0.00 0 0 1968 8 0.28 2 0.07 2010 25 0.35 3 0.04 1927 1 0.25 0 0 1969 4 0.13 2 0.07 2011 10 0.14 5 0.07 1928 2 0.40 0 0 1970 7 0.22 0 0.00 2012 15 0.21 6 0.08 1929 2 0.40 0 0 1971 8 0.24 1 0.03 2013 17 0.24 4 0.06 1930 0 0.00 0 0 1972 8 0.24 0 0.00 2014 16 0.22 3 0.04 1931 2 0.40 0 0 1973 9 0.26 0 0.00 2015 9 0.13 1 0.01 1932 0 0.00 0 0 1974 10 0.29 0 0.00 2016 24 0.33 2 0.03 1933 1 0.17 0 0 1975 10 0.26 0 0.00 2017 12 0.17 1 0.01 1934 1 0.14 0 0 1976 5 0.12 0 0.00 2018 15 0.21 1 0.01 1935 1 0.14 0 0 1977 11 0.26 0 0.00 2019 23 0.32 1 0.01 1936 2 0.29 0 0 1978 7 0.17 0 0.00 2020 17 0.24 2 0.03 1937 1 0.14 0 0 1979 11 0.24 0 0.00 2021 14 0.19 2 0.03
下载: 导出CSV
表 5 生长扰动强度的新分类阈值
Table 5. New classification of growth disturbances (GDs)
Ts Tm Tw Ti 年轮变化(%) 持续时间 年轮变化(%) 持续时间 年轮变化(%) 持续时间 生长抑制 ≥60% ≥5年 ≥60% < 5年且≥2年 ≥25%且 < 40% ≤2年 创伤及愈伤组织 存在创伤及愈伤组织
下载: 导出CSV
-
Corona, C., Lopez Saez, J., Stoffel, M., et al., 2012. How much of the Real Avalanche Activity Can Be Captured with Tree Rings? An Evaluation of Classic Dendrogeomorphic Approaches and Comparison with Historical Archives. Cold Regions Science and Technology, 74: 31-42. https://doi.org/10.1016/j.coldregions.2012.01.003 Ding, M., Bai, S. B., Wang, J., et al., 2016. Basic Procedures of Using Tree Rings to Reconstruct the Ttime of Landslide Reactivation. Mountain Research, 34(5): 545-554(in Chinese with English abstract). Franco-Ramos, O., Stoffel, M., Ballesteros-Cánovas, J. A., 2019. Reconstruction of Debris-Flow Activity in a Temperate Mountain Forest Catchment of Central Mexico. Journal of Mountain Science, 16(9): 2096-2109. https://doi.org/10.1007/s11629-019-5496-6 Huang, T., 2019. Study on the Response of Debris Flow Activity in the Upper Reaches of Minjiang River to Vertical Climate(Differentiation). Southwest University of Science and Technology, Mianyang (in Chinese with English abstract). Kogelnig-Mayer, B., Stoffel, M., Schneuwly-Bollschweiler, M., et al., 2011. Possibilities and Limitations of Dendrogeomorphic Time-Series Reconstructions on Sites Influenced by Debris Flows and Frequent Snow Avalanche Activity. Arctic, Antarctic, and Alpine Research, 43(4): 649-658. https://doi.org/10.1657/1938-4246-43.4.649 Lai, Z. P., Yang, A. N., Cong, L., et al., 2021. A Review on the Dating Techniques for Mountain Hazards-Induced Sediments. Earth Science Frontiers, 28(2): 1-18(in Chinese with English abstract). Li, Y., Cui, Y. F., Li, Z. H., et al., 2022. Evolution of Glacier Debris Flow and lts Monitoring System along Sichuan-Tibet Traffic Corridor. Earth Science, 47(6): 1969-1984(in Chinese with English abstract). Lundström, T., Stoffel, M., Stöckli, V., 2008. Fresh-Stem Bending of Silver Fir and Norway Spruce. Tree Physiology, 28(3): 355-366. https://doi.org/10.1093/treephys/28.3.355 Lyu, L. Q., 2017. Research on the Initiation and Motion of Gully Debris Flows in Tibetan Plateau(Dissertation). Tsinghua University, Beijing, 3-5(in Chinese with English abstract). Lyu, L. Q., Wang, Z. Y., Meng, Z., 2022. Reconstruction of Debris Flow Disasters in Polong Gully Based on Dendrochronology. Earth Science, 49(1): 335-346(in Chinese with English abstract). Mayer, B., Stoffel, M., Bollschweiler, M., et al., 2010. Frequency and Spread of Debris Floods on Fans: a Dendrogeomorphic Case Study from a Dolomite Catchment in the Austrian Alps. Geomorphology, 118(1/2): 199-206. https://doi.org/10.1016/j.geomorph.2009.12.019 Meng, Z., Lyu, L. Q., Yu, G. A., et al., 2022. Reconstruction of Glacial Debris Flow Disaster Based on Dendrochronology: A Case Study on Tianmo Gully, Tibet. Science Technology and Engineering, 22(32): 14124-14136(in Chinese with English abstract). Miao, X. Q., 2022. Research on the Movement Characteristics and Parameters of Debris Flow in Rangdaqu. Journal of Railway Engineering Society, 39(4): 20-25(in Chinese with English abstract). Peitzsch, E., Hendrikx, J., Stahle, D., et al., 2013. A Regional Spatiotemporal Analysis of Large Magnitude Snow Avalanches Using Tree Rings. Natural Hazards and Earth System Sciences, 21: 533-557. https://doi.org/10.5194/NHESS-21-533-2021 Schneuwly, D. M., Stoffel, M., Dorren, L. K. A., et al., 2009. Three-Dimensional Analysis of the Anatomical Growth Response of European Conifers to Mechanical Disturbance. Tree Physiology, 29(10): 1247-1257. https://doi.org/10.1093/treephys/tpp056 Schneuwly-Bollschweiler, M., Corona, C., Stoffel, M., 2013. How to Improve Dating Quality and Reduce Noise in Tree-Ring Based Debris-Flow Reconstructions. Quaternary Geochronology, 18: 110-118. https://doi.org/10.1016/j.quageo.2013.05.001 Šilhán, K., 2021. A New Tree-Ring-Based Index for the Expression of Spatial Landslide Activity and the Assessment of Landslide Hazards. Geomatics, Natural Hazards and Risk, 12(1): 3409-3428. https://doi.org/10.1080/19475705.2021.2011790 Stoffel, M., Bollschweiler, M., Hassler, G., 2006. Differentiating Past Events on a Cone Influenced by Debris‐Flow and Snow Avalanche Activity: a Dendrogeomorphological Approach. Earth Surface Processes and Landforms, 31(11): 1424-1437. Strunk, H., 1997. Dating of Geomorphological Processes Using Dendrogeomorphological Methods. CATENA, 31(1/2): 137-151. https://doi.org/10.1016/S0341-8162(97)00031-3 Tie, Y. B., Malik, I., Owczarek, P., 2014. Dendrochronological Dating of Debris Flow Historical Events in High Mountain Area: Take Daozao Debris Flow as an Example. Mountain Research, 32(2): 226-232(in Chinese with English abstract). Vădean, R., Arghiuş, V., Pop, O., 2015. Dendrogeomorphic Reconstruction of Past Debris-Flood Activity along a Torrential Channel: an Example from Negoiul Basin (Apuseni Mountains, Romanian Carpathians). Zeitschrift Fur Geomorphologie, 59(3): 319-335. https://doi.org/10.1127/zfg/2014/0156 Wang, K., 2016. Larixgmelini Tree-Ring Width Index with Responses to Climate Change In the West of Sichuan Subalpine Zone(Dissertation). Sichuan Agricultural University, Chengdu(in Chinese with English abstract). Wang, Z. L., Ma, C., Wu, J. L., et al., 2022. Debris Flow Event in Xiaoxitian Watershed of Miyun Based on Tree Ringreconstruction. Journal of Natural Disasters, 31(5): 183-192(in Chinese with English abstract). Wu, J. L., Ma, C., Wang, R., et al., 2021. Reconstruction of Torrent and Debris Flow Events Based Ondendro Geomorphology: A Case Study of Longtangou Basin in Miyun District, Beijing. Journal of Natural Disasters, 30(1): 183-190(in Chinese with English abstract). Zeng, X. Y., Zhang, J. J., Yang, D. X., et al., 2019. Characteristics and Geneses of Low Frequency Debris Flow along Parlongzangbo River Zone: Take Chaobulongba Gully as an Example. Science Technology and Engineering, (34): 103-107(in Chinese with English abstract). Zhang, J. S., Xie, H., Wang, X. D., et al., 2015. Debris-Flow of Jianmupuqu Ravine in Tibet. Journal of Catastrophology, 30(3): 99-103(in Chinese with English abstract). 曾宪阳, 张佳佳, 杨东旭, 等, 2019. 帕隆藏布流域低频泥石流的成因机制分析——以倾多镇抄布隆巴泥石流为例. 科学技术与工程, 19(34): 103-107. 丁苗, 白世彪, 王建, 等, 2016. 利用树木年轮重建滑坡复活时间的方法. 山地学报, 34(5): 545-554. 黄涛, 2019. 岷江上游泥石流活动对气候垂直分异的响应研究(硕士学位论文). 绵阳: 西南科技大学. 赖忠平, 杨安娜, 丛禄, 等, 2021. 山地灾害沉积物的测年综述. 地学前缘, 28(2): 1-18. 李尧, 崔一飞, 李振洪, 等, 2022. 川藏交通廊道林波段冰川泥石流发育动态演化分析及监测预警方案. 地球科学, 47(6): 1969-1984. doi: 10.3799/dqkx.2021.194 吕立群, 2017. 青藏高原泥石流的形成运动过程研究(博士学位论文). 北京: 清华大学, 3-5. 吕立群, 王兆印, 孟哲, 2024. 基于树木年代学的迫龙沟泥石流灾害历史重建. 地球科学, 49(1): 335-346. doi: 10.3799/dqkx.2022.142 孟哲, 吕立群, 余国安, 等, 2022. 基于树木年代学的冰川型泥石流灾害历史重建方法——以西藏天摩沟为例. 科学技术与工程, 22(32): 14124-14136. 苗晓岐, 2022. 瓤打曲泥石流运动特征及参数研究. 铁道工程学报, 39(4): 20-25. 铁永波, Ireneusz Malik, Piotr Owczarek, 2014. 树木年代学在高寒山区泥石流历史事件重建中的应用——以磨西河流域倒灶沟为例. 山地学报, 32(2): 226-232. 王珂, 2016. 川西亚高山落叶松径向生长与气候变化的关系(硕士学位论文). 成都: 四川农业大学. 王志兰, 马超, 吴佳亮, 等, 2022. 基于树木年轮重建密云小西天流域泥石流事件. 自然灾害学报, 31(5): 183-192. 吴佳亮, 马超, 王锐, 等, 2021. 基于树木地貌法重建山洪泥石流历史事件——以北京密云区龙潭沟流域为例. 自然灾害学报, 30(1): 183-190. 张金山, 谢洪, 王小丹, 等, 2015. 西藏尖姆普曲泥石流. 灾害学, 30(3): 99-103. 期刊类型引用(1)
1. 徐文君，马超，吕立群，杜翠. 西藏热玛沟流域冰川型泥石流历史的树木年代学重建. 山地学报. 2025(01): 122-131 . 
百度学术其他类型引用(0)
-
下载:
百度学术
点击查看大图
计量
- 文章访问数: 195
- HTML全文浏览量: 94
- PDF下载量: 17
- 被引次数: 1
