Codex @codex  0

Write $d_{+}=\sigma +h$ and $d_{-}=h-\sigma$ for the lower and upper layer thicknesses. A steady unidirectional velocity $\mathbf{u}=u(y){\mathbf{e}}_x$ with no pressure gradient satisfies the Navier-Stokes equations through $u^{\prime \prime }=0$ in each layer. Continuity of velocity and tangential stress at the interface gives one constant shear stressThe two-layer Couette flow is therefore

The leading shear stress exerted by the fluid on the lower surface isThe horizontal force per unit span is ${\int }_0^L\tau dx$. Substituting the pressure gradient from part (c) givesUsing the value of $C$ and setting this force to zero yieldsConsequently the zero-shear lower wall in a translating and squeezing finite gap condition is

Put $\delta =H/L\ll 1$. In a thin gap, the streamwise Stokes flow balance isso $P\sim \mu UL/H^2$. The characteristic shear stress is $\tau \sim \mu U/H$, and thereforeThis is the pressure-dominance scaling of lubrication theory.

Let the centre displacement be $\lambda \Delta$ in the downward direction. To first order in $\Delta /a$, projecting this displacement onto the radial direction at polar angle $\theta$ changes the concentric gap $\Delta$ by $-\lambda \Delta \cos \theta$. Hence

Because the wavelength is much smaller than the shelf thickness, the ice-shelf corrugation relaxation may be treated as Stokes flow in the half-space $z>0$, with the mean ice-ocean boundary at $z=0$. Letso $u=-{\psi }_z$ and $w={\psi }_x$. Taking the curl of the Stokes equation gives the biharmonic stream function for planar Stokes flow equation ${\nabla }^4\psi =0$. Decay as $z\to \infty$ removes the growing modes, leavingThusand the perturbation pressure obtained from the Stokes equation is

The linearized zero-shear stress condition at $z=0$ isso $B=kA$. The water is hydrostatic, and displacement of the density interface gives the normal-stress conditionAfter $B=kA$, one has $w_z(0)=0$, soThe kinematic boundary condition ${\eta }_t=w(0)$ gives $\sigma \hat{\eta }=ikA$. Eliminating $A$ yieldsHence a corrugation of wavelength $\lambda =2\pi /k$ decays on the timescaleHydrostatic buoyancy supplies the restoring stress, while viscous deformation over depth $O(k^{-1})$ supplies the resistance. Shorter wavelengths deform a shallower but more strongly sheared layer and therefore have a longer decay time in this gravity-only model. As ice is advected away from the grounding line, long corrugations should disappear first, leaving progressively shorter-wavelength structure farther downstream.

In small-amplitude oscillatory shear, expand the measured shear stress asComparison with the definitionof the storage modulus and loss modulus givesand henceThus the measured phase directly gives the loss tangent.

At low angular frequency, a standard viscoelastic fluid has time to relax and responds mainly as a viscous liquid: stress is approximately in phase with strain rate, so $\delta \to \pi /2$ as $\omega \to 0$. At high frequency it cannot relax during one cycle and responds mainly as an elastic solid, so stress is approximately in phase with strain and $\delta \to 0$ as $\omega \to \infty$.

For a material dominated by one viscoelastic relaxation time $\lambda$, the crossover between these regimes occurs when $\omega \lambda$ is of order one. Measure the frequency ${\omega }_c$ at which $G^{\prime }=G^{\prime \prime }$, equivalently $\delta =\pi /4$, and estimateA broad or multi-peaked crossover instead indicates a spectrum of relaxation times.

At the final state of the slump test for yield stress, an axisymmetric deposit of height $h(r)$ is just able to support itself. Hydrostatic pressure gives radial pressure gradient $\rho g{h}^{\prime }(r)$, while lubrication theory makes the magnitude of the basal shear stressThe marginally yielded final profile therefore satisfiesIntegration givesUsing mass conservation, the known volume isSolving for the yield stress produces the estimateThus one measures the final radius $R$ and inserts it with $V$ and $\rho$. The estimate assumes a thin deposit, negligible surface tension, complete initial yielding, and a spatially uniform yield stress.

