地震破裂带附近同震位移和变形模式对于深入理解地震破裂过程、断层行为及活动断层与地形地貌关系等至关重要。文中提出一种迭代最近点(ICP)算法, 利用地震前后地形点云进行差分确定断层近场三维同震地表位移。在川西大凉山断裂带交际河断层上选取2期SfM地形点云叠加同震位移场模拟震前-震后点云, 测试ICP方法获取同震位移场的精度。该方法在网格边长>50m的条件下可准确恢复同震位移场的方向和幅度, 水平和垂直精度分别为20~10cm, 这与地形点云定位精度相当。随着点云密度和网格窗口尺寸减小, 该方法恢复同震位移场等的精度将降低。通过分析树木生长、房屋建设、河流侵蚀等地形变化对恢复位移场的潜在影响, 研究结果表明扩大网格尺寸可使震前-震后点云具有足够地形结构进行匹配, 以减小局部地形变化对恢复位移场的影响, 网格窗口尺寸是在具有足够地形结构的大尺度和具有更精细分辨率的小尺度之间的权衡。文中所述的ICP方法利用震前-震后高精度点云可获取地震破裂带附近精细的三维地表位移场, 为浅层断层滑动和破裂带变形提供了新的约束, 有助于研究地震破裂过程和断层生长过程。
龙门山断裂带是青藏高原块体和华南块体的重要边界断裂带, 是一条重要的地震活动构造带。几何上自西向东可分为后山断裂、 中央断裂、 前山断裂和前陆盆地(成都平原)变形带。中央断裂带是龙门山断裂带的主干断裂, 映秀-北川断裂是中央断裂带中段活动最强烈的一支断裂。前人对龙门山断裂带映秀-北川断裂滑动速率进行了部分研究, 但结果差异很大, 范围从0.07mm/a到1.1mm/a不等。断层滑动速率是分析断层活动特征、 研究区域运动学和动力学机制的重要依据。文中依据遥感精细解译和野外调查资料, 通过10个阶地横剖面及阶地年龄对龙门山断裂带映秀-北川断裂白沙河段的滑动速率进行约束, 得到映秀-北川断裂白沙河段晚更新世以来的垂直滑动速率为(0.10±0.02)~(0.30±0.05)mm/a, 平均滑动速率为0.19mm/a。该结果可为理解汶川地震的发震构造、 龙门山断裂带的整体活动特征及评价龙门山断裂带的长期地震危险性提供基础资料。
随着断层面形貌测量、 宇宙成因核素测年等技术的发展, 基岩断层面成为在基岩区古地震研究的良好研究对象。但是, 受沟谷侵蚀、 沉积覆盖和人为改造等非构造因素影响, 并非所有出露的基岩断层面都能完整记录和保存古地震信息。因此, 必须在基岩断层面上选择合适的研究点。传统的野外测量方法耗时费力, 难以掌握基岩断层面的整体信息, 即使通过高分辨率卫星影像也无法获得基岩断层面的精细结构。小型无人机航测(sUAV)采用低空摄影测量技术, 可快速获取基岩断层面高精度的地表三维结构, 为筛选目标工作点提供便利。文中以山西地堑系六棱山北麓断裂的马家窑基岩断坎为例, 采用小型无人机航测结合SfM摄影测量技术, 获得了基岩断层面的高精度三维地形数据, 通过精细地貌解译和断层坎剖面分析识别了沿断层坎分布的崩塌破坏、 沉积埋藏和侵蚀剥露等非构造因素, 分析了侵蚀、 埋藏和构造这3种典型的基岩断层面的形态特征, 选取了适合开展古地震研究的基岩断层面点, 展现了小型无人机在基岩区古地震研究工作中发挥的重要应用潜力。
陡坎是一种自然界常见的台阶状地貌, 但其形成年龄通常很难直接测定。发育在松散堆积物中的陡坎经过初期短暂的重力垮塌之后, 将经历漫长的低能退化过程。如果陡坎剖面形态的演化可基于扩散方程来模拟, 且扩散系数可独立标定, 即可利用陡坎地形剖面估算其年龄, 这种方法被称为形貌测年。文中简要回顾了陡坎形貌测年的研究历史, 介绍并讨论了陡坎退化的概念模型与扩散模型, 特别是非线性扩散模型的建立及求解、 参数在扩散模型中所起的作用、 最佳陡坎形貌年龄的确定流程等, 分析了陡坎上、 下地貌面坡度对陡坎退化的影响, 编制了非线性扩散模型的年龄图版, 给出了图版的应用实例, 验证了形貌测年方法的有效性。线性扩散模型和非线性扩散模型均可用于单次事件陡坎的退化分析, 但对于年轻的单次事件陡坎推荐使用非线性扩散模型。断层重复活动形成的陡坎的退化分析则需要谨慎对待, 恒定滑动速率陡坎的非线性扩散模型适用于模拟年龄<10ka、 活动速率高的断层陡坎的演化; 多次事件陡坎模型(包括线性扩散和非线性扩散)需要仔细评估每次事件在陡坎剖面上的断错位置及其位移量。尽管陡坎形貌测年方法存在很多假设条件, 但目前快速获取一定范围内的高分辨率地形数据已成为现实, 从这些数据中可以沿着同一陡坎提取大量剖面进行分析, 继而得到具有统计意义的结果, 这为陡坎退化分析和形貌测年方法提供了广阔的应用前景。
同震地表破裂带的空间展布及形变特征是地球深部断层活动在地表的直观地貌表现, 不但记录着地震破裂和断层运动的信息, 还反映了区域应力和地壳运动状况。因此, 开展震后地震地表破裂带调查对于了解发震断层的构造活动尤为重要。高精度地形观测技术可以获取前所未有的高时空分辨率的地球表面特征, 为辨别历史地震地表破裂遗迹、 提取地表同震位移、 活动构造地质填图等提供高质量数据。文中选取富蕴1931年地震地表破裂带作为研究区, 利用SfM(Structural from Motion)摄影测量技术生成分辨率为1m的数字高程模型(DEM), 详细识别地表破裂并测量冲沟的右旋位移。基于地表破裂的几何及构造地貌特征, 将富蕴地震地表破裂带由北向南分为S1、 S2、 S3、 S4 4段, 其间以挤压隆起或拉分盆地相连接。沿破裂带共获得194组最新冲沟的右旋水平位移, 得到1931年同震位移的平均值为(5.06±0.13)m。同震位移局部缺失或突变的区域与几何阶区的位置也有良好的对应关系。以上结果填补了对富蕴地震地表破裂精细形态研究的空白, 也进一步展示高分辨率的地形数据在活动构造研究中良好的应用价值。
研究地震地表破裂的几何结构与同震位移分布的空间关联性, 对于理解断层面几何形貌与地震活动的相关性具有重要意义。 文中运用滑动窗口方法, 对28个震例的地表破裂非规则几何结构特征(走向变化、 阶区大小和粗糙度)和同震位移分别进行了量化和均值计算, 并计算了二者之间的相关性, 最后分析讨论了走滑型和倾滑型震例的相关性之间的差异性及成因。 研究结果表明: 1)相对于倾滑型, 走滑型地震地表破裂的几何结构与同震位移的相关性更高, 表明走滑型断层(破裂)的几何特征在分段中具有更高的参考价值; 2)在走滑型震例中, 地表破裂的几何结构与同震位移之间存在负相关关系, 即地表破裂上某段的阶区、 走向变化、 粗糙度的特征值越大, 其对应的同震位移量越小; 3)走滑型地震的地表破裂非规则几何结构特征中, 与同震位移分布的相关程度由高到低依次为阶区大小、 走向变化、 粗糙度。
构造地貌能够记录构造运动的长期积累,并能通过相关参数定量表达。但是,由于地貌还受到气候和岩性等因素的影响,所以如何利用相关参数来表达地貌与构造之间的关系是一个研究热点。文中利用水平分辨率为90m的SRTM数据,采用元分维模型,通过变差函数法计算了滇西南地区的二维分形参数。并对比了在不同岩性不同气候条件下分形参数的分布特征。结果显示,二维分形参数与岩性、气候因素相关性不强。而分维值与构造活动因素相关性很好,在构造活动强烈区域分维值较低,在构造活动不强烈区域分维值较高。这说明分维值能在一定程度上反映构造活动的强弱,为研究区域构造活动提供了一种新的手段。
大凉山断裂带是中国大陆大型强震断裂带鲜水河-小江断裂系的重要组成部分,其中段的普雄断裂是断裂带中最长的1条次级断裂。该次级断裂的古地震活动情况对评价该区域地震活动水平和建立地震灾害预防措施具有重要意义。近年开展的活动断层填图和古地震探槽研究表明普雄断裂是1条晚更新世以来活动强烈的略向W倾的高角度左旋走滑断层。沿断裂开挖的2个探槽分别揭露了2次和3次古地震事件,分别发生在8206 BC-1172 AD、1084-1549 AD和17434-7557 BC、1577-959 BC、927-1360 AD。结合探槽古地震结果和历史地震记载,进行建模分析得到断裂的离逝时间约为0.7ka,与距今倒数第2次事件的时间间隔约为2.3ka。同时,根据震级与地表破裂长度的经验公式估算古地震事件震级在7级以上。
探槽技术虽然已经发展成为古地震研究的重要手段, 结合定年技术可以识别出古地震事件与地震重复间隔, 但是仍然存在较大的不确定性和局限性。例如, 目视解译难以区分黄土沉积内部的细层理, 严重影响了古地震事件发生位置或时间的判定。如何提高古地震研究的精度和准确性, 降低古地震事件判定的不确定性, 是目前面临的一个迫切问题。山西洞峪沟剖面位于山西临汾盆地东北角, 横跨霍山山前断裂带, 不仅揭示了较好的黄土沉积序列, 还揭露了明显的地层错断事件。因此, 该剖面是一个开展高精度探槽古地震研究, 降低判定古地震事件不确定性的理想场所。根据高精度的粒度与磁化率变化曲线, 结合目视分层解译结果, 对洞峪沟黄土剖面进行了精细分层, 界定了各层的厚度和边界。依据细分层和断层两盘地层的对应关系, 将u6层沉积以来的3次断错事件的发生位置和时间作了再限定, 它们分别发生在u5-7、u4以及u2的顶部, 对应断层上盘埋深7.1m、4.7m与2.9m。根据释光测年结果以及断层上盘地层的平均沉积速率, 推测3次断错事件发生的时间分别在 (48.1±1.5)~(43.2±2.5)ka BP接近45.8ka BP、(35.0±2.4)~(30.6±1.3)ka BP接近32.8kaBP、(26.4±0.8)~(20.9±0.7)ka BP接近23.3ka BP。根据3次黄土-古土壤沉积旋回的厚度差, 判定3次地震的同震垂直位移分别为0.5m、0.4m和1.3m, 累计位移2.2m。依据高精度粒度与磁化率变化曲线的地层分层方法, 为有效降低黄土地区古地震研究的不确定性提供了一种很好的参考。
地形数据的质量(精度和分辨率)影响着地球科学的研究水平。LiDAR测量是目前获取高分辨率地形数据的有效技术方法之一, 但是其高昂的测量成本和相对复杂的后期数据处理限制了LiDAR技术的大众化应用。近年来, 一种被称为SfM(Structure from Motion)的适合大众化使用的新的高精度3维地形数据获取技术开始引起人们的注意。这种新型数字摄影测量技术可以利用高效的图像特征匹配算法从多视角照片中提取重叠区域的3维地形数据。由于SfM技术仅需要目标物体的照片, 而且对相机拍摄位置、图像尺度及拍摄焦距没有要求, 因此利用简单测量平台采集地面照片就可以获取高质量的3维地形数据。与LiDAR技术相比, 大大降低了获取高精度数据的成本, 使得高精度3维地形数据的使用大众化。文中介绍了SfM技术的基本原理和流程, 展示了SfM技术获取高精度3维地形数据的简单而有效的特性, 特别适合于植被稀少的区域。文中利用近千米高空拍摄的、具有约70%重叠度的一套随LiDAR飞行采集的数字航空照片生成具有真彩色的高密度SfM点云数据, 点密度高达25.5个/m2, 可生成分辨率0.2m的DEM(数字高程模型)。对比相同区域的LiDAR点云数据, 统计分析表明58.3%的LiDAR数据与SfM数据的垂直偏差<0.1m, 88.3%的LiDAR数据的垂直偏差<0.2m; 而且发现不同地貌的SfM数据精度存在差异, 平缓地形的SfM数据精度高于陡峭地形的SfM数据精度。文中还介绍了以氦气球作为拍摄平台的SfM测量系统, 可以快捷地获取高精度的3D地形数据和正射影像, 比目前常用的差分GPS测量具有更高的效率和数据精度。
大凉山断裂带是川滇菱形活动块体东边界的重要组成部分, 其活动习性和滑动速率对于了解青藏高原东南缘的块体运动和构造变形具有重要的理论意义, 同时对于该地区地震危险性评价和地震中长期预测有着重要的应用价值。大凉山断裂带由6条分支断裂构成, 竹马断裂是其北段的一条分支。文中在详细野外调查的基础上, 通过典型构造地貌的高精度GPS测量和断错时间的约束, 得出竹马断裂是一条以左旋走滑为主的全新世活动断裂, 全新世以来的水平滑动速率为1.5~3.1mm/a; 利用古地震探槽揭示出竹马断裂的2次古地震事件时间分别为(50.3±5.7)~30ka BP和30~(17.4±1.2)ka BP。考虑到断裂带北段还存在一条与竹马断裂平行的分支断裂——公益海断裂, 大凉山断裂带北段的滑动速率与南段基本一致。鲜水河-小江断裂系中段在大凉山断裂带上分配的滑动速率与并行的安宁河断裂带和则木河断裂带也大体相当, 断裂系中段的滑动速率之和与南北两段的速率大致吻合。大凉山断裂带的存在使断裂系在几何上成为一个完整的弧形构造, 也弥补了中段滑动速率的亏损, 使各段的滑动保持协调一致。此外, 沿竹马断裂大量发育的冲洪积台地, 根据沉积特征和测年结果, 认为是新仙女木事件后末次冰消期的冰川融水形成的冰水堆积, 类似的地貌面可能广泛分布于青藏高原东南缘。
基岩断层面形貌特征的定量分析方法是研究断层活动历史和识别古地震的有效手段, 是对探槽技术的有力补充, 可以弥补在基岩区开展古地震研究技术手段的严重不足。文中以山西霍山山前断裂上的断层面(崖)为例, 采用各向同性变差函数法计算了断层面形貌的2D分维值, 并通过沿断层走向逐条平均的方式获得了2D分维值在断层倾向上的分布特征, 结果显示断层面形貌在倾向上具有显著的分带性特征, 由此获得的每个分带的特征分维值随断层面高度的增加呈阶跃式增加。这种断层面形貌在倾向上的显著分带性特征、特征分维值的阶跃式增加反映了断层面的出露方式不是渐进式, 而是不连续的间歇式; 这种间歇式分带性出露方式与断层的周期性地震活动相关, 每个分带的特征值与每个分带的出露时间相关, 即与地震发生时间相关。因此: 1)可以通过每一个断层面分带的特征分维值估计该分带出露的时间; 2)断层面分带的宽度可以用来估计每次破裂地震的同震位移量。霍山山前断裂带上3个基岩断层面形貌分带揭示出3次地震事件; 3个2D特征分维值不仅显示3个分带自上而下由老变新的阶跃特征, 而且3个阶跃之间大致相等的特征分维值差显示了3次地震发生间隔是大致相等的; 断层面分带的宽度反映出破裂地震的倾向同震位移量大约为3.5m。此外, 在阶跃式分带之间还存在宽0.5~1m的较窄分带, 这些窄分带的特征分维值是随断层面高度增加而逐渐增大的, 显示在间震期基岩断层面被逐渐剥露的特征。
机载LiDAR技术为描绘活动构造相关构造地貌和最新的地表形变提供更精确的基础数据。如何将LiDAR新技术、新数据应用于活动构造填图和活动断层地震危险性评价等方面,是今后活动构造研究领域的一个重要的发展方向。文中以新疆天山北麓的独山子背斜-逆冲断裂带为试验区,开展了基于LiDAR数据的活动构造填图实验研究。首先,采用机载LiDAR技术进行数据采集,获得点云密度为6.6个/m2、平均点间距为0.39m的LiDAR原始数据;其次,利用试验区内12个测量精度可达mm级的GPS静态测量点评估LiDAR的相对垂直精度为0.12m、均方差值为0.078m;最后,对密度为6.4个/m2的地面点云数据进行DEM最佳分辨率评估,利用反距离权重算法获得0.5m分辨率的数据高程模型(DEM)。该分辨率的DEM数据足以完成独山子背斜-逆冲断裂带的精细构造地貌特征的确定以及高精度的空间解译。文中仅使用DEM可视化工具从不同虚拟的视角、不同色度或其他处理方式来识别微构造地貌、划分地貌面和确定断层位置等,宏观上获得与前人通过航片解译和野外调查一致的断裂分布特征,微观上较前者具有更高的精细程度。此外,数据采集、数据质量检验、数据处理及数据应用等技术和方法适用于其他能够获得LiDAR地形数据的活动断裂研究工作。
2013年4月20日的芦山 "4·20"MS7.0地震发生在龙门山断裂带西南段,震中地区分布多条NE向断裂,构造较为复杂。这次地震震源机制解显示为逆冲型地震,破裂面为NE走向,与龙门山断裂带的运动性质和走向一致。地表调查只在大川-双石断裂(前山断裂)和新开店断裂(大邑断裂南段)发现局部分布的NE向地表裂缝、沿地表裂缝分布的喷砂冒水和砂土液化,不规则的边坡开裂等地表变形,以及断裂沿线较严重的滑坡崩塌和房屋破坏。野外调查没有发现明显的地震地表破裂。GPS测量结果显示,此次地震的发震断裂位于芦山县城附近或其以东,而芦山西侧的断裂也可能参与了部分活动。根据野外地质调查、GPS观测、震源机制解、震源深度、余震分布等结果综合判定,芦山7.0级地震的主要发震构造是芦山之下、大川-双石断裂和新开店断裂之间的龙门山前缘滑脱带。此滑脱带在该段的运动导致了这次地震的发生,并可能带动了它上面的大川-双石和新开店等断裂的活动。
2013年4月20日,四川省芦山县发生了MS 7.0地震。文中简要介绍了芦山地震的基本情况与芦山地震区历史地震及其相关地震滑坡情况。依据2008年汶川地震滑坡与地震动峰值加速度(PGA)的空间关系,对芦山地震滑坡大体分布范围进行了推测。根据地震滑坡分类学,将芦山地震滑坡分为破坏型滑坡、连贯型滑坡、流滑型滑坡3大类。其中,破坏型滑坡包括岩质崩塌、岩质滑动、岩质崩滑、土质崩塌、土质滑动等5类; 连贯型滑坡包括土质坍塌与慢土流2类; 流滑型滑坡为快速流滑。破坏型滑坡如岩质崩塌、岩石滑动、土质崩塌这3类是芦山地震滑坡中最常见的类型。基于震后可利用的高分辨率航片,初步解译得到3 883处滑坡位置点数据。最后,从余震对滑坡的影响,芦山地震滑坡与邻区地震滑坡对比分析,对后续基于高分辨率遥感影像的滑坡精细解译的启示等3个方面开展了分析与讨论。
以霍山山前断裂为实验区,基于高分辨率IRS-P5 DEM数据提取了横穿断裂的冲沟,在23条冲沟纵剖面上识别出断裂活动诱发裂点。在其中5条冲沟上识别出1级裂点,裂点几乎位于断裂上; 在2条冲沟上识别出2级裂点,最新裂点也位于断裂上。这些位于断裂上的裂点,高度达4~9m,推测它们形成之后并未向上游明显迁移,并可能接受了多次事件的叠加。在其他16条冲沟上识别出2~3级裂点,最新一级裂点自断裂向上游迁移了40~70m,次新一级裂点主要分布在距断裂150~250m的河段内,第3级裂点距断裂300~500m。设定最新裂点形成于1303年洪洞8级地震(事件Ⅲ)发生时,再假定全新世中晚期以来研究区内河流溯源侵蚀速率基本不变,则推算出次新裂点对应的古地震事件(事件Ⅱ)发生在距今3 336~2 269a之间,最早一次事件(事件Ⅰ)发生在距今5 618~4 504a之间,霍山山前断裂强震重复周期为1 500~2 600a。这与前人通过探槽揭露的古地震事件序列和强震重复周期基本一致。