Polygonal deformation bands

We report for the first time the occurrence of polygonal faults in sandstone, which is compelling given that layer-bound polygonal fault systems have been observed so far only in fine-grained sediments such as clay and chalk. The polygonal faults are shear deformation bands that developed under shallow burial conditions via strain hardening in dm-wide zones. The edges of the polygons are 1–5 m long. The shear deformation bands are organized as conjugate faults along each edge of the polygon and form characteristic horst-like structures. The individual deformation bands have slip magnitudes ranging from a few mm to 1.5 cm; the cumulative average slip magnitude in a zone is up to 10 cm. The deformation bands heaves, in aggregate form, accommodate a small isotropic horizontal extension (strain <0.005). The individual shear deformation bands show abutting T-junctions, veering, curving, and merging where they mechanically interact. Crosscutting relationships are rare. The interactions of the deformation bands are similar to those of mode I opening fractures. The documented fault networks have important implications for evaluating the geometry of km-scale polygonal fault systems in the subsurface, top seal integrity, as well as constraining paleo-tectonic stress regimes.     

Strain variation with progressive deformation in basement-involved trishear structures

Basement-involved structures associated with reverse, vertical and normal faults commonly involve non-parallel shear within a triangular deformation (trishear) zone located on the front limbs of the structures. Deformation within the trishear zone is characterized by shear gradients and an associated decrease in the dips of the beds in stratigraphically higher units. Geometric models suggest that the layer-parallel strain within the trishear zone depends on the type of fault (normal, reverse, or vertical), the dip and throw of the fault, the dip of the anticlinal or synclinal axial surfaces, and the distance of any unit above the initial tip of the trishear zone, located at the basement-sediment contact. At any given location, reverse faults typically show increasing layer parallel shortening, followed by decreasing layer parallel shortening and a transition to extension, with increasing throw. The transition from contraction to extension occurs at lower values of throw for stratigraphically lower units and also for faults with smaller dips. Vertical and normal faults exhibit increasing layer-parallel extension of all units with increasing throw, with larger extension for stratigraphically lower units. Experimental models suggest that the trishear zone can expand with increasing fault throw. The strain within the trishear zones is accommodated largely by secondary faults, which are rotated with progressive deformation. The strain variations in the experiments closely mimic those predicted by the geometric models for reverse, vertical, and normal faults.     

Influence of host lithofacies on fault rock variation in carbonate fault zones: A case study from the Island of Malta

Relatively few studies have examined fault rock microstructures in carbonates. Understanding fault core production helps predict the hydraulic behaviour of faults and the potential for reservoir compartmentalisation. Normal faults on Malta, ranging from <1 m to 90 m displacement, cut two carbonate lithofacies, micrite-dominated and grain-dominated carbonates, allowing the investigation of fault rock evolution with increasing displacement in differing lithofacies. Lithological heterogeneity leads to a variety of deformation mechanisms. Nine different fault rock types have been identified, with a range of deformation microstructures along an individual slip surface. The deformation style, and hence type of fault rock produced, is a function of host rock texture, specifically grain size and sorting, porosity and uniaxial compressive strength. Homogeneously fine-grained micrtie-dominated carbonates are characterised by dispersed deformation with large fracture networks that develop into breccias. Alternatively, this lithofacies is commonly recrystallised. In contrast, in the coarse-grained, heterogeneous grain-dominated carbonates the development of faulting is characterised by localised deformation, creating protocataclasite and cataclasite fault rocks. Cementation also occurs within some grain-dominated carbonates close to and on slip surfaces. Fault rock variation is a function of displacement as well as juxtaposed lithofacies. An increase in fault rock variability is observed at higher displacements, potentially creating a more transmissible fault, which opposes what may be expected in siliciclastic and crystalline faults. Significant heterogeneity in the fault rock types formed is likely to create variable permeability along fault-strike, potentially allowing across-fault fluid flow. However, areas with homogeneous fault rocks may generate barriers to fluid flow.     

煤岩两体模型变形破坏数值模拟