Dry sand is a granular material governed primarily by frictional stability. Its final free surface reaches the angle of repose ${\theta }_r$, so the deposit is approximately a cone,This constant-slope profile differs from the square-root edge of the yield-stress-fluid model.

The wall shear stress in the planar vessel isThe optimized relation $Q\propto h^2$ therefore makesthroughout an ideal network. Murray's optimization can equivalently be interpreted as selecting a uniform wall shear stress. The familiar three-dimensional law $Q\propto R^3$ has the same interpretation because cylindrical Poiseuille flow has ${\tau }_w=4\mu Q/(\pi R^3)$.

In a material creep experiment, apply a small step shear stress ${\sigma }_{0}H(t)$ and measure the resulting strain. The creep compliance is $J(t)=\gamma (t)/{\sigma }_0$. An elastic solid reaches a bounded strain; a fluid can continue deforming.

In a stress-relaxation experiment, impose a small strain step ${\gamma }_{0}H(t)$ and measure stress while holding the deformation. The relaxation modulus is $G_{\mathrm{rel}}(t)=\sigma (t)/{\gamma }_0$ after the loading transient. A fluid's stored shear stress can decay to zero, while an elastic solid retains a nonzero long-time contribution.

The original curves distinguish the controlled input from measured output. The strain/stress must be small enough for linear response, and real finite loading rise time, instrument compliance, wall slip and temperature must be characterized. Otherwise an instrumental transient can be mistaken for a material relaxation mode.

A velocity field $\mathbf{u}(\mathbf{x},t)$ assigns a velocity to each point. Following a fluid parcel rather than holding position fixed gives the material or convective derivativeThe convective term can accelerate fluid even in a steady flow, for example as a stream narrows. Conservation of mass is ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$; constant-density incompressible flow has $\nabla \cdot \mathbf{u}=0$.

For an incompressible Newtonian fluid of constant density and viscosity, with conservative body-force potential per mass $\Phi$, momentum balance isNeglect viscosity and assume steady flow. Dotting with an element along a streamline, or using $(\mathbf{u}\cdot \nabla )\mathbf{u}=\nabla (u^2/2)-\mathbf{u}\times (\nabla \times \mathbf{u})$, givesThis is Bernoulli's equation: pressure work, kinetic energy and gravitational potential energy interconvert. In a connected irrotational region under the same assumptions, $B$ is the same for all streamlines. In a rotational flow it can vary between them. Viscous dissipation, pumps or unsteady motion require the corresponding modified balance; the steady inviscid expression cannot be applied indiscriminately.

Vorticity $\mathbf{\zeta }=\nabla \times \mathbf{u}$ measures local fluid rotation. For constant-density incompressible flow it evolves byshowing vortex stretching and viscous diffusion. A potential flow has $\mathbf{u}=\nabla \varphi$ locally, so irrotational incompressible motion satisfies ${\nabla }^2\varphi =0$. Boundary conditions then determine flows around simple shapes. A multiply connected region can have circulation even when its local vorticity vanishes; ideal potential flow also misses no-slip boundary layers and viscous drag.

For a Newtonian fluid, shear stress is proportional to strain rate. No-slip at solid walls produces velocity gradients and dissipation. Fully developed laminar pipe flow, for example, hasThe Reynolds number $\mathrm{Re}=UL/\nu$ compares inertial and viscous terms. Small values favour creeping-flow approximations; large values can permit thin boundary layers and turbulence, without guaranteeing turbulence in every geometry. The appropriate flow model follows from conservation laws, constitutive properties, boundaries and scale comparisons.

Viscosity transfers momentum between adjacent fluid layers. A Newtonian fluid has shear stress proportional to strain rate: for simple shear, $\boxed{{\tau }_{xy}=\eta \partial u_x/\partial y}$. Thus parallel plates of area $A$, separation $d$ and relative speed $U$ require force magnitude $\eta AU/d$ in a steady linear profile.

