Mean flow and turbulence structure over fixed, two-dimensional dunes: implications for sediment transport and bedform stability

Detailed measurements of flow velocity and its turbulent fluctuation were obtained over fixed, two-dimensional dunes in a laboratory channel. Laser Doppler anemometry was used to measure the downstream and vertical components of velocity at more than 1800 points over one dune wavelength. The density of the sampling grid allowed construction of a unique set of contour maps for all mean flow and turbulence parameters, which are assessed using higher moment measures and quadrant analysis. These flow field maps illustrate that: (1) the time-averaged downstream and vertical velocities agree well with previous studies of quasi-equilibrium flow over fixed and mobile bedforms and show a remarkable symmetry from crest to crest; (2) the maximum root-mean-square (RMS) of the downstream velocity values occur at and just downstream of flow reattachment and within the flow separation cell; (3) the maximum vertical RMS values occur within and above the zone of flow separation along the shear layer and this zone advects and diffuses downstream, extending almost to the next crest; (4) positive downstream skewness values occur within the separation cell, whereas positive vertical skewness values are restricted to the shear layer; (5) the highest Reynolds stresses are located within the zone of flow separation and along the shear layer; (6) high-magnitude, high-frequency quadrant-2 events (‘ejections’) are concentrated along the shear layer (Kelvin-Helmholtz instabilities) and dominate the contribution to the local Reynolds stress; and (7) high-magnitude, high-frequency quadrant-4 events occur bounding the separation zone, near reattachment and close to the dune crest, and are significant contributors to the local Reynolds stress at each location. These data demonstrate that the turbulence structure associated with dunes is controlled intrinsically by the formation, magnitude and downstream extent of the flow separation zone and resultant shear layer. Furthermore, the origin of dune-related macroturbulence lies in the dynamics of the shear layer rather than classical turbulent boundary layer bursting. The fluid dynamic distinction between dunes and ripples is reasoned to be linked to the velocity differential across the shear layer and hence the magnitude of the Kelvin-Helmholtz instabilities, which are both greater for dunes than ripples. These instabilities control the local flow and turbulence structure and dictate the modes of sediment entrainment and their transport rates.     

Numerical modeling of the flow and transport of radionuclides in heterogeneous porous media

This paper is concerned with numerical methods for the modeling of flow and transport of contaminant in porous media. The numerical methods feature the mixed finite element method over triangles as a solver to the Darcy flow equation and a conservative finite volume scheme for the concentration equation. The convective term is approximated with a Godunov scheme over the dual finite volume mesh, whereas the diffusion–dispersion term is discretized by piecewise linear conforming triangular finite elements. It is shown that the scheme satisfies a discrete maximum principle. Numerical examples demonstrate the effectiveness of the methodology for a coupled system that includes an elliptic equation and a diffusion–convection–reaction equation arising when modeling flow and transport in heterogeneous porous media. The proposed scheme is robust, conservative, efficient, and stable, as confirmed by numerical simulations.      

超高压变质岩的塑性流变：显微构造和变形机制

超高压变质岩是大陆深俯冲作用的产物。超高压变质岩在深俯,中和快速折返过程中,经历了长距离地构造搬运和构造力的作用。其构造变形主要集中在韧性剪切带中,并发生强烈地塑性流变。研究超高压变质构造岩的显微构造及其变形机制对于深入了解大陆壳岩石在深俯;中过程中的流变学行为有十分重要的意义。山东仰口的超高压韧性剪切带中榴辉岩质和花岗质糜棱岩记录了超高压变形的历史。在超高压条件下的稳定矿物绿辉石、多硅白云母、兰晶石和钾长石具有不规则波状消光、亚晶界、核幔构造和动态重结晶等显微构造特征,TEM研究揭示了大量的位错构造,表明位错蠕变是其主要的变形机制。在花岗质糜棱岩中,金红石在刚性矿物的压力影中沉积,细粒的石榴石条带平行片理延伸,都说明超高压变形过程中有流体存在,流体助力的物质扩散迁移是又一个重要的变形机制、依据现有的流变学定律估算的流变应力应该在几十兆帕以上。     

新疆塔什库尔干温泉地区花岗岩体侵入与新生代构造变形:对东北帕米尔中新世构造演化的启示