采用拉格朗日元法,在弹性岩石与弹性-应变软化煤体所构成的平面应变两体模型的上、下端面上不存在水平方向摩擦力条件下,模拟了模型的破坏过程、岩石高度对模型及煤体全程应力-应变曲线、煤体变形速率、煤体破坏模式及剪切应变增量分布的影响。结果表明,当模型的全程应力-应变曲线达到峰值时煤体内部的剪切带图案已经十分明显,在模型的应变硬化阶段,煤体中的应变局部化可视为模型失稳破坏的前兆,随岩石高度的增加,模型应力-应变曲线的软化段变得陡峭,这与单轴压缩条件下的解析解在定性上是一致的;煤体应力-应变曲线的软化段变得平缓,煤体消耗能量的能力增强;弹性阶段煤体的变形速率降低;煤体内部的剪切应变增量增加。煤体应力-应变曲线的软化段的斜率、弹性阶段煤体的变形速率、煤体内部的剪切应变增量及塑性耗散能都受岩石高度的影响,说明了岩石几何尺寸对煤体的影响(煤岩相互作用)是不容忽视的。     

Reverse fault slip through soft rock and sand strata by centrifuge modeling tests

The devastating damage after the 1999 Chi-Chi and 1999 Izmit earthquakes has greatly motivated soil–reverse fault interaction studies. However, most centrifuge modeling studies have employed a single homogeneous soil layer during testing, which does not represent in situ conditions. Indeed, while geological conditions vary spatially, engineering soils are often underlain by soft rocks. Therefore, four centrifuge models were developed to evaluate the effect of soft rock layers on the ground surface and subsurface deformation. Sand–cement mixtures of varying thicknesses with a uniaxial compressive strength of 0.975 MPa, simulating extremely soft rock, were overlain by pluviated sandy soil. The model thickness was 100 mm, corresponding to 8 m in the prototype scale when spun at 80 g. Every model was subjected to a vertical offset of 50 mm/4 m (0.5 H; H: total sedimentary deposit thickness) along a reverse fault with a 60° dip. The results indicate that the presence of a soft rock stratum results in the creation of a horst profile at the ground surface. Additionally, the thinner the soil layer on top of the soft rock stratum is, the longer and higher the horst created at the ground surface. Consequently, the fault deformation zone lengthens proportionally with the increasing thickness ratio of the soft rock. Furthermore, the presence of soft rock as an intermediary stratum between bedrock and soil causes the deformation zone boundary on the hanging wall side to move in the direction of fault movement.

     

Damage zone and slip-surface evolution over μm to km scales in high-porosity Navajo sandstone, Utah

Detailed mapping of throw variations and deformation along two km-scale normal faults in the high-porosity Navajo sandstone, Utah, has been used to investigate fault growth in this lithology. The faults consist of one or more through-going, striated, slip-surfaces, accommodating the greater part of the offset surrounded by a damage zone consisting of deformation band clusters and short, unconnected slip-surfaces. In contrast to previous models for deformation in this lithology, we find that the nucleation of slip-surfaces begins where measurable throw is negligible and deformation bands are forming and increasing in number. The microstructure and porosity of deformation bands and slip surfaces are distinct and independent of the amount of offset that they accommodate, i.e. they represent different and yet contemporaneous deformation mechanisms. The point where measurable throw begins to accumulate (the fault tip) is marked by the first through-going connected slip-surface. Increase in throw towards the centre of the fault results in a three-dimensional strain field, producing orthorhombic structural geometries within the damage zone. We find that the total width of the damage zone increases as offset is accumulated. For these faults, the damage zone width is approximately 2.5 times the total fault throw.     

Geometry and architecture of faults in a syn-rift normal fault array: The Nukhul half-graben, Suez rift, Egypt

