Rheology and deformation mechanisms

Continuous deformation or creep

Permanent (non-recoverable) deformation that is not associated with a loss of cohesion is variously termed continuous deformation, flow, ductile flow or creep.

In creep, it is the strain rate, not the finite strain, that is a function of the applied stress.

The exact form of the relationship between strain rate and stress is called the constitutive law for the material.

Ideal Viscous (Newtonian) behaviour

Liquids flow as soon as a shear stress is applied, and the rate of shear strain is directly proportional to the applied shear stress. The constant of proportionality is called the viscosity. Some solid rocks behave in this way at high temperatures, close to their melting points, though the actual viscosities are very high.

The equation for viscous flow is

strain rate is proportional to differential stress

or Viscous

where η is the viscosity

Visco-elastic material

Visco-elastic are essentially elastic, but show a delayed response to stress as the material flows into its elastically deformed state.

Elastico-viscous material

Elastico-viscous material deforms viscously at low strain rates and low stress, flowing like a liquid, but when the strain rate increases, elastic (and eventually brittle) behaviour takes over. Some rocks at high temperatures, and the silly putty used in the lecture demonstration, exhibit behaviour that it approximately elastico-viscous.

Ideal plastic behaviour

We may imagine a material which is rigid or elastic at low stress but which starts to flow easily when a fixed yield stress is reached.

differential stress is a constant (yield stress)

This is described as ideal plastic or Von Mises behaviour. The stress within the material can never exceed the yield stress because the material flows at whatever speed is necessary to get out of the way of the applied stress. Obviously this is an idealized model for flow, but it approximates some types of materials, especially if the strain rate is kept low.

Some geologic materials undergo a transition from brittle to plastic behaviour as depth increases. The failure envelope for such a material appears as in VDPM 06.21a

Power law creep

Many minerals do not show Newtonian behaviour. Instead, the apparent viscosity changes with strain rate. Some materials follow what is called a power law relationship.

or plc

where A is a constant, E* is an activation energy, T is the absolute temperature, R is the gas constant

Exponential creep

Under other conditions, it is found that the strain rate is proportional to an exponential function of the stress.

so it's better to write the exponential function thus

strain rate is proportional to exponential function of differential stress

or with constants added expcrp

where A is a constant, E* is an activation energy, T is the absolute temperature, R is the gas constant

Crystal defects

It's possible to calculate bond strengths in crystals from chemical data.  In most cases, the strengths indicate that it's impossible, with geologically reasonable stresses, to simultaneously break all the bonds in a crystal, so as to (for example) shear the crystal.

Nevertheless, we see clear evidence for shearing of crystals in deformed rocks.

It turns out that crystal defects - departures from the ideal lattice structure - are critically important in allowing crystals to deform.

Many of the illustrations in this section are from an excellent series of lecture notes posted by M Jessell, P Bons & P Rey 2002 at Microstructures Online http://www.earth.monash.edu.au/Teaching/mscourse/

Point defects

Point defects include:

Diagram showing types of point defects

http://www.virtualexplorer.com.au/special/meansvolume/contribs/jessell/

Point defects can move through crystals by diffusion (random motion of atoms). Vacancies are particularly important in allowing diffusion to take place

Animation of vacancy diffusion.

Line defects or dislocations

Edge dislocations

Line defects where there is a non-standard linking of atoms are also known as dislocations.

Diagram of edge dislocation (www.crystalmaker.com) [Equivalent diagram in VDPM 9.08a ]

Large numbers of dislocations produce a curved crystal lattice VDPM 09.24

 

This shows up in thin section as undulose extinction [with accessory plate]

Dislocations in a crystal under TEM

The nature of the defect is defined by theBurger's vector b.

To define the Burger's vector, describe a closed loop in an undeformed section of crystal.  Then follow a similar loop around the dislocation; the gap between the start and end point is the Burger's vector b.

An edge dislocation has b perpendicular to the dislocation line.  It marks the edge of an extra half-plane of atoms.

Screw dislocation

A screw dislocation has b parallel to the dislocation line.  Atoms are connected in a helix ("spiral staircase") around the screw dislocation.
In practice, a given dislocation can change orientation from parallel to perpendicular to b, and all orientations in between. VDPM 09.06

Planar defects

Twin planes

Typically, across a twin plane, the structure of a crystal is reflected. (Rotations are also possible).  Calcite and plagioclase feldspar are two minerals that typically show abundant twin planes.

Subgrain boundaries

Subgrain boundaries are surfaces across which the crystal lattice undergoes a subtle change in orientation.

A subgrain boundary is equivalent to a dislocation wall.

They are often visible under crossed polars in grains that show undulose extinction.

Microcracks

Microcracks are important in the initiation of brittle deformation, because their ends are regions of stress concentration that can nucleate faults.

Stacking faults

Stacking faults are planar sufaces across which unit cells are built ('stacked') in an order that does not match that of the perfect crystal.

Creep processes

Herring Nabarro Creep

Also known as volume diffusion - crystal changes shape by migration of point defects - particularly vacancies - through the crystal.

Typically Herring Nabarro creep follows a Newtonian or viscous constitutive law, but is also dependent on grain size, so that the viscosity is proportional to the cube of the grain size.

Coble Creep

Also known as grain boundary diffusion

Crystal changes shape by migration of atoms along grain boundaries. This works better at small grain sizes because the area of grain boundary is large.

The constitutive law for Coble creep is also viscous or Newtonian, but in this case the viscosity is proportional to the square of the grain size.

Solution Creep

Also known as pressure solution

Rate depends on amount, pressure, and chemistry of fluid present. Leads to development of stylolites and pressure solution cleavage

