Application of 3D Laser Scanning and Quantitative Morphology Analysis Method to Bedrock Fault Surface in Paleo-seismic Research in Bedrock Area
-
摘要: 由于对第四纪地层的严重依赖,传统古地震探槽研究方法在基岩区难以发挥作用,导致无法获取基岩区断层的强震活动历史。本研究以山西地堑系的交城断裂为目标断裂,以断裂北段2处基岩断层面为研究对象,通过三维激光扫描技术获取基岩断层面高精度形貌,基于变差函数法结合滑动窗口操作量化断层表面形貌特征,开展在基岩区提取断裂古地震信息的实例研究。结果显示,2处基岩断层面的形貌在高度上具有明显的分段特征,指示了断层面在地震事件作用下的分段出露过程。这种断层面形貌分段特征可以用来识别古地震事件和同震位移量。在思西村基岩断层面上,识别出由老到新的3次古地震事件,同震倾滑位移量依次为2.0 m、1.9 m和2.3 m,在上兰镇基岩断层面上,识别出由老到新的3次古地震事件,同震倾滑位移量依次为1.4 m、2.5 m和2.0 m,指示了交城断裂北、中段具有产生同震位移量大于2 m、震级大于7.5级的破裂型地震的能力。上述研究成果表明,基于三维激光扫描和形貌量化分析方法开展基岩断层面古地震研究,可以准确而高效地识别古地震事件次数和同震位移量,扩展古地震的研究对象,拓宽古地震的研究空间。在未来的研究中,可以适时地开展宇宙成因核素测年以测定断层面的暴露年龄,获得发震年代,给予地震序列年龄框架。Abstract: Due to the heavy reliance on Quaternary strata, traditional trench research methods are difficult to apply in the bedrock area, resulting in the inability to obtain strong seismic activity history on faults in the bedrock area. In this study, the Jiaocheng Fault in the Shanxi Rift is selected as the target fault, and the two bedrock fault surfaces in the northern segment are taken as the study object to carry out case studies. The variogram combined with the sliding window operation quantifies the morphologic features of the fault surface and extracts the paleo-seismic information. The results show that the morphology of the two bedrock fault surfaces has segmentation in height, indicating the exposure process of the fault surfaces under the action of periodic seismic events. Therefore, paleo-earthquake events and co-seismic slips can be identified by the morphological segments. Accordingly, three paleo-earthquake events were identified on the bedrock fault of Sixi village, and the co-seismic dip slips were 2.0 m, 1.9 m and 2.3 m in sequence. The co-seismic dip slips of the three paleoearthquakes are 1.4m, 2.5 m and 2.0 m in sequence. The results indicate that the northern segment of the northern and middle segments of Jiaocheng fault can produce ruptures with a co-seismic slip greater than 2 m and a magnitude greater than M7.5. The above results also show that the paleo-seismic research on bedrock fault surfaces based on t-LiDAR and morphological quantitative analysis can accurately and efficiently identify the number of paleo-earthquake and co-seismic slips, increase the objects of paleo-earthquakes, and expand research space of paleo-earthquakes. In future research, cosmogenic nuclide dating can be carried out timely to determine the exposure age of the fault surface, obtain the earthquake ages, and provide an age frame to the earthquake sequences.1) 2为区分变量,此处用D0表示同震位移,与刘静等(1996)中同震位移D含义一致。
-
引言
地质学家已经确定全球范围的全新世火山约有1500多座(许建东,2011)。火山类型主要为成层火山、复合式火山、熔岩锥、破火山口、盾状火山等(刘若新等,1999)。此外,还有众多火山位于深海海底,它们的喷发很少能到达水面,因此少有观察和记录。大部分火山的形成是板块运动的结果,是地球动力学过程的重要现象(Press等,1982)。全球现今活动构造划分为环太平洋、大洋中脊和大陆三大构造系统(马宗晋等,2003)。全球火山活动分西太平洋火山活动区、东太平洋火山活动区以及大西洋火山活动区(洪汉净等, 2003, 2009)。位于太平洋板块西南缘的印度尼西亚、菲律宾和日本等国家,拥有的火山数量为全球正在活动火山数量的1/3(Siebert等,2010),2017年该区域活动强度最大的火山为印度尼西亚巴厘岛的阿贡火山(王佳龙等,2018a)。
本文数据主要来自全球火山计划网站2,该网站由华盛顿国家自然历史博物馆矿物科学部门(Department of Mineral Sciences,National Museum of Natural History,Washington D.C)的史密森全球火山项目(Smithsonian Institution Global Volcanism Program)支撑,汇集了世界各地活动火山监测机构的监测信息及相应的火山监测网站链接,并于每周四发布监测周报。文中的活动火山指2016年12月30日—2017年12月31日发生活动并被该网站记录下来的火山。
2 http://volcano.si.edu/reports_weekly.cfm
通过整理该网站2017年发布的约858条全球火山监测活动信息,总结2017年度全球火山活动的时空关系,并针对阿贡火山的喷发过程进行了追踪,为监测火山喷发、研究火山活动机理、预测火山灾害等提供参考。
1. 2017年全球活动火山空间分布特征
火山警戒等级是衡量火山活动危险程度的标识,USGS火山灾害项目1将火山预警等级划分为正常(Ⅰ级)、咨询(Ⅱ级)、注意(Ⅲ级)和警戒(Ⅳ级)4个等级,其代表的危险程度依次上升,火山活动越强烈,引发的火山灾害也越严重。具体划分标准为:Ⅰ级,火山警戒级别为正常,航空颜色代码为绿色,火山有活动迹象,但活动处于正常范围内,主要表现为火山喷气及火山灰柱飘散;Ⅱ级,火山警戒级别为低度警戒,航空颜色代码为黄色,火山活动信号上升,高于正常范围,主要表现为火山出现热异常,火山喷气柱升高,出现火山爆炸声,火山由正常活动向危险活动过渡,该类型火山一般不稳定;Ⅲ级,火山警戒级别为中度警戒,航空颜色代码为橙色,火山活动信号进一步上升,具有喷发的前兆,主要表现为火山地震频率和强度增大,可以监测到熔岩流溢出,热异常和火山灰柱活动都有所加强;Ⅳ级,火山警戒级别为高度警戒,航空颜色代码为红色,表现为火山即将喷发或正在喷发。一般情况下,火山的警戒等级随火山活动的强弱而改变,可能会在数月内出现连续变化(张传杰等,2016)。
1 https://volcanoes.usgs.gov/vhp/about_alerts.html
全球活动火山约有80%分布在环太平洋地区(Siebert等,2010)。2017年度的活动火山约有90%分布在环太平洋地区,其中,太平洋北部阿留申群岛、堪察加半岛和北方四岛等分布12座活动火山,其警戒等级为Ⅱ级和Ⅲ级,警戒等级Ⅲ级的活动火山相对比较活跃;位于太平洋板块西缘的日本岛弧和东部海岛分布11座活动火山,警戒等级为Ⅰ级的火山相对比较活跃;位于澳大利亚板块北缘的印度尼西亚等一系列岛链,由澳大利亚板块向北俯冲到缅甸板块、欧亚板块和太平洋板块之下形成,构造复杂,火山活动能量大,2017年此区域有21座活动火山,包括3座警戒等级Ⅲ级的活动火山,可见其活动频度和强度都极高。
全球火山分布图(图 1)的底图来源于USGS,对其做罗宾逊投影变换,中央子午线设为155.0°E。若某火山2017年每周均有持续活动,则将该火山的活动频率设为1;若火山全年无任何活动,则将火山的活动频率设为0。火山的活动频率一般介于0和1之间,它反映了火山的活跃程度。
2. 2017年全球活动火山时间分布特征
2017年全球范围内,平均每月约有30座火山喷发的记录,最多达34座,最少为26座(图 2)。年度内火山活动数量和活动强度有逐渐增多、增强的趋势,表明俯冲板块活动的加剧。可将活动的火山分为3种类型(图 3):①“冒泡”型,一般活动强度较弱,持续时间较短,警戒等级一般为Ⅰ级或Ⅱ级,活动形式主要表现为火山喷气、火山震颤、火山地震等,通常不会造成地质灾害和影响周围居民的生活,典型代表有冰岛的卡特拉火山(Katla)、美国的希沙尔丁火山(Shishaldin)、印度尼西亚的赛梅鲁火山(Semeru)和俄罗斯的朱帕诺夫斯基火山(Zhupanovsky)等;②“持续”型,基本常年处于持续喷发状态,警戒等级一般为Ⅱ级或Ⅲ级,活动形式主要表现为火山喷气、火山地震、岩浆活动、火山碎屑流等,这类火山是全球火山监测的重点,对其进行的火山喷发监测和火山灾害预防都较为成熟,典型代表有美国夏威夷的基拉韦厄火山(Kilauea)(Fontijn等,2015;王佳龙,2018b)、印度尼西亚的锡纳朋火山(Sinabung)、秘鲁的萨班卡亚火山(Sabancaya)和俄罗斯勘察加半岛的希韦卢奇火山(Sheveluch),这几座火山几乎常年都在活动,火山区周围常引发地质灾害;③“剧烈”型,常年处于稳定状态,偶尔几周或几个月突然爆发,警戒等级一般为Ⅲ级或Ⅳ级,喷发的形式主要表现为喷发柱的突然上升,有时高达几千米至十几千米,随后火山碎屑物质喷出火口并造成坍塌,引发火山碎屑流,最后岩浆溢出形成熔岩流,典型代表有日本的新燃岳火山(Kirishimayama)、菲律宾的马荣火山(Mayon)以及印度尼西亚的阿贡火山(Agung)(王佳龙,2018a)。
3. 阿贡火山的喷发过程
全球的活火山主要分布在环太平洋地区,此区域也是著名的“火环”,而印度尼西亚岛链是西南太平洋最为活跃的1条“火链”(图 4,图中火山名称下方数字,如2003/46/2代表的含义为2003年第1次监测到该火山的活动,目前共计活动了46周,2017年度从2月开始变活跃),其居民数量为世界之最,阿贡火山则是该“火链”中最为活跃的火山之一(Self等1996)。阿贡火山的喷发周期约50年,最近1次喷发于1963年(Zen等1964;Marinelli等,1968),距今56年。据印尼国家灾害管理局1(Badan Nasional Penanggulangan Bencana)的监测数据,自2017年8月开始,阿贡火山的地震明显增多,有复苏和再次喷发的前兆,9月份火山警戒等级升至Ⅲ级,10月份开始有喷发柱喷出,11月喷发柱急剧上升至4km并开始有岩浆活动,警戒等级也升至最高级Ⅳ级。印度尼西亚岛链位于澳大利亚板块与欧亚板块之间,由于澳大利亚板块向北的俯冲,在此区域形成了1套完整的沟-弧-盆体系。5—6级地震的震中主要位于岛链与巽他海沟之间,与岛链大致平行,距离岛链约100—150km,地震震源深度约100km;活火山主要位于岛链内部,岩浆来源大致在深度150—200km板块的交汇处。由此可见,由澳大利亚板块俯冲而下的矿物在岛链下方150—200km处发生脱水作用,脱水诱导下覆地震的部分熔融,产生低盐度、钙碱性岩浆;岩浆上升侵入上覆板块的岩石圈中。
4. 阿贡火山喷发过程中的形变
在火山喷发过程中,前期地下岩浆活动,地表热异常,火山开始膨胀;当压力达到一定程度后,火山气体首先冲出,火山通道内部的压力减小,岩浆成分挥发,火山灰喷出,随后岩浆补给,岩浆爆炸喷出或溢出火口,后期岩浆冷却收缩,火口坍塌。目前,阿贡火山处在火山灰的喷发过程,未出现岩浆溢出。自2017年6月开始,阿贡火山有所膨胀,至2017年12月膨胀近15cm,随后开始逐渐收缩(图 5)。膨胀的中心位于阿贡火山北侧附近,并未在火山的正下方,说明岩浆通道向北倾斜。
采用“二通”法对收集到的SAR复数影像进行干涉差分处理。在处理过程中,使用SRTM4 DEM数据消除地形相位影响,同时,为保证较高的配准精度,采用结合轨道和地形数据的图像配准技术,配准精度优于0.001个像元。由于巴厘岛大部分地区为山地,且气候温和多雨,岛上植被茂盛,易造成干涉像对的失相干,同时SAR影像在山区易出现叠掩、透视收缩、阴影等,也将加剧失相干的程度。失相干区或低相干区将导致相位解缠误差,为此,使用迭代自适应滤波算法进行降噪处理,并采用基于狄洛尼三角剖分的最小费用流算法(MCF)。首先,对高质量的相位区进行解缠获得可靠的参考相位模型,再利用参考相位实现对低相干区域的解缠,从而得到全局的最优结果,最后,经过地理编码获取了阿贡火山喷发过程的InSAR形变场。
5. 阿贡火山喷发柱高度变化
阿贡火山自2017年9月开始活动以来,活动特征明显,浅源小震持续增多、热异常明显。2017年10月后,阿贡火山开始出现不同规模的爆炸,爆炸产生的火山灰柱上升至高空并随风向飘逸,爆炸的规模直接决定了喷发柱的高度。阿贡火山的整个喷发过程以及爆炸与喷发柱变化的关系,如图 6所示。由图可见,阿贡火山在2017年10月份主要表现为火山地震和震颤,是地下岩浆通过管道向上运移的过程,此时山体表面也开始不断膨胀,该过程一直持续到11月27日,当日发生1次剧烈爆炸,造成火山剧烈喷发,喷发柱升至海拔4km,大量火山碎屑流顺火山翼冲下,当地居民被迫迁移,航班被迫取消;在随后的3个月内,不断有地下岩浆补给,经常性地发生爆炸并产生喷发柱;在2018年2月27日之后,岩浆补给减弱,火山警戒等级下降,爆炸规模减小,喷发柱高度降低,火山表面逐渐收缩;自2017年9月末至2018年3月末,阿贡火山历经了整个喷发过程,从开始的地震活动,到喷发时大规模爆炸,再到最后阶段的平静收缩,此过程对研究其它“剧烈”型火山具有借鉴意义。
6. 结论
2017年全球火山喷发较往年强烈,年度内火山数量和强度均有增大的趋势(图 7)。其中“剧烈”型火山造成的火山灾害较为严重,应对其加强监测和防灾。活火山主要分布在环太平洋地区,位于太平洋西南的印度尼西亚是2017年受灾最为严重的地区,据印尼国家灾害管理局的统计,阿贡火山的喷发造成10万余人无家可归。阿贡火山历时半年多的活动过程从膨胀到收缩,从平静到爆炸再到平静,使我们对火山喷发的过程有了更全面的了解。
-
表 1 思西村基岩断层面的分形维数与分段高度
Table 1. Characteristic fractal value (D)and heights of morphological segments on Sixi village bedrock fault surface
断层面 滑动窗口尺寸 分形维数(D) 标准差(σ) 下段 中段 上段 下段 中段 上段 思西村
断层面66 mm×66 mm 2.4303 2.4118 2.3218 0.0506 0.0331 0.0485 130 mm×130 mm 2.4107 2.3876 2.2979 0.0500 0.0262 0.0458 258 mm×258 mm 2.3889 2.3771 2.2842 0.0441 0.0361 0.0339 分段高度H/m / 2.3 1.9 2.0 / / / 表 2 上兰镇基岩断层面的分形维数和分段高度
Table 2. Characteristic fractal value (D) and heights of morphological segments on Shanglanzhen bedrock fault surface
断层面 滑动窗口尺寸 分形维数(D) 标准差(σ) 下段 中段 上段 下段 中段 上段 上兰镇
断层面66 mm×66 mm 2.7105 2.5331 2.2361 0.0266 0.0815 0.0419 130 mm ×130 mm 2.6810 2.5204 2.2122 0.0441 0.0653 0.0458 258 mm×258 mm 2.6596 2.5025 2.1685 0.0165 0.0502 0.0174 分段高度H/m / 2.0 2.5 1.4 / / / -
[1] 郭慧, 江娃利, 谢新生, 2012. 山西交城断裂北端及中段3个大型探槽全新世断错现象分析. 地震地质, 34(1): 76—92 doi: 10.3969/j.issn.0253-4967.2012.01.008Guo H. , Jiang W. L. , Xie X. S. , 2012. Analysis of Holocene faulting phenomena revealed in the three trenches along the northern and central Jiaocheng Fault, Shanxi. Seismology and Geology, 34(1): 76—92. (in Chinese) doi: 10.3969/j.issn.0253-4967.2012.01.008 [2] 李峰, 徐锡伟, 陈桂华等, 2008. 高精度测量方法在汶川MS 8.0地震地表破裂带考察中的应用. 地震地质, 30(4): 1065—1075 doi: 10.3969/j.issn.0253-4967.2008.04.022Li F. , Xu X. W. , Chen G. H. , et al. , 2008. The application of different high-precision survey methods in the investigation of the MS 8.0 Wenchuan earthquake surface ruptures. Seismology and Geology, 30(4): 1065—1075. (in Chinese) doi: 10.3969/j.issn.0253-4967.2008.04.022 [3] 刘静, 汪良谋, 1996. 震级与地震破裂参数统计关系式系数的讨论. 地震地质, 18(3): 225—236Liu J. , Wang L. M. , 1996. Discussion on coefficients of statistical regressions between magnitude and earthquake rupture parameters. Seismology and Geology, 18(3): 225—236. (in Chinese) [4] 毛凤英, 张培震, 1995. 古地震研究的逐次限定方法与新疆北部主要断裂带的古地震研究. 活动断裂研究, (4): 153—164Mao F. Y. , Zhang P. Z. , 1995. Progressive constraining method in paleoseismic study and paleoearthquakes along the major active faults in northern Xinjiang. Research on Active Fault, (4): 153—164. (in Chinese) [5] 王林, 田勤俭, 李德文, 等, 2017. 蔚广盆地南缘断裂带唐山口段山前断层活动性研究. 震灾防御技术, 12(2): 319—328 doi: 10.11899/zzfy20170208Wang L. , Tian Q. J. , Li D. W. , et al. , 2017. The activity of the piedmont fault along the Tangshankou segment of Yuguang basin southern marginal fault. Technology for Earthquake Disaster Prevention, 12(2): 319—328. (in Chinese) doi: 10.11899/zzfy20170208 [6] 谢新生, 江娃利, 孙昌斌等, 2008. 山西交城断裂带多个大探槽全新世古地震活动对比研究. 地震地质, 30(2): 412—430 doi: 10.3969/j.issn.0253-4967.2008.02.007Xie X. S. , Jiang W. L. , Sun C. B. , et al. , 2008. Comparison study on Holocene paleoseismic activities among multi-trenches along the Jiaocheng Fault zone, Shanxi. Seismology and Geology, 30(2): 412—430. (in Chinese) doi: 10.3969/j.issn.0253-4967.2008.02.007 [7] Bi L. S., He H. L., Wei Z. Y., et al., 2012. Fractal properties of landforms in the ordos block and surrounding areas, China. Geomorphology, 175—176: 151—162. [8] Chase C. G. , 1992. Fluvial landsculpting and the fractal dimension of topography. Geomorphology, 5(1—2): 39—57. doi: 10.1016/0169-555X(92)90057-U [9] Davies S. , Hall P. , 1999. Fractal analysis of surface roughness by using spatial data. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 61(1): 3—37. doi: 10.1111/1467-9868.00160 [10] Elliot J. K. , 1989. An investigation of the change in surface roughness through time on the foreland of Austre Okstindbreen, North Norway. Computers & Geosciences, 15(2): 209—217. [11] Galli P. , Galadini F. , Pantosti D. , 2008. Twenty years of paleoseismology in Italy. Earth-Science Reviews, 88(1—2): 89—117. doi: 10.1016/j.earscirev.2008.01.001 [12] Giaccio B. , Galadini F. , Sposato A. , et al. , 2003. Image processing and roughness analysis of exposed bedrock fault planes as a tool for paleoseismological analysis: results from the Campo Felice fault (central Apennines, Italy). Geomorphology, 49(3—4): 281—301. doi: 10.1016/S0169-555X(02)00191-5 [13] Goodall H. J. , Gregory L. C. , Wedmore L. N. J. , et al. , 2021. Determining histories of slip on normal faults with bedrock scarps using cosmogenic nuclide exposure data. Tectonics, 40(3): e2020TC006457. [14] He H. L. , Wei Z. Y. , Densmore A. , 2016. Quantitative morphology of bedrock fault surfaces and identification of paleo-earthquakes. Tectonophysics, 693: 22—31. doi: 10.1016/j.tecto.2016.09.032 [15] Jones R. R. , Kokkalas S. , McCaffrey K. J. W. , 2009. Quantitative analysis and visualization of nonplanar fault surfaces using terrestrial laser scanning (LiDAR)—The Arkitsa fault, Central Greece, as a case study. Geosphere, 5(6): 465—482. [16] Klinkenberg B., 1992. Fractals and morphometric measures: is there a relationship? Geomorphology, 5(1—2): 5—20. [17] Mandelbrot B. , 1967. How long is the coast of Britain? Statistical self-similarity and fractional dimension. Science, 156(3775): 636—638. doi: 10.1126/science.156.3775.636 [18] McCalpin J., 2009. Paleoseismology. 2nd ed. Burlington: Academic Press. [19] Palumbo L. , Benedetti L. , Bourlès D. , et al. , 2004. Slip history of the Magnola fault (Apennines, Central Italy) from 36Cl surface exposure dating: evidence for strong earthquakes over the Holocene. Earth and Planetary Science Letters, 225(1—2): 163—176. doi: 10.1016/j.jpgl.2004.06.012 [20] Power W. L. , Tullis T. E. , 1991. Euclidean and fractal models for the description of rock surface roughness. Journal of Geophysical Research: Solid Earth, 96(B1): 415—424. doi: 10.1029/90JB02107 [21] Ran Y. K. , Chen L. C. , Chen J. , et al. , 2010. Paleoseismic evidence and repeat time of large earthquakes at three sites along the Longmenshan fault zone. Tectonophysics, 491(1—4): 141—153. doi: 10.1016/j.tecto.2010.01.009 [22] Sagy A. , Brodsky E. E. , Axen G. J. , 2007. Evolution of fault-surface roughness with slip. Geology, 35(3): 283—286. doi: 10.1130/G23235A.1 [23] Sagy A. , Brodsky E. E. , 2009. Geometric and rheological asperities in an exposed fault zone. Journal of Geophysical Research: Solid Earth, 114(B2): B02301. [24] Sieh K. E. , 1978. Prehistoric large earthquakes produced by slip on the San Andreas Fault at Pallett Creek, California. Journal of Geophysical Research: Solid Earth, 83(B8): 3907—3939. doi: 10.1029/JB083iB08p03907 [25] Sieh K. E. , 1984. Lateral offsets and revised dates of large prehistoric earthquakes at Pallett Creek, southern California. Journal of Geophysical Research: Solid Earth, 89(B9): 7641—7670. doi: 10.1029/JB089iB09p07641 [26] Sung Q. C., Chen Y. C., Chao P. C., 1998. Spatial variation of fractal parameters and its geological implications. Terrestrial, Atmospheric and Oceanic Sciences, 9(4): 655—672. [27] Sung Q. C. , Chen Y. C. , 2004. Self-affinity dimensions of topography and its implications in morphotectonics: an example from Taiwan. Geomorphology, 62(3—4): 181—198. doi: 10.1016/j.geomorph.2004.02.012 [28] Wells D. L. , Coppersmith K. J. , 1994. New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. Bulletin of the Seismological Society of America, 84(4): 974—1002. [29] Wiatr T. , Papanikolaou I. , Fernández-Steeger T. , et al. , 2015. Bedrock fault scarp history: insight from t-LiDAR backscatter behaviour and analysis of structure changes. Geomorphology, 228: 421—431. doi: 10.1016/j.geomorph.2014.09.021 [30] Young J. J. , Arrowsmith J. R. , Colini L. , et al. , 2002. Three-dimensional excavation and recent rupture history along the Cholame segment of the San Andreas fault. Bulletin of the Seismological Society of America, 92(7): 2670—2688. doi: 10.1785/0120000604 [31] Zou J. J. , He H. L. , Yokoyama Y. , et al. , 2020. Seismic history of a bedrock fault scarp using quantitative morphology together with multiple dating methods: a case study of the Luoyunshan piedmont fault, southwestern Shanxi Rift, China. Tectonophysics, 788: 228473. doi: 10.1016/j.tecto.2020.228473 [32] Zou J. J., He H. L., Yokoyama Y., et al., 2021. Identification of Paleoearthquakes and Coseismic Slips on a Normal Fault Using High-Precision Quantitative Morphology: application to the Jiaocheng Fault in the Shanxi Rift, China. Lithosphere, (Special 2): 2550879. [33] Zreda M. , Noll J. S. , 1998. Ages of prehistoric earthquakes revealed by cosmogenic chlorine-36 in a bedrock fault scarp at Hebgen Lake. Science, 282(5391): 1097—1099. doi: 10.1126/science.282.5391.1097 -