The geometry and architecture of a well exposed syn-rift normal fault array in the Suez rift is examined. At pre-rift level, the Nukhul fault consists of a single zone of intense deformation up to 10 m wide, with a significant monocline in the hanging wall and much more limited folding in the footwall. At syn-rift level, the fault zone is characterised by a single discrete fault zone less than 2 m wide, with damage zone faults up to approximately 200 m into the hanging wall, and with no significant monocline developed. The evolution of the fault from a buried structure with associated fault-propagation folding, to a surface-breaking structure with associated surface faulting, has led to enhanced bedding-parallel slip at lower levels that is absent at higher levels. Strain is enhanced at breached relay ramps and bends inherited from pre-existing structures that were reactivated during rifting. Damage zone faults observed within the pre-rift show ramp-flat geometries associated with contrast in competency of the layers cut and commonly contain zones of scaly shale or clay smear. Damage zone faults within the syn-rift are commonly very straight, and may be discrete fault planes with no visible fault rock at the scale of observation, or contain relatively thin and simple zones of scaly shale or gouge. The geometric and architectural evolution of the fault array is interpreted to be the result of (i) the evolution from distributed trishear deformation during upward propagation of buried fault tips to surface faulting after faults breach the surface; (ii) differences in deformation response between lithified pre-rift units that display high competence contrasts during deformation, and unlithified syn-rift units that display low competence contrasts during deformation, and; (iii) the history of segmentation, growth and linkage of the faults that make up the fault array. This has important implications for fluid flow in fault zones.     

The mechanics of a magma chamber-fault system in trans-tension with application to Coso

The mechanical interaction between an elliptically shaped magma chamber and a fault subject to transtension is investigated with particular reference to the Coso geothermal field. The geologic setting of the Coso field is interpreted as a releasing bend step-over structure formed by the Airport Lake and Owens Valley dextral strike-slip fault system. The role of the Coso volcano-magmatic center in the development of the “over-step” structure is examined by treating the magma chamber as a liquid inclusion in a viscoelastic crust containing a fault (Airport Lake). The problem is numerically solved using a 2D viscoelastic finite element model with thermally activated viscosity to account for thermal weakening of the rock. The temperature distribution around the magma body is calculated based on a 3D steady-state approach and using the mesh-less numerical method. The fault is modeled as a frictionless contact. The simulated distributions of stress and strain around the inclusion display a rotation caused by the shearing component of the applied transtension. The results indicate that the fault tends to overstep the chamber in a geometric pattern similar to a step-over. There is good correspondence between the computed distributions of the maximum shear stress in the vicinity of the magma chamber and the map of earthquake epicenters at a depth of 7–10 km in Coso.     

Kinematics-Based Mathematical Model for Deforming Thrust Wedges

Given the wealth of data concerning the kinematics of deforming fold-thrust belts (FTBs), first-order generalizations about how the major strain components vary within a deforming thrust wedges are considered. These generally observed strain patterns are used to constrain a general, kinematics-based, FTB-wedge model. We considered five strain components within a deforming thrust sheet: (1) thrust-parallel simple shear, (2) horizontal contractional strain, (3) thrust-normal reaction strain, (4) gravitational strain, and (5) a lateral confining boundary condition. After making assumptions about how these strain components vary within a model FTB-wedge, the incremental deformation matrix can be calculated for any given point within the deforming wedge. Thus, the material path of a given marker can be determined and an initially spherical marker’s strain path can be calculated as it moves through the deforming wedge. Furthermore, by illustrating various kinematic parameters of many initially spherical markers (for example, Flinn’s k-value, incremental octahedral shear strain, transport-perpendicular stretch), we have assembled representations of the kinematic properties of the entire model wedge. By including a flat-ramp-flat fault surface geometry for the model wedge, we are able to examine the kinematic effects of this relatively common structural geometry. Within the fault ramp segment there are greater incremental strain magnitudes, out-of-the-plane motion, and flattening strains. Additionally, data from this model suggests that gravitational strains potentially have a significant effect on the strain distribution within a deforming thrust wedge. M. Mookerjee is formerly Matthew Strine.     

On the lag time between internal strain and basement involved thrust induced exhumation: The case of the Colombian Eastern Cordillera