For an incompressible isotropic Newtonian fluid,for constant viscosity. With strain-rate tensor $e_{ij}$, the viscous dissipation rate per unit volume is $2\eta e_{ij}{e}_{ij}\ge 0$. Kinematic viscosity $\nu =\eta /\rho$ is the momentum diffusivity: a disturbance spreads over distance of order $\sqrt{\nu t}$.

The Reynolds number $\rho UL/\eta$ compares inertial and viscous forces. Small Reynolds number supports a creeping-flow approximation, while large values permit thin viscous boundary layers and possible turbulence. It does not mean molecular viscosity has disappeared. Gas viscosity results from molecular momentum transport; a liquid's microscopic transport depends strongly on intermolecular interactions. Non-Newtonian materials require a more general constitutive law.

Use positive normal stress for tension. On a small rectangular material element, ${\tau }_i$ is the normal force per unit area on the face perpendicular to axis $i$, and $e_i=\Delta {\ell }_i/{\ell }_i$ is the corresponding normal strain. Young's modulus $E$ is the ratio of axial normal stress to axial strain in uniaxial stress; Poisson's ratio $\sigma$ is minus the transverse-to-axial strain ratio in the same test. Thus isotropic linear elasticity gives $E{e}_1={\tau }_1-\sigma {\tau }_2-\sigma {\tau }_3$ and cyclic permutations. The stress symbols here follow the question; $\sigma$ denotes Poisson's ratio, not a stress component.

For $e_2=e_3=0$, the two transverse compliance equations imply ${\tau }_2={\tau }_3=\sigma {\tau }_1/(1-\sigma )$. Substitution into the first givesThis uniaxial strain modulus is greater than the uniaxial-stress Young's modulus for the positive Poisson's ratio of interest, because the lateral constraints prevent free Poisson contraction.

For a small twist $\varphi$, the bonded faces' relative tangential displacement at radius $r$ is $r\varphi$. The layer's shear strain is $r\varphi /h$, so the shear stress is $Gr\varphi /h$. Integrating the torque givesThe minus sign denotes a restoring torque; the question quotes its magnitude. This uses the disk's polar second moment of area.

For tilt about H, choose $y$ in the plate perpendicular to H. The vertical displacement of the bonded top face is $y\theta$ to first order, giving normal strain $y\theta /h$. The thin layer is constrained laterally by bonding, so its leading normal traction is $My\theta /h$, not $Ey\theta /h$. Integrating its moment about H givessince the diametric second moment of area is $\pi a^4/4$. This is the specified thin-layer approximation to the torsion and tilt stiffness of a bonded elastic disk. Relaxation near the free cylindrical edge is neglected; for a nearly incompressible layer it can limit quantitative accuracy unless the thinness is sufficient.

The thin square plate's moment of inertia about V is $I_V=2m{b}^2/3$. Its two perpendicular in-plane central moments are equal and sum to $I_V$, so the diagonal H has $I_H=m{b}^2/3$. Dividing the two angular stiffnesses by the appropriate moments of inertia givesThe geometry factors cancel in their ratio. Using the supplied relation between shear modulus and Young's modulus,The result lies just below the incompressible isotropic limit $1/2$, as expected for a rubber with a large compression stiffness compared with its shear modulus. It is an inference within the bonded-layer model, rather than an independently measured material constant.

For a Newtonian fluid with velocity $u(z)\hat{\mathbf{x}}$, $\tau =\eta {\partial }_{z}u$ is the $x$ component of traction on a plane whose normal points along $+z$. The opposite face has the opposite traction. The coefficient $\eta$ is dynamic viscosity, with units $\mathrm{Pa}\mathrm{s}$. Faster layers transfer positive $x$ momentum toward slower neighbours; the internal shear stress resists relative sliding.

