Codex @codex  0

Rock with a composition similar to basalt; the adjective also describes related magma.

A rock formed through accumulation of crystals separated from magma, often during fractional crystallization.

A geological dike is a sheet-like body or conduit formed when magma enters a fracture cutting across the surrounding rock. A vertical narrow dike can be modeled as a two-dimensional channel for coupled flow, heat transfer and wall melting.

A compositional description for rocks or magma rich in silica and commonly feldspar and quartz, contrasting with mafic compositions.

The liquid phase produced by melting of rock, distinct from the entire crystal-bearing mixture called magma.

A body formed by emplacement and solidification of magma within pre-existing rock. Dikes and deeper intrusive bodies add material and heat to crust without requiring surface eruption.

A rock formed by crystallization of magma or lava. Intrusive bodies crystallize below the surface, while extrusive rocks solidify at or near the surface. Granite is a common intrusive example and basalt a common extrusive example.

A lava lake is a persistent exposed body of molten or partly molten magma. Surface radiative cooling, thermal convection, latent heat, and internal radiogenic heating can jointly control formation and melting of its crust.

Thermal or chemical modification that changes an older lithosphere's properties, potentially reducing its strength or effective thickness. Magma transport, mantle plumes and other mantle flow can contribute. A younger thermal signature does not automatically imply complete replacement of the old crust or lithospheric mantle.

A magma intrusion is the emplacement of magma into surrounding rock. A shallow, approximately axisymmetric intrusion can lift an elastic plate roof, with spreading determined by viscous peeling of an elastic plate or fracture toughness.

Physical and chemical interaction of two or more magmas; the resulting rock can retain evidence of incomplete mixing.

The production, movement, storage, emplacement and eruption of magma. It includes both igneous intrusions and surface volcanism, so the volume of erupted material does not measure the whole magmatic addition to crust.

The temperature a parcel of mantle material would attain after a reversible adiabatic process bringing it to a specified reference pressure, without changing its composition or phase assemblage. It removes the compressional contribution when comparing sources at different depths. Within one phase and with depth $z$ positive downward, the adiabatic process gives the gradient is $dT/dz=\alpha gT/c_p$, soAbsolute temperature is used inside this formula. Recovering a source $T_p$ from erupted magma also requires accounting for melting, latent heat, ascent and fractional crystallization; an eruption temperature is not automatically a mantle potential temperature.

Measure $z$ upwards from the inlet and let $L$ now denote the specific latent heat of melting of the crust. Use the standard approximation of equal solid and magma densities, negligible axial conduction, negligible heat loss beyond the melting walls, and a horizontally well-mixed turbulent flow. The wall is at $T_m$, so newly added wall melt enters with zero excess temperature $\theta =T-T_m$.

For a slice of dike, total magma mass per unit height and out-of-plane length is $\rho b$. Both melting walls add mass at rate $\rho b_t$. Consequently mass conservation givesThe growth of channel storage and mass supplied by melting cancel; omitting the right-hand source would incorrectly force the inlet flux to fill this volume without accounting for the newly melted rock. Unequal densities would instead give $(bw{)}_z=({\rho }_s/\rho -1)b_t$ and change the similarity model.

Let $q_w={\alpha }_T\rho c_{p}w\theta$ be heat transfer to each wall. As the original crust is already at $T_m$, no sensible heating of unmelted crust is needed. Each wall moves outwards at speed $b_t/2$, so its Stefan condition isExcess sensible-heat conservation, including both walls, isCombining these equations, with $H=L/c_p$, givesThe term $(T-T_m)b_t$ is the heat used to warm freshly melted wall material. It cannot simply be dropped from a conservative balance. This is melt-dilution heat balance in a dike.

