Strain memory of 2D and 3D rigid inclusion populations in viscous flows — What is clast SPO telling us?

We model the development of shape preferred orientation (SPO) of a large population of two- and three-dimensional (2D and 3D) rigid clasts suspended in a linear viscous matrix deformed by superposed steady and continuously non-steady plane strain flows to investigate the sensitivity of clasts to changing boundary conditions during a single or superposed deformation events. Resultant clast SPOs are compared to one developed by an identical initial population that experienced a steady flow history of constant kinematic vorticity and reached an identical finite strain state, allowing examination of SPO sensitivity to deformation path. Rotation paths of individual triaxial inclusions are complex, even for steady plane strain flow histories. It has been suggested that the 3D nature of the system renders predictions based on 2D models inadequate for applied clast-based kinematic vorticity gauges. We demonstrate that for a large population of clasts, simplification to a 2D model does provide a good approximation to the SPO predicted by full 3D analysis for steady and non-steady plane strain deformation paths. Predictions of shape fabric development from 2D models are not only qualitatively similar to the more complex 3D analysis, but they display the same limitations of techniques based on clast SPO commonly used as a quantitative kinematic vorticity gauge. Our model results from steady, superposed, and non-steady flow histories with a significant pure shearing component at a wide range of finite strain resemble predictions for an identical initial population that experienced a single steady simple shearing deformation. We conclude that individual 2D and 3D clasts respond instantaneously to changes in boundary conditions, however, in aggregate, the SPO of a population of rigid inclusions does not reflect the late-stage kinematics of deformation, nor is it an indicator of the unique ‘mean’ kinematic vorticity experienced by a deformed rock volume.     

Analogue modelling of the influence of shape and particle/matrix interface lubrication on the rotational behaviour of rigid particles in simple shear

The rotational behaviour of a rigid particle embedded in a linear viscous matrix undergoing cylindrical simple shear (Couette) flow was studied in 2D rock-analogue experiments. The influence of particle shape (elliptical vs. monoclinic), aspect ratio and the nature of the matrix/particle interface (lubricated vs. unlubricated) was investigated. Both matrix (PDMS) and lubricant (liquid soap) were linear viscous, with a viscosity ratio of ca. 10⁴. Without lubricant, the rotational behaviour of all particles closely approximates the Jeffery theory. Lubricated monoclinic particles with the long diagonal initially parallel to the shear direction show back rotation and approach a stable position. Lubricated elliptical particles initially parallel to the shear direction also show back rotation but only transiently stabilize. Weak planar zones in the matrix adjacent to unlubricated elliptical particles do not induce backward rotation. In general for elliptical particles, rotation rate as a function of orientation depends on axial ratio and thickness of the lubricant mantle. For thick mantles (initially >10% of the volume of the particle), rotation rates are faster than Jeffery theory. For very thin mantles they are markedly slower compared with thick mantles, particularly when the long axis is nearly parallel to the shear direction. Rotation rates are never strictly zero, so true stabilization does not occur. However, for more elongate particles (axial RATIO=6) rotation rates are so slow that a very strong shape preferred orientation would develop in a lubricated elliptical particle population. In experiments, the volume of lubricant is constant and the thickness adjacent to the long side of the particle progressively decreases with increasing strain. In natural examples of porphyroclast systems, the weak mantle continually develops by recrystallization and/or cataclasis of the rigid clast core and a steady state between production and thinning could be attained, potentially leading to true stabilization for particles with a high axial ratio.     

2-D rotation behavior of a rigid ellipse in confined viscous simple shear: numerical experiments using FEM

The rotation behavior of rigid elliptical inclusions adherent to the viscous matrix in simple shear flow is investigated using a 2-D finite element numerical model. Several simulations were performed using different ratios (S) between shear zone width and inclusion's least principal axis. A computational strategy was devised to calculate pressure and viscous forces exerted on the inclusion and deduce its angular velocity. For large S values, results agree remarkably well with theoretical predictions, while for small S values results deviate significantly from theory but are in agreement with previous analogue experiments. The numerical model provided detailed and coherent information about the physical parameters involved in the process (e.g., pressure, strain rate and vorticity distributions within the model).     