Thrust sheets accumulate internal strain before they start moving along discrete fault planes. However, there are no previous studies evaluating the time difference between initiation of strain and fault displacement. In this paper we use observations from the Eastern Cordillera of Colombia to evaluate this interval. We utilize multiple thermochronometers and paleothermometers to refine the timing of deformation. Based on these new data we build time-temperature path estimates that together with geometric outcrop-based structural analysis and fluid inclusions allow us to assign relative timing to features associated with strain, such as cleavage, veins and certain types of fractures, and compare that with the timing of thrusting. We find that cleavage was only formed close to maximum paleotemperatures, almost coeval with the onset of thrust-induced denudation by the Late Oligocene. The corresponding structural level of fold-related veins suggest that they were formed later but still when the country rocks were at temperatures higher than 160 °C, mostly during the Early Miocene and still coexisted with the latest stages of cleavage formation. Our data show that the main period of strain hardening was short (probably a few million years) and occurred before first-order basement thrusting was dominant, but was associated with second-order folding.     

Shear localization,velocity weakening behavior,and development of cataclastic foliation in experimental granite gouge

Microstructural aspects of room-temperature deformation in experimental Westerly granite gouge were studied by a set of velocity stepping rotary-shear experiments at 25 MPa normal stress. The experiments were terminated at: (a) 44 mm, (b) 79 mm, and (c) 387 mm of sliding, all involving variable-amplitude fluctuations in friction. Microstructural attributes of the gouge were studied using scanning (SEM) and scanning transmission electron microscopy (STEM), image processing, and energy dispersive X-ray (EDX) analyses. The gouge was velocity weakening at sliding distances >10 mm as a core of cataclasites along a through-going shear zone developed within a mantle of less deformed gouge in all experiments. Unlike in experiment (a), the cataclasites in experiments (b) and (c) progressively developed a foliation defined by stacks of shear bands. The individual bands showed an asymmetric particle-size grading normal to shearing direction. These microstructures were subsequently disrupted and reworked by high-angle Riedel shears. While the microstructural evolution affected the effective thickness and frictional strength of the gouge, it did not affect its overall velocity dependence behavior. We suggest that the foliation resulted from competing shear localization and frictional slip hardening and that the velocity dependence of natural fault gouge depends upon compositional as well as microstructural evolution of the gouge.     

Strain compatibility and fault linkage in relay zones on normal faults

Relay zones on normal faults are unlikely to have tabular geometries as depicted in idealised models. Rotation of a relay ramp between non-parallel and non-planar relay-bounding faults will inevitably lead to strain compatibility problems causing open gaps or overlaps within the relay zone. Linkage of relay-bounding faults does not evolve from a single branch point. Rather, linkage occurs at multiple points along the fault tip lines giving rise to initially discontinuous branch lines. Where linkage occurs along a discontinuous slip-aligned branch line, displacement at different levels within the relay zone is partitioned between variable amounts of ramp rotation and slip across the branch line. The linking fault propagates when strain compatibility can no longer be maintained by continuous deformation processes, such as thickening or thinning of incompetent layers within the relay ramp. Step-like changes in vertical displacement vs. distance (d − x) profiles on horizons containing apparently intact relay ramps are probably indicative of incipient breaching and can be used predict the presence of a slip-aligned branch line in the sub-surface. Despite the complexity of the strain distribution within relay zones, the total vertical displacement across the relay remains geometrically coherent at all levels.     

The nature and origin of off-fault damage surrounding strike-slip fault zones with a wide range of displacements: A field study from the Atacama fault system, northern Chile

