Crustal structure of the Lofoten–Vesterålen continental margin, off Norway

By compiling wide-angle seismic velocity profiles along the 400-km-long Lofoten–Vesterålen continental margin off Norway, and integrating them with an extensive seismic reflection data set and crustal-scale two-dimensional gravity modelling, we outline the crustal margin structure. The structure is illustrated by across-margin regional transects and by contour maps of depth to Moho, thickness of the crystalline crust, and thickness of the 7+ km/s lower crustal body. The data reveal a normal thickness oceanic crust seaward of anomaly 23 and an increase in thickness towards the continent–ocean boundary associated with breakup magmatism. The southern boundary of the Lofoten–Vesterålen margin, the Bivrost Fracture Zone and its landward prolongation, appears as a major across-margin magmatic and structural crustal feature that governed the evolution of the margin. In particular, a steeply dipping and relatively narrow, 10–40-km-wide, Moho-gradient zone exists within a continent–ocean transition, which decreases in width northward along the Lofoten–Vesterålen margin. To the south, the zone continues along the Vøring margin, however it is offset 70–80 km to the northwest along the Bivrost Fracture Zone/Lineament. Here, the Moho-gradient zone corresponds to a distinct, 25-km-wide, zone of rapid landward increase in crustal thickness that defines the transition between the Lofoten platform and the Vøring Basin. The continental crust on the Lofoten–Vesterålen margin reaches a thickness of 26 km and appears to have experienced only moderate extension, contrasting with the greatly extended crust in the Vøring Basin farther south. There are also distinct differences between the Lofoten and Vesterålen margin segments as revealed by changes in structural style and crustal thickness as well as in the extent of elongate potential-field anomalies. These changes may be related to transfer zones. Gravity modelling shows that the prominent belt of shelf-edge gravity anomalies results from a shallow basement structural relief, while the elongate Lofoten Islands belt requires increased lower crustal densities along the entire area of crustal thinning beneath the islands. Furthermore, gravity modelling offers a robust diagnostic tool for the existence of the lower crustal body. From modelling results and previous studies on- and off-shore mid-Norway, we postulate that the development of a core complex in the middle to lower crust in the Lofoten Islands region, which has been exhumed along detachments during large-scale extension, brought high-grade, lower crustal rocks, possibly including accreted decompressional melts, to shallower levels.     

Investigation into the effect of fault properties on wave transmission

The estimation of wave transmission across the fractured rock masses is of great importance for rock engineers to assess the stability of rock slopes in open pit mines. Presence of fault, as a major discontinuity, in the jointed rock mass can significantly impact on the peak particle velocity and transmission of blast waves, particularly where a fault contains a thick infilling with weak mechanical properties. This paper aims to study the effect of fault properties on transmission of blasting waves using the distinct element method. First, a validation study was carried out on the wave transmission across a single joint and different rock mediums through undertaking a comparative study against analytical models. Then, the transmission of blast wave across a fault with thick infilling in the Golgohar iron mine, Iran, was numerically studied, and the results were compared with the field measurements. The blast wave was numerically simulated using a hybrid finite element and finite difference code which then the outcome was used as the input for the distinct element method analysis. The measured uplift of hanging wall, as a result of wave transmission across the fault, in the numerical model agrees well with the recorded field measurement. Finally, the validated numerical model was used to study the effect of fault properties on wave transmission. It was found that the fault inclination angle is the most effective parameter on the peak particle velocity and uplift. Copyright © 2017 John Wiley & Sons, Ltd.     

Constitutive and numerical framework for modeling joints and faults in rock

A three‐dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre‐existing gaps or the formation of gaps. Representative stress–strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright © 2015 John Wiley & Sons, Ltd.     

Inverse modelling using PS‐InSAR data for improved land subsidence simulation in Las Vegas Valley,Nevada

Las Vegas Valley has had a long history of groundwater development and subsequent surface deformation. InSAR interferograms have revealed detailed and complex spatial patterns of subsidence in the Las Vegas Valley area that do not coincide with major pumping regions. This research represents the first effort to use high spatial and temporal resolution subsidence observations from InSAR and hydraulic head data to inversely calibrate transmissivities (T), elastic and inelastic skeletal storage coefficients (S_ke and S_kv) of the developed‐zone aquifer and conductance (CR) of the basin‐fill faults for the entire Las Vegas basin. The results indicate that the subsidence observations from InSAR are extremely beneficial for accurately quantifying hydraulic parameters, and the model calibration results are far more accurate than when using only groundwater levels as observations, and just a limited number of subsidence observations. The discrepancy between distributions of pumping and greatest levels of subsidence is found to be attributed to spatial variations in clay thickness. The Eglington fault separates thicker interbeds to the northwest from thinner interbeds to the southeast and the fault may act as a groundwater‐flow barrier and/or subsidence boundary, although the influence of the groundwater barrier to this area is found to be insignificant. Copyright © 2016 John Wiley & Sons, Ltd.     