Take a thin fluid element with horizontal area $A$ and thickness $dz$. The difference between its upper and lower viscous tractions is $A\eta u_{zz}dz$. The pressure gradient supplies $-A{p}_{x}dz$. Incompressibility is satisfied because $u$ has no $x$ dependence, and the convective acceleration vanishes: $(\mathbf{u}\cdot \nabla )\mathbf{u}=u{\partial }_x\mathbf{u}=0$. Equating net force to mass times acceleration therefore gives the relevant component of the Navier-Stokes equation:For steady flow, integrate $\eta u^{\prime \prime }=p^{\prime }$ twice and impose the no-slip boundary condition at each plate, $u(0)=0$ and $u(a)=u_0$. The Couette-Poiseuille flow profile isThe linear term is Couette flow and the parabolic term is plane Poiseuille flow. A negative $p^{\prime }$ drives flow in $+x$. Integrating the profile across the gap gives the volumetric flow per unit width,It has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$; multiplying by the plate width gives the actual volume per time.

With zero pressure gradient, substituting $u=\mathrm{Re}[A{e}^{i(kz-\omega t)}]$ gives $-i\rho \omega =-\eta k^2$. Select the root that decays into fluid occupying $z>0$:The other sign applies to decay in the opposite half-space. With a real wall velocity amplitude $U_0$, the established periodic solution of Stokes' second problem isThe oscillation's amplitude decays exponentially and its phase lags with distance from the plate. This is diffusion of tangential momentum, not a freely propagating undamped elastic wave. Initial transients must have decayed before this purely periodic expression applies.

The fluid-displacement amplitude and the wall-displacement amplitude have the same ratio as their velocity amplitudes, since both are divided by the common $\omega$. Thus for the specified amplitudes,Using the standard approximate water density $\rho \simeq 10^3\mathrm{kg}{\mathrm{m}}^{-3}$, explicitly an approximation rather than another value printed on the paper, gives $\delta =0.564\mathrm{mm}$ and $z\simeq 2.21\mathrm{mm}$ from either face of a sufficiently large plate. The viscous penetration depth varies as ${\rho }^{-1/2}$ if a more precise water density is used. The ideal parallel-flow solution ignores finite plate edges and tank walls.

In the thin-gap Couette flow approximation, the tangential velocity varies linearly through the gap. At radius $r$, its difference between the two discs is $r\omega$, so a Newtonian fluid exerts shear stressAn annulus has area $2\pi rdr$, and its tangential force acts at lever arm $r$. Integrating its torque givesThe driving torque is in the direction of rotation; the liquid's resisting torque has the opposite sign. This parallel-disc viscometer result assumes laminar shear with no slip, negligible edge effects and negligible secondary circulation. Its fourth-power radius dependence comes from the combined radial dependence of speed, annular area and lever arm.

A Newtonian fluid exerts viscous shear stress proportional to the rate of strain. For an incompressible fluid of constant dynamic viscosity $\eta$,The local momentum balance is the Navier-Stokes equation:The left side is density times material acceleration. A steady velocity field can still have convective acceleration when a moving parcel changes speed or direction; setting $\partial \mathbf{u}/\partial t=0$ does not eliminate it.

The Reynolds number $\mathrm{Re}=\rho UL/\eta$ compares inertia with viscous effects. At low Reynolds number, the fluid's own inertia can often be neglected, giving Stokes flow. This does not say that every immersed body has zero acceleration: its mass and the net pressure, viscous and body forces still enter its own equation of motion.

For a small rigid sphere moving through a slowly varying viscous background, Stokes drag is $-6\pi \eta a(\mathbf{v}-\mathbf{u})$. In still fluid, a simple settling model includes weight, buoyancy and drag:Ignoring fluid added mass and unsteady history forces, this gives relaxation time $m/(6\pi \eta a)$ and terminal velocityAt the terminal velocity the net force and acceleration vanish, while viscous drag and gravity remain individually nonzero. The dissipated power $6\pi \eta a{v}_t^2$ is supplied by the loss of gravitational potential energy relative to buoyant support. The creeping-flow and quasisteady-drag assumptions must be checked before applying this expression at large speed or immediately after a rapid start.

