基于密集台阵的青藏高原东北缘地壳精细结构及九寨沟地震震源区结构特征分析

基于青藏高原东北缘密集宽频带野外流动观测台阵以及固定台站资料,利用双差层析成像方法对地震位置和研究区的地壳速度结构进行了反演.最终用于联合反演的地震事件合计9644个.结果显示青藏高原东北缘速度结构具有明显的横向不均匀性.从整体上看,青藏高原地区表现为低速异常,鄂尔多斯表现为高速异常,而扬子地块亦表现为高速异常.不同深度处速度结构表现不一致,同一深度处P波速度结构和S波速度结构也有明显差异.由西秦岭北缘断裂带、临潭-宕昌断裂以及礼县-罗家堡断裂围限的地震活动强烈的区域中,P波速度结构由深度0 km时呈现的低速异常,逐渐过渡到5 km时高低速相间分布的特征;而S波速度结构在此区域中,由近地表0 km时高低速相间分布的特征,逐渐过渡到30 km时几乎表现为低速异常.2017年8月8日九寨沟7级地震所在的塔藏断裂、岷江断裂和雪山断裂围限区域,在深度20 km处的P波速度结构和周围存在明显差异,九寨沟地震处于高速异常与低速异常的过渡带内.此外,2013年7月22日发生在青藏高原东北缘的岷漳县6.6级地震,震源区所在的临潭-宕昌断裂附近的P波速度结构在15 km深度处也有明显特征,震源位置所在区域也处于高低速过渡带.该区域这种地壳内部高低速过渡带可能是应力比较容易积累而发生中强地震的一个重要场所.     

地震槽波的数学-物理模拟初探

针对地震槽波在低速层的传播特性,开展了煤层内地震槽波勘探的数值模拟和物理模拟研究的初探工作.在数值模拟研究方面,采用交错网格有限差分法对煤层中的地震槽波进行三分量全波场模拟.基于波场快照和人工合成地震记录研究了不同模型中的波场特征和各种波型的传播规律.在物理模拟方面,通过选用不同配比的环氧树脂和硅橡胶类材料构建地震槽波物理模型,利用透射法和反射法观测系统获得了清晰的地震槽波记录以研究槽波的地震学特征.研究表明,在煤层内槽波的地震波场中,Love型槽波的能量小于Rayleigh型槽波的SV分量,大于Rayleigh型槽波的SH分量.相对于Love型槽波和Rayleigh型槽波的SH分量,Rayleigh型槽波的SV分量在围岩中的泄露能量较强.在煤层界面附近的围岩中,地震波仍以槽波形式传播,随着距离的增加能量逐渐衰减.随着煤层变薄,煤层槽波主频向高频方向移动,频散现象增强,传播速度增大.     

滇西南地区南汀河断裂带三维深部电性结构及其孕震环境

南汀河断裂带为滇西南地区活动断裂体系中规模最大的一条北东向断裂,其构造活动及地震危险性一直备受关注.本文基于覆盖云南境内南汀河断裂带的大地电磁测深宽频带阵列数据,利用大地电磁三维反演解释技术,首次获得了南汀河断裂带的精细三维深部电性结构.在上地壳深度,南汀河断裂带西南段与中段的电性结构表现出沿构造走向的高导条带特征,北东段表现为高阻结构.该高阻结构可能为临沧—勐海花岗岩体的电性反映,指示南汀河断裂可能未切穿该花岗岩带.在中下地壳深度,南汀河断裂带西南段存在大范围高导层,北东段则表现为整体性的高阻地壳,因此南汀河断裂北东段可能具有发生强震的介质结构背景.南汀河断裂带西南段的耿马地震区深部呈现北东向与北北西向的"X"型高导构造样式,该高导结构以南存在一个显著高阻异常体,1988年耿马M_S7.2地震以及2015年沧源M_S5.5级地震均发生于该高阻体与"X"型高导条带的电性边界.青藏高原东南缘绕东构造结流入滇西地区的中下地壳流可能受到南汀河断裂北东段中下地壳高阻体的阻挡而呈分流式分布于保山地块以及澜沧江断裂以东.     

The effect of lithology on valley width,terrace distribution,and bedload provenance in a tectonically stable catchment with flat‐lying stratigraphy

