Fabric stability in oblique convergence and divergence

Forward modeling of transpression–transtension, assuming homogeneous strain and a direct relationship between finite strain axes and foliation–lineation in tectonites, investigates fields of stability of foliation and lineation orientations in oblique convergence and divergence. Vertical foliation–horizontal lineation (VF–HL) develop for angles of convergence–divergence between 0 and 20°. With increasing finite strain, this narrow window of stability is further reduced; lineation switches to vertical in transpression and foliation switches to horizontal in transtension. If a shear zone contains VF–HL, it either developed as a zone very close to pure wrenching, or recorded low finite strain. The stability of VF–HL at high strain and higher angles of convergence is enhanced by lateral extrusion of material along transpression zones. VF–HL may be stabilized in magmatic bodies that progressively intrude transtension zones, if the wrench component of deformation partitions within them. Alternatively, if these bodies are dike-like, cool fast, and do not record large deformation, they take up the extension component of transtension through anisotropic volume addition, leaving a larger component of wrench deformation in the country rocks; this effect stabilizes VF–HL effectively at low strain, but only marginally so at high strain.     

Strike-slip faults parallel to crustal spreading axes: data from Iceland and the Afar Depression

Slickenside studies in regions of crustal spreading such as Iceland and the Afar Depression, East Africa, reveal that a significant number of faults parallel and close to rift axes are strike-slip rather than normal. Therefore, the pattern of brittle deformation in these regions does not conform to the classic two-dimensional schemes of oceanic tectonics and pre-oceanic rifting. Dip-slip and strike-slip faulting presumably alternated along or in the vicinity of spreading axes, indicate a varying stress field and a combination of transverse and longitudinal movements. In Iceland, strike-slip faults parallel to rifts are observed both west and east of the rift system as well as in a median area between overlapping rifts; the mechanisms proposed for their origin include accommodation of oblique convergence or divergence of crustal sections due to variations of spreading directions along axis and the interaction of overlapping rifts. In the Afar Depression this kind of fault is recorded west of the rift of Asal and can be imputed to reflect an interaction among rifts in the vicinity of the Afar triple junction. Rift-parallel strike-slip faults cannot however be assumed to be a feature of all crustal spreading axes due to the peculiarity of the examined regions: both of them are hot-spot areas and the Afar Depression lies at a triple junction.     

The Role of Pressure as a Source of Gravity and Wormholes

The role of pressure as a source of gravity in general relativity is considered. In homogeneous, anisotropic cosmological models, the anisotropic pressure components along different coordinate axes are included as terms in the equations determining the accelerations along these axes; in this sense, this situation is similar to how the mass density is included in these equations. A change in the pressure component along one of the coordinates leads to an anisotropic change in the accelerations along all three coordinates. However, the pressure components are not present in the equation determining the initial conditions. On the other hand, with an inhomogeneous matter density distribution, the mass density and pressure play fundamentally different roles as sources of gravity. The consequences of this are considered.     

Axial directions of folds in rocks with linear/planar fabrics

This paper examines the geological implications of an analysis of constraints on the orientation of fold axes in orthotropic materials. We argue that rocks with penetrative linear and planar shape fabrics may have orthotropic (anisotropic) properties during deformation. Two forms of anisotropy (rheological and structural) could be potentially important in the control of fold axial directions. We discuss a model of deformation of rocks with linear/planar fabrics where, in a single deformation event, major fold axes need not be parallel to minor fold axes and neither need be perpendicular to the principal compressive stress direction. Geological and model examples of anistropic control on fold axial directions are given.     

青藏高原东北缘柴达木盆地路乐河地区新生代构造变形的古地磁证据