Damage surrounding the core of faults is represented by deformation on a range of scales from microfracturing of the rock matrix to macroscopic fracture networks. The spatial distribution and geometric characterization of damage at various scales can help to predict fault growth processes, subsequent mechanics, bulk hydraulic and seismological properties of a fault zone. Within the excellently exposed Atacama fault system, northern Chile, micro- and macroscale fracture densities and orientation surrounding strike-slip faults with well-constrained displacements ranging over nearly 5 orders of magnitude (0.12 m–5000 m) have been analyzed. Faults have been studied that cut granodiorite and have been passively exhumed from 6 to 10 km depth. This allows direct comparison of the damage surrounding faults of different displacements. The faults consist of a fault core and associated damage zone. Macrofractures in the damage zone are predominantly shear fractures orientated at high angles to the faults studied. They have a reasonably well-defined exponential decrease with distance from the fault core. Microfractures are a combination of open, healed, partially healed and fluid inclusion planes (FIPs). FIPs are the earliest set of fractures and show an exponential decrease in fracture density with perpendicular distance from the fault core. Later microfractures do not show a clear relationship of microfracture density with perpendicular distance from the fault core. Damage zone widths defined by the density of FIPs scale with fault displacement but appear to reach a maximum at a few km displacement. One fault, where damage was characterized on both sides of the fault core shows no damage asymmetry. All faults appear to have a critical microfracture density at the fault core/damage zone boundary that is independent of displacement. An empirical relationship for microfracture density distribution with displacement is presented. Preferred FIP orientations have a high angle to the fault close to the fault core and become more diffuse with distance. Models that predict off-fault damage such as a migrating process zone during fault formation, wear from geometrical irregularities and dynamic rupture are all consistent with our data. We conclude it is very difficult to distinguish between them on the basis of field data alone, at least within the limits of this study.     

Distinctive diamagnetic fabrics in dolostones evolved at fault cores,the Dead Sea Transform

We resolve the anisotropy of magnetic susceptibility (AMS) axes along fault planes, cores and damage zones in rocks that crop out next to the Dead Sea Transform (DST) plate boundary. We measured 261 samples of mainly diamagnetic dolostones that were collected from 15 stations. To test the possible effect of the iron content on the AMS we analyzed the Fe concentrations of the samples in different rock phases. Dolostones with mean magnetic susceptibility value lower than −4 × 10⁻⁶ SI and iron content less than ∼1000 ppm are suitable for diamagnetic AMS-based strain analysis. The dolostones along fault planes display AMS fabrics that significantly deviate from the primary “sedimentary fabric”. The characteristics of these fabrics include well-grouped, sub-horizontal, minimum principal AMS axes (k₃) and sub-vertical magnetic foliations commonly defined by maximum and intermediate principal AMS axes (k₁ and k₂ axes, respectively). These fabrics are distinctive along fault planes located tens of kilometers apart, with strikes ranging between NNW-SSE and NNE-SSW and different senses of motion. The obtained magnetic foliations (k₁–k₂) are sub-parallel (within ∼20°) to the fault planes. Based on rock magnetic and geochemical analyses, we interpret the AMS fabrics as the product of both shape and crystallographic anisotropy of the dolostones. Preferred shape alignment evolves due to mechanical rotation of subordinate particles and rock fragments at the fault core. Preferred crystallographic orientation results from elevated frictional heating (>300 °C) during faulting, which enhances c-axes alignment in the cement-supported dolomite breccia due to crystal-plastic processes. The penetrative deformation within fault zones resulted from the local, fault-related strain field and does not reflect the regional strain field. The analyzed AMS fabrics together with fault-plane kinematics provide valuable information on faulting characteristics in the uppermost crust.     

The initiation of brittle faults in crystalline rock

Faults in the upper crust initiate from pre-existing (inherited) or precursory (early-formed) structures and typically grow by the mechanical interaction and linkage of these structures. In crystalline rock, rock architecture, composition, cooling, and exhumation influence the initiation of faults, with contrasting styles observed in plutonic rocks, extrusive igneous rocks, and foliated metamorphic rocks. Brittle fault growth in granitic rock is commonly controlled by the architecture of inherited joints or preexisting dikes. In basalt, abundant joints control the surface expression of faulting, and enhanced compliance due to abundant joints leads to folding and deformation asymmetry in the fault zone. Highly reactive mafic minerals likely become rapidly evolving fault rocks. In foliated metamorphic rocks, fault initiation style is strongly influenced by strength anisotropy relative to the principal stress directions, with fracturing favored when the foliation is aligned with the directions of principal stress. The continuity of micas within the foliation also influences the micromechanics of fault initiation. Brittle kink bands are an example of a strain-hardening precursory structure unique to foliated rock. Each of these fault initiation processes produces different initial fault geometry and spatial heterogeneity that influence such properties as fault permeability and seismogenesis.     

