利用跨断层资料研究首都圈地区断层活动特征及地震危险性

本文针对首都圈的跨断层形变资料（1984-2016年）,利用11个场地的基线和水准观测资料,计算了断层三维活动量,以此来分析首都圈断层现今的活动特征。同时计算了断层形变协调比序列,分析其与首都圈中强地震的关系。结果表明:首都圈的断层活动有明显的区域特征,受张家口-渤海断裂带影响的西北地区的断层以左旋张性活动为主,而中部和东部断层表现为右旋压性特征。断层协调比结果显示在首都圈及其邻区出现中强地震前后,协调比会出现发散现象,这种现象可以作为今后首都圈及邻近地区中强地震危险性分析的一种依据。     

Effect of the time-weakening friction law during the nucleation process

Based on dynamic rupture simulations on a planar fault in a homogeneous half-space, we investigated the nucleation processes using the time-weakening friction law. Both the characteristic time and the rupture speed in the nucleation asperity play an important role in determining rupture behaviors on a fault plane following the time-weakening friction law, with which rupture starts from a single point in the nucleation asperity and propagates at a given speed toward the boundary of the nucleation area. Rupture with a small characteristic time or a large rupture speed in the nucleation asperity propagates earlier from the hypocenter. Rupture following the slip-weakening friction law requires a smaller radius of nucleation patch to have similar rupture front contours of the time-weakening friction law. Even if the rupture velocity in the nucleation patch of the time-weakening friction law increases to infinity, the peak slip rate in the nucleation asperity is smaller than that of the slip-weakening law. The peak ground velocity distributions of ruptures following the two friction laws are also compared.     

2015年尼泊尔Gorkha地震强地面运动记录分析

2015年4月25日在尼泊尔Gorkha地区发生M_W7.8地震,距离发震断层约11 km的KATNP台站完整记录了主震的加速度时程.本文根据KATNP台站记录的加速度数据分析了Gorkha地震的地震动特征.结果表明Gorkha地震在KATNP台站处产生的水平向峰值加速度为0.17 g,竖直向峰值加速度为0.19 g,该数值小于科学家们对如此大规模地震产生的地震动的预期,初步推测这可能是由加德满都山谷产生的非线性响应造成的（Dixit et al.,2015）;地震在KATNP台站处产生了地表永久位移,其中竖向永久位移为131.9 cm,水平向永久位移的绝对值为159.2 cm,方向为南偏西19°（199°）,据此可简单推算出断层走向约为289°（109°）.地震产生了脉冲型地震动,影响因素有盆地效应、地震破裂的向前的方向性效应以及滑冲效应,其中盆地效应的周期约为5 s左右,方向性效应产生的速度脉冲的周期约为8 s左右.加速度反应谱显示在0.5 s和5.0 s左右各有一个峰值,前者是由地震破裂的脉冲式滑移产生的大量高频地震动造成的,后者是由于盆地效应和地震破裂的方向性效应造成的.基于阿里亚斯烈度计算的地震动持时约在36~46 s之间,小于与其规模相当的地震产生的地震动持时,并且不同方向上的地震动持时可能与地震破裂方向有关.阿里亚斯烈度随时间的变化比较简单,也反映了Gorkha地震是一次连续的、能量释放相对简单的地震事件.     

Using UAV acquired photography and structure from motion techniques for studying glacier landforms: application to the glacial flutes at Isfallsglaciären

Glacier and ice sheet retreat exposes freshly deglaciated terrain which often contains small‐scale fragile geomorphological features which could provide insight into subglacial or submarginal processes. Subaerial exposure results in potentially rapid landscape modification or even disappearance of the minor‐relief landforms as wind, weather, water and vegetation impact on the newly exposed surface. Ongoing retreat of many ice masses means there is a growing opportunity to obtain high resolution geospatial data from glacier forelands to aid in the understanding of recent subglacial and submarginal processes. Here we used an unmanned aerial vehicle to capture close‐range aerial photography of the foreland of Isfallsglaciären, a small polythermal glacier situated in Swedish Lapland. An orthophoto and a digital elevation model with ~2 cm horizontal resolution were created from this photography using structure from motion software. These geospatial data was used to create a geomorphological map of the foreland, documenting moraines, fans, channels and flutes. The unprecedented resolution of the data enabled us to derive morphological metrics (length, width and relief) of the smallest flutes, which is not possible with other data products normally used for glacial landform metrics mapping. The map and flute metrics compare well with previous studies, highlighting the potential of this technique for rapidly documenting glacier foreland geomorphology at an unprecedented scale and resolution. The vast majority of flutes were found to have an associated stoss‐side boulder, with the remainder having a likely explanation for boulder absence (burial or erosion). Furthermore, the size of this boulder was found to strongly correlate with the width and relief of the lee‐side flute. This is consistent with the lee‐side cavity infill model of flute formation. Whether this model is applicable to all flutes, or multiple mechanisms are required, awaits further study. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Subaerial gravel size measurement using topographic data derived from a UAV‐SfM approach