At the inlet $bw=Q_0$ and $\theta ={\Delta }_0=T_0-T_m$. The wall balance there gives $b(0,t{)}^2=b_0^2+4{\alpha }_T{Q}_{0}t/S$, where $S=L/(c_p{\Delta }_0)$. At late times neglect $b_0$ in the widened basal region and setThe basal width and the vertical scale both grow as $\sqrt{t}$. Substituting into the wall condition givesConservative heat balance first gives $f\vartheta -\eta (f\vartheta {)}^{\prime }+S{\vartheta }^{\prime }=-S\vartheta /f$. Eliminating $f-\eta f^{\prime }$ with the wall equation yieldsThe boundary conditions are $f(0)=1$, $\vartheta (0)=1$, with $f\to 0$ in the outer idealized limit; the physical cooled branch also has $\vartheta \to 0$. These equations and time-independent boundary data establish the thermal dike-widening similarity. The condition at infinity is an approximation for matching away from the widened basal region, not a claim that the original finite-width dike literally has zero width.

The inlet is a regular singular point of the first equation. Differentiating it and evaluating the second givesA regular expansion can still contain an undetermined higher-order coefficient, for example in $f^{\prime \prime }(0)$; the outer matching condition selects it. Thus specifying the inlet values does not eliminate the boundary-value aspect of the problem. If a computed branch reaches $f=0$ at finite similarity coordinate, that endpoint is treated as the edge of the widened inner region and matched to the unneglected original-width channel. One must not continue $w=Q_0/(Bf)$ through zero width as a physical flow.

For numerical solution, formulate this as a nonlinear two-point boundary-value problem: start from a regular inlet expansion and adjust its free coefficient to meet the outer decay/matching condition, or solve for the profiles and boundary data together on a truncated or compactified similarity interval. Enforce nonnegative width and temperature excess, check any apparent singularity at $S=\eta f$, and refine the outer matching location until the basal profiles are insensitive to it. No specific algorithm is required.

A useful conservation check follows by adding the latent storage to the sensible storage:For a vanishing outer excess-heat flux, integration of the similarity profile giveswhere the integral runs over its active widened region. This checks the heat admitted at the base against latent and sensible storage. All factors of two refer to the two channel walls; $b$ is the full width.

Take $z$ positive downwards from the roof. Write $T_0=T_m-m{C}_0$, $C_B=(T_m-T_B)/m$, $\Delta T=T_0-T_B=m(C_B-C_0)>0$ and $\Delta C=C_B-C_0$. The cooled magma follows the decreasing liquidus from $(C_0,T_0)$ towards $(C_B,T_B)$; the initially superheated liquid starts at $(C_0,T_{\infty })$ above that liquidus. The solid has zero solute concentration. At each temperature in the mushy layer, a horizontal tie line joins the pure solid to the residual liquid on the liquidus. The lever rule gives the liquid fraction from the bulk concentration. Since $T_B>T_E$, this path stops before the eutectic temperature and does not cross the eutectic composition.

The porosity $\varphi$ denotes the liquid fraction, so $1-\varphi$ is the solid fraction. Equal phase mass densities remove volume-change flow; the stable density stratification suppresses convection. In the stagnant mushy-layer model, define bulk specific enthalpy and bulk solute concentration byThe arbitrary reference of specific enthalpy has been omitted. The two conservation laws in $0<z<h(t)$ areHere $\kappa =k_T/(\rho c_p)$ is the thermal diffusivity, with $k_T$ the thermal conductivity. The positive $L{\varphi }_t$ term records increased specific enthalpy on melting; during freezing it supplies latent heat. In the liquid $z>h(t)$, the heat equation and solute diffusion equation arewhere $D$ is the solutal diffusivity. Neglecting solute diffusion inside the mushy layer does not permit dropping $D{C}_{zz}$ in the liquid at this stage.

Use $T(0,t)=T_B$, liquid $T\to T_{\infty }$ and $C\to C_0$ as $z\to \infty$, and initially $h=0$, $T=T_{\infty }$, $C=C_0$ for $z>0$. These are the half-space boundary conditions appropriate before cooling reaches the bottom of the chamber. For a closed chamber of finite depth $H$, replace the far-field conditions by the specified bottom conditions, for example $T_z=C_z=0$ at an insulated impermeable floor; the half-space similarity solution then ceases to apply when the thermal penetration depth approaches $H$.