Influence of a low-viscosity layer between rigid inclusion and viscous matrix on inclusion rotation and matrix flow: A numerical study

We have used 2-D finite element modelling to investigate the influence of a permanent low-viscosity layer between matrix and inclusion on matrix flow and inclusion rotation under viscous simple shear flow. Rigid inclusions of different shape (circle, square, ellipse, lozenge, rectangle and skewed rectangles) and aspect ratio (R) were used. The calculated matrix flow pattern is neither bow tie nor eye-shaped. It is a new flow pattern that we call cat eyes-shaped, which is characterized by: (i) straight streamlines that slightly bend inwards at the inclusion's crests; (ii) elongate eye-shaped streamlines on each side of the inclusion; (iii) stagnation points in the centre of the eyes; (iv) absence of closed streamlines surrounding the inclusion; (v) changes in flow configuration with inclusion orientation; the lines of flow reversal bend and tilt, closed streamline circuits may disappear, and streamlines may bend outwards at the inclusion's crests.Concerning inclusion rotation, the numerical results show that: (i) a low-viscosity layer (LVL) makes inclusions with R = 1 rotate synthetically, but the rotation rate depends upon shape (circle or square) and orientation. Therefore, shape matters in the slipping mode. (ii) All studied shapes with R > 1 rotate antithetically when starting with the greatest principal axis (e₁) parallel to the shear direction ( = 0°); (iii) rotation is limited because there is a stable equilibrium orientation (_se) for all studied shapes with R > 1. (iii) There is also an unstable equilibrium orientation (_ue), and both _se and _ue depend upon inclusion's R and shape.The present numerical results closely agree with previous results of analogue experiments with a permanent low viscosity interface. Only minor deviations related with small shape differences were detected.     

Pulsating oblate and prolate three-dimensional strains

The progressive ductile deformation of competent spherical inclusions is modeled analytically. Results of this study may help to understand better the limitations connected to geological field methods using competent inclusions for strain analyses. Parameters studied and quantified here are the strain magnitude, the progressive change in inclusion shape, the orientation of the finite strain axes, the frequency of pulsation, and the coupling between the strain ellipticity and viscosity contrast. Competent inclusions develop pulsating apparent strains if the host material is subjected to a component of simple shear and provided time or strain rate is sufficient to complete the strain cycle. The disparity between the strain magnitude inferred from competent viscous inclusions and that undergone by the host rock, increases for larger viscosity between them. The pulsation of the inclusion may suggest zero strain after a strain cycle has been completed, even though strain in the host rock is extremely large. The inclusion will develop pulsating oblate strains if a shortening rate is superposed normal to the plane of pulsation. Conversely, pulsating prolate strains occur if an extension rate is superposed instead of shortening. Stretching lineations outlined by deformed competent inclusions within shear zones beneath collapsing nappe sheets may even point perpendicular to the direction of nappe transport. This finding offers an explanation for the occurrence of mutually perpendicular pebble elongations in nearby locations within the Bygdin conglomerate beneath the Jotun nappe, Norwegian Caledonides.     

Influence of a permanent low-friction boundary on rotation and flow in rigid inclusion/viscous matrix systems from an analogue perspective