帕米尔高原是受到印度-亚洲大陆碰撞、持续汇聚影响最显著的地区之一,以强烈地壳增厚和缩短、大量断裂和片麻岩穹窿的形成以及广泛的陆内岩浆活动为特征。以往有关帕米尔陆内岩浆岩的讨论多集中于对其地球化学成分及其所指示的构造背景的研究方面,而对岩浆形成与大型新生代构造之间的联系研究较少。本文通过对东北帕米尔塔什库尔干温泉地区新生代花岗岩及其围岩中锆石U-Pb和Ar-Ar年代学研究,结合该地区新生代构造变形分析,揭示岩浆侵入与区域构造变形之间的关系。U-Pb及Ar-Ar测年结果显示温泉地区花岗岩脉形成时代为中新世(11.8±0.2Ma),其及其围岩在10.8±0.1Ma冷却到300℃左右;中新世花岗岩脉中继承锆石及围岩片岩中碎屑锆石U-Pb的年龄分布特征迥异,它们分别具有类似松潘-甘孜地体以及中帕米尔地体的物质来源。花岗岩内部几乎无变形,围岩片岩变形主要体现为近E-W向的伸展构造,反映其形成于拉张为主的构造应力环境。结合区域构造背景,推测温泉中新世花岗岩脉为公格尔山、慕士塔格峰构造单元岩石部分熔融产物,是在东北帕米尔地壳从挤压增厚向局部伸展转换的过程中形成的,此时,公格尔伸展断裂系可能已经开始发育;此后,直到6~4Ma,公格尔伸展断裂系开始快速运动,与之相伴公格尔-慕士塔格片麻岩穹窿快速折返。

     

Subcritical compaction and yielding of granular quartz sand

Cylindrical samples of water-saturated, initially loose, St. Peter quartz sand were consolidated using triaxial deformation apparatus at room temperature, constant fluid pressure (12.5 MPa), and elevated confining pressures (up to 262.5 MPa). The samples were deformed along four loading paths: (1) hydrostatic stressing tests in which confining pressure was monotonically increased; (2) hydrostatic stress cycling similar to (1) except that effective pressure was periodically decreased to initial conditions; (3) triaxial deformation at constant effective pressure in which differential stress was applied after raising effective pressure to an elevated level; and (4) triaxial stress cycling similar to (3) except that the axial differential stress was periodically decreased to zero. Hydrostatic stressing at a constant rate results in a complex nonlinear consolidation response. At low pressures, large strains occur without significant acoustic emission (AE) activity. With increased pressure, the stress versus strain curve becomes quasi-linear with a corresponding nonlinear increase in AE rates. At elevated pressures, macroscopic yielding is marked by the onset of large strains, high AE rates, and significant grain failure. Stress cycling experiments show that measurable inelastic strain occurs at all stages of hydrostatic loading. The reload portions of stress cycles are characterized by a poro-elastic response and lower AE rates than during constant rate hydrostatic stressing. As the stress nears and exceeds the level that was applied during previous loading cycles, strain and AE rates increase in a manner consistent with yielding. Triaxial stressing cycles achieve greater consolidation and AE rates than hydrostatic loading at similar mean stress levels. By comparing our results with previously published studies, we construct a three-component model to describe elastic and inelastic compaction of granular sand. This model involves acoustically silent grain rearrangement that contributes significant inelastic strain at low pressures, poro-elastic (Hertzian) deformation at all pressures, and inelastic strain related to granular cracking and particle failure which increases in significance at greater pressures.     

Strain analysis in Jurassic argillites of the Monte Sirino area (Lagonegro Zone,southern Apennines,Italy) and implications for deformation paths in pelitic rocks

