The importance of the effective intermediate principal stress (σ′2) to fault slip patterns

In normal faulting regimes, the magnitudes and orientations of the maximum and minimum principal compressive stresses may be known with some confidence. However, the magnitude of the intermediate principal compressive stress is generally much more difficult to constrain and is often not considered to be an important factor. In this paper, we show that the slip characteristics of faults and fractures with complex or nonoptimal geometry are highly sensitive to variation or uncertainty in the ambient effective intermediate principal stress (σ′₂). Optimally oriented faults and fractures may be less sensitive to such variations or uncertainties. Slip tendency (T_s) analysis provides a basis for quantifying the effects of uncertainty in the magnitudes and orientations of all principal stresses and in any stress regime, thereby focusing efforts on the most important components of the system. We also show, for a normal faulting stress regime, that the proportion of potential surfaces experiencing high slip tendency (e.g., T_s ≥ 0.6) decreases from a maximum of about 38% where σ′₂ = σ′₃, to a minimum of approximately 14% where σ′₂ is halfway between σ′₃ and σ′₁, and increases to another high of approximately 29% where σ′₂ = σ′₁. This analysis illustrates the influence of the magnitude of σ′₂ on rock mass strength, an observation previously documented by experimental rock deformation studies. Because of the link between fault and fracture slip characteristics and transmissivity in critically stressed rock, this analysis can provide new insights into stress-controlled fault transmissivity.     

Geometry of the fault zone of the 2003 Ms = 7.5 Chuya earthquake and associated stress fields, Gorny Altai

The co-seismic deformations produced during the September 27, 2003 Chuya earthquake (Ms = 7.5) that affected the Gorny Altai, Russia, are described and discussed along a 30 km long segment. The co-seismic deformations have manifested themselves both in unconsolidated sediments as R- and R′-shears, extension fractures and contraction structures, and in bedrock as the reactivation of preexisting schistosity zones and individual fractures, as well as development of new ruptures and coarse crushing zones. It has been established that the pattern of earthquake ruptures represents a typical fault zone trending NW–SE with a width reaching 4–5 km and a dextral strike–slip kinematics. The initial stress field that produced the whole structural pattern of co-seismic deformations during the Chuya earthquake, is associated with a transcurrent regime with a NNW–SSE, almost N–S, trending of compressional stress axis (σ₁), and a ENE–WSW, almost E–W, trending of tensional stress axis (σ₃). The state of stress in the newly-formed fault zone is relatively uniform. The local stress variations are expressed in insignificant deviation of σ₁ from N–S to NW–SE or NE–SW, in short-term fluctuations of relative stress values in keeping their spatial orientations, or in a local increase of the plunge angle of the σ₁. The geometry of the fault zone associated with the Chuya earthquake has been compared with the mechanical model of fracturing in large continental fault zones with dextral strike–slip kinematics. It is apparent that the observed fracture pattern corresponds to the late disjunctive stage of faulting when the master fault is not fully developed but its segments are already clearly defined. It has been shown that fracturing in widely different rocks follows the common laws of the deformation of solid bodies, even close to the Earth surface, and with high rates of movements.     

Investigations on the mechanical significance of crenulation cleavage

Measurements of total, incremental and progressive strains associated with the development of small scale crenulation cleavage in some low-grade metamorphic rocks from Australia and Switzerland are applied to a discussion of the mechanical significance of the cleavage.Limits are placed on the amount of incremental and total slip or simple shear possible along the cleavage by the observation that the XY principal plane trace of bulk total crenulation strain coincides within 4° of the crenulation cleavage trace in all cases where this strain has been measured or estimated. The measurements are made on eight specimens using deformed porphyroblasts, crystal fibres in pressure-shadows around pyrite and flattened folds and include deformations with coaxial and non-coaxial histories.Further measurements derived from pressure-shadow fibres (eight specimens) show that the style and orientation of incremental deformation are essentially independent of the crenulation cleavage, except for a limit (43°) to the obliquity of the principal incremental extension axis during a given cleavage episode. The only special deformation related to the cleavage is the coaxial one. An indication of passive cleavage behaviour at high strain is shown by the progressive strain history of one specimen. Evidence for passive rotation of a transected axial plane is shown by another. A model is proposed to account for these observations, especially the conditions necessary for initiation and continued development of a new cleavage fabric.Some further applications of existing strain measurement techniques are described: of the R_f/Ø_f method to heterogeneously superposed tectonic strains and of an improved procedure of t′_α/α flattening analysis.     