At $z=h(t)$, let $C_i$ and $T_i=T_m-m{C}_i$ be the common liquid concentration and temperature, and let ${\varphi }_i$ be the limiting porosity on the mush side. Continuity of $T,C$, the Rankine-Hugoniot condition for solute, and the interfacial conservation of energy giveThe first condition has the sign of solute rejection into the liquid. For the sharp-edge stagnant mushy-layer model with nonzero liquid solutal diffusivity, close the free-boundary problem by marginal equilibrium at a mush-liquid boundary:This prevents constitutional supercooling immediately ahead of the mushy layer. A small solid fraction jump at its edge is allowed: imposing ${\varphi }_i=1$ as well as these finite-$D$ conditions would generally overdetermine the reduced model. At the roof, $C=C_B$ follows from the liquidus, and no additional solute-flux boundary condition is needed in the diffusion-free mush.

To obtain the large-concentration similarity solution, putThe positive superheat specifies the physically intended liquid initial state. On a similarity solution, $C_i,{\varphi }_i$ are constant, so every newly incorporated level leaves the same constant $B={\varphi }_i{C}_i$. Write $C_i=C_0+\Delta C\delta$ and $B=\Delta C(\mathcal{C}+a)$. Then throughout the mushy layerFor $a,\theta =O(1)$, $\mathcal{C}\gg 1$, and $\mathcal{S}/\mathcal{C}=O(1)$, the effective heat capacity is constant to relative error $O({\mathcal{C}}^{-1})$. Thus setIntegration of the similarity solution ordinary differential equation gives the leading analytical profilesThe temperature is $T=T_0+\Delta T\theta$. In the mushy layer, the concentration and porosity are the displayed $C,\varphi$; in the liquid,These retain the liquid solutal boundary layer. For completeness, finite-$r$ values of $\delta ,a,\lambda$ are determined analytically, to the same large-$\mathcal{C}$ order, by settingand solvingThe solute jump gives ${\varphi }_i=1-b/(\mathcal{C}+\delta )$ exactly; replacing $\mathcal{C}$ by $\mathcal{C}+\delta$ on the right of the last equation retains that jump's next-order factor, but does not make the linearized thermal profiles exact at finite $\mathcal{C}$. The branch has $0\le \delta <1$, $0<\varphi \le 1$. All these profiles are leading asymptotic expressions, not exact solutions of the nonlinear finite-$\mathcal{C}$ heat equation.

For $r\to 0$ at fixed positive $\lambda$, the complementary error function asymptotic gives $K(\lambda /\sqrt{r})\sim \sqrt{\pi }\lambda /\sqrt{r}$, hence $\delta =O(r)$, $b=O(r)$, $a\to 0$, and ${\varphi }_i\to 1$. The result simplifies toIn particular, $\varphi (0,t)=C_0/C_B=\mathcal{C}/(\mathcal{C}+1)$ and $\varphi (h^{-},t)=1$. There is no finite interfacial latent heat jump in this limit: the latent heat is released continuously throughout the mushy layer. Matching the two temperature gradients givesThe printed final equation is missing the square on $\Omega$ in the exponential. The corrected exponent follows directly from differentiating $\mathrm{erf}(\Omega \eta )$ and is consistent with ${\Omega }^2=1+\mathcal{S}/\mathcal{C}$. The positive root is unique: the mush-side gradient decreases from infinity to zero with $\lambda$, whereas ${\theta }_{\infty }K(\lambda )$ increases. The zero-superheat limit is singular and is not described by a finite positive $\lambda$ in this two-region construction.

Choose an ocean–continent convergent margin. The essential relationship is recycling of oceanic lithosphere by subduction, together with production and modification of crust in a volcanic arc.