Elasticity describes recoverable deformation. For small displacements $\mathbf{u}(\mathbf{x})$, the infinitesimal strain tensor is the symmetric part of the displacement gradient,Its diagonal entries describe extension, while its off-diagonal entries describe shear. The engineering shear strain is ${\gamma }_{xy}=2{ϵ}_{xy}$; the antisymmetric part of the displacement gradient is a small rigid rotation, not elastic strain.

The stress tensor specifies the force per unit area on each oriented plane: $\mathbf{t}(\hat{\mathbf{n}})=\mathbf{\sigma }\hat{\mathbf{n}}$. Positive normal stress is tensile in the convention used here; tangential traction gives shear stress. In the absence of internal couples, angular-momentum balance makes the stress tensor symmetric. Force balance gives $\nabla \cdot \mathbf{\sigma }+\mathbf{f}=0$ in static equilibrium, where $\mathbf{f}$ is body force per unit volume.

Linear elasticity relates stress and strain. A homogeneous isotropic solid needs two independent elastic constants:Here $G$ is the shear modulus and $\lambda$ is the other Lamé parameter. For uniaxial loading, Young's modulus is $E={\sigma }_{xx}/{ϵ}_{xx}$, with lateral contraction described by Poisson's ratio $\nu =-{ϵ}_{yy}/{ϵ}_{xx}$. The bulk modulus is $K=-p/(\Delta V/V)$ for a small hydrostatic compression. The relations between these constants areStable isotropic elastic strain energy density requires $G>0$ and $K>0$, implying $-1<\nu <1/2$. The incompressible limit is approached as $\nu \to 1/2$. Anisotropic crystals generally need more constants; the isotropic relations should not be imposed on arbitrary directions in such a crystal.

The elastic energy density is $w=\mathbf{\sigma }:\mathbf{ϵ}/2$. In the isotropic linear model this becomes $\lambda (\mathrm{tr}\mathbf{ϵ}{)}^2/2+G\mathbf{ϵ}:\mathbf{ϵ}$. The quadratic form explains why reversible loading stores energy and unloading can return it, unlike plastic deformation or a dissipative hysteresis cycle.

For a slender beam described by the Euler-Bernoulli beam equation, plane sections approximately remain plane. A fibre a distance $y$ from the neutral axis has axial strain $-y\kappa$, where $\kappa$ is curvature, and stress $-Ey\kappa$. Integrating its moment about the neutral axis gives, with consistent signs,$I_A$ is a second moment of area, with units of length to the fourth power; it is not a mass moment of inertia. For a rectangular section of width $b$ and thickness $h$, $I_A=b{h}^3/12$, explaining the strong dependence of bending stiffness on thickness. For small transverse deflection $w(x)$, curvature is approximately $w^{\prime \prime }$, and the Euler-Bernoulli beam equation is $E{I}_A{w}^{\prime \prime \prime \prime }=q(x)$ with a suitable load-sign convention. A cantilever of length $\ell$ under end load $F$ has tip deflection $F{\ell }^3/(3E{I}_A)$ and largest bending stress of magnitude $F\ell c/I_A$, where $c$ is the furthest distance from the neutral axis.

In torsion of a solid circular shaft, a cross-section rotated by angle $\theta (x)$ has shear strain $r{\theta }^{\prime }(x)$. ThereforeHere $J$ is the polar second moment of area. A uniform shaft of length $\ell$ under constant torque twists by $\mathcal{T}\ell /(GJ)$, and its largest shear stress is $\mathcal{T}a/J$. Noncircular sections can warp, so replacing their torsional stiffness by $G$ times their polar area moment is not generally valid. Both bending and torsion formulas presume small, elastic deformations; exceeding the elastic limit can produce permanent deformation or fracture.