Brittle-plastic deformation in initially dry rocks at fluid-present conditions: transient behaviour of feldspar at mid-crustal levels

We present detailed microstructural and chemical analyses from an initially dry anorthositic rock deformed during wet amphibolite facies conditions. Three different domains representing the microstructural variation of the deformed samples are investigated in detail in terms of fracture morphology and mode, grain characteristics and chemistry of present phases. Results show transient deformational behaviour where a close interaction between brittle, plastic and fluid-assisted deformation mechanisms can be observed. Our analysis allows us to describe the succession, interrelationships and effects of active mechanisms with progressively increasing strain in three so-called stages. In Stage 1, initial fracturing along cleavage planes promoted fluid influx that caused fragmentation and chemical reactions, producing fine-grained mineral assemblages in the fractures. Deformation twins and dislocations developed in clast pieces due to stress relaxation. Passive rotation of conjugate fracture sets and interconnection of intracrystalline fractures formed micro-shear-zones, constituting Stage 2. Microstructures and grain relationships indicate the activity and fluctuation between fracturing, dissolution-precipitation creep, grain boundary sliding and locally dislocation creep, reflecting the transient behaviour of brittle and plastic deformation mechanisms. Further rotation and widening of fractures into overall foliation parallel shear-bands (Stage 3) promoted strain partitioning into these areas through increased fluid influx, influence of fluid-assisted grain boundary sliding, phase mixing and presence of weak phases such as white mica. We suggest that local differences in fluid availability, volume fraction of weak phases produced by fluid present metamorphic reactions coupled with volume increase and local variations in stress concentration induced transient brittle-plastic behaviour. The studied shear-zone represents an example of the transformation of a rigid dry rock to a soft wet rock during deformation through syntectonic fracturing.     

Determination of Shear Strength and Three-dimensional Yield Strength for the Hoek-Brown Criterion

Summary. Most currently used techniques for analysing the stability of near surface structures, such as rock slopes, are based on the application of the effective Coulomb shear strength parameters cohesion c′, and the angle of friction φ′ on some known or anticipated shear surface subjected to an effective normal stress σ′_n. The most widely used of these techniques are the variants of the method of slices and related upper bound techniques. If the Hoek-Brown criterion is to be used to model the strength of near surface fractured rocks, it is necessary to determine equivalent Coulomb shear strength parameters for the specified level of effective normal stress. Calculation of the equivalent Coulomb parameters for the Hoek-Brown criterion for cases when a ≠ 0.5 is not a straightforward matter. A simple procedure for calculating instantaneous values of c′_i and φ′_i has been developed based on spreadsheet calculations and the application of a numerical optimisation routine. This procedure can also be applied to calculating the Hoek-Brown envelope plotted in shear stress/normal stress space. A simple closed form solution for c′_i and tan φ′_i has also been developed for the special case when a = 1. A three-dimensional version of the Hoek-Brown criterion has been developed by combining it with the Drucker-Prager criterion. This new yield criterion has been implemented by numerical solution of the governing equations. A simplification of this three-dimensional yield criterion has been developed by introducing an intermediate principal stress weighting factor. Comparison with published results demonstrates that this simplified criterion has the capacity to model the results of true triaxial tests for a range of different rock types over a wide range of stress levels. The new three-dimensional yield criterion has the advantage that its input parameters can be determined from routine uniaxial compression tests and mineralogical examination.     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

European continuous active tectonic strain–stress map

This paper shows a new continuous strain–stress map for Europe obtained from the direct inversion of earthquake focal mechanisms calculated from the centroid tensor method. A total of 1608 focal mechanisms have been selected with several quality criteria from different catalogues (CMT Harvard, ETH, Med-Net, I.G.N. and I.A.G.) from 1973 to the present day. Values for the maximum horizontal shortening direction and brittle strain–stress regime defined by the k′ ratio (e_y/e_z, horizontal maximum/vertical strain) have been calculated following in Europe and Pannonian Basin the slip model of tri-axial deformation. The individual results including Dey and the shape of the active brittle strain ellipsoid have been interpolated to a final 15′ regular grid taking into account the relationship between the tectonic horizontal strain–stress value and the vertical load. Both continuous strain regime and maximum horizontal shortening (Dey) maps show a good correlation with the primary tectonic forces generated along the plate boundaries, plate kinematics and also some local perturbations related with main crustal heterogeneities and topography, as well as significant spatial variations in integrated crustal strength.     