青藏高原东北缘是高原由西南向东北方向扩展的前缘位置，其新生代构造变形对揭示青藏高原隆升、扩展的过程与动力学机制具有重要的意义。柴达木盆地是青藏高原东北缘最大的新生代沉积盆地，发育巨厚的新生代地层，这些地层所记录的古地磁极旋转信息是定量约束柴达木盆地新生代以来构造变形发生的时间、方式与幅度的载体。本文以柴达木盆地北缘新生代地层出露良好、具有精确地层年代控制的路乐河剖面为研究对象，开展了古地磁极旋转研究，统计分析路乐河剖面24. 6~5. 2 Ma之间1477个可靠古地磁样品的特征剩磁方向(ChRM)，发现柴达木盆地北缘路乐河地区在24. 6~16. 4 Ma发生小幅度(不显著)的逆时针旋转，旋转角度约为8. 4°±6. 1°；16. 4~13. 9 Ma路乐河地区发生显著的顺时针旋转，旋转角度可达36. 1°±6. 0°；13. 9~5. 2 Ma 该地区未发生明显的构造旋转；5. 2 Ma以后路乐河地区逆时针旋转了~6°。结合柴达木盆地北缘区域构造变形的分析，我们提出柴达木盆地北缘路乐河地区在16. 4~13. 9 Ma 之间发生强烈的顺时针旋转构造变形(~36°)可能代表了盆地北缘中中新世遭受强烈的地壳差异缩短变形，从而成为高原最新形成的部分。     

Analysis of Long-Term Anisotropic Convergence in Drifts Excavated in Callovo-Oxfordian Claystone

The main purpose of this paper is to analyze the convergence measurements in drifts of the Underground Research Laboratory (URL) of the French National Radioactive Waste Management Agency (Andra), excavated in Callovo-Oxfordian claystone. These measurements show an anisotropic closure, which depends on the drift orientations with respect to the horizontal in situ stresses. This anisotropic character of the deformation is taken into account by assuming that the drifts section evolves following an elliptical shape. The characteristics of the deformed elliptical section are evaluated following the methodology proposed by Vu et al. (Rock Mech Rock Eng 46:231–246, 2013). Then, using the semi-empirical law proposed by Sulem et al. (Int J Rock Mech Min Sci Geomech Abstr 24:145–154, 1987), the convergence evolution is fitted independently for each axis of the ellipse. This method allows to distinguish two effects: the face advance and the time-dependent behavior of the ground. The results for the two drift orientations (along the major horizontal stress and perpendicular to it) show very close values for the parameters describing the time-dependent properties of the ground, the distance of influence of the face, and the extent of the decompressed zone around the drift. Finally, the model is validated by keeping these parameters as constants and by simulating the convergence data on a new drift. It is shown that with a period of about 40 days of convergence monitoring, the model can provide valuable insights for predictions of the convergence evolution in the long term.     

Coupled modeling of Excavation Damaged Zone in Boom clay: Strain localization in rock and distribution of contact pressure on the gallery’s lining

Around galleries excavated at depth in geological media, the creation of a damaged zone with significant irreversible deformation is generally unavoidable. In the case of a geological disposal system for high-level radioactive waste, the resulting change in the host rock properties in this damaged zone may potentially be important with respect to the long-term evolution and the performance of that system. In this context, predicting the extent of the so-called Excavation Damaged Zone (EDZ) and, possibly, the fractures’ network topology remains a challenge. This study is aimed to simulate numerically the extension of this zone at the large scale’s excavation, around the Connecting gallery (HADES URL, Mol, Belgium), in Boom clay host rock through analyzing the evolution of strain localization in shear bands mode. To realistically model the involved phenomena, the concrete lining is considered on the gallery wall highlighting its impacts on the evolution of convergence and EDZ around the gallery. The focus of the current paper is made on analyzing the coupled hydro-mechanical behavior of Boom clay host rock during and after the gallery excavation with respect to the evolution of localized shear bands around the gallery. This study is accompanied by the analysis of the contact mechanism on the interface between the clay massive and the lining. The obtained results reveal some interesting features regarding the contact phenomenon relatively to the evolution pattern of shear bands within the clay around the gallery. To assess the reliability of the proposed approach, a discussion on some in-situ observations during the gallery’s construction is also performed based on which a good agreement is found between the in-situ evidence and simulated results.     

Structural analysis of sheath folds and geochronology in the Cuonadong Dome,southern Tibet,China: New constraints on the timing of the South Tibetan detachment system and its relationship to North Himalayan Gneiss Domes