Analysis of strain in Jurassic argillites forming part of the folded and thrusted sedimentary succession of the Lagonegro basin (southern Italian Apennines) has been carried out using ellipsoid-shaped reduction spots as strain markers. Most of the determined finite strain ellipsoids are of oblate type and show a peculiar distribution of the maximum extension direction (X), with maxima either subparallel or subperpendicular to the local fold axes. Using the strain matrix method, two different deformation histories have been considered to assist the interpretation of the observed finite strain pattern. A first deformation history involved vertical compaction followed by horizontal shortening (occurring by a combination of true tectonic strain and volume loss), whereby all strain is coaxial and there is no change in the intermediate axis of the strain ellipsoid. By this type of deformation sequence, which produces a deformation path where total strain moves from the oblate to the prolate strain field and back to the oblate field, prolate strain ellipsoids can be generated and may be recorded where tectonic deformation has not been large enough to reverse pretectonic compaction. This type of deformation history may be of local importance within the study area (i.e. it may characterize some fold hinge regions) and, more generally, is probably of limited occurrence in deformed pelitic rocks. A second deformation sequence considered the superposition of pre-tectonic compaction and tectonic strain consisting of initial layer-parallel shortening followed by layer-parallel shear (related to flexural folding). Also in this instance, volume change during tectonic deformation and tectonic plane strain have been assumed. For geologically reasonable amounts of volume loss due to compaction and of initial layer-parallel shortening, this type of deformation history is capable of producing a deformation path entirely lying within the oblate strain field, but still characterized by a changeover, during deformation, of the maximum extension axis (X) from a position parallel to the fold axis to one perpendicular to it. This type of deformation sequence may explain the main strain features observed in the study area, where most of the measured finite strain ellipsoids, determined from the limb regions of flexural folds, display an oblate shape, irrespective of the orientation of their maximum extension direction (X) with respect to the local structural trends. More generally, this type of deformation history provides a mechanism to account for the predominance of oblate strains in deformed pelitic rocks.     

Orebody geometry in lode gold deposits from Zimbabwe: implications for fluid flow,deformation and mineralization

The geometry of some orebodies can be described simply and accurately by three orthogonal axes, U≥V≥W. The ratios between these axes can be expressed as a parameter j=(U/V−1)/(V/W−1), and represented by a graph of U/V plotted against V/W, analogous to the treatment of strain ellipsoids. The orientations of orebodies can be plotted simply on projections using the UVW axes. Measurements of ore bodies from two examples of lode gold deposits from the Zimbabwe craton show that most of these orebodies are oblate. However, orebodies can have significant U/V ratios, implying a component of pipe-like fluid flow during mineralization. Pipe flow is demonstrated to be orders of magnitude more conductive than flow in planar veins and faults. There are significant variations in orebody geometry between deposits and within different sections of a single deposit. W values appear to be influenced by host rock: more permeable rocks have higher W. A negative trend of j value with orebody volume indicates that orebodies do not evolve in a self-similar way, but tend to more oblate shapes with increasing volume.     

Unstable flow,magma mixing and magma-rock deformation in a deep-seated conduit: the Gil-Márquez Complex,south-west Spain

The Gil-Marquez Complex is an exceptional outcrop of plutonic rocks ranging in composition from diorites to granites emplaced into Devonian terrigenous metasediments of the southernmost part of the Hercynian basement of Iberia. A combined study of this complex, including field geology, petrology, structural geology and geochemistry, reveals that it represents an ancient conduit of magma transport through the continental crust. This conduit allowed the intrusion of magmas of contrasted compositions. Two end-members and several hybrids are identified. The first end-member is a biotite granite and the second is a basaltic magma generated by partial melting of a depletedmantle source. Both magmas rose through a common channel in which favorable conditions for unstable flow and magma mixing occurred. The observed relations in the Gil-Márquez Complex show that mixing in conduits may be an important mechanism for producing homogeneous hybrid magmas.     

可可西里西段羊湖盆地沉积、构造特征及其动力学意义

对可可西里西段新生代盆地缺乏了解是导致该区新生代地质演化存在争议的重要原因.本文以沉积学和构造变形分析为主要手段,对可可西里西段羊湖盆地时代、充填序列、物源区和变形特征进行了分析,结果表明,羊湖盆地新生界沉积厚度大于1302m,主要由下部雅西错组冲积扇相碎屑岩和上部五道梁群湖泊相碳酸盐岩组成,其岩石组合和充填序列与可可西里东段具有一致性,同时古流向和碎屑锆石U-Pb年代学分析显示盆地物源来南部的羌塘地块,盆地形成演化受南部褶皱冲断带制约,盆地构造变形强烈,沿褶皱冲断带和羊湖盆地地壳分别发生51％和41％的缩短.沉积充填结构和变形特征表明,羊湖盆地与东段可可西里盆地具有相同的演化历史和性质,预示青藏高原中部在渐新世-中新世在存在一个大的、统一的可可西里盆地.     