Very high strains recorded in mylonites along the Alpine Fault, New Zealand: implications for the deep structure of plate boundary faults

Oblique displacement on the Alpine Fault, which forms the principal structure along the Australian–Pacific plate boundary in South Island, New Zealand, has resulted in exhumation of a kilometre-wide mylonite zone in the hanging wall adjacent to the current brittle fault trace. The mylonites formed under amphibolite facies conditions at depths of ca. 25 km and have been uplifted during the past 5 Ma. A suite of 65–70 Ma pegmatite veins in the hanging wall Alpine schists has been progressively deformed within the mylonite zone and sheared out over a strike length of ca. 100 km. Measurements of the thickness distribution of the pegmatite veins within the non-mylonitised schists and at three localities within the progressively strained mylonites have been used to estimate strain values within the mylonites. The thicknesses approximate a log-normal distribution, with a mean value that is progressively reduced through the protomylonites, mylonites and ultramylonites. By assuming that the thickness distribution currently observed in the schists was the same for the pegmatites within the mylonites before strain, a model of deformation incorporating simple shear and simultaneous pure shear is used to strain the undeformed veins until a fit is obtained with the strained distributions. Shear strains calculated range from 12 to 22 for the protomylonites, 120 to 200 for the mylonites and 180 to 300 for the ultramylonites, corresponding to pure shear values of 1–3 in each case. These values are compatible with the strains predicted if most of the surface displacement on the fault over the past 5 Ma were accommodated within a 1–2-km-wide mylonite zone through the middle and lower crusts. The results suggest that processes such as erosional focussing of deformation and thermal weakening may cause intense strain localisation within the lower crust, with plate boundary deformation restricted to narrow zones rather than becoming increasingly distributed over a widening shear zone with depth.     

Fault geometry and mechanics of marly carbonate multilayers: An integrated field and laboratory study from the Northern Apennines,Italy

Sealing layers are often represented by sedimentary sequences characterized by alternating strong and weak lithologies. When involved in faulting processes, these mechanically heterogeneous multilayers develop complex fault geometries. Here we investigate fault initiation and evolution within a mechanical multilayer by integrating field observations and rock deformation experiments. Faults initiate with a staircase trajectory that partially reflects the mechanical properties of the involved lithologies, as suggested by our deformation experiments. However, some faults initiating at low angles in calcite-rich layers (θ_i = 5°–20°) and at high angles in clay-rich layers (θ_i = 45°–86°) indicate the important role of structural inheritance at the onset of faulting. With increasing displacement, faults develop well-organized fault cores characterized by a marly, foliated matrix embedding fragments of limestone. The angles of fault reactivation, which concentrate between 30° and 60°, are consistent with the low friction coefficient measured during our experiments on marls (μ_s = 0.39), indicating that clay minerals exert a main control on fault mechanics. Moreover, our integrated analysis suggests that fracturing and faulting are the main mechanisms allowing fluid circulation within the low-permeability multilayer, and that its sealing integrity can be compromised only by the activity of larger faults cutting across its entire thickness.     

High pressure deformation in two-phase aggregates