Structural analysis and new zircon U/Pb geochronology were performed in the Cuonadong Dome in order to test the hypothesis that the top‐to‐North deformation in the North Himalayan Gneiss Domes (NHGD) has been correlated with deformation along the South Tibetan detachment system (STDS). Structural analysis suggests that the X‐axes of sheath folds at the different flanks of the dome are consistent, and show average a trend and plunge of 350° and 18°. Many shear sense indicators are recorded in the ductile shear zone during the top‐to‐north deformation, including the augen structure, σ‐type and δ‐type rotating porphyroclast, S–C fabric and pressure shadow structure. The development of the sheath folds first provide a direct‐field structural evidence that represents a deeper ductile manifestation of the STDS, and further suggests that the top‐to‐north deformation in the NHGD have been correlated with deformation along the STDS. The emplacement of the synkinematic granites at ca. 32 Ma exposed at the core of the sheath folds is also interpreted as the initiation of slip along the STDS.     

大渡河丹巴水电站长大深埋软岩各向异性及宏观力学参数研究

研究了层状岩体各向异性以及宏观力学参数分析方法,提出基于大型岩土分析软件ADINA各向异性材料模块,采用分区分级宏观力学参数模拟分析方法,结合“代表单元体”理论思想,对丹巴水电站二云片岩进行了各向异性模拟和宏观力学参数分析,研究得到丹巴片岩不同片理倾角宏观力学参数以及各向异性特征变化规律。研究结果显示,丹巴二云片岩尺度增加,平行片理向、垂直片理向变形模量 、 收敛幅度不同步, 变化幅度随倾角增大而增大,收敛速率更高,水平向变形模量 趋势较平缓,泊松比的变化规律恰好与之相反;随片理面倾角增大,丹巴二云片岩力学参数趋于极限值规律明显。     

Geodynamic evolution of the South Variscan Iberian Suture as recorded by mineral transformations

New structural, petrographic, mineralogical and geochemical data from the Beja-Acebuches Ophiolite Complex (BAOC) are presented, and reviewed together with data published elsewhere. The new data obtained shed light on questions such as: 1) the relative importance of the obduction event; 2) its geological record in the deep levels of BAOC; 3) the nature and intensity of the Variscan metamorphism and deformation during subsequent continental (arc) collision; 4) the age relationships between BAOC and the Beja Igneous Complex; and 5) by means of numerical modelling, the thermal metamorphism of the Ossa-Morena autochthonous terranes induced by the ophiolite obduction. The emerging picture is that of a fairly simple overall geological evolution for BAOC, seamlessly integrated within the evolution of the southern branch of the Iberian Variscides. Obduction of BAOC is a relatively minor early event in the general NE–SW convergence that gave rise to the orogen as seen regionally and is recorded by an anisotropic, high-temperature, metamorphic fabric at the gabbro levels and by subtle features of the chemical composition of primary minerals at the underlying peridotite level; it caused chilling of the obducted ophiolitic slab and no significant metamorphism on the autochtonous rocks of the Ossa-Morena Zone. BAOC underwent most of its deformation and (amphibolite facies) metamorphism during a later collisional event, that took place as the most primitive rocks of the Beja Igneous complex were being intruded, and whose waning stages are responsible for extensive serpentinisation of peridotites and for important aquocarbonic fluid discharges along the semibrittle–brittle shear zones meanwhile developed.     

Interaction of chemical and physical processes during deformation at fluid-present conditions: a case study from an anorthosite–leucogabbro deformed at amphibolite facies conditions

