地球与空间科学 · 第 5 册 · 本科三年级

University Earth and Space Sciences — Year 3

书籍暂未推出简体中文版——当前显示英文版。

University Earth and Space Sciences — Year 3 · 本科三年级

7Deformation Mechanisms and Shear Zones

At the surface, the granite of a mountain shatters along faults. Twenty kilometres down, the same granite flows like toffee in a shear zone a few hundred metres wide, its feldspars rolled into eyes and its quartz drawn out into ribbons. A crystal can flow without breaking because it is never perfect: it moves by its defects, one row of atoms at a time, or by atoms diffusing from one face of a grain to another. This chapter follows those mechanisms from the crystal lattice to the laws that geodynamics uses, to the maps that say which mechanism wins, and to the shear zones in which the rock records how it moved.

You already know

Point defects, the Arrhenius law and the activation energy are in the chemistry books’ university volumes; the diffusion coefficient in the physics books’. Stress and strain, the strain ellipse, simple shear and its ellipse (S1,2=1+γ2/4±γ/2S_{1,2} = \sqrt{1 + \gamma^2/4} \pm \gamma/2, long axis at θ′\theta' from the shear direction with tan⁡2θ′=2/γ\tan2\theta' = 2/\gamma), power-law creep with its constants for wet quartzite and dry olivine, the strength envelope and the brittle–ductile transition are in the Year 1 volume, in its chapter on stress, strain and rheology; foliation, lineation and the descriptive mylonite in its chapter on faults, folds and fabrics. The effective viscosity of the mantle is in Chapter 6.

7.1 Crystal defects

To shear a perfect crystal, every bond across a lattice plane would have to break at once, and that needs a stress of the order of a tenth of the shear modulus, which in the mantle is about 8×1010 Pa8 \times 10^{10}\,\mathrm{Pa} (Chapter 6): thousands of megapascals. Real crystals yield at far lower stresses, because they slip one line of atoms at a time.

Definition 7.1 (Dislocation)

A dislocation is a line defect in a crystal: the boundary, inside the crystal, between a part of a lattice plane that has slipped and a part that has not. Its Burgers vector b\mathbf b is the slip that it carries, one lattice spacing in a lattice direction. In an edge dislocation b\mathbf b is perpendicular to the line, which is the edge of an extra half-plane of atoms; in a screw dislocation it is parallel to the line. A dislocation moves by glide on its slip plane, and an edge dislocation leaves that plane by climb, as atoms or vacancies diffuse to or from the edge of its half-plane.

(a) An edge dislocation: an extra half-plane of atoms (red) ends on the slip plane, distorting the lattice around its edge (the symbol ). Under the shear stress  it glides along the slip plane by breaking and re-forming one row of bonds at a time. (b) When it has crossed the crystal, the upper half has slipped by one Burgers vector b, and the lattice is perfect again.
(a) An edge dislocation: an extra half-plane of atoms (red) ends on the slip plane, distorting the lattice around its edge (the symbol ⊥\perp). Under the shear stress τ\tau it glides along the slip plane by breaking and re-forming one row of bonds at a time. (b) When it has crossed the crystal, the upper half has slipped by one Burgers vector b\mathbf b, and the lattice is perfect again.

The lattice resists the glide of a dislocation with a small stress, the Peierls stress, which depends on the strength and the direction of the bonds: that is why some slip systems of a mineral are easy and others hard. As a crystal deforms, dislocations multiply and tangle; they block one another, and the crystal hardens, unless they can climb round obstacles or annihilate, which needs diffusion, hence heat.

7.2 Mechanisms of creep

Definition 7.2 (Mechanisms of creep)

In dislocation creep crystals deform by the glide of dislocations, with climb releasing them from obstacles. In diffusion creep they change shape by the diffusion of atoms, through the lattice (Nabarro–Herring creep) or along the grain boundaries (Coble creep), from the faces under the highest normal stress to those under the lowest. In grain-boundary sliding grains slide past one another along their boundaries, the gaps being closed by diffusion or dislocation creep. In pressure solution the matter dissolves in a fluid film on the faces under high stress and is deposited where the stress is low, or carried away.

Proposition 7.3 (Grain-size dependence of diffusion creep)

At low stress σ\sigma, a rock of grain size dd creeps by diffusion at the rate

ε˙≈BLDLΩkBTσd2(through the lattice),ε˙≈BBδDBΩkBTσd3(along the boundaries),\dot\varepsilon \approx B_L\frac{D_L\Omega}{k_BT}\frac{\sigma}{d^2} \quad\text{(through the lattice)},\qquad \dot\varepsilon \approx B_B\frac{\delta D_B\Omega}{k_BT}\frac{\sigma}{d^3} \quad\text{(along the boundaries)},

where DLD_L and DBD_B are the diffusion coefficients in the lattice and in the boundaries, δ\delta the width of a boundary, Ω\Omega the volume of an atom, kBk_B Boltzmann’s constant and BLB_L, BBB_B numbers of order ten that depend on the shape of the grains.