Dislocation generation, slip systems, and dynamic recrystallization in experimentally deformed plagioclase single crystals

Three samples of gem quality plagioclase crystals of An60 were experimentally deformed at 900 °C, 1 GPa confining pressure and strain rates of 7.5–8.7×10⁻⁷ s⁻¹. The starting material is effectively dislocation-free so that all observed defects were introduced during the experiments. Two samples were shortened normal to one of the principal slip planes (010), corresponding to a “hard” orientation, and one sample was deformed with a Schmid factor of 0.45 for the principal slip system [001](010), corresponding to a “soft” orientation. Several slip systems were activated in the “soft” sample: dislocations of the [001](010) and 110(001) system are about equally abundant, whereas 110{111} and [101] in ( 31) to ( 42) are less common. In the “soft” sample plastic deformation is pervasive and deformation bands are abundant. In the “hard” samples the plastic deformation is concentrated in rims along the sample boundaries. Deformation bands and shear fractures are common. Twinning occurs in close association with fracturing, and the processes are clearly interrelated. Glissile dislocations of all observed slip systems are associated with fractures and deformation bands indicating that deformation bands and fractures are important sites of dislocation generation. Grain boundaries of tiny, defect-free grains in healed fracture zones have migrated subsequent to fracturing. These grains represent former fragments of the fracture process and may act as nuclei for new grains during dynamic recrystallization. Nucleation via small fragments can explain a non-host-controlled orientation of recrystallized grains in plagioclase and possibly in other silicate materials which have been plastically deformed near the semi-brittle to plastic transition.     

The strain history of a dalradian slide: Using pebbles with low fluctuations in axis orientation