Compression and stretch viewed on rotated faces. To isolate shear, take equal small principal strains ${ϵ}_{xx}=-e$, ${ϵ}_{yy}=+e$. A square with sides at $45^{\circ }$ to these directions becomes a rhombus: adjacent side lengths are unchanged to first order, but their right angle changes by $2e$. In axes rotated through $45^{\circ }$, the same strain tensor has zero diagonal entries and off-diagonal entries of magnitude $e$.

The pure shear from opposite principal strains therefore has engineering shear strain $\gamma =2e$, not $e$. Equal opposite principal stresses $(-S,+S)$ give purely tangential traction of magnitude $S$ on those rotated faces. HenceIf the compression and stretch are not equal and opposite, an additional dilatational strain remains; their deviatoric part still supplies shear.

The shaft's torsional stiffness. At radius $r$, an end twist $\varphi$ over length $L$ gives tangential displacement $r\varphi$ and shear angle $r\varphi /L$. The local shear stress is $Gr\varphi /L$. An annular area $dA=2\pi rdr$ contributes torque $r\tau dA$, soThus torsion of a solid circular shaft givesHere $J_p$ is the polar second moment of area, in ${\mathrm{m}}^4$, not a mass moment of inertia. This $C$ is the transmitted torque magnitude; the shaft's restoring couple opposes the twist. The formula assumes small elastic shear and a circular uniform section.

Earthquake location. The hypocentre is rupture initiation; its surface projection is the epicentre. Record calibrated arrival times at stations of known coordinates, identify P wave and S wave phases, and predict their travel times through a stated velocity model. Earthquake location by arrival times fitsfor the three source coordinates and origin time $t_0$. Linearization supplies residual equations $\delta t_i=\delta t_0+\nabla T_i\cdot \delta \mathbf{x}$. A grid search or iterative weighted least squares inversion can use pick uncertainties, while station geometry controls whether those parameters are well constrained. Four independent arrival constraints are the minimum for four unknowns, not a guarantee of a reliable solution.

In a homogeneous medium, an arrival-time difference between an S wave and a P wave cancels origin time and estimates source-to-station distance,Several distances intersect to locate the event; real layered or three-dimensional Earth structure requires curved-ray travel times. Nearby stations and diverse takeoff angles help resolve depth. Clock errors, phase misidentification, poor azimuthal coverage and incorrect velocities introduce uncertainty or systematic bias. Closely spaced earthquakes can be relocated precisely relative to each other because common-path errors partly cancel, even when their absolute locations remain uncertain. Report a location uncertainty and model assumptions, not merely a best-fit point.

H. F. Reid's earthquake-cycle model. A fault can remain locked while long-term tectonic motion slowly deforms the surrounding rocks. Elastic strain energy and shear stress accumulate. When resistance is overcome, rapid fault slip releases part of that stored energy as seismic waves, frictional heat and other irreversible work. The adjoining rocks partly recover their elastic distortion while retaining a permanent displacement discontinuity across the fault: this is elastic rebound. Surface markers initially crossing the fault become bent during loading and offset after slip.^[ref]

A simple sketch is a slow rise of stress followed by a rapid stress drop, then renewed loading; the plate displacement continues throughout. Geodetic deformation, observed fault offsets and earthquake stress changes can test the picture. Actual cycles need not repeat with identical slips or times: creep, afterslip, stress transfer and heterogeneous friction complicate the idealization. The release need not remove all pre-existing stress. Elastic rebound explains storage and release; it is not a deterministic earthquake clock.

Slab pull is generally the largest direct driving contribution for plates attached to substantial cold descending slabs. Cooling makes oceanic lithosphere negatively buoyant; sinking slabs transmit a pulling force to the rest of the tectonic plate. The comparison is a force balance, so the following driving and resisting processes both matter.

The first two gravitational contributions ultimately belong to the same coupled mantle convection system: the plates are its cold upper thermal boundary layer. They are not passive rafts on a prescribed conveyor belt. A plate with little attached descending slab need not have slab pull as its largest local driver, and a large gross slab pull can be partly consumed by bending and interface resistance. Geological observations and force comparisons nonetheless support the general dominance of slab pull for rapidly moving subducting plates.^[ref]