We investigate the rheological behavior of multi-phase aggregates at high pressure and high temperature. Using synchrotron X-ray radiation as the probing tool, we are able to quantify the stress state of individual phases within the aggregates. This method provides fundamental information in interpreting the behavior of two phase/multi-phase mixtures, which contribute to our understanding of the deformation process at deep earth conditions. We choose MgAl₂O₄ spinel and MgO periclase as our model materials. Mixtures of various volume proportions were deformed in a multi-anvil high pressure deformation apparatus at pressure of 5 GPa and elevated temperatures. Stress is determined from X-ray diffraction, providing a measure of stress in each individual phase of the mixture in situ during the deformation. Macroscopic strain is determined from X-ray imaging. We compare the steady state strength of various mixtures at 1000 °C and 800 °C and at the strain rate in the range of 1.8 to 8.8 × 10^− 5 s^− 1. Our data indicate that the weak phase (MgO) is responsible for most of the accumulated strains while the strong phase (spinel) is supporting most of the stress when the volume proportion is 75% spinel and 25% MgO. The intermediate compositions (40/60) are much weaker than either of the end members, while the grain sizes for the intermediate compositions (submicrons) are much smaller than the end members (5–10 μm). We conclude that a change in flow mechanism resulting from these smaller grains is responsible for the low strength of the intermediate composition mixtures. This study demonstrates an approach of using synchrotron X-rays to study the deformation behaviors of multi-phase aggregates at high pressure and high temperature.     

Deformation inside a paleosubduction channel – Insights from microstructures and crystallographic preferred orientations of eclogites and metasediments from the Tauern Window,Austria

The Eclogite Zone, of the Tauern Window is an exhumed subduction channel comprising eclogites with different grades of retrogression in a matrix of high-pressure metasediments. The rocks were exposed to 600 °C and 20–25 kbars, and then retrogressed during their exhumation, first under blueschist facies and later under amphibolite facies metamorphism. To gain insights into the deformation within the subduction channel during subduction and exhumation, both fresh and retrogressed eclogites, as well as the surrounding metasediments were investigated with respect to their deformation microstructures and crystallographic preferred orientations (CPOs). Pristine and retrogressed eclogites show grain boundary migration and subgrain rotation recrystallization microstructures in omphacite. A misorientation axes analysis reveals the activity of complementary deformation mechanisms including grain boundary sliding and dislocation creep. The omphacite CPOs of the eclogites correspond to dominant SL-fabrics characteristic of plane strain deformation, though there are local variations towards flattening or constriction within the paleosubduction channel. The glaucophane CPOs in retrogressed eclogites match those of omphacite, suggesting that a constant strain geometry persisted during exhumation at blueschist facies conditions. Plastic deformation of the host high-pressure metasediments outlasted that of the eclogites, as indicated by white mica fabrics and quartz CPO. The latter is consistently asymmetric, pointing to the operation of non-coaxial deformation. The microstructures and CPO data indicate a continuous plastic deformation cycle with eclogite and blueschist facies metamorphism related to subduction and exhumation of the different rock units.     

Do cataclastic deformation bands form parallel to lines of no finite elongation (LNFE) or zero extension directions?

Conjugate cataclastic deformation bands cut unconsolidated sand and gravel at McKinleyville, California, and dip shallowly towards the north-northeast and south-southwest. The acute dihedral angle between the two sets of deformation bands is 47° and is bisected by the sub-horizontal, north-northeast directed incremental and finite shortening directions. Trishear models of fault propagation folding above the McKinleyville fault predict two sets of LNFE (lines of no finite elongation) that plunge steeply and shallowly to the south and north. These predictions are inconsistent with deformation band orientations and suggest that deformation bands did not form parallel to these LNFE. During plane strain, zero extension directions with acute dihedral angles of 47° develop when the dilatancy rate (dV/dε₁) is −4.3. Experimental dilatancy rates for Vosges sandstone (cohesion > 0) and unconsolidated Hostun sand suggest the deformation bands either developed parallel to zero extension directions or in accordance with the Mohr–Coulomb criterion, assuming initial porosities of 22% and 39%, respectively. An empirical relationship between dV/dε₁, relative density and mean stress suggests that dilatancy rates for Vosges sandstone overestimate dV/dε₁ at McKinleyville. Deformation bands at McKinleyville likely developed either in a Mohr–Coulomb orientation, or an intermediate orientation bounded by the Mohr–Coulomb (θ_C) and Roscoe (θ_R) angles.