Fracture initiation, growth and effect on stress field: a numerical investigation

A numerical investigation was made of the relationships between fracture initiation, growth, stress field and boundary conditions. Two-dimensional plane strain continuum models were used in which fractures appeared as zones of strain localization developed through application of a strain softening Mohr–Coulomb constitutive model. R and R′ fractures developed first, followed by Y fractures at larger strains. The models showed that equal development of conjugate R and R′ fractures is easily changed to favor one or the other set by minor variations in model initial conditions. Strength loss in fractures caused stress field rotations in regions bounded by fractures, altering the orientation of subsequent fractures. The amount and sense of stress field rotation is dependent on the strength loss during displacement on the fractures, the orientation of fractures, and on the boundary conditions. Y oriented fractures could be explained on the basis of a Mohr–Coulomb failure criterion provided that stress field rotation is accounted for. Monitoring of fracture slip activity showed that, under conditions of constant boundary velocity, slip was discontinuous in time, alternating on fractures throughout the model.     

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.     

Creep of saturated materials as a chemically enhanced rate‐dependent damage process

Material behaviour that exhibits characteristics of creep induced by a spontaneous mineral dissolution enhanced by material damage is studied. It is believed that the characteristic rates of the chemical processes involved determine the time‐rate dependence of the resulting strain. A basic model of a combined chemo‐plastic softening and chemically enhanced deviatoric strain hardening for saturated geomaterials is presented. Chemical softening is postulated to occur as a consequence of the net mass removal resulting from dissolution and precipitation of specific minerals occurring at the damage‐generated inter‐phase interfaces. Closed and open systems are discussed. In the former case, deformation at constant stress results entirely from a local compensation mechanism between the chemical softening and strain hardening. The classical three stages of creep are interpreted in terms of mechanisms of dissolution and precipitation, as well as the variation in the reaction surface areas involved in the mass exchange. In an open system, the above local mechanism is enhanced by the removal of mass via diffusion of species affecting the mass balance. Such a system is addressed via a boundary value problem as shown in an example. Copyright © 2007 John Wiley & Sons, Ltd.     

Evolution of a calcite marble shear zone complex on Thassos Island, Greece: microstructural and textural fabrics and their kinematic significance

The deformation history of a monophase calcite marble shear zone complex on Thassos Island, Northern Greece, is reconstructed by detailed geometric studies of the textural and microstructural patterns relative to a fixed reference system (shear zone boundary, SZB). Strain localization within the massive marble complex is linked to decreasing P–T conditions during the exhumation process of the metamorphic core complex. Solvus thermometry indicates that temperatures of 300–350°C prevailed during part of the shear zone deformation history. The coarse-grained marble protolith outside the shear zone is characterized by symmetrically oriented twin sets due to early coaxial deformation. A component of heterogeneous non-coaxial deformation is first recorded within the adjacent protomylonite. Enhanced strain weakening by dynamic recrystallization promoted strong localization of plastic deformation in the ultramylonite of the calcite shear zone, where high strain was accommodated by non-coaxial flow. This study demonstrates that both a pure shear and a simple shear strain path can result in similar crystallographic preferred orientations (single c-axis maximum perpendicular to the SZB) by different dominant deformation mechanisms. Separated a-axis pole figures (+a- and −a-axis) show different density distributions with orthorhombic texture symmetry in the protolith marble and monoclinic symmetry in the ultramylonite marble consistently with the observed grain fabric symmetry.     

Stresses and velocities in orogenic wedges with power-law rheology and linearly varying longitudinal strain rate