Proof. A normal stress σ\sigma on a face changes the work needed to form a vacancy there by σΩ\sigma\Omega, so the concentration CC of vacancies differs between faces under tension and under compression by ΔC≈CσΩ/(kBT)\Delta C \approx C\sigma\Omega/ (k_BT) when σΩ≪kBT\sigma\Omega \ll k_BT. Across a distance of order dd, the diffusive flux of vacancies, hence of atoms the other way, is of order DΔC/dD\Delta C/d per unit area. Through the lattice the flux crosses an area of order d2d^2, along the boundaries an area of order δd\delta d. The volume moved per unit time, divided by the grain’s volume d3d^3, is the strain rate: of order (DLσΩ/kBT) d2/(d⋅d3)=DLσΩ/(kBTd2)(D_L\sigma\Omega/k_BT)\,d^2/(d\cdot d^3) = D_L\sigma\Omega/(k_BTd^2) in the lattice, and (DBσΩ/kBT) δd/(d⋅d3)=δDBσΩ/(kBTd3)(D_B\sigma\Omega/k_BT)\,\delta d/(d\cdot d^3) = \delta D_B \sigma\Omega/(k_BTd^3) in the boundaries. The product C DC\,D is the self-diffusion coefficient written DD in the statement. ∎

Diffusion creep is linear in the stress, a true viscous flow, and it is fast in fine-grained rock: halving the grain size makes Coble creep eight times faster (figure below). Dislocation creep does not depend on the grain size, but grows with a high power of the stress.

Diffusion creep of a grain shortened along . Atoms leave the faces under the highest compression and are added to the faces under the lowest (blue), so the grain flattens and lengthens (grey arrows): (a) through the grain’s interior, Nabarro–Herring creep; (b) along its boundaries, Coble creep.
Diffusion creep of a grain shortened along σ\sigma. Atoms leave the faces under the highest compression and are added to the faces under the lowest (blue), so the grain flattens and lengthens (grey arrows): (a) through the grain’s interior, Nabarro–Herring creep; (b) along its boundaries, Coble creep.

Deformed crystals also mend themselves while they deform.

Definition 7.4 (Recovery, subgrain, dynamic recrystallisation)

Recovery is the rearrangement of dislocations, by climb and cross-slip, into lower-energy arrays, and their annihilation. Dislocations arranged in walls divide a grain into subgrains, parts of the grain whose lattices are tilted by a few degrees from one another. Dynamic recrystallisation is the replacement, during deformation, of strained grains by new grains with fewer dislocations, either by the rotation of subgrains until they become grains, or by the migration of grain boundaries into the more strained neighbours.

Definition 7.5 (Crystallographic preferred orientation)

A rock has a crystallographic preferred orientation (CPO) when the lattices of its grains are not oriented at random but cluster about directions set by the deformation; dislocation creep produces it, by rotating each grain’s slip systems towards the flow.

A rock with a CPO inherits the anisotropy of its crystals: its seismic speeds, for instance, depend on direction, which lets seismology see the fabric of the deep mantle. Diffusion creep and grain-boundary sliding, which let grains rotate freely, do not build one.

7.3 Flow laws

Definition 7.6 (Flow law, stress exponent, homologous temperature)

A flow law gives the strain rate of a rock deforming by one mechanism as a function of the differential stress σ\sigma, the grain size dd, the temperature TT and the pressure PP:

ε˙=A σnd−mexp⁡(−E+PVRT),\dot\varepsilon = A\,\sigma^nd^{-m}\exp\Bigl(-\frac{E + PV}{RT}\Bigr),

with AA, nn, mm, the activation energy EE and the activation volume VV measured in the laboratory. The stress exponent is nn, the grain-size exponent mm. The homologous temperature of a mineral is T/TmT/T_m, its absolute temperature over its absolute melting temperature.

In the laws used here, dislocation creep has nn of 3.5 or 4 and m=0m = 0; diffusion creep has n=1n = 1 and m=2m = 2 or 3 (Proposition 7.3). Creep grows as the homologous temperature approaches 1. Forsterite melts at 1890 ∘C1890\,{}^{\circ}\mathrm{C}, so the mantle at 1300 ∘C1300\,{}^{\circ}\mathrm{C} stands at T/Tm=0.73T/T_m = 0.73; ice at −10 ∘C-10\,{}^{\circ}\mathrm{C} at 0.96, which is why glaciers flow. The olivine values used here are those of a widely used laboratory compilation: dry olivine in dislocation creep, A=1.1×105A = 1.1 \times 10^{5}, n=3.5n = 3.5, E=530 kJ mol−1E = 530\,\mathrm{kJ}\,\mathrm{mol}^{-1}, V=14 cm3 mol−1V = 14\,\mathrm{cm}^{3}\,\mathrm{mol}^{-1}; in diffusion creep, A=109.2A = 10^{9.2}, n=1n = 1, m=3m = 3, E=375 kJ mol−1E = 375\,\mathrm{kJ}\,\mathrm{mol}^{-1}, V=5 cm3 mol−1V = 5\,\mathrm{cm}^{3}\,\mathrm{mol}^{-1} (stress in megapascals, grain size in micrometres, time in seconds); wet quartzite is the Year 1 volume’s. Water weakens olivine: at a water fugacity of about 300 MPa300\,\mathrm{MPa}, olivine aggregates creep five to six times faster than dry ones in dislocation creep, five times faster in diffusion creep.