Finite deformation analysis of liquefaction-induced flow failure in soil embankments

Summary A finite element formulation is proposed for finite deformation dynamic analysis of saturated soil systems. The formulation is based on an updated Lagrangian approach and specifically considers the finite deformation effects on the flow of water through a soil element which undergoes a large deformation or rotation. A two-surface plasticity model is used to model the stress-strain behaviour of the soil skeleton. The proposed formulation has been implemented and is applied to simulate the response of a centrifuge model embankment. The calculated response is in good agreement with the observed behaviour of the soil embankment in the centrifuge test.     

The role of dewatering in the progressive deformation of a sandy accretionary wedge: constraints from direct imagings of fluid flow and void structure

Geological investigation of the deformation structures and sedimentary setting of the Emi Group, a Miocene sand-rich accretionary complex, central Japan, revealed a six stage-structural evolution during shallow level accretion in a subduction zone. The early deformation (stage 1) is characterized by independent particulate flow in layer parallel faults, scaly cleavages and web structures, and upward dewatering in dish-and-pillar structures and breccia injections, while later deformation (stages 2–6) involve mappable scale folding, meso- to macro-scopic thrusts and web structures with cataclastic flow. Based on microscopic analyses of these structures, the early faulting with independent particulate flow (stage 1 deformation) is associated with dilatancy and preferred orientation of void space, whereas the later faulting with cataclastic flow (stage 2 deformation) occurs with compaction and crude preferred orientation. The former features imply more permeable fluid migration pathways, supported by the permeability measurements and direct imaging of fluid flow by X-ray CT. On the other hand, the later fault zone has lower permeability and porosity than intact rock, and plays as fluid sealing. Thus, in the early stage (stages 1), fluid flow occurs as focused flow through dilatant fault zones with independent particulate flow or fluid migration by upward dewatering forming dish-and-pillar structures and breccia injections, whereas no evidence of fluid flow is recognized at the later stages (stages 2–6). Namely the fault zones focus fluid flow during primary accretion in shallow levels, and the fluid flow is strongly controlled by the deformation mechanism. Furthermore, the change of the deformation mechanism could be effected by progressive increment of the confining pressure, accompanied with accretion and lithification in the accretionary prism. In the shallow, dilatant-faulting regime where the deformation mechanism is independent particulate flow, focused flow dominates, whereas in the deep, cataclastic regime distributed flow may play a main conduit rather than the focused flow.     

Characterisation of deformation and flow mechanics around porphyroclasts in a calcite marble ultramylonite by means of EBSD analysis

This paper studies the flow heterogeneity around porphyroclasts associated with greenschist facies deformation of a calcite marble shear zone. Microstructural data from electron backscatter diffraction analyses (EBSD) are used to constrain the flow mechanics of this dominantly non-coaxial type of deformation. The microstructure of the undisturbed ultramylonite (grain-size range 5–100 μm, mean 40 μm) is interpreted to represent steady-state (time-independent) flow conditions with flow planes parallel to the shear zone boundary. Single calcite porphyroclasts (grain-size 1–3 mm) caused flow perturbation in the fine-grained marble ultramylonite. It is the shape, in particular, of these rigid porphyroclasts that controls their rotational behaviour during deformation and, therefore, the development of specific flow fabrics. The flow planes around elongated-rhomboidal, stable porphyroclasts change the orientation to become roughly parallel to the porphyroclast margin, whereas the geometry of flow planes around nearly equant, rotating porphyroclasts describes a δ-type flow pattern. We infer that to some extent decoupling at the clast–matrix interface has occurred to guarantee a stable orientation of elongated porphyroclasts, but was not sufficient to reduce the rotation rate of equant clasts to zero. According to the flow deflection, the general crystallographic preferred orientation (CPO) with its single c-axis maximum perpendicular to the flow plane is rotated about an axis which is (sub)parallel to the kinematic rotation axis of the shear zone. Ultramylonite microstructures, CPOs and misorientation data are best explained by the dual operation of grain-size-insensitive (dislocation creep with recovery and recrystallization) and grain-size-sensitive (diffusion creep) mechanisms. The limited grain-size reduction around porphyroclasts suggests that the grain-size-insensitive mechanisms controlled rheology.     