A two-dimensional analytical solution for stress, strain rate, and velocity is obtained for parallel-sided and wedge-shaped blocks with generalized viscous rheology (linearly viscous and power-law) deforming in plane strain. The main assumptions used in the derivation of the solution are that the material is incompressible, the longitudinal gradient in shear stress is much less than the vertical gradient of vertical normal stress, and the longitudinal strain rate varies linearly in the horizontal direction. Velocity boundary conditions are specified at the top of the block, and shear stress boundary conditions at the base of the block. In the one-dimensional case (where stress and strain rate do not vary in the longitudinal direction), the solution reduces to a well-known solution for the deformation of parallel-sided ice sheets [Nye, J. F. (1957) The distribution of stress and velocity in glaciers and ice sheets. Proceedings of the Royal Society of London A-239, 113–133]. The stress equilibrium for tapered wedges [Platt, J. P. (1986) Dynamics of orogenic wedges and the uplift of high-pressure metamorphic rocks. Geological Society of America Bulletin 97, 1037–1053] is a special case of the present stress solution. Implementation of the solution requires the subdivision of the wedge into vertical segments, and yields the tectonic normal and shear stresses that must be applied to the rear of a block with specified rheology in order to maintain a given longitudinal strain rate. The solution makes it possible to model deformation patterns analytically with longitudinally varying strain rate (including coeval compression and extension) and with vertical components of velocity reflecting the effects of underplating.     

The stream function and Gauss' principle of least constraint: Two useful concepts for structural geology

To the extent that rock deformation can be approximated by a two-dimensional Newtonian model, a powerful stream-function simulation method is applicable. The significance of stream functions is that velocity, strain, stress and energy derived from the same stream function satisfy automatically three basic conditions of dynamics:

1. (1) the condition of continuity.

2. (2) the Navier-Stokes equations.

3. (3) conservation of energy.

Hence we state with Jaeger: “If a stream function can be found which satisfies the boundary conditions of a dynamic model the complete solution follows.” All pertinent bits of dynamic information are implied in the stream function from which they can be directly derived, guaranteed—so to speak—not to violate the basic conditions of dynamics. Stream functions useful in structural geology are solutions of: A double-polynomial solution of max. degree 14 is developed, in which the coefficients are related controlled by the ⁴ψ = 0 constraint, and their absolute values are determined by the boundary conditions of specific models and by the condition of maximum rate of energy dissipation or maximum rate of decline of potential energy. The polynomial stream function is applied to a collapsing viscous “nappe” consisting of a thin basal layer with low viscosity on which a thicker layer with high viscosity slides due to gravitational spreading. The velocity of forward movement depends upon absolute and relative values of the following parameters: viscosity, thickness, the aspect ratio and density. The velocity of a variety of nappes with different thicknesses, aspect ratios, viscosities and densities is determined.     

Rheological and kinematical responses to flow of two-phase rocks

Numerical simulations have been performed to investigate the strain-dependent behaviour of rheological and kinematical responses to flow of two-phase rocks using the commercial finite-difference program FLAC2D. It was assumed that the two phases have Maxwell rheology. Plane strain and velocity boundary condition, which produces a simple shear deformation, were also assumed. Two types of geometries were considered: strong phase supported (SPS) and weak phase supported (WPS). We calculated strain-dependent variations of effective viscosity and partitioning of strain rate, vorticity and kinematic vorticity number during deformation in both SPS and WPS structure models.The results show that the strain-dependent behaviour is largely influenced by the geometry of the composite. SPS models show both strain hardening and strain softening during the simulations, with strain hardening preceding strain softening. A critical shear strain is necessary to begin the strain softening behaviour. Strain hardening and strain softening are accompanied by a reduction and an increase of the partition of strain rate into the weak phase, respectively. On the other hand, WPS models show only weak strain hardening and strain softening, being the strain-dependent behaviour close to a steady state flow. In addition, the following results are obtained on vorticity and kinematic vorticity number; (1) in both SPS and WPS models the partition of vorticity into weak phase increases with progressive shear strain, i.e. the strong phase becomes less rotational, (2) in SPS models weak inclusions changes from sub-simple shear to super-simple shear with progressive strain, whereas the strong matrix changes from super-simple shear to sub-simple shear, (3) in WPS models the strong inclusions with high viscosity contrasts are less rotational but can be in super-simple shear condition to high strains.The observed strain-dependent behaviours have been compared with previous proposed analytical models. The degree of agreement is variable. Balshin and Ryshkewitch–Duckworth models are only applicable to SPS models. Ji-generalized mixture rule model is applicable to both models.The results suggest that polyphase rocks with SPS structure during ductile shear deformation respond as strain softening materials, after an initial strain hardening stage that may drive to the strain localization into the material.     

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.     