How rock resistance or erodibility affects fluvial landforms and processes is an outstanding question in geomorphology that has recently garnered attention owing to the recognition that the erosion rates of bedrock channels largely set the pace of landscape evolution. In this work, we evaluate valley width, terrace distribution, and bedload provenance in terms of reach scale variation in lithology in the study reach and discuss the implications for landscape evolution in a catchment with relatively flat‐lying stratigraphy and very little uplift. A reach of the Buffalo National River in Arkansas was partitioned into lithologic reaches and the mechanical and chemical resistance of the main lithologies making up the catchment was measured. Valley width and the spatial distribution of terraces were compared among the different lithologic reaches. The surface grain size and provenance of coarse (2–90 mm) sediment of both modern gravel bars and older terrace deposits that make up the former bedload were measured and defined. The results demonstrate a strong impact of lithology upon valley width, terrace distribution, and bedload provenance and therefore, upon landscape evolution processes. Channel down‐cutting through different lithologies creates variable patterns of resistance across catchments and continents. Particularly in post‐tectonic and non‐tectonic landscapes, the variation in resistance that arises from the exhumation of different rocks in channel longitudinal profiles can impact local base levels, initiating responses that can be propagated through channel networks. The rate at which that response is transmitted through channels is potentially amplified and/or mitigated by differences between the resistance of channel beds and bedload sediment loads. In the study reach, variation in lithologic resistance influences the prevalence of lateral and vertical processes, thus producing a spatial pattern of terraces that reflects rock type rather than climate, regional base level change, or hydrologic variability. Copyright © 2017 John Wiley & Sons, Ltd.     

Predicting gravel bed river response to environmental change: the strengths and limitations of a regime‐based approach

Rivers respond to environmental changes such as climate shifts, land use changes and the construction of hydro‐power dams in a variety of ways. Often there are multiple potential responses to any given change. Traditionally, potential stream channel response has been assessed using simple, qualitative frameworks based largely on professional judgement and field experience, or using some form of regime theory. Regime theory represents an attempt to use a physically based approach to predict the configuration of stable channels that can transport the imposed sediment supply with the available discharge. We review the development of regime theory, and then present a specific regime model that we have created as a stand‐alone computer program, called the UBC Regime Model (UBCRM). UBCRM differs from other regime models in that it constrains its predictions using a bank stability criterion, as well as a pattern stability criterion; it predicts both the stable channel cross‐sectional dimensions as well as the number of anabranches that the stream must have in order to establish a stable channel pattern. UBCRM also differs from other models in that it can be used in a stochastic modelling mode that translates uncertainty in the input variables into uncertainty in the predicted channel characteristics. However, since regime models are fundamentally based on the concept of grade, there are circumstances in which the model does not perform well. We explore the strengths and weaknesses of the UBCRM in this paper, and we attempt to illustrate how the UBCRM can be used to augment the existing qualitative frameworks, and to help guide professionals in their assessments. Copyright © 2016 John Wiley & Sons, Ltd.     

Morphodynamics of a bedrock‐alluvial meander bend that incises as it migrates outward: approximate solution of permanent form

In meandering rivers cut into bedrock, erosion across a channel cross‐section can be strongly asymmetric. At a meander apex, deep undercutting of the outer bank can result in the formation of a hanging cliff (which may drive hillslope failure), whereas the inner bank adjoins a slip‐off slope that connects to the hillslope itself. Here we propose a physically‐based model for predicting channel planform migration and incision, point bar and slip‐off slope formation, bedrock abrasion, the spatial distribution of alluvial cover, and adaptation of channel width in a mixed bedrock‐alluvial channel. We simplify the analysis by considering a numerical model of steady, uniform bend flow satisfying cyclic boundary conditions. Thus in our analysis, ‘sediment supply’, i.e. the total volume of alluvium in the system, is conserved. In our numerical simulations, the migration rate of the outer bank is a specified parameter. Our simulations demonstrate the existence of an approximate state of dynamic equilibrium corresponding to a near‐solution of permanent form in which a bend of constant curvature, width, cross‐sectional shape and alluvial cover distribution migrates diagonally downward at constant speed, leaving a bedrock equivalent of a point bar on the inside of the bend. Channel width is set internally by the processes of migration and incision. We find that equilibrium width increases with increasing sediment supply, but is insensitive to outer bank migration rate. The slope of the bedrock point bar varies inversely with both outer bank migration rate and sediment supply. Although the migration rate of the outer bank is externally imposed here, we discuss a model modification that would allow lateral side‐wall abrasion to be treated in a manner similar to the process of bedrock incision. Copyright © 2016 John Wiley & Sons, Ltd.     

Morphodynamic assessment of side channel systems using a simple one‐dimensional bifurcation model and a comparison with aerial images