Paleocollapse structure as a passageway for groundwater flow and contaminant transport

 Paleocollapse structure is a rock collapse, resulting from the failure in the geological history of the bedrock overlying karstified limestone. Depending on the present hydrogeological conditions within the area of paleocollapse and the internal properties of these structures, they can provide a means to facilitate groundwater flow and contaminant transport. Inactive paleocollapse structures can be reactivated by human activities such as dam construction, mining underground minerals, pumping groundwater, and development of landfills. They can also be reactivated by natural events such as earthquakes and neotectonic movements. In the mines of northern China, sudden inflow of karst water from Ordovician limestone into drifts and mining stopes through paleocollapse structures has caused significant economic loss. Water pumping tests and accompanied dye traces are effective approaches of locating water-conducting paleocollapse structures. Grouting is probably the best means of preventing them from becoming geohazards. Received: 26 November 1996 · Accepted: 17 June 1997     

Flow laws of multiphase materials and rocks from end-member flow laws

Rheological properties of polyphase rocks play an important role in the dynamics of the lithosphere and asthenosphere. However, flow laws for large portions of the polyphase rocks in the Earth's crust and mantle have not been well determined. An analysis based on the theory of mixtures has been made to calculate the general flow laws of coarse (≥5 μm), nearly equant-grained (aspect-ratio ≤3), and massive polyphase rocks and materials for which only the flow laws and volume fractions of the constituents are taken into consideration in the modeling of the bulk rheology and effects of microstructure could be ignored. The theoretical analysis is based on three assumptions: (1) the polyphase composite and the monophase aggregates of its constituents obey the same kind of flow laws (linear, power or exponential), (2) there is no change in the operative deformation mechanism of each phase when it is in the composite as compared to when it is in a monophase aggregate, and (3) neither chemical (metamorphic) reactions take place among the constituent phases nor eutectic melting occurs due to the phase mixing. The proposed iterative process allows to predict, to the first approximation, the flow laws for a large number of polyphase rocks in terms of the experimentally determined flow laws of a relatively small number of monomineralic aggregates. Applications of this approach to typical polyphase rocks such as granite, diorite, diabase, aplite and websterite as well as to synthetic two-phase materials such as forsterite–enstatite mixtures and water ice–ammonia dehydrate aggregates yield quite accurate approximations to the experimental values.     

Fracture step structure: geometrical characterization and effects on fluid flow and breakthrough curve

Groundwater flow and solute transport through fractured rock is highly responsive to the hydraulic anisotropy and heterogeneity that are specific to every major fracture. A major fracture is modeled as the combination of some primal master fractures and several splay fractures that branch out from primal master fractures: step structures (or jog parts). Step structures are commonly observed along a major fracture on various scales. Master fractures were formed and developed by shear movement while some splay fractures were formed by extension normal to their wall. This difference in fracturing process may lead to a permeability difference between master fractures and splay fractures which seems to be one of the major factors controlling flow and solute transport through the fracture networks due to its hydraulic anisotropic and heterogeneous features. This study is composed of two major components: (1) identification and characterization of a step structure from borehole data; (2) evaluation of effect of some idealized step structures on breakthrough curve by numerical simulations. The fracture data of four 1000-m boreholes were used to make clear fracture patterns in the Tono area of Japan. Some major fractures were identified using stereographic projection technique. On the basis of these results, several idealized models of a major fracture having a step was constructed for the numerical study. The obtained results from numerical simulations clearly imply that geometry of step structure plays an important role in flow and transport through the fracture networks.     

帕米尔东北缘地区构造变形特征与盆山结构

帕米尔东北缘乌泊尔地区是正确认识帕米尔北缘盆山结构和构造变形特征非常关键的地区,本文利用连续电磁剖面(CEMP)资料和地震资料,并结合野外地质调查资料和钻井资料,对帕米尔东北缘乌泊尔地区的盆山结构和构造变形特征进行了研究。认为帕米尔东北缘及其以北地区的盆山结构表现为帕米尔造山带向北冲断和南天山向南冲断所形成的对冲结构;帕米尔山前为基底卷入式构造,古生界—中生界沿高角度的逆冲断层推覆到新近系和第四系之上,形成山前的古生界—中生界逆冲推覆带;北侧由受乌泊尔断裂控制的深部隐伏冲断体系和浅部的第四纪背驮盆地所构成。研究区的新生代构造变形时间开始于上新世晚期,并持续变形至今,形成了下更新统西域组(Q_1x)与下伏上新统、Q_2与Q_1和Q_(3-4)与Q_2之间的不整合。研究区最小构造缩短量为48.6 km,缩短率为48.1%。     