Stress‐dependent hardening and failure surfaces of dry sand

During several triaxial compression experiments on plastic hardening, softening, and failure properties of dense sand specimens, it was found on various stress paths that the size of the failure surface was not constant. Instead, it changed depending on the current state of hydrostatic pressure. This finding is in contrast to the standard opinion consisting of the fact that the failure surface remains constant, once it has been reached during an experiment or in situ. In general, the behaviour of cohesionless granular‐material‐like sand is somehow characterised in between fluid and solid, where the solid behaviour results from the angle of internal friction and the confining pressure. Although the friction angle is an intrinsic material property, the confining pressure varies with the boundary conditions, thus defining different solid properties like plastic hardening, softening, and also failure. Based on our findings, it was the goal of the present contribution to introduce an improved setting for the plastic strain hardening and softening behaviour including the newly found yield properties at the limit state. For the identification of the material parameters, a complete triaxial experimental analysis of the tested sand is given. The overall elasto‐plasticity concept is validated by numerical computations of several laboratory foundation‐ and slope‐failure experiments. The performance of the proposed approach is compared with the standard concept of a constant failure surface, where the corresponding yield surfaces are understood as contours of equivalent plastic work or plastic strain. Copyright © 2012 John Wiley & Sons, Ltd.     

Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model

A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd.     

Senescence- and death-related alteration of chlorophylls and carotenoids in marine phytoplankton