We present microstructural and chemical analyses of chemically zoned and recrystallized plagioclase grains in variably strained samples of a naturally deformed anorthosite–leucogabbro, southern West Greenland. The recorded microstructures formed in the presence of fluids at mid-crustal conditions (620–640 °C, 7.4–8.6 kbar). Recrystallized plagioclase grains (average grain size 342 μm) with a random crystallographic orientation are volumetrically dominant in high-strain areas. They are characterized by asymmetric chemical zoning (An₈₀ cores and An₆₄ rims) that are directly associated with areas exhibiting high amphibole content and phase mixing. Analyses of zoning indicate anisotropic behaviour of bytownite plagioclase with a preferred replacement in the $ \left\langle {0 10} \right\rangle $ direction and along the (001) plane. In areas of high finite strain, recrystallization of plagioclase dominantly occurred by bulging recrystallization and is intimately linked to the chemical zoning. The lack of CPO as well as the developed asymmetric zoning can be explained by the activity of grain boundary sliding accommodated by dissolution and precipitation creep (DPC). In low-strain domains, grain size is on average larger and the rim distribution is not related to the inferred stress axes indicating chemically induced grain replacement instead of stress-related DPC. We suggest that during deformation, in high-strain areas, pre-existing phase mixture and stress induced DPC-caused grain rotations that allowed a deformation-enhanced heterogeneous fluid influx. This resulted in local plagioclase replacement through interface-coupled dissolution and precipitation and chemically induced grain boundary migration, accompanied by bulging recrystallization, along with neocrystallization of other phases. This study illustrates a strong interaction and feedback between physical and chemical processes where the amount of stress and fluids dictates the dominant active process. The interaction is a cause of deformation and external fluid infiltration with a result of strain localization and chemical re-equilibration at amphibolite facies conditions.     

Prediction of fracture trajectory in anisotropic rocks using modified maximum tangential stress criterion

A new criterion to predict crack propagation trajectory in anisotropic rocks with incorporating the concept of T-stress in formulating stress field near the crack tip was developed. The developed criterion along with enrichment functions and interaction integral in the extended finite element method (XFEM) framework made a sophisticated tool in modeling fracturing process in anisotropic media. Numerical results indicated that stress intensity factors considerably depend on orientation of anisotropy axes and ratio of the elastic modulus. The proposed formulation for anisotropic media provides a more accurate prediction of crack propagation trajectory compared with conventional methods, especially in mixed mode conditions.     

Transpressive Structures in the Ghadir Shear Belt,Eastern Desert,Egypt: Evidence for Partitioning of Oblique Convergence in the Arabian‐Nubian Shield during Gondwana Agglutination

Transpressional deformation has played an important role in the late Neoproterozoic evolution of the ArabianNubian Shield including the Central Eastern Desert of Egypt. The Ghadir Shear Belt is a 35 km-long, NW-oriented brittleductile shear zone that underwent overall sinistral transpression during the Late Neoproterozoic. Within this shear belt, strain is highly partitioned into shortening, oblique, extensional and strike-slip structures at multiple scales. Moreover, strain partitioning is heterogeneous along-strike giving rise to three distinct structural domains. In the East Ghadir and Ambaut shear belts, the strain is pure-shear dominated whereas the narrow sectors parallel to the shear walls in the West Ghadir Shear Zone are simple-shear dominated. These domains are comparable to splay-dominated and thrust-dominated strike-slip shear zones. The kinematic transition along the Ghadir shear belt is consistent with separate strike-slip and thrustsense shear zones. The earlier fabric(S1), is locally recognized in low strain areas and SW-ward thrusts. S2 is associated with a shallowly plunging stretching lineation(L2), and defines ～NW-SE major upright macroscopic folds in the East Ghadir shear belt. F2 folds are superimposed by ～NNW–SSE tight-minor and major F3 folds that are kinematically compatible with sinistral transpressional deformation along the West Ghadir Shear Zone and may represent strain partitioning during deformation. F2 and F3 folds are superimposed by ENE–WSW gentle F4 folds in the Ambaut shear belt. The sub-parallelism of F3 and F4 fold axes with the shear zones may have resulted from strain partitioning associated with simple shear deformation along narrow mylonite zones and pure shear-dominant deformation in fold zones. Dextral ENEstriking shear zones were subsequently active at ca. 595 Ma, coeval with sinistral shearing along NW-to NNW-striking shear zones. The occurrence of upright folds and folds with vertical axes suggests that transpression plays a significant role in the tectonic evolution of the Ghadir shear belt. Oblique convergence may have been provoked by the buckling of the Hafafit gneiss-cored domes and relative rotations between its segments. Upright folds, fold with vertical axes and sinistral strike-slip shear zones developed in response to strain partitioning. The West Ghadir Shear Zone contains thrusts and strikeslip shear zones that resulted from lateral escape tectonics associated with lateral imbrication and transpression in response to oblique squeezing of the Arabian-Nubian Shield during agglutination of East and West Gondwana.     