Extended finite element framework for fault rupture dynamics including bulk plasticity

We present an explicit extended finite element framework for fault rupture dynamics accommodating bulk plasticity near the fault. The technique is more robust than the standard split‐node method because it can accommodate a fault propagating freely through the interior of finite elements. To fully exploit the explicit algorithmic framework, we perform mass lumping on the enriched finite elements that preserve the kinetic energy of the rigid body and enrichment modes. We show that with this technique, the extended FE solution reproduces the standard split‐node solution, but with the added advantage that it can also accommodate randomly propagating faults. We use different elastoplastic constitutive models appropriate for geomaterials, including the Mohr–Coulomb, Drucker–Prager, modified Cam‐Clay, and a conical plasticity model with a compression cap, to capture off‐fault bulk plasticity. More specifically, the cap model adds robustness to the framework because it can accommodate various modes of deformation, including compaction, dilatation, and shearing. Copyright © 2013 John Wiley & Sons, Ltd.     

基于卫星遥感影像和DEM融合的地质断裂带研究

应用遥感影像进行地表地质分析,具有数据量大、耗时短、费用低廉等优势,可以快速获取大范围区域的相关地质信息,因此,近年来,卫星遥感影像广泛应用于地质分析工作中。但是由于卫星影像在不同时期内的特征不同,且有大量不同的数据源,在不同影像上获取地质特征困难,因此,如何利用影像处技术增强影像中的地质特征,成为遥感影像地质判读的研究重点。地质断裂带调查一般借助各种直接与间接的调查方法以确认断裂带的存在与否和位置,了解断裂带的活动性,并提供点状或剖面数据材料。卫星影像在追踪研究贯穿地表且可为影像所解析的断裂带上,提供真实的影像证据,有助于直接分析影像中断裂带走向和延展。利用遥感影像处理计算,将卫星影像数据叠加到DEM上编制三维地形图,并利用计算机仿真,通过不同方向的投射光源,找出最适合的观察角及最能凸显出重点区域线状地类的状态。对于地质现象的分析及解释上,可提供较平面影像细致、准确的参考依据。     

Weakness and mechanical anisotropy of phyllosilicate-rich cataclasites developed after mylonites of a low-angle normal fault (Simplon Line,Western Alps)

The Simplon Fault Zone is a late-collisional low-angle normal fault (LANF) of the Western Alps. The hanging wall shows evidence of brittle deformation only, while the footwall is characterized by a c. 1 km-thick shear zone (the Simplon Fault Zone), which continuously evolved, during exhumation and cooling, from amphibolite facies conditions to brittle-cataclastic deformations. Due to progressive localization of the active section of the shear zone, the thermal-rheological evolution of the footwall resulted in a layered structure, with higher temperature mylonites preserved at the periphery of the shear zone, and cataclasites occurring at the core (indicated as the Simplon Line). In order to investigate the weakness of the Simplon Line, we studied the evolution of brittle/cataclastic fault rocks, from nucleation to the most mature ones. Cataclasites are superposed on greenschist facies mylonites, and their nucleation can be studied at the periphery of the brittle fault zone. This is characterized by fractures, micro-faults and foliated ultracataclasite seams that develop along the mylonitic SCC′ fabric, exploiting the weak phases mainly represented by muscovite and chlorite. Approaching the fault core, both the thickness and frequency of cataclasite horizons increase, and, as their thickness increases, they become less and less foliated. The fault core itself is represented by a thicker non-foliated cataclasite horizon. No Andersonian faults or fractures can be found in the footwall damage zone and core zone, whilst they are present in the hanging wall and in the footwall further from the fault. Applying a stress model based on slip tendency, we have been able to calculate that the friction coefficient of the Simplon Line cataclasites was <0.25, hence this fault zone is absolutely weak. In contrast with other fault zones, the weakening effect of fluids was of secondary importance, since they accessed the fault zone only after an interconnected fracture network developed exploiting the cataclasite network.     

Smoothing and re-roughening processes: The geometric evolution of a single fault zone

