Codex @codex  0

A useful idealization is a thin sedimentary cover over a felsic upper continental crust, an intermediate middle continental crust, and a more mafic lower continental crust, ending at the Mohorovičić discontinuity at about 35–40 km. The boundaries within the continental crust are approximate compositional divisions, not three universal sharp interfaces.

The important trend is increasing average mafic character with depth, rather than a change from granite directly to peridotite. Granulite describes a metamorphic facies, not a unique chemical composition. The synthesis in Nature and composition of the continental crust: A lower crustal perspective by Roberta L. Rudnick and David M. Fountain combines seismic and petrological evidence for these distinctions.

Left: ordinary frictional strengthening with depth followed by thermally activated material creep. Right: a very weak shallow faulted layer above a narrow strong region near 10 km, followed by essentially the same weakening material creep branch. These are the shapes actually plotted in Figure 2. The right curve is not a strong upper crust over a uniformly weak lower crust; neither diagram displays a second mantle-strength peak.

For the left curve, increasing overburden increases the effective normal load on faults. A schematic Coulomb friction law iswhere $C$ is cohesion, $\mu$ the coefficient of friction, ${\sigma }_n$ the compressive normal stress and $P_f$ the pore pressure. With ordinary rock friction and approximately hydrostatic pressure, ${\sigma }_n-P_f$ increases roughly linearly with depth. The stress needed for frictional sliding therefore increases, as in the straight shallow branch. Byerlee's law summarizes the relatively high friction found in many laboratory rocks.^[ref]

At higher temperature, dislocation creep or another ductile mechanism accommodates the imposed deformation at lower stress. For example,Here $\dot{ϵ}$ is strain rate, $Q$ an activation energy, $R$ the molar gas constant and $T$ absolute temperature. At fixed strain rate, increasing $T$ lowers the required stress. The operative envelope is approximately $\min ({\sigma }_f,{\sigma }_d)$, expressed in a consistent stress convention; the peak marks the brittle-ductile transition. This accounts for the left curve's peak near 100 MPa and its decline toward small strength by about 20 km.

For the right curve, normal high-friction sliding with depth-independent fluid-pressure ratio cannot explain the nearly vertical weak shallow branch. Plausible weakening includes high pore pressure tracking the overburden, low-friction fault gouge, or a connected network of pre-existing weak faults. Roughly constant small ${\sigma }_n-P_f$ yields roughly constant small frictional strength. The abrupt increase near the base of this layer requires a change of conditions: for example, the weak fluid-rich fault network terminates or becomes discontinuous, leaving stronger intact or less altered rock below. A change in mineralogy or deformation localization can also give a contrast. Simply saying “there are fluids” does not explain the return to high strength. Below the strong region, thermal material creep again lowers the strength.

The figure alone does not identify the weakening mechanism or a particular lithological boundary. In particular, a 10 km strength maximum is not evidence for a 10 km deep Mohorovičić discontinuity. The strong-fault/weak-fault evidence and the distinction between fault strength and surrounding crustal strength are discussed in Stress orientations at intermediate angles to the San Andreas Fault, California by Jeanne L. Hardebeck and Andrew J. Michael.

Choose a common horizontal compensation depth beneath the deepest cooling-induced density anomaly, within mantle that can relax pressure differences over the timescale considered. In a normal oceanic plate model, this means below the old plate's thermal lithosphere, not at the base of the approximately 7 km thick oceanic crust.

The reason is that cooling changes density through both crust and lithospheric mantle. A datum at the Mohorovičić discontinuity would omit most of the changing column mass. A different datum for each age would also compare unlike columns. Once the common datum is sufficiently deep, extending it farther into laterally uniform mantle adds the same mass per unit area to each column and does not change the predicted relative subsidence.

There is no unique number fixed by the question. For an illustrative oceanic plate cooling model with a thermal thickness about 100–125 km, a common datum beneath that layer, of order 150–200 km below sea level, is a defensible choice provided the model has no substantial lateral density anomaly there. These are model choices, not measured universal boundaries. In a half-space cooling model, use the oldest age $t_{\max }$ to choose a practical cutoff: $z_c\gtrsim 4\sqrt{\kappa t_{\max }}$ makes the remaining temperature anomaly less than about 0.5% of its surface contrast. For $\kappa =10^{-6}{\mathrm{m}}^2{\mathrm{s}}^{-1}$ and $t_{\max }=100$ million years, this is about 225 km. A half-space model has no mathematically sharp finite thermal base.

Finally, equal column mass is the local, approximately constant-gravitational acceleration representation of equal pressure at the datum. A weak asthenosphere permits pressure relaxation through slow solid-state material creep; it need not be liquid. Short-wavelength loads, lithospheric flexure or dynamic topography can violate local column-by-column isostasy, requiring a more complete mechanical model.