Smectite to Illite: A Thermal Clock

Part 4, Part 4: Pore Pressure

Learning objectives

Describe smectite and illite as the same 2:1 sheet differing in the interlayer: hydrated exchangeable cations at 15 angstrom against bare potassium keyed at 10
Read mixed-layer illite-smectite as the proportion of individually collapsed interlayers, and connect that proportion to the basal spacing a diffractogram measures
Separate the thermodynamics, which decides that illite wins, from the kinetics, which decides when, and write the Arrhenius rate that gates the reaction
Explain why the reaction half-life falls through seven orders of magnitude between 40 and 180 degrees, making the transition a window rather than a switch
Predict the depth of the transition from geothermal gradient and burial rate, and state why it tracks an isotherm near 100 degrees rather than a fixed depth

One Sheet, Two Minerals

Smectite and illite are built from the same unit. Both are 2:1 phyllosilicates: a sheet of aluminium-oxygen octahedra sandwiched between two sheets of silicon-oxygen tetrahedra, the T-O-T layer that Part 3 met as the reason clay is mechanically weak. Stack those sandwiches and the mineral is defined by what sits between them, and that is where the two clays part company. Substitutions inside the lattice leave each layer with a net negative layer charge, which some cation must balance. In smectite the charge is low, roughly 0.4 per half cell, and it is balanced by exchangeable cations that arrive with their hydration shells intact. Water comes into the interlayer with them, props the sheets apart to a basal spacing $d_{001}$ 001 near 15 angstrom, and can leave or return: this is why smectite swells, why it has the high cation exchange capacity that makes it a drilling nuisance, and why it holds a reservoir of water that is neither free porosity nor mineral. In illite the layer charge is higher, near 0.9, and it is balanced by potassium, which fits the ditrigonal cavity in the tetrahedral sheet so precisely that it sheds its hydration shell and keys the two sheets together. The interlayer water is gone, the repeat collapses to 10 angstrom, and the structure no longer swells. Same sandwich, different filling, and every property that matters to a geomechanicist follows from that one difference.

The layer panel makes the point the usual textbook figure misses. Real mudstone almost never holds pure smectite or pure illite; it holds mixed-layer illite-smectite, a single crystal stack in which individual interlayers are either collapsed or expanded. The illite fraction is not a blend of two minerals sitting side by side, it is the proportion of interlayers in one stack that have taken their potassium and closed. That is why the reaction reads as a continuous variable on a log, and why X-ray diffraction returns a basal spacing between 10 and 15 angstrom rather than two separate peaks: the diffractogram is averaging over a stack that is partly shut. Watch an interlayer close: the water leaves the crystal in the moment the potassium keys in, the stack settles, and the bracket beside it measures the loss. That lost height is water the mudstone must now find room for, and the next section is entirely about where it goes.

Thermodynamics Decides Whether, Kinetics Decides When

Two different questions get confused here, and separating them is the whole framework. Thermodynamics answers whether the reaction should happen: at the temperatures and pressures of a buried basin, illite plus quartz has the lower Gibbs free energy, so smectite is metastable, holding a structure it has no thermodynamic right to keep. The reaction, written schematically, spends potassium and returns silica: $\text{smectite} + \text{K}^{+} + \text{Al}^{3+} \rightarrow \text{illite} + \text{Si}^{4+} + \text{H}_2\text{O}$ , with the potassium supplied by dissolving K-feldspar and mica, and the released silica precipitating as the quartz cement that stiffens deep shale. If thermodynamics were the whole story, every buried smectite would already be illite.

It is not, because the conversion must dismantle and rebuild a crystal lattice, and that path runs over an energy barrier. Kinetics answers when, through the Arrhenius law: the rate carries an activation energy $E_a$ a near 33 kilocalories per mole, so $k = A\,e^{-E_a/RT}$ a/RT, and the smectite fraction $S$ in the mixed-layer stack follows $\frac{dS}{dt} = -k\,S^{5}$ . The exponent is empirical and the fifth order is the classical basin-modelling form, but the exponential is the physics, and the clock panel shows what it costs. Hold a mudstone at 40 degrees and half its smectite converts in about 6800 million years, longer than the Earth has existed: the reaction is thermodynamically favoured and utterly forbidden in practice. Raise it to 100 degrees and the same half-conversion takes a little over a million years, brisk on basin time. At 180 degrees it takes a thousandth of that. Seven orders of magnitude across a hundred and forty degrees, and that steepness is exactly why the transition occupies a window: it is not a switch that trips at a threshold, it is a clock whose hands accelerate.

So Where Does It Happen?

Now the question the pore-pressure analyst actually asks, and the reason this section exists before the pressure one. A rock does not experience a temperature, it experiences a history of temperature, so conversion is the Arrhenius rate integrated along the burial path. For a layer buried at rate $v$ through a gradient $G$ , every metre of burial contributes $dz/v$ of time at the temperature of that depth, and the burial rate enters as $1/v$ : bury faster and the rock arrives deep having spent less time hot. The profile panel integrates exactly that, and its right edge carries the basin's own temperature scale so both coordinates stay in view at once. Two consequences follow, and they are different in size. Change the geothermal gradient from 20 to 45 degrees per kilometre and the half-converted depth moves from about 3.8 km to 1.9 km, a swing of two kilometres, while the temperature at that depth barely stirs, 97 to 103 degrees. Change the burial rate sixteenfold and the half-conversion temperature shifts by only about twenty degrees, because an exponential rate needs very little extra heat to compensate for far less time.

The lesson is worth stating plainly, because it is the thing practitioners carry and beginners get backwards. The reaction tracks an isotherm, not a depth. Somewhere near 100 degrees, in most basins, the mudstone is about half converted; the depth of the transition is simply wherever the basin happens to put that isotherm, which is why the same reaction is a shallow drilling problem on a hot margin and an irrelevance in a cold one. Burial rate is the second-order correction that trades time against temperature, and potassium supply is the third: starve the reaction of K-feldspar and it stalls, pushing the transition about half a kilometre deeper even with the heat unchanged. Carry those three levers and you can predict, before a well is drilled, the interval where the clay is rebuilding itself. The next section takes that interval and asks the question this course cares about: the reaction is expelling water and weakening the frame, so what happens to the pressure?

References

Pytte, A. M., & Reynolds, R. C. (1989). The thermal transformation of smectite to illite. In Thermal History of Sedimentary Basins, Springer, 133-140.
Hower, J., Eslinger, E. V., Hower, M. E., & Perry, E. A. (1976). Mechanism of burial metamorphism of argillaceous sediment. GSA Bulletin, 87(5), 725-737.
Huang, W.-L., Longo, J. M., & Pevear, D. R. (1993). An experimentally derived kinetic model for smectite-to-illite conversion. Clays and Clay Minerals, 41(2), 162-177.
Lahann, R. W., & Swarbrick, R. E. (2011). Overpressure generation by load transfer following shale framework weakening due to smectite diagenesis. Geofluids, 11(4), 362-375.