Accurate and reliable methods for quantifying grain size are important for river science, management and in various other sedimentological settings. Remote sensing offers methods of quantifying grain size, typically providing; (a) coarse outputs (c. 1 m) at the catchment scale where individual grains are at subpixel level, or; (b) fine resolution outputs (c. 1 mm) at the patch scale. Recently, approaches using unmanned aerial vehicles (UAVs) have started to fill the gap between these scales, providing hyperspatial resolution data (< 10 cm) over reaches a few hundred metres in length, where individual grains are at suprapixel level. This ‘mesoscale’ is critical to habitat assessments. Most existing UAV‐based approaches use two‐dimensional (2D) textural variables to predict grain size. Validation of results is largely absent however, despite significant differences in platform stability and image quality obtained by manned aircraft versus UAVs. Here, we provide the first quantitative assessment of the accuracy and precision of grain size estimates produced from a 2D image texture approach. Furthermore, we present a new method which predicts subaerial gravel size using three‐dimensional (3D) topographic data derived from UAV imagery. Data is collected from a small gravel‐bed river in Cumbria, UK. Results indicate that our new topographic method gives more accurate measures of grain size (mean residual error ‐0.0001 m). Better results for the image texture method may be precluded by our choice of texture measure, the scale of analysis or the effects of image blur resulting from an inadequate camera gimbal. We suggest that at our scale of assessment, grain size is more strongly related to 3D variation in elevation than to the 2D textural patterns expressed within the imagery. With on‐going improvements, our novel method has potential as the first grain size quantification approach where a trade‐off between coverage and resolution is not necessary or inherent. Copyright © 2017 John Wiley & Sons, Ltd.     

3‐D uncertainty‐based topographic change detection with structure‐from‐motion photogrammetry: precision maps for ground control and directly georeferenced surveys

Structure‐from‐motion (SfM) photogrammetry is revolutionising the collection of detailed topographic data, but insight into geomorphological processes is currently restricted by our limited understanding of SfM survey uncertainties. Here, we present an approach that, for the first time, specifically accounts for the spatially variable precision inherent to photo‐based surveys, and enables confidence‐bounded quantification of 3D topographic change. The method uses novel 3D precision maps that describe the 3D photogrammetric and georeferencing uncertainty, and determines change through an adapted state‐of‐the‐art fully 3D point‐cloud comparison (M3C2), which is particularly valuable for complex topography. We introduce this method by: (1) using simulated UAV surveys, processed in photogrammetric software, to illustrate the spatial variability of precision and the relative influences of photogrammetric (e.g. image network geometry, tie point quality) and georeferencing (e.g. control measurement) considerations; (2) we then present a new Monte Carlo procedure for deriving this information using standard SfM software and integrate it into confidence‐bounded change detection; before (3) demonstrating geomorphological application in which we use benchmark TLS data for validation and then estimate sediment budgets through differencing annual SfM surveys of an eroding badland. We show how 3D precision maps enable more probable erosion patterns to be identified than existing analyses, and how a similar overall survey precision could have been achieved with direct survey georeferencing for camera position data with precision half as good as the GCPs'. Where precision is limited by weak georeferencing (e.g. camera positions with multi‐metre precision, such as from a consumer UAV), then overall survey precision can scale as n^‐½ of the control precision (n = number of images). Our method also provides variance–covariance information for all parameters. Thus, we now open the door for SfM practitioners to use the comprehensive analyses that have underpinned rigorous photogrammetric approaches over the last half‐century. Copyright © 2017 John Wiley & Sons, Ltd.     

ANALYSIS ON THE CHARACTERISTICS OF DURATION AND PERIOD OF GROUND MOTION OF THE LUSHAN EARTHQUAKE BASED ON THE STATION RECORDS