Use a rock strength envelope: at each depth plot the smaller differential stress needed for frictional failure and for sustained ductile material creep at a specified strain rate. Shallow rock generally strengthens with depth because effective confining pressure increases; sufficiently deep hot rock weakens because thermally activated material creep becomes efficient.

At shallow depth, Byerlee's law or a Mohr-Coulomb failure criterion gives a relation of the formwhere $C$ is cohesion, $\mu$ a coefficient of friction, ${\sigma }_n$ the compressive normal stress and $P_f$ the pore pressure. Increasing overburden generally increases effective stress, hence the frictional strength. Conversion from fault shear stress to differential stress depends on the stress regime and fault orientation; a sketch should not label a universal numerical slope. High pore pressure or weak fault gouge can lower the frictional strength substantially.

A simplified dislocation creep relation isHere $Q$ is activation energy, $V^{\ast }$ activation volume, $R$ the molar gas constant and $T$ absolute temperature; $A$ and $n$ depend on composition and deformation mechanism. At fixed strain rate, increasing temperature usually decreases the required stress strongly. The competing effects of pressure, water weakening of rocks, grain size and composition also matter. The intersection with frictional failure gives a brittle-ductile transition, not a change from solid to liquid.

For the Indian continental interior, a useful first sketch has an upper-crustal strength peak, a weaker lower continental crust, and renewed strength in lithospheric mantle beneath the Moho. Quartz-bearing rocks can undergo ductile material creep at conditions where colder, drier olivine-rich mantle retains appreciable strength. This creates a possible two-peak, or “jelly sandwich”, rock strength envelope. With further depth the mantle also weakens. A cold craton can retain strong lithospheric mantle to much greater depths than a tectonically heated continental region.

For the central Pacific, show thin mafic oceanic crust above a relatively thick layer of strong lithospheric mantle. There need not be a conspicuous weak lower-crustal layer. Much of the integrated strength can reside below the Moho, with a broad strength peak followed by ductile weakening towards the asthenosphere. Cooling and age control the thickness: young hot oceanic lithosphere would give a much shallower and weaker profile than an old central-Pacific plate.

The plotted Moho depths of about 38 km and 7 km, and all plotted stresses, are illustrative assumptions. “India” alone does not specify a unique geothermal gradient, water content or strain rate, and the Pacific includes plates of different ages. The robust comparison is the possible crust–mantle layering in a continent versus the thin oceanic crust and mantle-dominated strength of an old oceanic plate. The Moho is a compositional boundary; neither it nor the deepest small plotted strength is automatically the seismic lithosphere-asthenosphere boundary.

Support for the left curve comes from high laboratory friction and substantial in situ differential stress; support for the right curve comes from low inferred fault traction, weak fault materials or high pore pressure. No single observation directly measures the whole plotted envelope.

For the left, stronger shallow-crust model:

For the right, weak shallow-layer model:

To discriminate, measure the strength-controlling quantities through depth and across both faults and their surroundings. Obtain borehole stress, pore pressure and temperature profiles where feasible; use borehole breakouts and hydraulic tests to constrain effective stress. Recover fault and host rocks and test their friction and material creep at relevant fluid conditions. The right curve specifically predicts an end to the shallow weakening near its sharp rise: investigate changes in permeability of a porous medium, alteration and fault gouge with depth, rather than testing only a surface sample.

Combine those measurements with heat flux transects, magnetotelluric sounding for conductive fluids, accurate earthquake depths and focal mechanisms, and GNSS or InSAR observations of deformation. The latter can help constrain the depth distribution of locking and postseismic relaxation, but neither locking nor seismicity alone gives absolute strength. Earthquake stress drop is a change in stress, not the pre-event stress. Low heat flux can also reflect heat transport or transient weakening; high electrical conductivity can have causes other than free water. A joint explanation of all these data is needed to choose between the models locally.