Side channel construction is a common intervention applied to increase a river's conveyance capacity and to increase its ecological value. Past modelling efforts suggest two mechanisms affecting the morphodynamic change of a side channel: (1) a difference in channel slope between the side channel and the main channel and (2) bend flow just upstream of the bifurcation. The objective of this paper was to assess the conditions under which side channels generally aggrade or degrade and to assess the characteristic timescales of the associated morphological change. We use a one‐dimensional bifurcation model to predict the development of side channel systems and the characteristic timescale for a wide range of conditions. We then compare these results to multitemporal aerial images of four side channel systems. We consider the following mechanisms at the bifurcation to be important for side channel development: sediment diversion due to the bifurcation angle, sediment diversion due to the transverse bed slope, partitioning of suspended load, mixed sediment processes such as sorting at the bifurcation, bank erosion, deposition due to vegetation, and floodplain sedimentation. There are limitations to using a one‐dimensional numerical model as it can only account for these mechanisms in a parametrized manner, but the model reproduces general behaviour of the natural side channels until floodplain‐forming processes become important. The main result is a set of stability diagrams with key model parameters that can be used to assess the development of a side channel system and the associated timescale, which will aid in the future design and maintenance of side channel systems. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Stream channel erosion in a rapidly urbanizing region of the US–Mexico border: documenting the importance of channel hardpoints with Structure‐from‐Motion photogrammetry

Urbanization can lead to accelerated stream channel erosion, especially in areas experiencing rapid population growth, unregulated urban development on erodible soils, and variable enforcement of environmental regulations. A combination of field surveys and Structure‐from‐Motion (SfM) photogrammetry techniques was used to document spatial patterns in stream channel geometry in a rapidly urbanizing watershed, Los Laureles Canyon (LLCW), in Tijuana, Mexico. Ground‐based SfM photogrammetry was used to map channel dimensions with 1 to 2 cm vertical mean error for four stream reaches (100–300 m long) that were highly variable and difficult to survey with a differential GPS. Regional channel geometry curves for LLCW had statistically larger slopes and intercepts compared with regional curves developed for comparable, undisturbed reference channels. Cross‐sectional areas of channels downstream of hardpoints, such as concrete reaches or culverts, were up to 64 times greater than reference channels, with enlargement persisting, in some cases, up to 230 m downstream. Percentage impervious cover was not a good predictor of channel enlargement. Proximity to upstream hardpoint, and lack of riparian and bank vegetation paired with highly erodible bed and bank materials may account for the instability of the highly enlarged and unstable cross‐sections. Channel erosion due to urbanization accounts for approximately 25–40% of the total sediment budget for the watershed, and channel erosion downstream of hardpoints accounts for one‐third of all channel erosion. Channels downstream of hardpoints should be stabilized to prevent increased inputs of sediment to the Tijuana Estuary and local hazards near the structures, especially in areas with urban settlements near the stream channel. Copyright © 2017 John Wiley & Sons, Ltd.     

3D P-WAVE VELOCITY STRUCTURE AT THE NORTHEASTERN MARGIN OF ORDOS BLOCK

Using the 7 100 absolute first arrivals of P waves and 91 513 relative P arrival times of 726 events at the northeastern margin of the Ordos block since 2009, the 3D fine structure of P wave velocity within the depth of 15km in the crust was inverted by the double difference seismic tomography method. The results show that there exist obvious high-speed continuous bodies in the northwest of the study area, and their lateral areas increase gradually with depth, while the velocity of east and south is relatively low. The velocity inhomogeneity exists and differs at different depths. The lateral differences of velocity are related to seismicity and faults. The 5~15km depth profile shows that earthquakes tend to occur in the area with relatively high velocity or high speed transition zones, which to some extent reflects the fragility of regional crustal media and the strong differential movement of faults in vertical and horizontal directions where the crust body is easy to absorb and store strain energy and generate major earthquakes. A "Y"-shape low-velocity channel is present in the lower crust around Liangcheng, corresponding to the NW-trending Heilaoyao-Shahukou fault set, which may reveal the migration path of the Late Tertiary-Quaternary basalt eruption. The Helingeer M6.2 earthquake in 1976 was related to the formation of the locking section of the thermal welding in this area. The three-dimensional fine structure of P wave velocity presented in this paper provides intuitive seismological evidence for physical and chemical properties of crustal media and the deep tectonic environment of earthquake preparation.     

双稀疏字典和FISTA的地震数据去噪

地震数据的随机噪声去除是地震数据处理中的一项重要步骤,双稀疏字典提供了两层稀疏模型,比单层稀疏模型可以更好地去除噪声.该方法首先利用contourlet变换对地震数据进行稀疏表示,然后在contourlet域中使用快速迭代收缩阈值算法（fast iterative shrinkage-thresholding algorithm,FISTA）对初始字典系数进行更新,接着采用数据驱动紧标架（data-driven tight frame,DDTF）在contourlet域中得到DDTF字典并通过FISTA得到更新后的字典系数,最后通过DDTF字典和更新后的字典系数获得新的contourlet系数,并对新的contourlet系数进行硬阈值和contourlet反变换得到去噪后的数据.通过模拟数据和实际数据的实验证明：与固定基变换去噪方法相比,该方法可以自适应地对地震数据进行稀疏表示,在地震数据较为复杂时得到更高的信噪比;与字典学习去噪方法相比,该方法不仅拥有较快的去噪速度,而且克服了字典学习因为缺少先验约束造成瑕疵的缺点.