Based on 60 records from the 20 stations within 100km to the epicenter of Lushan earthquake, the predominant period, period of peak response spectrum, duration of ground motion, and source duration are investigated. By the study, we conclude that within 100km to the epicenter, the scope of predominant period is 0.013~0. 275s in EW, NS and UD direction; the scope of period of peak response spectrum for 5% damping ratio is 0.03~0.65s; the scope of 90% ground motion durations is 5. 1~35. 9s; the scope of averaging source duration is 6.41 (EW), 6.05 (NS) and 5.47s (UD). Furthermore, based on the ground motion duration calculated by 20 stations dada, the predictive curve and equation of ground motion duration is regressed and compared with the recent equation of ground motion duration by Bommer (2009). We find that the ground motion durations of most stations are larger than predictively mean value by Bommer (2009), which means that the source effect of Lushan earthquake is relatively larger. Lastly, by the contour figure of ground motion duration and source duration, we conclude that the directivity character of duration is relatively apparent in NE direction. The relatively larger source duration and ground motion duration in NE direction indicate more energy release in the main shock of Lushan earthquake, which perhaps causes the relatively less aftershocks in this direction. The duration has no hanging wall effect, which perhaps results from the blind-reverse fault structure of Lushan earthquake.     

Development of duration‐dependent damage‐based inelastic response spectra

A simple relationship is proposed in this paper to construct damage‐based inelastic response spectra including the effect of ground motion duration that it can be used for damage control in seismic design of structures. This relation is established for three groups of ground motions with short‐duration, moderate‐duration, and long‐duration ranges. To develop the model, the duration effect is included in the cyclic ductility of structures by an energy‐based method, and then strength reduction factors are computed based on this modified ductility (named ). The strength reduction factors were calculated for 44 stiffness‐degrading oscillators having vibration periods between 0.05 and 4.0 s, four ultimate ductility capacities, and five damage levels subjected to 296 earthquake records. The results showed that ductility capacity, damage level, and ground motion duration are effective parameters in the energy dissipation of structures, which affect the spectra. The values of short‐period oscillators (e.g., low‐rise structures) under short‐duration records are generally greater than those under moderate‐duration and long‐duration records. Residual analysis has been made in terms of magnitude and distance to examine the validity of the proposed simple expression. Finally, the introduced spectra were compared with three previously published proposals. Copyright © 2016 John Wiley & Sons, Ltd.     

基于NNSIM随机有限断层法的7.0级九寨沟地震强地面运动场重建

北京时间2017年8月8日21时19分在四川阿坝州九寨沟县（北纬33.2°、东经103.8°）发生了7.0级地震,作者利用新提出的震源破裂过程控制的NNSIM随机有限断层方法,对震中约100 km范围内强地面运动进行了情景构建,并与研究区范围内的8个强震测震台观测记录进行了加速度、速度、反应谱对比,对比的一致性说明模拟过程中所应用的各种震源参数适合九寨沟7.0地震的近震源强地面运动重建.论文的后半部分基于相同参数模拟了九寨沟7.0级地震的近震源强地面运动场的空间分布,结果显示此次地震强地面运动的PGA、PGV以及反应谱分布均呈现出以震中为圆心的圆形分布,不具有大地震沿断层走向的分布特点.论文最后关于反应谱的讨论,揭示了此次九寨沟县城灾害不重的直接原因.     

基于强震动资料的破裂过程快速反演及其自动化的可行性

破裂过程快速反演是目前快速获取地震灾害特征的主要手段之一,是震后应急工作的重要内容.近十年来基于远震资料开展的手动快速反演工作取得了长足进步,但在响应时间方面存在固有的局限,阻碍了反演效率的持续提升.我们根据新近发展的IDS （Iterative Deconvolution and Stacking）自动反演方法,尝试反演近场强震动资料确定破裂过程,探讨破裂过程反演自动化的可行性.对近几年国内发生的强震——包括2013年芦山M_W6.6地震、2016年青海门源M_W5.9地震和2016年新疆阿克陶M_W6.6地震——的应用结果表明,采用IDS方法反演强震数据可以得到稳定可靠的破裂模型,且反演计算时间都控制在几十秒内.此外,以2008年汶川M_W7.9地震为例,测试了不同子断层尺度、截止频率和地壳速度结构模型对反演结果的影响,发现滑动分布主要特征不强烈依赖于反演参数和地壳模型,证实了自动反演的稳定性和很强的适应能力.这一研究表明,基于强震动资料的自动反演可能是破裂过程快速反演的主要发展方向.特别地,在未来强震动台网持续发展、强震动数据的质量和共享速度都得到进一步提高之后,这一工作可望纳入到地震参数的常规自动测定工作中,为震后应急和海域地震的海啸预警提供急需的震源模型.