The initiation and development of metamorphic foliation in the Otago Schist, Part 1: competitive oriented growth of white mica

The 3D shape, size and orientation data for white mica grains sampled along two transects of increasing metamorphic grade in the Otago Schist, New Zealand, reveal that metamorphic foliation, as defined by mica shape‐preferred orientation (SPO), developed rapidly at sub‐greenschist facies conditions early in the deformation history. The onset of penetrative strain metamorphism is marked by the rapid elimination of poorly oriented large clastic mica in favour of numerous new smaller grains of contrasting composition, higher aspect ratios and a strong preferred orientation. The metamorphic mica is blade shaped with long axes defining the linear aspect of the foliation and intermediate axes a partial girdle about the lineation. Once initiated, foliation progressively intensified by an increase in the aspect ratio, size and alignment of grains, although highest grade samples within the chlorite zone record a decrease in aspect ratio and reduction in SPO strength despite continued increase in grain size. These trends are interpreted in terms of progressive competitive anisotropic growth of blade‐shaped grains so that the fastest growth directions and blade lengths tend to parallel the extension direction during deformation. The competitive nature of mica growth is indicated by the progressive increase in size and resultant decrease in number of metamorphic mica with increasing grade, from c. 1000 relatively small mica grains per square millimetre of thin section at lower grades, to c. 100 relatively large grains per square millimetre in higher grade samples. Reversal of SPO intensity and grain aspect ratio trends in higher grade samples may reflect a reduction in the strain rate or reduction in the deviatoric component of the stress field.     

隧洞围岩收敛损失位移的求取方法及应用

隧洞围岩变形监测断面设置与开挖施工工序上的先后导致了损失位移问题,它对于判断围岩的变形状态、评价其稳定性具有重要的意义。复杂多变的围岩条件使基于收敛曲线的累积位移和增量位移参数反演方法的应用变得非常谨慎。针对这些问题,通过对一致性问题的讨论后认为,虽然变形演化规律上的一致性尚存在问题,但三维计算和现场围岩的总位移之间的一致性是可以保证的。以此问题为前提,以Hoek经验公式的通式为基础,依据规范规定,基于现场监测收敛位移曲线,建立了损失位移的求取步骤,并将其在锦屏II级水电站辅助洞A、B围岩变形规律分析中进行了应用,得出了两个监测断面各测线的最终监测位移收敛值、损失位移、总收敛位移和相对收敛值等,拟合而得的全收敛曲线与监测曲线相吻合,能够很好地描述现场围岩随空间效应变化的演化规律,为围岩稳定性评价和基于全位移的参数反演分析奠定了基础。     

定轴排水剪切试验中各向异性砂土的力学响应

南京云母砂主要由片状颗粒组成,具有较为明显的各向异性特性。采用香港科技大学的自动控制空心圆柱扭剪仪进行了应力主轴方向角分别为10°、30°、60°、80°的密实云母砂定轴剪切试验,研究云母砂的初始各向异性规律。为研究中主应力的影响,又对同一密度的试样在中主应力系数b＝0.1、0.5、1.0三种不同情况下进行定轴剪切试验。试验结果表明,在大主应力方向角从0°～90°的变化过程中,强度先减小后增大,且随着b值的增加强度增大;变形表现出明显的各向异性,随着b值的增加体变增大     

Convergence analysis of fixed stress split iterative scheme for anisotropic poroelasticity with tensor Biot parameter

We perform a convergence analysis of the fixed stress split iterative scheme for the Biot system modeling coupled flow and deformation in anisotropic poroelastic media with tensor Biot parameter. The fixed stress split iterative scheme solves the flow subproblem with all components of the stress tensor frozen using a multipoint flux mixed finite element method, followed by the poromechanics subproblem using a conforming Galerkin method in every coupling iteration at each time step. The coupling iterations are repeated until convergence and Backward Euler is employed for time marching. The convergence analysis is based on studying the equations satisfied by the difference of iterates to show that the fixed stress split iterative scheme for anisotropic poroelasticity with Biot tensor is contractive. We also demonstrate that the scheme is numerically convergent using the classical Mandel’s problem solution for transverse isotropy.     

Present-day crustal deformation along the Philippine Fault in Luzon,Philippines

