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.     

Spatial correlation of deformation and mineral reaction in experimentally deformed plagioclase–olivine aggregates

Shear deformation of hot pressed plagioclase–olivine aggregates was studied in the presence and absence of mineral reaction. Experiments were performed at 900 °C, 1500 MPa, and a constant shear strain rate of 5×10⁻⁵ s⁻¹ in a solid medium apparatus. Whether the mineral reaction between plagioclase and olivine takes place or not is controlled by choosing the appropriate plagioclase composition; labradorite (An₆₀) does not react, anorthite (An₉₂) does. Labradorite–olivine aggregates deformed without reaction are very strong and show strain hardening throughout the experiment. Syndeformational reaction between olivine and anorthite causes a pronounced strain weakening. The reaction produces fine-grained opx–cpx–spinel aggregates, which accommodate a large fraction of the finite strain. Deformation and reaction are localised within a 0.5-mm-wide sample. Three representative samples were analysed for their fabric anisotropy R* and shape-preferred orientation α* (fabric angle with the shear plane) using the autocorrelation function (ACF). Fabric anisotropy can be calibrated to quantify strain variations across the sheared samples. In the deformed and reacted anorthite–olivine aggregate, there is a strong correlation between reaction progress and strain; regions of large shear strain correspond to regions of maximum reaction progress. Within the sample, the derived strain rate variations range up to almost one order of magnitude.     

The frictional behavior of lizardite and antigorite serpentinites: Experiments,constitutive models,and implications for natural faults

Laboratory studies of the frictional behavior of rocks can provide important information about the strength and sliding stability of natural faults. We have conducted friction experiments on antigorite and lizardite serpentinites, rocks common to both continental and oceanic crustal faults. We conducted both velocity-step tests and timed-hold tests on bare surfaces and gouge layers of serpentinite at room temperature. We find that the coefficient of friction of lizardite serpentinite is quite low (0.15–0.35) and could explain the apparent low stresses observed on crustal transform faults, while that of antigorite serpentinite is comparable to other crustal rocks (0.50–0.85). The frictional behavior of both types of serpentinite is well described by a two-mechanism model combining state-variable-dominated behavior at high slip velocities and flow-dominated behavior at low velocities. The two-mechanism model is supported by data from velocity-step tests and timed-hold tests. The low velocity behavior of serpentinite is strongly rate strengthening and should result in stable fault creep on natural faults containing either antigorite or lizardite serpentinite.     

Compositional changes of minerals associated with dynamic recrystallizatin

The rate of compositional and isotopic exchange between minerals may be enhanced significantly if the rock is deformed simultaneously. The enhanced exchange rate may result from a reduction in grain size (shorter distance for volume diffusion), dissolution and growth of grains by diffusion creep (pressure solution), or the movement of high-angle grain boundaries through strained grains during recrystallization in the dislocation creep regime. The migration of high-angle grain boundaries provides high diffusivity paths for the rapid exchange of components during recrystallization. The operation of the latter process has been demonstrated by deforming aggregates consisting of two plagioclases (An₁ and An₇₉) at 900°C, 1 GPa confining pressure, and a strain rate of ∼2x10^-6s^-1. The polygonal, recrystallized grains were analyzed using an analytical transmission electron microscope and have a variable but often intermediate composition. At the conditions of these experiments, the volume interdiffusion rate of NaSi/CaAl is too slow to produce any observable chemical change, and microstructural-chemical relations indicate that the contribution from diffusion creep was insignificant except for initially fine-grained (2–10 μm) aggregates. These results indicate that strain-induced recrystallization can be an effective mechanism for enhancing the kinetics of metamorphic reactions and for resetting the isotope systematics of minerals such as feldspars, pyroxenes, and amphiboles.     

Reflections on measurement of residual stress in rock

Consideration of the behaviour of elastic bodies shows that it is not possible for a set of overcoring measurements that are made within isolated blocks to show residual strains or stresses that have a non-zero average, unless the size of the equilibrium volume over which the residual stresses balance is both considerably larger than the volume of the overcore and smaller than the size of the isolated block. Since some results have been reported that do not match these constraints, non-elastic behaviour must have occurred during overcoring. A possible explanation, is that oriented microcracks are somehow opened by the overcoring. In some cases, stress fields induced by microcracking near the overcoring cut may explain measured strain changes far from the cut. Consideration of various reported measurements in terms of relative sizes of overcoring and equilibrium volumes shows that care is necessary when interpreting residual stress measurements.     

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.     

An Independent Assessment of the Load/Unload Response Ratio (LURR) Proposed Method of Earthquake Prediction

The Load/Unload Response Ratio (LURR) method is a proposed technique to predict earthquakes that was first put forward by Yin (1987). LURR is based on the idea that when an area enters the damage regime, the rate of seismic activity during loading of the tidal cycle increases relative to the rate of seismic activity during unloading in the months to one year preceding a large earthquake. Since earth tides generally contribute the largest temporal variations in crustal stress, it seems plausible that earth tides would trigger earthquakes in areas that are close to failure (e.g., Vidale et al., 1998). However, the vast majority of studies have shown that earth tides do not trigger earthquakes (e.g., Vidale et al., 1998; Heaton, 1982; Rydelek et al., 1992). In this study, we conduct an independent test of the LURR method, since there would be important scientific and social implications if it were proven to be a robust method of earthquake prediction. Smith and Sammis (2004) undertook a similar study and found no evidence that there was predictive significance to the LURR method. We have repeated calculations of LURR for the Northridge earthquake in California, following both the parameters of X.C. Yin (personal communication) and the somewhat different ones of Smith and Sammis (2004). Though we have followed both sets of parameters closely, we have been unable to reproduce either set of results. Our examinations have shown that the LURR method is very sensitive to certain parameters. Thus it seems likely that the discrepancies between our results and those of previous studies are due to unaccounted for differences in the calculation parameters. A general agreement was made at the 2004 ACES Workshop in China between research groups studying LURR to work cooperatively to resolve the differences in methods and results, and thus permit more definitive conclusions on the potential usefulness of the LURR method in earthquake prediction.     

Hydrogeology,Chemical and Microbial Activity Measurement Through Deep Permafrost

Little is known about hydrogeochemical conditions beneath thick permafrost, particularly in fractured crystalline rock, due to difficulty in accessing this environment. The purpose of this investigation was to develop methods to obtain physical, chemical, and microbial information about the subpermafrost environment from a surface‐drilled borehole. Using a U‐tube, gas and water samples were collected, along with temperature, pressure, and hydraulic conductivity measurements, 420 m below ground surface, within a 535 m long, angled borehole at High Lake, Nunavut, Canada, in an area with 460‐m‐thick permafrost. Piezometric head was well above the base of the permafrost, near land surface. Initial water samples were contaminated with drill fluid, with later samples <40% drill fluid. The salinity of the non‐drill fluid component was <20,000 mg/L, had a Ca/Na ratio above 1, with δ¹⁸O values ～5‰ lower than the local surface water. The fluid isotopic composition was affected by the permafrost‐formation process. Nonbacteriogenic CH₄ was present and the sample location was within methane hydrate stability field. Sampling lines froze before uncontaminated samples from the subpermafrost environment could be obtained, yet the available time to obtain water samples was extended compared to previous studies. Temperature measurements collected from a distributed temperature sensor indicated that this issue can be overcome easily in the future. The lack of methanogenic CH₄ is consistent with the high sulfate concentrations observed in cores. The combined surface‐drilled borehole/U‐tube approach can provide a large amount of physical, chemical, and microbial data from the subpermafrost environment with few, controllable, sources of contamination.