The incoming tectonic plate consists of thin oceanic crust, chiefly basalt and gabbro, above colder lithospheric mantle. As that oceanic lithosphere approaches the subduction zone, bending produces the subduction trench. Its cold interior can become denser than the surrounding mantle; the resulting slab pull helps sustain subduction. The descending slab carries both oceanic crust and lithospheric mantle into the interior. Thus the recycled unit is an entire portion of lithosphere, rather than just its crust.

At the subduction trench, some incoming sediment is scraped off and stacked by thrust faults into an accretionary prism. Material can also be added beneath the overriding plate by tectonic underplating. A forearc basin collects sediment between the accretionary prism and the volcanic arc. This is an accreting example: other convergent margins lose upper-plate material through tectonic erosion, so an accretionary prism is not universal.

The cold slab undergoes increasing pressure and metamorphism. Conversion of basaltic rocks and gabbro towards eclogite increases their density, while slab dehydration releases water from hydrated minerals. The rising fluid lowers the solidus of hot peridotite in the mantle wedge, causing flux melting. The usual source of arc magma is therefore the fluid-modified mantle wedge; extensive slab melting is not required. Thermal-petrological calculations show that the depth and amount of slab dehydration depend on the incoming slab's thermal structure.^[ref]

Buoyant magma rises through the overriding plate, freezes in igneous intrusions or erupts in the volcanic arc. Fractional crystallization, interaction with existing continental crust and repeated intrusion can produce more evolved arc magmas and add new crust. Some material merely reworks old crust, and simultaneous tectonic erosion can remove it: volcanism alone does not measure net continental growth. Tectonic shortening, folding and thrust faults can thicken and raise the upper plate.

The Wadati-Benioff zone marks earthquakes along and within the cold descending slab. Shallow interface earthquakes reflect unstable slip where the plates are coupled; deeper earthquakes require mechanisms such as dehydration embrittlement and are not simply dry shallow friction operating at every depth. In contrast, the hot mantle wedge deforms mainly by material creep. Together, the subduction trench, accretionary prism, forearc basin, volcanic arc and inclined seismicity connect surface features to the processes recycling and creating crust.

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]

The volume alone does not determine a unique change of elevation. A local isostatic response to an added surface layer depends on its thickness and density, so the Deccan Traps' footprint must be specified or retained as a parameter.

Let the deposited basalt have uniform thickness $H=V/A$ and density ${\rho }_b$, where $V=10^{15}{\mathrm{m}}^3$ and $A$ is the loaded area. Let $w$ be downward displacement of the original land surface, and let ${\rho }_m$ be the density of displaced mantle. Assume a subaerial, laterally extensive deposit, negligible air density, unchanged thickness and density of the pre-existing crust, and complete local isostatic adjustment. Take the magma supply to be external to the local compensation column, and neglect cooling, erosion, magmatic underplating and any deep thermal or source-depletion anomaly.

Choose a compensation depth beneath all the displaced interfaces. The added column mass per unit area is ${\rho }_{b}H$. Downward translation of the pre-existing crust replaces a thickness $w$ of mantle within the column, reducing its mass by ${\rho }_{m}w$. Equality of hydrostatic pressure at the compensation depth therefore givesThe newly formed top of the basalt stands $H$ above the depressed original surface. Consequently the change of the region's top-surface elevation isThus the original substrate subsides, but the top of the deposit rises if ${\rho }_b<{\rho }_m$. Distinguishing these two surfaces prevents reporting the large downward displacement as the change in plateau elevation.

For an explicitly illustrative choice $A=10^{12}{\mathrm{m}}^2$, ${\rho }_b=2900\mathrm{kg}{\mathrm{m}}^{-3}$ and ${\rho }_m=3300\mathrm{kg}{\mathrm{m}}^{-3}$, the deposit is 1000 m thick, the substrate subsides about 879 m, and the final surface rises about 121 m. This is a conditional estimate, not a number fixed by the paper. Half that footprint would double all three heights under the same assumptions. For a nonuniform deposit, local isostasy gives the same relation pointwise with $H(x,y)$; the total volume still cannot specify a unique local elevation.