The Philippine Fault results from the oblique convergence between the Philippine Sea Plate and the Sunda Block/Eurasian Plate. The fault exhibits left-lateral slip and transects the Philippine archipelago from the northwest corner of Luzon to the southeast end of Mindanao for about 1200 km. To better understand fault slip behavior along the Philippine Fault, eight GPS surveys were conducted from 1996 to 2008 in the Luzon region. We combine the 12-yr survey-mode GPS data in the Luzon region and continuous GPS data in Taiwan, along with additional 15 International GNSS Service sites in the Asia-Pacific region, and use the GAMIT/GLOBK software to calculate site coordinates. We then estimate the site velocity from position time series by linear regression. Our results show that the horizontal velocities with respect to the Sunda Block gradually decrease from north to south along the western Luzon at rates of 85–49 mm/yr in the west–northwest direction. This feature also implies a southward decrease of convergence rate along the Manila Trench. Significant internal deformation is observed near the Philippine Fault. Using a two dimensional elastic dislocation model and GPS velocities, we invert for fault geometries and back-slip rates of the Philippine Fault. The results indicate that the back-slip rates on the Philippine Fault increase from north to south, with the rates of 22, 37 and 40 mm/yr, respectively, on the northern, central, and southern segments. The inferred long-term fault slip rates of 24–40 mm/yr are very close to back-slip rates on locked fault segments, suggesting the Philippine Fault is fully locked. The stress tensor inversions from earthquake focal mechanisms indicate a transpressional regime in the Luzon area. Directions of σ₁ axes and maximum horizontal compressive axes are between 90° and 110°, consistent with major tectonic features in the Philippines. The high angle between σ₁ axes and the Philippine Fault in central Luzon suggests a weak fault zone possibly associated with fluid pressure.     

石榴石Lu-Hf、原位独居石U-Pb定年对多期变形的时代制约——以北祁连托勒牧场为例

多期变质变形事件的精确年代限定是造山构造年代学研究的热点问题之一。本文尝试运用面理弯切轴测量技术,结合石榴石Lu-Hf和原位独居石U-Pb定年,厘定北祁连托勒牧场地区石榴石和斜长石斑晶记录的两期构造变形事件:石榴石斑晶生长记录的早期构造变形事件年代为512.3±2.7Ma;斜长石斑晶生长记录的晚期构造变形事件年代不早于481.0±2.3Ma,并推断该期构造变形水平挤压主应力方向为北东-南西。斜长石斑晶内未发现独居石,用于年代学测试的独居石颗粒均位于斜长石斑晶外基质中。显微构造分析认为,独居石生长不早于斜长石斑晶。481.0±2.3Ma的独居石U-Pb年龄,应为斜长石斑晶所记录构造变形的时代下限。结合前人锆石U-Pb定年和Hf同位素研究结果分析认为,获得的512.3±2.7Ma石榴石-全岩Lu-Hf等时线年龄,代表了祁连洋俯冲过程中石榴石的生长时间,后期变质变形作用未对石榴石的Lu-Hf同位素体系产生明显影响。结合显微构造分析,石榴石Lu-Hf定年可为早期构造变形提供有效年代学制约。     

Constitutive modeling of principal stress rotation by considering inherent and induced anisotropy of soils

In order to simulate the soil response during principal stress rotation, anisotropic unified hardening (UH) model is developed within the framework of elastoplastic theory. Without introducing any additional mechanism to display the role of stress rotation specifically, this model achieves the simulation by considering the material anisotropy. The effect of inherent anisotropy is reflected using the anisotropic transformed stress method, but a new formula for the stress mapping is adopted to keep the mean stress unchanged. Analysis indicates that from the view of the transformed stress tensor, the anisotropic soil is subjected to loading during pure rotation of principal stress axes, so that plastic strains can be calculated. To represent the induced anisotropy, a fabric evolution law is proposed based on laboratory and numerical test results. At the critical state, the fabric tensor reaches a stable value determined by the stress state, while the critical state line is unique in the plane of void ratio versus mean stress. The anisotropic UH model has concise formulation and explicit elastoplastic flexibility matrix and can provide reasonable predictions for the deformation of anisotropic soils when principal stresses rotate.