Natural occurrences and recent experimental work show that a low-friction inclusion/matrix boundary can be responsible for antithetical rotation and development of stable shape preferred orientations in simple shear. The flow of a viscous matrix around a rigid inclusion to which it may or may not be adherent is still not well studied, but it is relevant to the understanding of the behaviour of structural elements in mylonite zones. We used two-dimensional (2-D) analogue experiments to address these issues. The experimental results with a permanent low-friction inclusion/matrix boundary (nonadherent mode) show the following. (1) The rotation behaviour of inclusions in this mode is markedly different from the theoretical predictions and experimental results for the adherent mode. (2) Inclusions with aspect ratio equal to 1 rotate synthetically at a rate that depends upon inclusion shape and orientation. (3) Ellipse-, rectangle- and lozenge-shaped inclusions rotate antithetically when starting with their greatest axis parallel to the shear direction. (4) Back rotation is limited in all cases studied, and the angle between the inclusion greatest axis and the shear direction represents a stable orientation, which depends on inclusion aspect ratio and shape. (5) A metastable orientation exists, which is strongly dependent upon inclusion shape and aspect ratio, and separates fields of antithetic and synthetic rotation.Our experimental results show that the overall flow pattern is bow tie-shaped in adherent and nonadherent modes. However, there are major differences in the way the matrix flows near the inclusion. (i) In the nonadherent mode, the nearby matrix flows past the inclusion, whereas in the adherent mode, the nearby matrix flows with and follows the inclusion. Therefore, in the adherent mode, passive marker lines parallel to the shear direction and streamlines show considerable deflections at the inclusion crests, in marked contrast with their straight character in the nonadherent mode. (ii) Stagnation points or closed flow lines near the inclusion were not observed in the nonadherent mode, which means that there is no closed separatrix around the inclusion in this mode, despite the fact that the overall flow shape is bow tie. (iii) In the adherent mode, the line of flow reversal is stable throughout deformation, but in the nonadherent mode, it changes position and orientation with progressive shearing. This shifting of the line of flow reversal can be an important factor controlling rotation behaviour in the nonadherent mode. (iv) In the nonadherent mode, the inclusion behaviour is similar to that observed in confined flow. (v) The flow pattern in the nonadherent mode provides an explanation for the observed lack of drag folds associated with small-scale rigid inclusions in mylonites.     

Pulsating strains

Pulsating strains are modelled analytically using a continuum model with a well defined heterogeneous rheology. The mechanical continuum comprises single circular cylindrical inclusions of isotropic homogeneous viscosity hosted in a less viscous matrix. The analysis concentrates on the behaviour of the competent inclusion and demonstrates that the progressive deformation of the inclusion is entirely controlled by its viscosity ratio to the host and the orientation of the far field stress. The rate of deformation is determined by the stress magnitude and the viscosity of the inclusion, but does not affect the deformation path followed. An extreme case is provided by pure shear deformations where, by definition, the shear rate vanishes from the deformation tensor. Vorticity will then be zero, and strains will develop coaxially without pulsation.     

The influence of matrix rheology and vorticity on fabric development of populations of rigid objects during plane strain deformation

The influence of vorticity and rheology of matrix material on the development of shape-preferred orientation (SPO) of populations of rigid objects was experimentally studied. Experiments in plane strain monoclinic flow were performed to model the fabric development of two populations of rectangular rigid objects with object aspect ratios (R_ob) 2 and 3. The density of the rigid object populations was 14% of the total area. Objects were dispersed in a Newtonian and a non-Newtonian, power law matrix material with a power law exponent n of 1.2. The kinematic vorticity number (W_n) of the plane strain monoclinic flow was 1, 0.8 and 0.6 with finite simple shear strain of 4.6, 3.0 and 0.9, respectively. In experiments with R_ob=3, the SPO is strongly influenced by W_n and the material properties of the matrix. Deformation of a power law matrix material and low W_n resulted in a stronger SPO than deformation of a linear viscous matrix and high W_n. Strain localization coupled with particle interaction plays a significant role in the development of a shape-preferred orientation. High strain simple shear zones separate trains of rigid objects that are surrounded by low strain zones with W_n lower than the bulk W_n. In fabrics involving populations of objects with R_ob=2, rheology of the matrix materials does not systematically influence the intensity of the SPO.     

A general deformation matrix for three-dimensions