中国阿尔泰南缘构造变形层次与变形序列

中国阿尔泰造山带发育了一套被认为是早古生代的变质岩系,保留了原始成分层理S₀。岩石在下-中地壳层次经历了塑性变形及最深达到角闪岩相变质作用的改造。通过对岩石构造层次及变形序列分析,在造山带南缘识别出由深至浅、由早到晚发育的4个构造层次的变形样式,分别为流变构造、韧性剪切构造、褶皱劈理构造和断裂节理构造。流变构造以轴面NW走向、枢纽近直立的流动褶曲和NW走向近直立片麻岩面理为特征,局部发育长英质和镁铁质矿物分异条带,相当于下构造层次变形（约15～18 km）。韧性剪切构造以NW-NNW向和NNE-NEE向近直立的糜棱岩面理和片理面为特征,相当于下-中构造层次变形（约8～12 km）。褶皱劈理构造对应轴向NW、枢纽近水平的等厚褶皱、NWW和NNW向近直立劈理的产出,为中-上构造层次变形（约5～8 km）。上构造层次的断裂、节理构造在造山带南缘广泛发育,变形深度约浅于3～ 4 km。根据变形几何学、运动学关系及变形次序,确认它们是3期构造应力场作用的结果。在二叠—三叠纪,阿尔泰经历NE向（47°～55°）挤压应力场,在下-中构造层次发育317°～325°走向的流动褶皱、片麻理构造S_1-1,306°～308°走向左行和328°～332°走向右行近直立韧性剪切面理S_1-2,290°～295°走向左行和350°～355°走向右行、近直立片理面S_1-3。在古近纪NW向挤压应力场下,中-上构造层次发育20°～22°走向左行和78°走向右行片理面S_2-1,10°～20°和70°～78°近直立脆—韧性剪切面S_2-2和5°～8°压扭性断裂构造S_2-3。新近纪近N-S向应力场下,在浅地壳层次产出320°～325°和340°～345°走向右行走滑断裂,NWW向破劈理和NE向节理构造S₃。下构造层次变形被中-上构造层次变形叠加改造,体现变形深度层次与变形序列的一致性,反映了同构造挤压期的地壳抬升现象。     

Structural-event framework for the eastern Yilgarn Craton,Western Australia,and its implications for orogenic gold

The NNW-trending tectonic grain of the eastern Yilgarn Craton (EYC) was established as a result of predominantly ENE–WSW directed extension (D1 and D3) and E(ENE)–W(WSW) (D2, D4) to NE–SW directed (D5) contraction. The result has been a succession of NNW-striking temporally discrete fabric elements, which can be difficult to interpret reliably at any single location. Despite this, many past workers interpreted the NNW-striking fabric as the result of only one regional contractional event, and used it as a marker for correlating structural events across the region. In order to unravel the complexity, this paper presents a new sixfold (D1–D6) deformation nomenclature based on >10,000 new mesoscale structural observations, including their kinematic analysis and cross-cutting relationships. These mesoscale data were referenced with regional 3D map patterns, stratigraphic-magmatic-metallogenic considerations, and deep seismic reflection images. This integrated geodynamic-architectural approach is applicable to solving structural-event histories in other polydeformed terrains. Gold mineralisation occurred during the first five events, but was particularly concentrated from D3 onwards. The D3 event marked the most profound change in the tectonic evolution of the EYC, with changes in greenstones, granites and tectonic mode (lithospheric extension and core complexes), with the first significant gold deposited within extensional shear zones that dissect the crust. Later contraction (D4) was imposed at a high angle to the previously established anisotropic architecture. The outcome was the creation of a new dynamic permeability framework, which resulted in gold mineralisation during NNW-striking sinistral strike-slip faulting and associated thrusting. A further stress switch (D5) further modified the architecture resulting in N- to NNE-striking dextral strike-slip faulting, and the final period of gold mineralisation, before late-stage extension (D6).