The strain associated with the Horn Head Slide, a major tectonic break in the Dalradian rocks of NW Ireland, is recorded by pebbles in an adjacent quartzite horizon. Mean X/Y ratios of the deformed pebble shapes in excess of 8.0 are seen closest to the slide and the field of three-dimensional shapes lies along the K = 1 line. The usual methods of separating initial shape ratio (R_i) and tectonic strain ratio (R_s) from the deformed shape ratio (R_f) of ellipsoidal markers are based on measurements of variation in fluctuation (e.g., the (R_f/φ technique). However, due to the high X/Y strains in this situation and since the pebbles initially lay parallel to bedding and to a principal plane of the finite strain ellipsoid, fluctuation is generally very low. Thus, except for the least deformed X/Y data, the R_f/φ technique is inapplicable and other methods are used. For X/Y data with mean (R_f > 4.0: R_s is calculated as the harmonic mean of R_f; maximum R_i values only are obtained from the range of R_f data. For all Y/Z and X/Z data: R_i is calculated from the logarithmic range (ω_log) of R_f; R_s is simply obtained from the geometric mean of R_f modified by R_i. It is concluded from this that a varying prolate tectonic strain (K - 1.5) reaching X/Y values in excess of 8.00 was coaxially superimposed on an initial oblate shape fabric to produce the present field of deformed pebbles in the quartzite near the slide.     

Steeply-dipping extension fractures in the Newark basin, New Jersey

Late Triassic and Early Jurassic bedrock in the Newark basin is pervasively fractured as a result of Mesozoic rifting of the east–central North American continental margin. Tectonic rifting imparted systematic sets of steeply-dipping, en échelon, Mode I, extension fractures in basin strata including ordinary joints and veins. These fractures are arranged in transitional-tensional arrays resembling normal dip-slip shear zones. They contributed to crustal stretching, sagging, and eventual faulting of basin rift deposits. Extension fractures display progressive linkage and spatial clustering that probably controlled incipient fault growth. They cluster into three prominent strike groups correlated to early, intermediate, and late-stage tectonic events reflecting about 50– 60° of counterclockwise rotation of incremental stretching directions. Finite strain analyses show that extension fractures allowed the stretching of basin strata by a few percent, and these fractures impart stratigraphic dips up to a few degrees in directions opposing fracture dips. Fracture groups display three-dimensional spatial variability but consistent geometric relations. Younger fractures locally cut across and terminate against older fractures having more complex vein-cement morphologies and bed-normal folds from stratigraphic compaction. A fourth, youngest group of extension fractures occur sporadically and strike about E–W in obliquely inverted crustal blocks. A geometric analysis of overlapping fracture sets shows how fracture groups result from incremental rotation of an extending tectonic plate, and that old fractures can reactivate with oblique slip components in the contemporary, compressive stress regime.     

Dynamic recrystallization near the brittle-plastic transition in naturally and experimentally deformed quartz aggregates

The electron backscattering diffraction technique (EBSD) was used to analyze bulging recrystallization microstructures from naturally and experimentally deformed quartz aggregates, both of which are characterized by porphyroclasts with finely serrated grain boundaries and grain boundary bulges set in a matrix of very fine recrystallized grains. For the Tonale mylonites we investigated, a temperature range of 300–380 °C, 0.25 GPa confining pressure, a flow stress range of ~ 0.1–0.2 GPa, and a strain rate of ~ 10^− 13 s^− 1 were estimated. Experimental samples of Black Hills quartzite were analyzed, which had been deformed in axial compression at 700 °C, 1.2–1.5 GPa confining pressure, a flow stress of ~ 0.3–0.4 GPa, a strain rate of ~ 10^− 6 s^− 1, and to 44% to 73% axial shortening. Using orientation imaging we investigated the dynamic recrystallization microstructures and discuss which processes may contribute to their development. Our results suggest that several deformation processes are important for the dismantling of the porphyroclasts and the formation of recrystallized grains. Grain boundary bulges are not only formed by local grain boundary migration, but they also display a lattice misorientation indicative of subgrain rotation. Dynamic recrystallization affects especially the rims of host porphyroclasts with a hard orientation, i.e. with an orientation unsuitable for easy basal slip. In addition, Dauphiné twins within porphyroclasts are preferred sites for recrystallization. We interpret large misorientation angles in the experimental samples, which increase with increasing strain, as formed by the activity of fluid-assisted grain boundary sliding.     

The initiation and linkage of surface fractures above a buried strike-slip fault: An experimental approach

Surface fractures in the overburdened sedimentary rocks, formed above a deep-seated basement fault, often provide important information about the kinematics of the underlying master fault. It has already been established that these surface fractures dynamically evolve and link one another with progressive displacement on the master fault below. In the present study, two different series of riedel-type experiments were carried out with clay analogue models under different boundary conditions viz., (i) heterogeneous simple shear of the cover rocks above a buried strike slip fault (wrench system) and (ii) heterogeneous simple shear with a component of shear-normal compression of the overburden package above a basement fault (transpressional system), to observe the initiation and linkage of surface fractures with varying T′ (where T′ = thickness of the overburden normalized with respect to the width of the master fault). In the wrench system, Riedel (R) shears were linked by principal displacement (Y) shears at early stages (shear strain of 0.8 to 1) in thin (2 < T′ < 18) models and finally (at a minimum shear strain of 1.4) gave rise to a through-going fault parallel to the basement fault without development of any other fracture. Conjugate Riedel (R′) shears develop only within the thicker (T′ > 18) clay models at a minimum shear strain of 0.7. With increasing deformation (at a minimum shear strain of 1.2) two R′ shears were joined by an R shear and finally opened up to make a sigmoidal vein with an asymmetry antithetic to the major fault-movement sense. Under transpression, the results were similar to that of heterogeneous simple shear for layers 2 < T′ < 15. In layers of intermediate thickness (15 < T′ < 25) early formed high angle R shears were cross cut by low angle R shears (at a minimum shear strain of 0.5 and shortening of 0.028) and “Riedel-within-Riedel” shears were formed within thick (T′ > 25) models (at minimum shear strain of 0.7 and shortening of 0.1), with marked angularity of secondary fault zone with the master fault at depth.     

Deformation pattern in a heterogeneous material: Folded and cleaved sedimentary cover immediately overlying a crystalline basement (Oisans, French Alps)

Particularly well exposed structures in folded and cleaved sedimentary cover immediately overlying a crystalline basement have been studied. Chemical analysis (X.R.F. and microprobe) reveal pressure solution process and give the possibilities of measurement of mass transfer. Study of fluid inclusion veins has determined the temperature pressure conditions: thermal effect of the basement and decrease of temperature and pressure with the age of various synkinematic veins.Characteristic examples of the behaviour of a heterogeneous material during coaxial and non-coaxial deformation are shown:

1. (1) Successive different asymmetrical folds, various cleavages and fractures appear in a shear zone parallel to the main fabric with variations of thickness and rock behaviour.

2. (2) Evolution of cleavage in such a shear zone (with or without slipping) is linked to the relations between the rotation of contraction direction and the rate of the cleavage process.

3. (3) Fold axes changed from the horizontal y direction to the vertical (or E—W transversal to the crystalline massif) X direction, with increase of the (X/Z) and (X/Y) ratios (obtained by fossils and reconstructed fold shape). This strain is always heterogeneous and the most deformed zone frequently evolves to discontinuities with slip.

4. (4) Indentation exists on all scale: from hard object (100 μ, with parenthesis form of pressure solution cleavage apparent on map distribution of various element) to basement block (with variation of strain value in the indented cover).

A model of the evolution of the deformation of sedimentary cover immediately overlying a crystalline basement is given in conclusion.     

Talc friction in the temperature range 25°–400 °C: Relevance for Fault-Zone Weakening

Talc is one of the weakest minerals that is associated with fault zones. Triaxial friction experiments conducted on water-saturated talc gouge at room temperature yield values of the coefficient of friction, μ (shear stress, τ/effective normal stress, σ′_N) in the range 0.16–0.23, and μ increases with increasing σ′_N. Talc gouge heated to temperatures of 100°–400 °C is consistently weaker than at room temperature, and μ < 0.1 at slow strain rates in some heated experiments. Talc also is characterized by inherently stable, velocity-strengthening behavior (strength increases with increasing shear rate) at all conditions tested. The low strength of talc is a consequence of its layered crystal structure and, in particular, its very weak interlayer bond. Its hydrophobic character may be responsible for the relatively small increase in μ with increasing σ′_N at room temperature compared to other sheet silicates.Talc has a temperature–pressure range of stability that extends from surficial to eclogite-facies conditions, making it of potential significance in a variety of faulting environments. Talc has been identified in exhumed subduction zone thrusts, in fault gouge collected from oceanic transform and detachment faults associated with rift systems, and recently in serpentinite from the central creeping section of the San Andreas fault. Typically, talc crystallized in the active fault zones as a result of the reaction of ultramafic rocks with silica-saturated hydrothermal fluids. This mode of formation of talc is a prime example of a fault-zone weakening process. Because of its velocity-strengthening behavior, talc may play a role in stabilizing slip at depth in subduction zones and in the creeping faults of central and northern California that are associated with ophiolitic rocks.     

The relationship between regional stress field,fracture orientation and depth of weathering and implications for groundwater prospecting in crystalline rocks

In a uniform granite gneiss study area in central Zimbabwe, lineaments oriented parallel to the maximum regional compressive stress orientation exhibit the thickest regolith development, while lineaments oriented perpendicular to the maximum compressive stress show the shallowest development of weathered regolith. The principal fracture set orientations were mapped using aerial imagery. The regional stress field, estimated from global stress maps, was used to determine the stresses acting on each principal lineament orientation. Multi-electrode resistivity profiling was carried out across fractures with different orientations to determine their subsurface regolith conditions. The results indicate that the 360 and 060° lineaments, which are sub-parallel to the principal compressive stress orientation (σ₁) exhibit maximum development of the regolith, while 130° lineaments perpendicular to σ₁ do not exhibit significant regolith development. Since regolith thickness has been positively correlated with groundwater resources, it is suggested that fractures with orientations sub-parallel to the principal compressive stress direction constitute favourable groundwater targets. Knowledge of the regional stress field and fracture set orientations can be used as an effective low cost tool for locating potentially higher yielding boreholes in crystalline rock terrains.     

Probabilistic Analysis of Fracture Reactivation Associated with Deep Underground CO2 Injection

In the context of carbon capture and storage, deep underground injection of CO₂ induces the geomechanical changes within and around the injection zone and their impact on CO₂ storage security should be evaluated. In this study, we conduct coupled multiphase fluid flow and geomechanical modeling to investigate such geomechanical changes, focusing on probabilistic analysis of injection-induced fracture reactivation (such as shear slip) that could lead to enhanced permeability and CO₂ migration across otherwise low-permeability caprock formations. Fracture reactivation in terms of shear slip was analyzed by implicitly considering the fracture orientations generated using the Latin hypercube sampling method, in one case using published fracture statistics from a CO₂ storage site. The analysis was conducted by a coupled multiphase fluid flow and geomechanical simulation to first calculate the three-dimensional stress evolution during a hypothetical CO₂ injection operation and then evaluate the probability of shear slip considering the statistical fracture distribution and a Coulomb failure analysis. We evaluate the probability of shear slip at different points within the injection zone and in the caprock just above the injection zone and relate this to the potential for opening of new flow paths through the caprock. Our analysis showed that a reverse faulting stress field would be most favorable for avoiding fracture shear reactivation, but site-specific analyses will be required because of strong dependency of the local stress field and fracture orientations.     

Compressional- and transpressional-stress pattern for Pliocene and Quaternary brittle deformation in fore arc and intra-arc zones (Andes of Central and Southern Chile)

Kinematic analysis of fault slip data for stress determination was carried out on Late Miocene to Quaternary rocks from the fore arc and intra-arc regions of the Chilean Andes, between 33° and 46° south latitudes. Studies of Neogene and Quaternary infilling (the Central Depression), as well as plutonic rocks of the North Patagonian Batholith along the Liquiñe–Ofqui Fault Zone, have revealed various compressional and/or transpressional states of stress. In the Pliocene, the maximum compressional stress (σ₁) was generally oriented east–west. During the Quaternary, the deformation was partitioned into two coeval distinctive states of stress. In the fore arc zone, the state of stress was compressional, with σ₁ oriented in a N–S to NNE–SSW direction. In the intra-arc zone the state of stress was transpressional with σ₁ striking NE–SW. Along the coast, in one site (37°30′S) the Quaternary strain deformation is extensional, with an E–W direction, which can be explained by a co-seismic crustal bending readjustment.     

Double fracture method ofin situ stress measurement in brittle rock

Summary A new technique, the Double Fracture Method ofin situ stress measurement, is described with regard to its theoretical principles and field instrumentation. The method differs from the hydraulic fracture method in that the loading fluid is prevented from penetrating into the induced fractures. In the new method, the loading pressure is raised to create a set of two mutually perpendicular microfractures on the borehole boundary. The principal stresses and their orientation in the plane normal to the borehole can be deduced from the diametral deformation of the borehole. Laboratory and field measurements validating the method are discussed.     

From geometry to dynamics of microstructure: using boundary lengths to quantify boundary misorientations and anisotropy

The microstructure of a quartzite experimentally deformed and partially recrystallised at 900 °C, 1.2 GPa confining pressure and strain rate 10⁻⁶/s was investigated using orientation contrast and electron backscatter diffraction (EBSD). Boundaries between misoriented domains (grains or subgrains) were determined by image analysis of orientation contrast images. In each domain, EBSD measurements gave the complete quartz lattice orientation and enabled calculation of misorientation angles across every domain boundary. Results are analysed in terms of the boundary density, which for any range of misorientations is the boundary length for that range divided by image area. This allows a more direct comparison of misorientation statistics between different parts of a sample than does a treatment in terms of boundary number.The strain in the quartzite sample is heterogeneous. A 100×150 μm low-strain partially recrystallised subarea C was compared with a high-strain completely recrystallised subarea E. The density of high-angle (>10°) boundaries in E is roughly double that in C, reflecting the greater degree of recrystallisation. Low-angle boundaries in C and E are produced by subgrain rotation. In the low-angle range 0–10° boundary densities in both C and E show an exponential decrease with increasing misorientation. The densities scale with exp(−θ/λ) where λ is approximately 2° in C and 1° in E; in other words, E has a comparative dearth of boundaries in the 8–10° range. We explain this dearth in terms of mobile high-angle boundaries sweeping through and consuming low-angle boundaries as the latter increase misorientation through time. In E, the density of high-angle boundaries is larger than in C, so this sweeping would have been more efficient and could explain the relative paucity of 8–10° boundaries.The boundary density can be generalised to a directional property that gives the degree of anisotropy of the boundary network and its preferred orientation. Despite the imposed strain, the analysed samples show that boundaries are not, on average, strongly aligned. This is a function of the strong sinuosity of high-angle boundaries, caused by grain boundary migration. Low-angle boundaries might be expected, on average, to be aligned in relation to imposed strain but this is not found.Boundary densities and their generalisation in terms of directional properties provide objective measures of microstructure. In this study the patterns they show are interpreted in terms of combined subgrain rotation and migration recrystallisation, but it may be that other microstructural processes give distinctive patterns when analysed in this fashion.