A deformation that is obtained by any simultaneous combination of two steady-state progressive deformations: simple shearing and a coaxial progressive deformation, involving or not a volume change, can be expressed by a single transformation, or deformation matrix. In the general situation of simple shearing in a direction non-orthogonal with the principal strains of the coaxial progressive deformation, this deformation matrix is a function of the strain components and the orientation of shearing. In this example, two coordinate systems are defined: one for the coaxial progressive deformation (x_i system), where the principal and intermediate strains are two horizontal coordinate axes, and another for the simple shear (x _i ^t’ system), with any orientation in space. For steady-state progressive deformations, from the direction cosines matrix that defines the orientation of shear strains in the x_i coordinate system, an asymmetric finite-deformation matrix is derived. From this deformation matrix, the orientation and ellipticity of the strain ellipse, or the strain ellipsoid for three-dimensional deformations, can be determined. This deformation matrix also can be described as a combination of a rigid-body rotation and a stretching represented by a general coaxial progressive deformation. The kinematic vorticity number (W _k is derived for the general deformation matrix to characterize the non-coaxiality of the three-dimensional deformation. An application of the deformation matrix concept is given as an example, analyzing the changes in orientation and stretching that variously-oriented passive linear markers undergo after a general two-dimensional deformation. The influence of the kinematic vorticity number, the simple and pure shear strains, and the obliquity between the two deformation components, on the linear marker distribution after deformation is discussed.     

Shape preferred orientation of object populations: automatic analysis of digitized images

An automatic technique is presented for systematic shape preferred orientation (SPO) analysis of crystal fabrics in rocks. This technique is based on digitization of an image with a CCD camera and separation of an object population from its matrix. The image is analysed globally using a rotating grid that is superimposed on the square pixel grid of the computer screen. The parametric method is based on counting intercepts in all directions. The interceps rose provides information on both the orientation and strength of the SPO. Problems of hardware-related counting noise anisotropy are avoided using an original method of filtering the intercept function with a linear, sliding mask on the counting grid. The technique was tested both on single objects of different shapes and on different populations of identical objects. A very high degree of accuracy was obtained for orientation measurements. Fabric shape parameters are introduced which are specific to this intercept technique; they compare favourably with shape parameters used in axial fabric analysis.     

Numerical evaluation of three non‐coaxial kinematic models using the distinct element method for elliptical granular materials

This paper presents a numerical evaluation of three non‐coaxial kinematic models by performing Distinct Element Method (DEM) simple shear tests on specimens composed of elliptical particles with different aspect ratios of 1.4 and 1.7. The models evaluated are the double‐shearing model, the double‐sliding free‐rotating model and the double slip and rotation rate model (DSR² model). Two modes of monotonic and cyclic simple shear tests were simulated to evaluate the role played by the inherent anisotropy of the specimens. The main findings are supported by all the DEM simple shear tests, irrespective of particle shape, specimen density or shear mode. The evaluation demonstrates that the assumption in the double‐shearing model is inconsistent with the DEM results and that the energy dissipation requirements in the double‐sliding free‐rotating model appear to be too restrictive to describe the kinematic flow of elliptical particle systems. In contrast, the predictions made by the DSR² model agree reasonably well with the DEM data, which demonstrates that the DSR² model can effectively predict the non‐coaxial kinematic behavior of elliptical particle systems. Copyright © 2016 John Wiley & Sons, Ltd.     

Using AMS combined with mineral shape preferred orientation analysis to understand the emplacement fabrics of the Apiaí gabbro-norite (Ribeira Belt,SE Brazil)

The Apiaí gabbro-norite is a massive fine-grained Neoproterozoic intrusion emplaced in a core of synformal structure that deforms low-grade marine metasedimentary rocks of the Ribeira Belt of south-eastern Brazil. The lack of visible magmatic layering or any internal fabric has been a major limitation in deciding whether the emplacement occurred before or after the regional folding. To assist in the tectonic interpretations, we combine low-field anisotropy of magnetic susceptibility (AMS) and silicate shape preferred orientation (SPO) to reveal the internal structure of the mafic intrusion. Magnetic data indicate a mean susceptibility of about 10⁻² SI and a mean anisotropy degree (P) of about 1.08, essentially yielded by titanomagnetite. The magnetic and silicate foliations for P ≥ 1.10 are parallel to each other, while the lineations tend to scatter on the foliation plane, in agreement with the dominant oblate symmetry of the AMS and SPO ellipsoids. For lower P values, the magnetic and silicate fabrics vary from coaxial to oblique, and for P ≤ 1.05, their shapes and orientations can be quite distinct. The crystal size distribution (CSD) of plagioclase for P > 1.05 is log linear, in agreement with a bulk simple crystallisation history. These results combined show that for a strong SPO, corresponding to a magnetic anisotropy above 1.10, AMS is a reliable indicator of the magmatic fabric. They indicate that the Apiaí gabbro-norite consists of sill-like body that was inclined gently to the north by the regional folding.     

富钙长石-橄榄石包体与其他部分典型包体W-L边的成因

球粒陨石中的富Ca、Al包体(简称CAI)形成于星云演化的最初始阶段,保存了大量星云形成和演化的各种信息。研究认为,包体的成因主要包括星云直接凝聚和熔融结晶,少部分甚至经历过高温蒸发过程。部分CAI最外层具有由一种或几种矿物组成的Warking-Lovering边(简称为W-L边),CAI和其W-L边对于认识早期星云环境和界定CAI的形成时间等均具有重要意义。目前,对于W-L边的形成过程研究并不深入,且一直存在争议。本文主要介绍了三个典型包体:C#1(富钙长石-橄榄石包体)、GRV 022459-2RI5(A型包体)和GRV 021579-3RI5(富尖晶石球粒状包体)及其W-L边的矿物岩石学和氧同位素组成特征。C#1包体明显经历过熔融结晶过程,W-L边氧同位素组成具有与包体内部矿物相似的富~(16)O同位素特征,表明W-L边的成因与包体的形成过程密切相关,形成于同一富~(16)O同位素组成区域,且W-L边属于包体熔融结晶过程后期的产物。矿物岩石学特征表明,GRV 022459-2RI5属于星云直接凝聚形成,其W-L边为包体形成过程最晚期星云凝聚产物。GRV021579-3RI5经历过熔融结晶过程,其W-L边为包体结晶最后阶段的产物。     

Refractory inclusions in the Murchison meteorite

Mineralogical and petrographic studies of a wide variety of refractory objects from the Murchison C2 chondrite have revealed for the first time melilite-rich and feldspathoid-bearing inclusions in this meteorite, but none of these is identical to any inclusion yet found in Allende. Blue spinel-hibonite spherules have textures indicating that they were once molten, and thus their SiO₂-poor bulk composition requires that they were exposed to higher temperatures (>1550°C) than those deduced so far from any Allende inclusion. Melilite-rich inclusions are similar to Allende compact Type A's, but are more Al-, Ti-rich. One inclusion (MUCH-1) consists of a delicate radial aggregate of hibonite crystals surrounded by alteration products, and probably originated by direct condensation of hibonite from the solar nebular vapor. The sinuous, nodular and layered structures of another group of inclusions, spinel-pyroxene aggregates, suggest that these also originated by direct condensation from the solar nebular gas. Each type of inclusion is characterized by a different suite of alteration products and/or rim layers from all the other types, indicating modification of the inclusions in a wide range of different physico-chemical environments after their primary crystallization. All of these inclusions contain some iron-free rim phases. These could not have formed by reaction of the inclusions with fluids in the Murchison parent body because the latter would presumably have been very rich in oxidized iron. Other rim phases and alteration products could have formed at relatively low temperatures in the parent body, but some inclusions were not in the locations in which they were discovered when this took place. Some of these inclusions are too fragile to have been transported from one region to another in the parent body, indicating that low temperature alteration of these may have occurred in the solar nebula.     

Influence of shape fabric and crystal texture on marble degradation phenomena: simulations

Microstructure-based finite element simulations were used to study the influence of grain shape fabric and crystal texture on thermoelastic responses related to marble degradation phenomena. Calcite was used as an illustrative example for studying extremes of shape preferred orientation (SPO) in shape fabric and lattice preferred orientation (LPO) in crystal texture. Three SPOs were analyzed: equiaxed grains, elongated grains, and a mixture of equiaxed and elongated grains. Three LPOs were considered: a random orientation distribution function and two degrees of strong directional crystal texture. Finally, the correlation between the direction of the LPO with respect to that of the SPO was examined. Results show that certain combinations of SPO, LPO, and their directional relationship have significant influence on the thermomechanical behavior of marble. For instance, while there is no major dependence of the elastic strain energy density and the maximum principal stress on SPO for randomly textured microstructures, there is a strong synergy between LPO and its directional relationship with respect to the SPO direction. Microcracking precursors, elastic strain energy density, and maximum principal stress, decrease when the crystalline c-axes have fiber texture perpendicular to the SPO direction, but increase significantly when the c-axes have fiber texture parallel to the SPO direction. Moreover, the microstructural variability increases dramatically for these latter configurations. In general, the influence of LPO was as expected, namely, the strain energy density and the maximum principal stress decreased with more crystal texture, apart from for the exception noted above. Spatial variations of these precursors indicated regions in the microstructure with a propensity for microcracking. Unexpectedly, important variables were the microstructural standard deviations of the spatial distributions of the microcracking indicators. These microstructural standard deviations were as large as or larger than the variables themselves. The elastic misfit-strain contributions to the coefficients of thermal expansion were also calculated, but their dependence was as expected.     

2D rotation of rigid inclusions in confined bulk simple shear flow: a numerical study

We have used incompressible Navier–Stokes in 2D finite element modelling to investigate rigid inclusion rotation under confined bulk simple shear flow. Confinement is defined as the ratio (S) between the channel width (H) and the inclusion's least axis (e₂)(S=H/e₂). The numerical results show that (i) inclusion rotation is strongly influenced by S and, when the confinement is effective, aspect ratio (R) and shape also play an important role. (ii) Back rotation is limited because inclusions reach a stable equilibrium orientation (_se). (iii) There is also an unstable equilibrium orientation (_ue), which defines an antithetic rotation field with _se, and both _se and _ue depend on S, and inclusion R and shape.     

The initiation and development of metamorphic foliation in the Otago Schist, Part 2: evidence from quartz grain-shape data

Shape, size and orientation measurements of quartz grains sampled along two transects that cross zones of increasing metamorphic grade in the Otago Schist, New Zealand, reveal the role of quartz in the progressive development of metamorphic foliation. Sedimentary compaction and diagenesis contributed little to the formation of a shape‐preferred orientation (SPO) within the analysed samples. Metamorphic foliation was initiated at sub‐greenschist facies conditions as part of a composite S1‐bedding structure parallel to the axial planes of tight to isoclinal F1 folds. An important component of this foliation is a pronounced quartz SPO that formed dominantly by the effect of dissolution–precipitation creep on detrital grains in association with F1 strain. With increasing grade, the following trends are evident from the SPO data: (i) a progressive increase in the aspect ratio of grains in sections parallel to lineation, and the development of blade‐shaped grains; (ii) the early development of a strong shape preferred orientation so that blade lengths define the linear aspect of the foliation (lineation) and the intermediate axes of the blades define a partial girdle about the lineation; (iii) a slight thinning and reduction in volume of grains in the one transect; and (iv) an actual increase in thickness and volume in the survivor grains of the second transect. The highest‐grade samples, within the chlorite zone of the greenschist facies, record segregation into quartz‐ and mica‐rich layers. This segregation resulted largely from F2 crenulation and marks a key change in the distribution, deformation and SPO of the quartz grains. The contribution of quartz SPO to defining the foliation lessens as the previously discrete and aligned detrital quartz grains are replaced by aggregates and layers of dynamically recrystallized quartz grains of reduced aspect ratio and reduced alignment. Pressure solution now affects the margins of quartz‐rich layers rather than individual grains. In higher‐grade samples, therefore, the rock structure is characterized increasingly by segregation layering parallel to a foliation defined predominantly by mica SPO.     

Synthetic fluid inclusions - VI. Quantitative evaluation of the decrepitation behaviour of fluid inclusions in quartz at one atmosphere confining pressure

ABSTRACT The decrepitation behaviour of fluid inclusions in quartz at one atmosphere confining pressure has been evaluated using pure H₂O synthetic inclusions formed by healing fractures in natural quartz. Three different modes of non-elastic deformation, referred to as stretching, leakage or partial decrepitation, and total decrepitation have been observed. The internal pressure required to initiate non-elastic deformation is inversely related to inclusion size according to the equation: internal pressure (kbar) = 4.26 D^-0.423 where D is the inclusion diameter in microns. Regularly shaped inclusions require a higher internal pressure to initiate non-elastic deformation than do irregularly shaped inclusions of similar size. Heating inclusions through the α/β quartz inversion results in mechanical instability in the quartz crystal and leads to mass decrepitation of inclusions owing to structural mismatches generated by pressure gradients in the quartz around each inclusion. Long-term heating experiments (∼2 years) suggest that the internal pressure required to initiate non-elastic deformation does not decrease significantly with time and indicates that short-lived thermal fluctuations in natural systems should not alter the inclusion density and homogenization temperature. Inclusions that do exhibit decreased density (higher homogenization temperature) are, however, always accompanied by a change in shape from irregular to that of a negative crystal. Observations of this study are consistent with elasticity theory related to fracture generation and propagation around inclusions in minerals. These results indicate that an inclusion will not be influenced by a neighbouring inclusion, or other defect in the host phase, as long as the distance between the two is >2–4 diameters of the larger of the two inclusions.     

Fabric analysis of quartzites with negative magnetic susceptibility – Does AMS provide information of SPO or CPO of quartz?

We ask the question whether petrofabric data from anisotropy of magnetic susceptibility (AMS) analysis of deformed quartzites gives information about shape preferred orientation (SPO) or crystallographic preferred orientation (CPO) of quartz. Since quartz is diamagnetic and has a negative magnetic susceptibility, 11 samples of nearly pure quartzites with a negative magnetic susceptibility were chosen for this study. After performing AMS analysis, electron backscatter diffraction (EBSD) analysis was done in thin sections prepared parallel to the K₁K₃ plane of the AMS ellipsoid. Results show that in all the samples quartz SPO is sub-parallel to the orientation of the magnetic foliation. However, in most samples no clear correspondance is observed between quartz CPO and K₁ (magnetic lineation) direction. This is contrary to the parallelism observed between K₁ direction and orientation of quartz c-axis in the case of undeformed single quartz crystal. Pole figures of quartz indicate that quartz c-axis tends to be parallel to K₁ direction only in the case where intracrystalline deformation of quartz is accommodated by prism <c> slip. It is therefore established that AMS investigation of quartz from deformed rocks gives information of SPO. Thus, it is concluded that petrofabric information of quartzite obtained from AMS is a manifestation of its shape anisotropy and not crystallographic preferred orientation.     

Flow kinematics in the deeper part of a subduction channel, as inferred from inclusion trails in plagioclase porphyroblasts from the Sambagawa metamorphic rocks, south-west Japan

The subduction and exhumation of accretionary prism metasedimentary rocks are accompanied by large‐strain ductile deformations which may be recorded in microstructures. Porphyroblast microstructures have been a key to unravel the kinematics in such deformed belts. Shape‐preferred orientation (SPO) of epidote and amphibole inclusions that define S‐shaped trails in prograde cores of plagioclase porphyroblasts were analysed from the high‐P/T Sambagawa metamorphic rocks. Inclusions are found to be elongate parallel to the [010] and [001] directions, respectively, and their long‐axis orientations define an internal foliation S_i (best‐fit great circle) and lineation L_i (maximum on the S_i). S‐shaped inclusion trails in the orthogonal sections do not exhibit the same geometries, but rather are grouped into two types, where the foliation intersection axes (FIAs) are nearly perpendicular and parallel to L_i, respectively. These two types of S‐shaped inclusion trails are seen in the sections inclined at low and high angles to the L_i, respectively. However, the latter type commonly consists of composite trails, where the S_i is first rotated about an FIA perpendicular to the L_i (i.e. unique axis), and then about an FIA parallel to the L_i. The S‐shaped inclusion trails are interpreted to have formed by the successive overgrowth of matrix minerals and rotation of the plagioclase porphyroblast cores about a unique axis in non‐coaxial deformation. The rotation of S_i about an FIA nearly parallel to the L_i is perhaps an apparent rotation, caused by the deflection of foliation around the growing prismatic plagioclase grain prior to inclusion into the porphyroblast. This study has for the first time documented the 3‐D geometry of S‐shaped inclusion trails in porphyroblasts from accretionary prism metasedimentary rocks and identified their origin, which helps to understand the flow kinematics in the deeper part of a subduction channel.