Proposition 7.7 (Effective viscosity)

A rock obeying a flow law ε˙=A′σn\dot\varepsilon = A'\sigma^n at given TT, PP and dd has the effective viscosity

ηeff=σ2ε˙=12A′−1/n ε˙(1−n)/n:\eta_{\mathrm{eff}} = \frac{\sigma}{2\dot\varepsilon} = \tfrac12A'^{-1/n}\,\dot\varepsilon^{(1-n)/n}:

for n>1n > 1 it falls as the rock deforms faster.

Proof. For a Newtonian fluid, σ=2ηε˙\sigma = 2\eta\dot\varepsilon in a uniaxial flow, by the definition of viscosity; ηeff\eta_{\mathrm{eff}} is the viscosity of the Newtonian fluid that would flow at the same rate under the same stress. Substituting σ=(ε˙/A′)1/n\sigma = (\dot\varepsilon/A')^{1/n} gives the second form. ∎

Wet quartzite at 400 ∘C400\,{}^{\circ}\mathrm{C} flowing at 10−14 s−110^{-14}\,\mathrm{s}^{-1} needs 65 MPa65\,\mathrm{MPa} and has an effective viscosity of 3.3×1021 Pa s3.3 \times 10^{21}\,\mathrm{Pa}\,\mathrm{s}; at 500 ∘C500\,{}^{\circ}\mathrm{C}, 18 MPa18\,\mathrm{MPa} and 9×1020 Pa s9 \times 10^{20}\,\mathrm{Pa}\,\mathrm{s}: a hundred degrees divides the strength by more than three (figure below). Dry olivine at 1100 ∘C1100\,{}^{\circ}\mathrm{C} and 1 GPa1\,\mathrm{GPa} needs only 3 MPa3\,\mathrm{MPa}. The strength envelopes of the Year 1 volume are built from exactly these laws.

Dislocation-creep flow laws: strain rate against stress at three temperatures. The slope on these axes is the stress exponent, 4 for quartzite and 3.5 for olivine; the grey line marks a geological strain rate, 10-14\, s-1. Olivine needs temperatures some 600\, C higher than quartz to flow as easily.
Dislocation-creep flow laws: strain rate against stress at three temperatures. The slope on these axes is the stress exponent, 4 for quartzite and 3.5 for olivine; the grey line marks a geological strain rate, 10−14 s−110^{-14}\,\mathrm{s}^{-1}. Olivine needs temperatures some 600 ∘C600\,{}^{\circ}\mathrm{C} higher than quartz to flow as easily.

7.4 Deformation mechanism maps

Several mechanisms run in parallel, and the fastest carries the deformation.

Definition 7.8 (Deformation mechanism map)

A deformation mechanism map of a mineral or rock shows, on axes of two of the variables stress, grain size and temperature, the fields in which each mechanism is the fastest, with contours of the total strain rate.

Proposition 7.9 (Mechanism boundary)

At given TT and PP, dislocation creep (A1A_1, n1n_1, E1E_1, V1V_1, no grain-size dependence) and diffusion creep (A2A_2, n2=1n_2 = 1, mm, E2E_2, V2V_2) run equally fast at the stress

σb=(A2d−mA1exp⁡(E1−E2)+P(V1−V2)RT)1/(n1−1),\sigma_b = \biggl(\frac{A_2d^{-m}}{A_1} \exp\frac{(E_1 - E_2) + P(V_1 - V_2)}{RT}\biggr)^{1/(n_1 - 1)} ,

so the boundary is a straight line of slope −m/(n1−1)-m/(n_1 - 1) on logarithmic axes of grain size and stress; dislocation creep wins above it.

Proof. Equate the two flow laws and solve for σ\sigma. Above σb\sigma_b, the ratio of the two rates, proportional to σn1−1\sigma^{n_1 - 1}, exceeds 1. ∎

For dry olivine at 1100 ∘C1100\,{}^{\circ}\mathrm{C} and 2 GPa2\,\mathrm{GPa}, the boundary lies at 5 MPa5\,\mathrm{MPa} for a grain size of 1 mm1\,\mathrm{mm} and at 79 MPa79\,\mathrm{MPa} for 100 µm100\,\text{µ}\mathrm{m} (figure below). Fine-grained mantle creeps by diffusion, coarse mantle by dislocations; laboratory and geophysical evidence suggests that both regimes occur in the upper mantle, the hot, shallow mantle flowing by dislocation creep and the colder or deeper by diffusion creep, and that where the stress rises and the grains shrink, the mantle near the boundary weakens and the deformation localises.

Deformation mechanism map of dry olivine at 1100\, C and 2\, GPa. The thick line divides the fields where dislocation creep (above) and diffusion creep (below) is faster; the thin lines are contours of the total strain rate, flat where dislocation creep, which ignores the grain size, dominates.
Deformation mechanism map of dry olivine at 1100 ∘C1100\,{}^{\circ}\mathrm{C} and 2 GPa2\,\mathrm{GPa}. The thick line divides the fields where dislocation creep (above) and diffusion creep (below) is faster; the thin lines are contours of the total strain rate, flat where dislocation creep, which ignores the grain size, dominates.

The grain size that dynamic recrystallisation produces depends on the stress, and this makes the grains a stress gauge.

Definition 7.10 (Palaeopiezometer)

A palaeopiezometer is a feature of a deformed rock whose size or density is set by the differential stress under which it formed, calibrated in the laboratory, from which the stress of a past deformation can be read.

For quartz, experiments at 1.5 GPa1.5\,\mathrm{GPa} and 800 to 1100 ∘C1100\,{}^{\circ}\mathrm{C}, at flow stresses from 34 to 268 MPa268\,\mathrm{MPa}, gave recrystallised grains from 46 down to 3.2 µm3.2\,\text{µ}\mathrm{m}, fitted by

D=103.56 σ−1.26(D in µm, σ in MPa),D = 10^{3.56}\,\sigma^{-1.26}\qquad(D\text{ in $\text{µ}\mathrm{m}$},\ \sigma\text{ in MPa}),

with no detectable effect of the temperature. A quartz mylonite with recrystallised grains 20 µm20\,\text{µ}\mathrm{m} across deformed under about 62 MPa62\,\mathrm{MPa}.

7.5 Shear zones

Definition 7.11 (Shear zone, strain softening)

A shear zone is a tabular zone, between walls that are little deformed, in which the strain is concentrated and is close to a simple shear parallel to the zone. A rock undergoes strain softening when its strength falls as it deforms.

Strain softening is what keeps a shear zone narrow: once a band is weaker than its walls, the deformation stays in it. Grain-size reduction by dynamic recrystallisation, which may tip the rock into diffusion creep, the growth of a CPO that aligns easy slip systems, the entry of water and reactions that make weak minerals such as micas all soften the rock. The mylonites of the Year 1 volume are their product. In a shear zone of simple shear, a marker that crossed the zone at right angles is deflected into it: if inside the zone it makes the angle ψ\psi with its original direction, the shear strain is γ=tan⁡ψ\gamma = \tan\psi, and the displacement across a zone of width ww is γw\gamma w.

Definition 7.12 (Kinematic indicators)

A kinematic indicator is a structure whose asymmetry records the sense of shear, the sense in which one wall of a shear zone moved relative to the other, read in the plane that contains the stretching lineation and is perpendicular to the foliation. An S–C fabric is a pair of foliations: S, the foliation of flattened grains, oblique to the zone and turned towards the direction of stretching, and C, shear bands parallel to the zone. A porphyroclast system is a strong grain with tails of recrystallised material: in a σ\sigma-type the tails step from one side of the clast’s centre line to the other without crossing it, in a δ\delta-type, dragged round by the clast’s rotation, they cross it. A mica fish is a lozenge-shaped mica grain inclined to the shear plane, with thin tails, stepping like those of a σ\sigma-type clast.

Three kinematic indicators of the same dextral shear (top to the right, arrows), seen in the plane of the lineation, normal to the foliation. (a) The S foliation (red) curves into the C shear bands (blue), leaning in the direction of stretching. (b) The tails of a -type clast leave its upper side to the right and its lower side to the left, stepping up in the direction of shear without crossing the centre line (dashed). (c) A mica fish leans like S, its tails stepping like a -type clast’s. Not to scale.
Three kinematic indicators of the same dextral shear (top to the right, arrows), seen in the plane of the lineation, normal to the foliation. (a) The S foliation (red) curves into the C shear bands (blue), leaning in the direction of stretching. (b) The tails of a σ\sigma-type clast leave its upper side to the right and its lower side to the left, stepping up in the direction of shear without crossing the centre line (dashed). (c) A mica fish leans like S, its tails stepping like a σ\sigma-type clast’s. Not to scale.

Read together, a few indicators that agree give the sense of shear with confidence; a single one can mislead, since back-rotated clasts and asymmetric folds of other origins exist. The photograph in the box below shows porphyroclasts in a mylonite.

In the lab — Squeezing rock at mantle conditions

A rock-deformation press with a jacketed sample assembly (AI-generated illustration).
A sawn augen mylonite from Norway: large feldspar porphyroclasts in a fine matrix of recrystallised quartz, feldspar and micas. Photo Mikenorton, CC BY-SA 4.0.

In a gas-medium deformation apparatus, a cylinder of rock or of hot-pressed olivine powder, a centimetre or two long, is sealed in a thin metal jacket and placed with its pistons inside a furnace, in a steel vessel filled with a compressed inert gas that sets the confining pressure, a few hundred megapascals. An actuator drives the piston at a constant rate, or under a constant load, while an internal load cell measures the stress on the sample itself, free of the friction of the seals. Olivine is deformed this way at temperatures near 1200 ∘C1200\,{}^{\circ}\mathrm{C}; at the end of a run the sample is quenched, cut, polished and examined under the microscope and the electron microscope, where its grain size, subgrains and dislocations are measured.

7.6 Exercises

Exercise 7.1 ★

Forsterite melts at 1890 ∘C1890\,{}^{\circ}\mathrm{C} and ice at 0 ∘C0\,{}^{\circ}\mathrm{C}. Compute the homologous temperature of forsterite at 1300 ∘C1300\,{}^{\circ}\mathrm{C} and of ice at −10 ∘C-10\,{}^{\circ}\mathrm{C}. Which creeps more readily?

Solution

Solution of Exercise 7.1.

Forsterite: 1573/2163=0.731573/2163 = 0.73; ice: 263/273=0.96263/273 = 0.96. Ice is far closer to its melting point and creeps more readily: glaciers flow at the Earth’s surface.

Exercise 7.2 ★

A dyke crossing a shear zone at right angles is deflected inside the zone by 70∘70{}^{\circ}, and another by 80∘80{}^{\circ}. Compute the shear strains. Why are large strains hard to measure this way?

Solution

Solution of Exercise 7.2.

γ=tan⁡70∘=2.75\gamma = \tan70{}^{\circ} = 2.75 and tan⁡80∘=5.67\tan80{}^{\circ} = 5.67. Near 90∘90{}^{\circ} the tangent changes very fast with the angle: one degree of error at ψ=84∘\psi = 84{}^{\circ} changes γ\gamma by a fifth.

Exercise 7.3 ★

With the quartz piezometer, what stress does a recrystallised grain size of 10 µm10\,\text{µ}\mathrm{m} indicate?

Solution

Solution of Exercise 7.3.

σ=(103.56/10)1/1.26=363.10.794=108 MPa\sigma = (10^{3.56}/10)^{1/1.26} = 363.1^{0.794} = 108\,\mathrm{MPa}.

Exercise 7.4 ★

With the flow law of wet quartzite (A=1.1×10−4A = 1.1 \times 10^{-4}, n=4n = 4, E=223 kJ mol−1E = 223\,\mathrm{kJ}\,\mathrm{mol}^{-1}, R=8.314 J mol−1 K−1R = 8.314\,\mathrm{J}\,\mathrm{mol}^{-1}\,\mathrm{K}^{-1}), compute the strain rate at 500 ∘C500\,{}^{\circ}\mathrm{C} under 20 MPa20\,\mathrm{MPa}.

Solution

Solution of Exercise 7.4.

E/RT=223000/(8.314×773.15)=34.69E/RT = 223000/(8.314\times773.15) = 34.69, so ε˙=1.1×10−4×204×e−34.69=1.1×10−4×1.6×105×8.6×10−16=1.5×10−14 s−1\dot\varepsilon = 1.1 \times 10^{-4}\times20^4\times e^{-34.69} = 1.1 \times 10^{-4}\times1.6 \times 10^{5}\times 8.6 \times 10^{-16} = 1.5 \times 10^{-14}\,\mathrm{s}^{-1}.

Exercise 7.5 ★★

Compute the effective viscosity of wet quartzite flowing at 10−14 s−110^{-14}\,\mathrm{s}^{-1} at 400 ∘C400\,{}^{\circ}\mathrm{C} (65.4 MPa65.4\,\mathrm{MPa}) and at 500 ∘C500\,{}^{\circ}\mathrm{C} (18.0 MPa18.0\,\mathrm{MPa}). Compare with the mantle viscosities of Chapter 6.

Solution

Solution of Exercise 7.5.

η=σ/(2ε˙)\eta = \sigma/(2\dot\varepsilon): 65.4×106/2×10−14=3.3×1021 Pa s65.4 \times 10^{6}/2 \times 10^{-14} = 3.3 \times 10^{21}\,\mathrm{Pa}\,\mathrm{s} at 400 ∘C400\,{}^{\circ}\mathrm{C} and 18.0×106/2×10−14=9.0×1020 Pa s18.0 \times 10^{6}/2 \times 10^{-14} = 9.0 \times 10^{20}\,\mathrm{Pa}\,\mathrm{s} at 500 ∘C500\,{}^{\circ}\mathrm{C}: the middle crust at these temperatures is about as viscous as the mantle beneath the ice sheets.

Exercise 7.6 ★★

Two experiments at the same temperature give 10−5 s−110^{-5}\,\mathrm{s}^{-1} under 100 MPa100\,\mathrm{MPa} and 1.6×10−4 s−11.6 \times 10^{-4}\,\mathrm{s}^{-1} under 200 MPa200\,\mathrm{MPa}. Compute the stress exponent and name the likely mechanism.

Solution

Solution of Exercise 7.6.

n=log⁡(1.6×10−4/10−5)/log⁡2=log⁡16/log⁡2=4n = \log(1.6 \times 10^{-4}/10^{-5})/\log2 = \log16/\log2 = 4: dislocation creep, whose stress exponent is well above 1.

Exercise 7.7 ★★

At the same stress, a sample creeps ten times faster at 1100 K1100\,\mathrm{K} than at 1000 K1000\,\mathrm{K}. Compute the activation energy.

Solution

Solution of Exercise 7.7.

ln⁡10=(E/R)(1/1000−1/1100)\ln10 = (E/R)(1/1000 - 1/1100), so E=8.314×2.3026/9.091×10−5=211 kJ mol−1E = 8.314\times2.3026/9.091 \times 10^{-5} = 211\,\mathrm{kJ}\,\mathrm{mol}^{-1}.

Exercise 7.8 ★★

By what factor does the strain rate of diffusion creep change when the grain size is divided by ten, for Nabarro–Herring and for Coble creep? Why is Coble creep the more sensitive?

Solution

Solution of Exercise 7.8.

Nabarro–Herring (d−2d^{-2}): 100 times faster; Coble (d−3d^{-3}): 1000 times. Along the boundaries the diffusing atoms use a path whose cross-section, δd\delta d, shrinks with the grain, on top of the shorter distance and the smaller volume to move.

Exercise 7.9 ★★

For dry olivine at 1100 ∘C1100\,{}^{\circ}\mathrm{C} and 2 GPa2\,\mathrm{GPa}, the two flow laws give ε˙dis=6.54×10−17 σ3.5\dot\varepsilon_{\mathrm{dis}} = 6.54 \times 10^{-17}\,\sigma^{3.5} and ε˙dif=3.59×10−6 σd−3\dot\varepsilon_{\mathrm{dif}} = 3.59 \times 10^{-6}\,\sigma d^{-3} (units of the text). Compute the boundary stress for d=1 mmd = 1\,\mathrm{mm} and d=100 µmd = 100\,\text{µ}\mathrm{m}, and the slope of the boundary.

Solution

Solution of Exercise 7.9.

σb2.5=3.59×10−6 d−3/6.54×10−17\sigma_b^{2.5} = 3.59 \times 10^{-6}\,d^{-3}/6.54 \times 10^{-17}: for d=1000d = 1000, σb=54.90.4=5.0 MPa\sigma_b = 54.9^{0.4} = 5.0\,\mathrm{MPa}; for d=100d = 100, σb=54 9000.4=79 MPa\sigma_b = 54\,900^{0.4} = 79\,\mathrm{MPa}. The slope is −3/2.5=−1.2-3/2.5 = -1.2: a tenfold finer grain raises the boundary by 101.2=1610^{1.2} = 16.

Exercise 7.10 ★★★

Following the proof of Proposition 7.3, explain why the lattice path gives d−2d^{-2} and the boundary path d−3d^{-3}, and show that boundary diffusion wins below the grain size d∗≈(BB/BL) δDB/DLd^* \approx (B_B/B_L)\,\delta D_B/D_L.

Solution

Solution of Exercise 7.10.

In both paths the driving difference of concentration is the same and the distance is dd; the volume to move per unit strain is d3d^3. The lattice offers a cross-section ∝d2\propto d^2, so the rate goes as d2/(d d3)=d−2d^2/(d\,d^3) = d^{-2}; the boundaries offer ∝δd\propto\delta d, so δd/(d d3)=δd−3\delta d/(d\,d^3) = \delta d^{-3}. The ratio of boundary to lattice rates is (BBδDB)/(BLDLd)(B_B\delta D_B)/(B_LD_Ld), larger than 1 when d<(BB/BL) δDB/DLd < (B_B/B_L)\,\delta D_B/D_L.

Exercise 7.11 ★★★

For a flow law with n=3.5n = 3.5, by what factor does the effective viscosity fall when the strain rate is multiplied by ten? Why does this favour the localisation of strain into shear zones?

Solution

Solution of Exercise 7.11.

η∝ε˙(1−n)/n=ε˙−2.5/3.5\eta \propto \dot\varepsilon^{(1-n)/n} = \dot\varepsilon^{-2.5/3.5}: ten times faster divides the viscosity by 100.714=5.210^{0.714} = 5.2. A band that happens to deform faster becomes weaker, so it takes more of the deformation: the strain gathers into it.

Exercise 7.12 ★★★

In the dry olivine of exercise 9, a mantle shear zone under 10 MPa10\,\mathrm{MPa} recrystallises from 3 mm3\,\mathrm{mm} to 30 µm30\,\text{µ}\mathrm{m} grains. Compute the total strain rate before and after, and the factor of softening. Which mechanism carries the deformation in each case?

Solution

Solution of Exercise 7.12.

Before: ε˙dis=6.54×10−17×103.5=2.07×10−13 s−1\dot\varepsilon_{\mathrm{dis}} = 6.54 \times 10^{-17}\times10^{3.5} = 2.07 \times 10^{-13}\,\mathrm{s}^{-1} and ε˙dif=3.59×10−6×10/30003=1.3×10−15 s−1\dot\varepsilon_{\mathrm{dif}} = 3.59 \times 10^{-6} \times10/3000^3 = 1.3 \times 10^{-15}\,\mathrm{s}^{-1}: 2.08×10−13 s−12.08 \times 10^{-13}\,\mathrm{s}^{-1}, dislocation creep. After: ε˙dif=3.59×10−6×10/303=1.33×10−9 s−1\dot\varepsilon_{\mathrm{dif}} = 3.59 \times 10^{-6}\times10/30^3 = 1.33 \times 10^{-9}\,\mathrm{s}^{-1}, diffusion creep, 64006400 times faster at the same stress: recrystallisation has softened the zone enormously.

7.7 Problem: The Varden Shear Zone

Problem 7.1

Weekend problem — the strain, stress, strain rate and age of an invented shear zone

Beneath the granite of the Varden Hills, the series’ invented region, a borehole and a deep road cut (both invented) reveal a ductile shear zone 400 m400\,\mathrm{m} wide, formed at about 15 km15\,\mathrm{km} depth and 450 ∘C450\,{}^{\circ}\mathrm{C}. A dolerite dyke that crosses the walls at right angles makes, inside the zone, an angle of only 6∘6{}^{\circ} with the zone’s boundaries. Its quartz is a mylonite of recrystallised grains 20 µm20\,\text{µ}\mathrm{m} across (from 15 to 25 µm25\,\text{µ}\mathrm{m}). Take the quartz piezometer D=103.56σ−1.26D = 10^{3.56}\sigma^{-1.26}, the flow law of wet quartzite (A=1.1×10−4A = 1.1 \times 10^{-4}, n=4n = 4, E=223 kJ mol−1E = 223\,\mathrm{kJ}\,\mathrm{mol}^{-1}), R=8.314 J mol−1 K−1R = 8.314\,\mathrm{J}\,\mathrm{mol}^{-1}\,\mathrm{K}^{-1}, 1 yr1\,\mathrm{yr} =3.156×107 s= 3.156 \times 10^{7}\,\mathrm{s}, and, in simple shear, γ˙=2ε˙\dot\gamma = 2\dot\varepsilon.

Part I — The strain.

  1. Compute the deflection ψ\psi of the dyke and the shear strain.
  2. Compute the displacement of one wall relative to the other.
  3. Compute the principal stretches of the strain ellipse and its axial ratio.
  4. Compute the angle of the long axis with the zone. What does the foliation of the zone do as γ\gamma grows?
  5. Seen from above, the dyke, followed from the southern wall, bends to the east as it enters the east–west zone, then returns to its northward course at the northern wall. Which way did the northern wall move, and is the zone dextral or sinistral?
  6. In which plane should a thin section be cut to read the kinematic indicators, and how should its S foliation lean?

Part II — The stress.

  1. Compute the stress indicated by the 20 µm20\,\text{µ}\mathrm{m} grains.
  2. Compute the stresses for 15 and 25 µm25\,\text{µ}\mathrm{m}.
  3. Is the stress within the range of the calibration?
  4. The piezometer showed no effect of temperature. Why does this matter here?
  5. What grain size would a stress of 100 MPa100\,\mathrm{MPa} have produced?

Part III — The strain rate.

  1. Compute E/RTE/RT at 450 ∘C450\,{}^{\circ}\mathrm{C} and the Arrhenius factor.
  2. Compute the strain rate of the quartz under the stress of question 7 (take 62 MPa62\,\mathrm{MPa}).
  3. Compute the shear strain rate.
  4. Compute the effective viscosity of the zone and compare it with the upper mantle of Chapter 6.
  5. By what factors would the strain rate change at 400 and 500 ∘C500\,{}^{\circ}\mathrm{C}, at the same stress?
  6. By what factor if the stress were the 78 MPa78\,\mathrm{MPa} of the finest grains?

Part IV — The time.

  1. Compute the time needed to accumulate the shear strain of question 1, in seconds.
  2. Convert it to thousands of years.
  3. Compute the mean slip rate across the zone, in mm/yr\mathrm{mm}/\mathrm{yr}.
  4. Compare it with plate-boundary rates of tens of millimetres a year.
  5. If the zone had carried 20 mm/yr20\,\mathrm{mm}/\mathrm{yr}, what shear strain rate, stress and recrystallised grain size would it have had?
  6. How do the uncertainties of the grain size and of the piezometer’s exponent affect the time?
  7. Give the time the Varden shear zone took to accumulate its strain, in thousands of years.
Solution

Solution of Problem 7.1.

1. The dyke turns from 90∘90{}^{\circ} to 6∘6{}^{\circ} from the zone: ψ=84∘\psi = 84{}^{\circ} and γ=tan⁡84∘=9.51\gamma = \tan84{}^{\circ} = 9.51.

2. γw=9.51×400=3804 m\gamma w = 9.51\times400 = 3804\,\mathrm{m}, about 3.8 km3.8\,\mathrm{km}.

3. S1=1+9.512/4+9.51/2=4.859+4.755=9.614S_1 = \sqrt{1 + 9.51^2/4} + 9.51/2 = 4.859 + 4.755 = 9.614, S2=1/S1=0.104S_2 = 1/S_1 = 0.104, and S1/S2=92.4S_1/S_2 = 92.4.

4. tan⁡2θ′=2/9.51=0.210\tan2\theta' = 2/9.51 = 0.210, θ′=5.9∘\theta' = 5.9{}^{\circ}: the foliation, close to the long axis, turns ever closer to the zone as γ\gamma grows.

5. Followed northwards, the dyke is carried east inside the zone and leaves it further east: the northern wall moved east relative to the southern. Seen from the southern wall, the far side moved to the right: the zone is dextral.

6. In the plane containing the stretching lineation, horizontal here for a strike-slip zone, and perpendicular to the vertical foliation: a horizontal section. Its S foliation leans anticlockwise from the zone, towards east-north-east, like the long axis of the strain ellipse.

7. σ=(103.56/20)1/1.26=181.50.794=62 MPa\sigma = (10^{3.56}/20)^{1/1.26} = 181.5^{0.794} = 62\,\mathrm{MPa}.

8. 242.10.794=78 MPa242.1^{0.794} = 78\,\mathrm{MPa} for 15 µm15\,\text{µ}\mathrm{m} and 145.20.794=52 MPa145.2^{0.794} = 52\,\mathrm{MPa} for 25 µm25\,\text{µ}\mathrm{m}.

9. Yes: the experiments spanned 34 to 268 MPa268\,\mathrm{MPa}.

10. The zone deformed at 450 ∘C450\,{}^{\circ}\mathrm{C}, far below the 800 to 1100 ∘C1100\,{}^{\circ}\mathrm{C} of the experiments; a relation that does not depend on the temperature can be carried down to it, though it is still an extrapolation.

11. D=3631/1001.26=3631/331=11 µmD = 3631/100^{1.26} = 3631/331 = 11\,\text{µ}\mathrm{m}.

12. E/RT=223000/(8.314×723.15)=37.09E/RT = 223000/(8.314\times723.15) = 37.09 and e−37.09=7.79×10−17e^{-37.09} = 7.79 \times 10^{-17}.

13. ε˙=1.1×10−4×624×7.79×10−17=1.1×10−4×1.478×107×7.79×10−17=1.27×10−13 s−1\dot\varepsilon = 1.1 \times 10^{-4}\times62^4\times7.79 \times 10^{-17} = 1.1 \times 10^{-4}\times1.478 \times 10^{7}\times7.79 \times 10^{-17} = 1.27 \times 10^{-13}\,\mathrm{s}^{-1}.

14. γ˙=2ε˙=2.53×10−13 s−1\dot\gamma = 2\dot\varepsilon = 2.53 \times 10^{-13}\,\mathrm{s}^{-1}.

15. η=62×106/(2×1.27×10−13)=2.4×1020 Pa s\eta = 62 \times 10^{6}/(2\times1.27 \times 10^{-13}) = 2.4 \times 10^{20}\,\mathrm{Pa}\,\mathrm{s}, close to the 3×1020 Pa s3 \times 10^{20}\,\mathrm{Pa}\,\mathrm{s} of the upper mantle beneath Lake Ardo: the zone was as weak as the mantle.

16. e−E/RTe^{-E/RT} changes by the factor 0.064 at 400 ∘C400\,{}^{\circ}\mathrm{C} (sixteen times slower) and 11.0 at 500 ∘C500\,{}^{\circ}\mathrm{C}.

17. (78/62)4=2.5(78/62)^4 = 2.5 times faster.

18. t=9.51/2.53×10−13=3.75×1013 st = 9.51/2.53 \times 10^{-13} = 3.75 \times 10^{13}\,\mathrm{s}.

19. 3.75×1013/3.156×1010=1190 kyr3.75 \times 10^{13}/3.156 \times 10^{10} = 1190\,\mathrm{kyr}.

20. 3804 m/1.19×106 yr=3.2 mm/yr3804\,\mathrm{m}/1.19 \times 10^{6}\,\mathrm{yr} = 3.2\,\mathrm{mm}/\mathrm{yr}.

21. Ten or more times slower than the main plate boundaries: the zone was a minor strand, or one of several that shared a boundary’s motion.

22. γ˙=0.02/(3.156×107×400)=1.58×10−12 s−1\dot\gamma = 0.02/(3.156 \times 10^{7}\times400) = 1.58 \times 10^{-12}\,\mathrm{s}^{-1}, 6.2 times faster; the stress rises by 6.21/4=1.586.2^{1/4} = 1.58, to 98 MPa98\,\mathrm{MPa}, and the grains would be 3631/981.26=11 µm3631/98^{1.26} = 11\,\text{µ}\mathrm{m}.

23. The time goes as σ−4\sigma^{-4}: the grain-size range alone gives (62/78)4=0.40(62/78)^4 = 0.40 to (62/52)4=2.0(62/52)^4 = 2.0 times the estimate, and the uncertainty of the exponent 1.26±0.131.26 \pm 0.13 and of the flow law’s constants widens it further: the result is an order of magnitude, not a date.

24. The Varden shear zone took about 1190 kyr1190\,\mathrm{kyr} to accumulate its strain.

本章定义的术语

在术语表中查看全部 957 个术语