The geometry of a fault zone exerts a major control on earthquake rupture processes and source parameters. Observations previously compiled from multiple faults suggest that fault surface shape evolves with displacement, but the specific processes driving the evolution of fault geometry within a single fault zone are not well understood. Here, we characterize the deformation history and geometry of an extraordinarily well-exposed fault using maps of cross-sectional exposures constructed with the Structure from Motion photogrammetric method. The La Quinta Fault, located in southern California, experienced at least three phases of deformation. Multiple layers of ultracataclasite formed during the most recent phase. Crosscutting relations between the layers define the evolution of the structures and demonstrate that new layers formed successively during the deformation history. Wear processes such as grain plucking from one layer into a younger layer and truncation of asperities at layer edges indicate that the layers were slip zones and the contacts between them slip surfaces. Slip surfaces that were not reactivated or modified after they were abandoned exhibit self-affine geometry, preserving the fault roughness from different stages of faulting. Roughness varies little between surfaces, except the last slip zone to form in the fault, which is the smoothest. This layer contains a distinct mineral assemblage, indicating that the composition of the fault rock exerts a control on roughness. In contrast, the similar roughness of the older slip zones, which have comparable mineralogy but clearly crosscut one another, suggests that as the fault matured the roughness of the active slip surface stayed approximately constant. Wear processes affected these layers, so for roughness to stay constant the roughening and smoothing effects of fault slip must have been approximately balanced. These observations suggest fault surface evolution occurs by nucleation of new surfaces and wear by competing smoothing and re-roughening processes.     

广西油麻坡钨钼矿床成岩成矿年代学研究及其地质意义

广西博白_岑溪断裂带位于钦杭成矿带南段,既是一条长期活动的地体边界断裂带,也是一条岩浆活动频繁的W_Mo_Sn_Pb_Zn_Sb_Au_Ag多金属成矿带。文章以该带中的广西博白县油麻坡矽卡岩型钨钼矿床为研究对象,在详细的岩性鉴别和划分的基础上,对矿区内花岗岩类进行了单颗粒锆石LA_ICP_MS锆石U_Pb测年,获得深灰色花岗斑岩的等时线年龄为(479.7±3.9)Ma,灰色黑云母花岗闪长岩为(109.7±1.1)Ma,浅灰色细粒白云母花岗岩为(103.3±1.2)Ma。辉钼矿Re_Os测年获得加权平均年龄为(97.4±1.9)Ma。由此可见,油麻坡岩体为一个复式岩体,由早奥陶世(加里东期)的花岗斑岩和早白垩世(燕山晚期)的黑云母花岗闪长岩和细粒白云母花岗岩组成,其中,黑云母花岗闪长岩构成油麻坡岩体的主体。钨钼矿形成于早白垩世晚期(即燕山晚期),与黑云母花岗闪长岩和白云母花岗岩密切相关。这些高精度测年数据的获得,进一步表明博白断裂带是一条长期活动的断裂,其在加里东期已经活动,并伴有岩浆岩的侵位。燕山晚期,该断裂重新复活,并导致大量的岩浆活动和一定强度的W_Mo_Sn_Pb_Zn_Sb_Au_Ag成矿作用,形成由多个大_中型矿床组成的多金属成矿带。博白_岑溪成矿带属于华南100~80 Ma大规模成矿的一部分,成矿背景与华南地区白垩纪的地壳伸展、钦杭成矿带的再次裂陷有关。     

依兰-伊通断裂新构造活动规律分析

作为郯庐断裂带北段主干的依兰-伊通断裂, 其新构造活动性与活动规律仍然存在不同的认识.本次工作通过详细的野外调查, 发现该断裂内活断层广泛存在, 由东、西两支北东走向的主干活断层构成, 沿着古近纪地堑边界断层发育.这些活断层主要呈破碎型结构, 多为逆右行平移活动.通过对这些活断层一系列实测擦痕反演应力场, 显示它们多是在东西向挤压中活动的, 而现今应力场转变为北东东-南西西向区域性挤压.依据本次野外观察与¹⁴ C定年, 并结合前人定年结果与近代地震分布, 表明依兰-伊通西支活断层的最新活动时代为全新世与晚更新世相间, 而东支活断层的最新活动时代主要为早-中更新世.依兰-伊通断裂内活断层显示了明显的差异性活动, 表现为西支的活动强度明显大于东支, 西支的最新活动时代皆晚于东支, 沿走向上活动性强、弱相间与最新活动时代不断变化, 以及近代地震活动不均一分布.它们沿走向上的分段性、差异性活动主要是因为被一系列北西向断层切断所致.