This report extends previous work ([Louda et al., 1998a] and [Louda et al., 1998b]. Chlorophyll degradation during senescence and death. Organic Geochemistry 29, 1233–1251.) in which we detailed type-I (alteration) and -II (destruction) degradation of chlorophyll with representative fresh water phytoplankton. The present study covers similar experiments with marine phytoplankton, namely, a cyanobacterium (“ANA” Anacystis sp), a coccolithophore (“COC” Coccolithophora sp.), a dinoflagellate (“GYM” Gymnodinium sp.) and two diatoms (“CYC” Cyclotella meneghiniana and “THAL” Thalassiosira sp.). Mg loss (‘pheophytinization') was rapid and continuous in all species under room-oxic conditions and slow or sporadic under anoxic conditions. The proportion of dephytylated pigments (pheophorbides-a, chlorophyllides-a), relative to the phytylated forms (chlorophyll-a, pheophytins-a), increased over the first year under room-oxic conditions and in room-anoxic conditions only in “CYC”. Pheophorbide-a was converted to pyropheophorbide-a within 15 months only in “THAL” and “ANA”, and slightly in “COC”. After 9–15 months of oxic incubation, “COC” was found to contain traces of purpurin-18 phytyl ester. Consideration of carotenoid pigments is also included herein. All fucoxanthin containing species, except “THAL”, exhibited conversion of fucoxanthin to fucoxanthinol in room-oxic conditions. Diadinoxanthin was rapidly de-epoxidized to give diatoxanthin within the first 2–4 weeks. Diatoxanthin then disappeared from all species by 15 months with a concurrent increase in a pigment which we tentatively identify as a cis-zeaxanthin. Incubations of pure cultures are found to be an effective way by which to model the early type-I reactions for both chlorophylls and carotenoids. The influence of oxygen during senescence-death and the onset of early diagenesis is of paramount importance. The absence of oxygen and, by inference, aerobic microbiota, retards the breakdown of these pigments dramatically.     

Dynamic recrystallization of wet synthetic polycrystalline halite: dependence of grain size distribution on flow stress, temperature and strain

It is often observed that dynamic recrystallization results in a recrystallized grain size distribution with a mean grain size that is inversely related to the flow stress. However, it is still open to discussion if theoretical models that underpin recrystallized grain size–stress relations offer a satisfactorily microphysical basis. The temperature dependence of recrystallized grain size, predicted by most of these models, is rarely observed, possibly because it is usually not systematically investigated. In this study, samples of wet halite containing >10 ppm water (by weight) were deformed in axial compression at 50 MPa confining pressure. The evolution of the recrystallized grain size distribution with strain was investigated using experiments achieving natural strains of 0.07, 0.12 and 0.25 at a strain rate of 5×10⁻⁷ s⁻¹ and a temperature of 125 °C. The stress and temperature dependence of recrystallized grain size was systematically investigated using experiments achieving fixed strains of 0.29–0.46 (and one to a strain of 0.68) at constant strain rates of 5×10⁻⁷–1×10⁻⁴ s⁻¹ and temperatures of 75–240 °C, yielding stresses of 7–22 MPa. The microstructures and full grain size distributions of all samples were analyzed. The results showed that deformation occurred by a combination of dislocation creep and solution-precipitation creep. Dynamic recrystallization occurred in all samples and was dominated by fluid assisted grain boundary migration. During deformation, grain boundary migration results in a competition between grain growth due to the removal of grains with high internal strain energy and grain size reduction due to grain dissection (i.e. moving boundaries that crosscut or consume parts of neighbouring grains). At steady state, grain growth and grain size reduction processes balance, yielding constant flow stress and recrystallized grain size that is inversely related to stress and temperature. Evaluation of the recrystallized grain size data against the different models for the development of mean steady state recrystallized grain size revealed that the data are best described by a model based on the hypothesis that recrystallized grain size organizes itself in the boundary between the (grain size sensitive) solution-precipitation and (grain size insensitive) dislocation creep fields. Application of a piezometer, calibrated using the recrystallized grain size data, to natural halite rock revealed that paleostresses can vary significantly with temperature (up to a factor of 2.5 for T=50–200 °C) and that the existing temperature independent recrystallized grain size–stress piezometer may significantly underestimate flow stresses in natural halite rock.     

The optical features and hydrocarbon-generating model of “barkinite” from Late Permian coals in South China

Leping coal is known for its high content of “barkinite”, which is a unique liptinite maceral apparently found only in the Late Permian coals of South China. “Barkinite” has previously identified as suberinite, but on the basis of further investigations, most coal petrologists conclude that “barkinite” is not suberinite, but a distinct maceral. The term “barkinite” was introduced by (State Bureau of Technical Supervision of the People's Republic of China, 1991, GB 12937-91 (in Chinese)), but it has not been recognized by ICCP and has not been accepted internationally.In this paper, elemental analyses (EA), pyrolysis-gas chromatography, Rock-Eval pyrolysis and optical techniques were used to study the optical features and the hydrocarbon-generating model of “barkinite”. The results show that “barkinite” with imbricate structure usually occurs in single or multiple layers or in a circular form, and no definite border exists between the cell walls and fillings, but there exist clear aperture among the cells.“Barkinite” is characterized by fluorescing in relatively high rank coals. At low maturity of 0.60–0.80%R_o, “barkinite” shows strong bright orange–yellow fluorescence, and the fluorescent colors of different cells are inhomogeneous in one sample. As vitrinite reflectance increases up to 0.90%R_o, “barkinite” also displays strong yellow or yellow–brown fluorescence; and most of “barkinite” lose fluorescence at the maturity of 1.20–1.30%R_o. However, most of suberinite types lose fluorescence at a vitrinite reflectance of 0.50% Ro, or at the stage of high volatile C bituminous coal. In particular, the cell walls of “barkinite” usually show red color, whereas the cell fillings show yellow color under transmitted light. This character is contrary to suberinite.“Barkinite” is also characterized by late generation of large amounts of liquid oil, which is different from the early generation of large amounts of liquid hydrocarbon. In addition, “barkinite” with high hydrocarbon generation potential, high elemental hydrogen, and low carbon content. The pyrolysis products of “barkinite” are dominated by aliphatic compounds, followed by low molecular-weight aromatic compounds (benzene, toluene, xylene and naphthalene), and a few isoprenoids. The pyrolysis hydrocarbons of “barkinite” are mostly composed of light oil (C₆–C₁₄) and wet gas (C₂–C₅), and that heavy oil (C₁₅⁺) and methane (C₁) are the minor hydrocarbon.In addition, suberinite is defined only as suberinized cell walls—it does not include the cell fillings, and the cell lumens were empty or filled by corpocollinites, which do not show any fluorescence. Whereas, “barkinite” not only includes the cell walls, but also includes the cell fillings, and the cell fillings show bright yellow fluorescence.Since the optical features and the hydrocarbon-generating model of “barkinite” are quite different from suberinite. We suggest that “barkinite” is a new type of maceral.     

Tectonomagmatism in continental arcs: evidence from the Sark arc complex

The island of Sark (Channel Islands, UK) exposes syntectonic plutons and country rock gneisses within a Precambrian (Cadomian) continental arc. This Sark arc complex records sequential pulses of magmatism over a period of 7 Ma (ca. 616–609 Ma). The earliest intrusion (ca. 616 Ma) was a composite sill that shows an ultramafic base overlain by a magma-mingled net vein complex subsequently deformed at near-solidus temperatures into the amphibolitic and tonalitic Tintageu banded gneisses. The deformation was synchronous with D₂ deformation of the paragneissic envelope, with both intrusion and country rock showing flat, top-to-the-south LS fabrics. Later plutonism injected three homogeneous quartz diorite–granodiorite sheets: the Creux–Moulin pluton (150–250 m; ca. 614 Ma), the Little Sark pluton (>700 m; 611 Ma), and the Northern pluton (>500 m; 609 Ma). Similar but thinner sheets in the south (Derrible–Hogsback–Dixcart) and west (Port es Saies–Brecqhou) are interpreted as offshoots from the Creux–Moulin pluton and Little Sark pluton, respectively. All these plutons show the same LS fabric seen in the older gneisses, with rare magmatic fabrics and common solid state fabrics recording syntectonic crystallisation and cooling. The cooling rate increased rapidly with decreasing crystallisation age: >9 Ma for the oldest intrusion to cool to lower amphibolite conditions, 7–8 Ma for the Creux Moulin pluton, 5–6 Ma for the Little Sark pluton, and <3 Ma for the Northern pluton. This cooling pattern is interpreted as recording extensional exhumation during D₂. The initiation of the D₂ event is suggested to have been a response to the intrusion of the Tintageu magma which promoted a rapid increase in strain rate (>10⁻¹⁴ s⁻¹) that focussed extensional deformation into the Sark area. The increased rates of extension allowed ingress of the subsequent quartz diorite–granodiorite sheets, although strain rate slowly declined as the whole complex cooled during exhumation. The regional architecture of syntectonic Cadomian arc complexes includes flat-lying “Sark-type” and steep “Guernsey-type” domains produced synchronously in shear zone networks induced by oblique subduction: a pattern seen in other continental arcs such as that running from Alaska to California.     

Large-strain deformation and strain partitioning in polyphase rocks: Dislocation creep of olivine–magnesiowüstite aggregates

Aggregates composed of olivine and magnesiowüstite have been deformed to large strains at high pressure and temperature to investigate stress and strain partitioning, phase segregation and possible localization of deformation in a polyphase material. Samples with 20 vol.% of natural olivine and 80 vol.% of (Mg_0.7Fe_0.3)O were synthesized and deformed in a gas-medium torsion apparatus at temperatures of 1127 °C and 1250 °C, a confining pressure of 300 MPa and constant angular displacement rates equivalent to constant shear strain rates of 1–3.3 × 10^− 4 s^− 1. The samples deformed homogeneously to total shear strains of up to γ  15. During constant strain rate measurements the flow stress remained approximately stable at 1250 °C while it progressively decreased after the initial yield stress at the lower temperature. Mechanical data, microstructures and textures indicate that both phases were deforming in the dislocation creep regime. The weaker component, magnesiowüstite, controlled the rheological behavior of the bulk material and accommodated most of the strain. Deformation and dynamic recrystallization lead to grain refinement and to textures that were not previously observed in pure magnesiowüstite and may have developed due to the presence of the second phase. At 1127 °C, olivine grains behaved as semi-rigid inclusions rotating in a viscous matrix. At 1250 °C, some olivine grains remained largely undeformed while deformation and recrystallization of other grains oriented for a-slip on (010) resulted in a weak foliation and a texture typical for pure dry olivine aggregates. Both a-slip and c-slip on (010) were activated in olivine even though the nominal stresses were up to 2 orders of magnitude lower than those needed to activate these slip systems in pure olivine at the same conditions.     

A comparison of intracrystalline deformation in naturally and experimentally deformed quartzites

A relatively undeformed quartzite sample from the Weverton formation was experimentally deformed in plane strain at a temperature of 700° C, a confining pressure of 15 kb and a constant strain rate of 10⁻⁶/sec, in a modified Griggs apparatus. A comparison of the known experimental strain for the sample with that measured from deformed rutile needles within the quartz grains shows fairly close agreement between the two values. This confirms the validity of using the needles as intracrystalline strain markers. A comparison has been made of the microstructures and preferred orientations in the experimentally deformed sample and a naturally deformed sample of the same quartzite which has undergone the same strain. The experimentally deformed sample exhibits more inhomogeneous intragranular deformation and a “double funnel” pattern of c axes, while the naturally deformed sample exhibits more homogeneous intragranular deformation and a broad great circle girdle of c axes normal to the foliation and lineation.     

Deformation mechanism maps for feldspar rocks

Deformation mechanism maps for feldspar rocks were constructed based on recently published constitutive laws for dislocation and grain boundary diffusion creep of wet and dry plagioclase aggregates. The maps display constant temperature contours in stress-grain size space for strain rates ranging from 10⁻¹⁶ to 10⁻¹² s⁻¹.Two fields of dominance of grain boundary diffusion-controlled creep and dislocation creep are separated by a strongly grain size-sensitive transition zone. For wet rocks, diffusion-controlled creep dominates below a grain size of about 0.1–1 mm, depending on temperature, stress, strain rate and feldspar composition. Plagioclase aggregates containing up to 0.3 wt.% water as often found in natural feldspars are more than 2 orders of magnitude weaker than dry rocks. The strength of water-bearing feldspar rocks is moderately dependent on composition and water fugacity.For a grain size range of about 10–50 μm commonly observed in natural ultramylonites, the deformation maps predict that diffusion-controlled creep is dominant at greenschist to granulite facies conditions. Low viscosity estimates of 10¹⁸–10¹⁹ Pa·s from modeling postseismic stress relaxation and channel flow of the continental lower crust can only be reconciled with laboratory experiments assuming dislocation creep at high temperatures >900 °C or, at lower temperatures, diffusion creep of fine-grained rocks possibly localized in abundant high strain shear zones. For similar thermodynamic conditions and grain size, lower crustal rocks are predicted to be less than order of magnitude weaker than upper mantle rocks.     

New observational information on the precursory accelerating and decelerating strain energy release

Recent reliable data are used to study the behavior of seismic activity before 46 strong shallow earthquakes (M ≥ 6.0), which correspond to five complete samples of mainshocks. These samples include 6 mainshocks (M = 6.0–7.1) that occurred in western Mediterranean since 1980, 17 mainshocks (M = 6.0–7.2) which occurred in the Aegean (Greece and surrounding area) since 1980, 5 mainshocks (M = 6.4–7.5) that occurred in Anatolia since 1980, 12 mainshocks (M = 6.0–7.3) that occurred in California since 1980 and 6 mainshocks (M = 7.0–8.3) that occurred in Japan since 1990. In all 46 cases, a similar precursory seismicity pattern is observed. Specifically, it is observed that accelerating Benioff strain (square root of seismic energy) release caused by preshocks occurs in a broad circular region (critical region), with a radius about eight times larger than the fault length of the mainshock, in agreement with results obtained by various research groups during the last two decades. However, in a much smaller circular region (seismogenic region), with a radius about four times the fault length, the corresponding preshock strain decelerates with the time to the mainshock. The time variation of the strain follows in both cases a power law but the exponent power is smaller than unit (m ¯ = 0.3) in the case of the accelerating preshock strain and larger than unit (m ¯ = 3.0) in the case of the decelerating preshock strain. Predictive properties of this “Decelerating In–Accelerating Out Strain” model are expressed by empirical relations. The possibility of using this model for intermediate-term earthquake prediction is discussed and the relative model uncertainties are estimated.     

Microstructure mechanism map of dynamically recrystallized marble

The relative nucleus density (RND) model of dynamically recrystallized grain size [Sakai, T., Jonas, J.J. 1984. Dynamic recrystallization: mechanical and microctructutal consideration. Acta metallurgica, 32, 198–209] was applied to experimentally and to naturally deformed marbles that have undergone dynamic recrystallization. The model shows that a relationship between initial grain size (D₀) and stable dynamically recrystallized grain size (D_S) for a given value of temperature-corrected strain-rate (Z) controls grain size evolution during dynamic recrystallization. New microstructural mechanism maps (MM-maps) for experimentally and naturally deformed marbles (based on previously published data) were defined in log grain size–log Z space and show two distinct regions of grain reduction and grain coarsening. The boundary between these two regions corresponds to an equation relating dynamically recrystallized grain size and temperature corrected strain rate, as proposed in this work. The new MM-map was used to trace semi-quantitatively microstructural and grain size evolution in naturally deformed marbles that underwent dynamic recrystallization at different thermal conditions. The boundary between grain coarsening and grain reduction does not necessarily coincide with the boundary between rotation and migration recrystallization mechanisms. Assessment of available natural data shows that the boundary condition D₀ = 2D_S between grain-coarsening and grain-reduction introduced by Sakai and Jonas [Sakai, T., Jonas, J.J. 1984. Dynamic recrystallization: mechanical and microctructutal consideration. Acta metallurgica, 32, 198–209] is not required for naturally deformed marble.     

Experimental study of strain localization in sensitive clays

For evaluation of slope stability in materials displaying strain-softening behavior, knowledge concerning the failed state material response is of importance. Here, soft sensitive clay is studied. Such clays behave contractant at failure, which for undrained conditions yields a strain-softening behavior governed by the generation of excess pore water pressure. Strain softening is further linked with material instability and the phenomenon of strain localization. In the case of shear band formation, internal pore pressure gradients are then expected to be present for globally undrained conditions in the sensitive clay due to its low permeability. In the present study, this hypothesis and its implications on the global response and shear band properties are investigated. Utilizing an experimental setup with a modified triaxial cell allowing for shear band formation, the effect of varying the displacement rate is studied. Onset of strain localization is interpreted to occur just before or at the peak shear strength. A strong rate dependency of the softening response is observed. Increasing displacement rates give raised brittleness in terms of the slope of the global softening curve due to accumulating pore pressure. Also, reduced shear band thickness and a shear band inclination approaching 45° are obtained for increasing rates. In the context of slope failure in such materials, the rate dependency in the post-peak state opens up for a large variation in behavior, all depending on time as an important factor.