Strain rates for pressure solution are typically highly dependent on chemical parameters and cannot easily be fitted into one of the constitutive laws for creep.

Superplasticity

In some highly deformed materials, grains can move past each other as the material deforms, leaving behind a fine-grained aggregate that does not show strong CPO. How does this happen? It seems that corners and other irregularities that would normally prevent grain-boundary-sliding can be diffused out of the way- the crystals remould themselves to allow grains to slide past each other. This phenomenon is called superplasticity. It's occurrence is a little controversial but it has been cited in some examples of fine grained rocks that appear to have sheared but which do not show the characteristic CPO of mylonites.

VDPM 09.32

The flow law for superplasticity is viscous or Newtonian.

Dislocation-based mechanisms

Many of the illustrations in this section are from an excellent series of lecture notes posted by M Jessell, P Bons & P Rey 2002 at Microstructures Online http://www.earth.monash.edu.au/Teaching/mscourse/

Dislocation glide

Dislocation glide is the fundamental process of deformation at high strain rates, by which mylonites flow

The constitutive law for dislocation glide is typically exponential creep.

Dislocation glide leads to huge densities of dislocations within crystals. Eventually leads to strain hardening because the dislocations get tangled - they need good, intact crystal to migrate, and so are held up by other dislocations.

Dislocation creep

However, at slower strain rates and/or higher temperatures, there are other processes involving diffusion that can assist with dislocation motion.

Deformation involving a combination of diffusion with dislocation glide are known as dislocation creep. In metallurgy, diffusion-assisted reorganization of dislocations is known as recovery.

Several processes are involved in dislocation creep and recovery:

Dislocation creep typically follows a power law relationship in which strain is proportional to stress raised to some power between 1 and 5, typically around 3

Dynamic recrystallization

Notice that the crystal lattice is distorted in the vicinity of a dislocation wherever one occurs. This results in availability of energy to break bonds, and thus there is a tendency for a crystal containing a lot of dislocations to rebuild itself into numerous small crystals. This is the process of dynamic recrystallization that produces mylonite.

Dynamic recrystallization produces a strong crystal preferred orientation (CPO) or lattice preferred orientation (LPO).

Thin section of mylonite

Thin section of mylonite xpl

Stereographic projection of quartz C-axes from mylonite.

Mechanisms based on plane defects

Planar defects are also sites of high energy, and can play a role in deformation

Twin gliding

Strain twinning (twin gliding) is a major mechanism of deformation in carbonates

Grain boundary migration

Grain boundaries can be thought of as major planar defects, and are associated with extra energy (surface energy) that can be released by the transfer of atoms across grain boundaries

Deformation mechanism maps <this section under construction>

Behaviour of common minerals

For most of these types of deformation it's now possible to determine a theoretical relationship between stress and strain, based on knowledge of the kinetic behaviour of atoms at various temperatures, and on the bond strengths for different types of bonds (from mineral chemistry).

Hence we can determine which mechanism gives the fastest strain rate at a given set of conditions (temperature and differential stress).

This leads to the idea of a deformation map - typically a plot with temperature and differential stress on the axes, divided into areas according to which mechanism is fastest. Here are some deformation maps for common minerals.

These maps show the fastest deformation mechanism under given conditions of T and differential stress. Contours are drawn within the diagrams showing the strain rate predicted under those given conditions. (Note: the lines are numbered with the log of the strain rate: where it says '-9' this means 10-9).

In general, only the highest strain rates can be experimentally verified: the lower strain-rate parts of the diagram are predicted based on chemically determined diffusion constants and bond strengths.

Strength of the lithosphere

We can use these diagrams to make predictions about the strength of the lithosphere at various depths, by assuming a geothermal gradient (rate of temperature increase with depth).

In order to do this we have to bring brittle deformation into the picture. Brittle deformation is strongly dependent on pressure. The Mohr-Coulomb failure relationship or 'Byerlee's Law' allows us to predict the increasing differential stress required for fracture with increasing pressure. Hence the strength depth relationship shows an initial more or less linear increase of strength with depth.

(The relationship is a little different depending on whether the lithosphere is being extended or shortened. In general if the maximum compressive stress is vertical, a rather lower differential stress is required than if the maximum compressive stress is horizontal. This follows from the Mohr construction for brittle failure.)

Diagram

At some depth, the brittle deformation curve intersects one the strength-depth curve for one of the ductile mechanisms. This point marks the brittle-ductile transistion. Below this point, creep requires less differential stress than fracture, so rocks flow.

If we are going to make deductions about the strength of the lithosphere we need to assume some compositions: in the following diagram we assume that quartz is the strength controlling-mineral in the crust, and that olivine controls strength in the mantle. For an average continental geothermal gradient this predicts a strong mid-upper crust, a weak lower crust, and a strong upper mantle.

VDPM14.08c

We can use somewhat more sophisticated models. The diagram below shows strength profiles for a crust that has a wet-quartz-dominated 'granitic' upper layer, a dry-quartz-dominated mid crust, a feldspar-dominated 'granulite' lower crustal layer, and an olivine-dominated peridotite lithospheric mantle.

Diagram

In the lower part of the diagram is a drawing based on a seismic cross-section of part of a rift zone on the Atlantic margin in NW Scotland. It shows the variation in structural style based on the inferred strength distribution. In the west (left) heat flow is low, so geothermal gradients are low. Fault-like behaviour is seen to great depth (note that the deformation map tells us that below a few km the 'faults' must be shear zones). In the east, heat-flow is higher, the crust is softer, and 'faults' extend to shallow depth only.