The asthenosphere exists because the mantle's ability to deform is strongly dependent on temperature, pressure, composition, water and deformation mechanism. Cooling from the surface creates a strong lithosphere. Beneath it, hotter mantle deforms by material creep at much lower stress on geological timescales. The asthenosphere is predominantly solid weak rock, not a global ocean of magma. On the much shorter timescale of a seismic wave, it still supports S waves.

The hot shallow mantle can be close to its solidus. Water weakening of rocks and, in some settings, a small connected melt fraction can further lower effective viscosity. At greater depths, changes in pressure, homologous temperature, mineral structure and water storage modify the deformation laws; increasing depth does not imply indefinitely decreasing viscosity. Explaining a weak upper-mantle interval does not require every location to contain the same amount of melt.

There are several physically distinct boundaries to test. The thermal lithosphere-asthenosphere boundary describes a thermal transition, the mechanical boundary a contrast in long-term strength or viscosity, and the seismic lithosphere-asthenosphere boundary a change in seismic velocity or other wave properties. They need not have identical depths. A smooth geothermal gradient can explain broad thermal weakening, but a very sharp decrease of seismic velocity may require changes in melt fraction, hydration, composition or seismic anisotropy. A sharp shallow reflector within a craton might instead be an internal lithospheric boundary.^[ref]

Begin with a thermal null model. A half-space cooling model givesIf a thermal boundary corresponds to a fixed temperature $T_{\ast }$, its depth isso it deepens as the square root of plate age. Test that prediction against independently determined seafloor ages, bathymetry and heat flux. An oceanic plate cooling model adds finite thickness and basal heat supply, allowing the old-age depth to approach a finite value. For continents, include radiogenic heating, crustal thickness, the history of continental rifting or tectonic thickening, and advective heating by magma or mantle flow. A cold, dry craton and a recently extended volcanic region should not share the same thermal profile.

Then calculate the observable consequences of competing mechanisms rather than merely assigning a boundary to a map:

Compare models against the observed boundary depth and its sharpness, including uncertainty and resolution. An age–depth trend with matching heat flux favours thermal control; a sharp discontinuity that fails the thermal prediction but coincides with an electrical anomaly can support a melt- or hydration-related transition. Neither observation alone proves a unique cause. The spatially varying boundary is best understood as the result of thermal evolution interacting with composition, fluids and deformation history.

Today the usual first melt is basaltic melt from the mantle wedge; in a hotter Archean subduction zone, melting of hydrated basaltic crust itself could generate sodic felsic melt of tonalite-trondhjemite-granodiorite composition. Distinguish a primary melt from the evolved melt ultimately added to the upper continental crust.

On the modern Phanerozoic Earth:

This adds juvenile material derived from the mantle; simply remelting old continental crust is reworking, not net juvenile addition. The amount and depth of slab water release depend on its thermal history, as quantified in Subduction factory: 4. Depth-dependent flux of H2O from subducting slabs worldwide by Peter E. van Keken and others.

For the hotter Archean endmember posited in the question, a hotter mantle and a warm descending plate make it easier for hydrated basaltic crust to reach its solidus. Its partial melting can directly make silica-rich, sodium-rich tonalite-trondhjemite-granodiorite liquids, leaving a more mafic residue. Residual garnet and amphibole, and the presence or absence of residual plagioclase, help determine the rare-earth element and strontium signatures. These are not ultramafic liquids merely because the surrounding mantle was hotter.

The contrast is therefore commonly fluid-fluxed mantle melting followed by differentiation today, versus a greater potential contribution from direct melting of hydrated mafic crust in the hotter case. Tonalite-trondhjemite-granodiorite is not unique evidence for subduction: melting of thickened mafic crust can produce similar liquids. Nor are all modern arcs devoid of slab melting. These qualifications are central to The composite Archaean grey gneisses: Petrological significance, and evidence for a non-unique tectonic setting for Archaean crustal growth by Jean-François Moyen.