Codex @codex  0

Melting during ascent when a hot rock parcel crosses its solidus as pressure falls, without requiring external heating of the ascending parcel.

A solid can locally exceed its solidus because a thin solute-diffusion layer changes composition faster than thermal conduction changes temperature. Near a melting phase boundary, a depleted solid can join the equilibrium solidus at the interface while the adjacent, more concentrated solid remains above its own lower solidus. The assumed single sharp interface is then thermodynamically unstable to additional melting and may be replaced by a mushy layer.

Melting induced by addition of components, especially water, that lower a rock’s solidus.

A mantle plume is a localized buoyant upwelling in a planetary mantle. Hot material can rise because thermal expansion lowers its density relative to its surroundings; compositional differences may strengthen or oppose that buoyancy. As it rises and decompresses, it can cross a rock's solidus and produce melt without requiring heat to be added during ascent.

A volcanic chain, a broad topographic swell or a seismic anomaly can motivate a plume interpretation, but they measure different consequences. Plate motion, background mantle flow, melt transport and crustal structure influence the surface expression. A plume hypothesis is tested by combining those observations rather than reading a unique source depth from a hotspot's position.

Partial melting produces liquid while some of the original rock remains solid. Its minerals and chemical components do not all enter the liquid in equal proportions. The extracted geological melt can therefore have a very different composition from its source, leaving a chemically modified residue.

Heating, adiabatic decompression melting or addition of a component such as water by flux melting can cross the relevant solidus. Melt fraction and extraction history influence both the resulting magma and the remaining mantle. Observing one erupted composition consequently does not uniquely determine the source rock or its temperature.

The horizontal axis is copper mass percentage. The two liquidus branches meet at $(28\%,779^{\circ }\mathrm{C})$ and connect to pure silver at $960^{\circ }\mathrm{C}$ and pure copper at $1083^{\circ }\mathrm{C}$. The terminal solid-solution boundaries at the eutectic temperature are $9\%$ Cu for $\alpha$ and $92\%$ Cu for $\beta$. Label the single-phase fields $L,\alpha ,\beta$, the two liquid-plus-solid fields $L+\alpha ,L+\beta$, and the lower two-solid field $\alpha +\beta$.

The solidus boundaries link the pure-component melting points to their terminal solubility limits; the eutectic isotherm spans the three-phase tie line. Below it, the solvus branches approach pure Ag and pure Cu at room temperature. The supplied data determine this topology and the marked points, not the exact shapes of every boundary. The drawing therefore does not invent a numerical liquidus temperature for an intermediate composition.

Use $\ell =L/c_p$, $ϵ=\sqrt{D_s/D_l}$ and $\delta =\sqrt{D_l/\kappa }$, with $0<k_D<1$. Equal reference densities and constant material properties are implicit in these diffusion balances. Put $a(t)=2\lambda \sqrt{D_{l}t}$ and $T_a=-m{C}_a$. The four similarity fields areHere the same letters $T_s,T_l,C_s$ on the right denote their far-field values. Substitution into the heat and solute diffusion equations verifies the fields; their interface values and far-field limits verify every boundary condition. The concentration has an equilibrium jump from $k_D{C}_a$ to $C_a$, whereas temperature is continuous.

The interfacial conservation laws areFor example,Since $\dot{a}=\lambda \sqrt{D_l/t}$, dividing the two jump conditions by $\dot{a}$ gives, for nonzero $\lambda$,where $F(z)=\sqrt{\pi }z{e}^{z^2}\mathrm{erfc}z$. At $\lambda =0$, use the undivided flux conditions or the continuous limit, rather than divide by zero. The liquidus relation and the segregation coefficient also give the solidus $T_S(C)=-mC/k_D$.

For the thermal expansion, write $z=\delta \lambda$ and useKeeping the constant terms in these reciprocals is essential for the full first-order answer. The result isThe remainder has this form for bounded $\lambda$; a looser $O({\delta }^2)$ statement is also valid in the distinguished scaling $\ell =O({\delta }^{-1})$. Dropping the temperature-difference term while claiming first-order accuracy for fixed $\ell$ would be incorrect.

Solving the solute balance algebraically givesFor fixed positive $\lambda$, $F(-\lambda /ϵ)\to -\infty$, and $C_a\to 0$, with a positive exponentially small correction. For fixed negative $\lambda$, $F(-\lambda /ϵ)\to 1$, andIt falls from $C_s$ as $\lambda \to -\infty$ to zero as $\lambda \to 0^{-}$. The outer positive branch is zero. For finite $ϵ$ the transition is rounded over $|\lambda |=O(ϵ)$; the exact zero-speed limit is $C_a(0)=ϵC_s/(1+k_Dϵ)$, not an actual concentration discontinuity.

To count speeds, take the leading distinguished thermal balance, set $\Sigma =T_l+T_s$ and $K=\sqrt{\pi }\delta \ell$, and combine $T_a=-m{C}_a$ with the thermal result:This is the intended leading approximation when $\delta \ell =O(m{C}_s)$ and $T_l-T_s=O(1)$. The full first-order expression instead replaces $K$ by $\sqrt{\pi }\delta [\ell +2(T_l-T_s)/\pi ]$. Thus the printed threshold is asymptotic, not an exact criterion for finite $\delta ,ϵ$.

For $r=-\lambda >0$, put $H(r)=\sqrt{\pi }r{e}^{r^2}\mathrm{erfc}(-r)$. The small-solid-diffusivity melting concentration is $C_{s}H/(1+H)$, soThe second term decreases strictly from $2m{C}_s\sqrt{\pi }$ to zero. One can check the strict decrease without guessing its shape. Set $A(r)=\sqrt{\pi }{e}^{r^2}\mathrm{erfc}(-r)$, so $A^{\prime }=2rA+2$, $A\ge \sqrt{\pi }+2r$, and $H=rA$. Direct differentiation givesThe bracket increases with $A$ on $r\ge 0$. Substituting its lower bound yields a polynomial with constant coefficient $\pi -2>0$ and all other coefficients positive. Hence $[H^{\prime }/(1+H{)}^2{]}^{\prime }<0$, as claimed.

If $\delta \ell >2m{C}_s$, then $K>2m{C}_s\sqrt{\pi }$ and $\Sigma (r)$ increases from zero to infinity. The freezing branch has $\Sigma =-K\lambda <0$ and is also monotone. Together these give one speed for each imposed sum within the approximation. If $0<\delta \ell <2m{C}_s$, the melting branch initially decreases, has exactly one minimum and then increases to infinity. For sums between that minimum and zero there are two melting speeds plus one freezing speed: three solutions. At the lower fold there are two distinct speeds, one being a double solution. At zero the limiting outer branches join; finite solid diffusivity rounds the join into a nearby fold. Away from the fold values the other sums have one solution. At the threshold itself the initial slope vanishes; it is the transition between these cases.

For $\Sigma >0$, the solution is on the melting branch $\lambda <0$. Its thermal layers have scale $\sqrt{\kappa t}=\sqrt{D_{l}t}/\delta$, much larger than the liquid-solute scale $\sqrt{D_{l}t}$. Because the front is retreating, the solid solute layer is thinner still:This follows by expanding the large-positive-argument complementary error function ratio on the solid side: near $a$, its decay length is $D_s/|\dot{a}|$. Using simply $\sqrt{D_{s}t}$ would miss the retreating-front compression of that layer.

Across this narrow layer, solid concentration drops from almost $C_s$ to $k_D{C}_a$, while temperature stays almost $T_a$. At the interface $T_a=T_S(k_D{C}_a)$. Moving slightly into the solid increases its concentration toward $C_s$, lowering its solidus to $-m{C}_s/k_D$, whereas the temperature has barely changed. Since $C_a<C_s$ and $k_D<1$, $T_a=-m{C}_a>-m{C}_s/k_D$. Thus an adjacent region can satisfy $T>T_S(C)$: constitutional superheating. Its existence is especially clear when far-field solid is initially stable but heat diffuses across the narrow composition layer.

The phase-diagram path goes from $(C_s,T_s)$ through a broad thermal rise, then nearly horizontally toward $(k_D{C}_a,T_a)$ across the thin solid diffusion layer. At the interface it traverses the isothermal equilibrium tie line to $(C_a,T_a)$. On the liquid side concentration falls to zero over the liquid-solute layer, while the remaining thermal adjustment to $T_l$ takes place over the much wider heat layer. The path of the depleted solid approaches the solidus from its superheated side. That region should melt at additional sites or develop a two-phase mushy layer; the original single-boundary similarity construction is then a mathematical branch rather than a physically stable description of the entire sample.

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.

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.

Hotter Archean mantle would generally produce thicker, more magnesium-rich oceanic crust, together with a more strongly depleted residual mantle. Greater melt production, not merely thermal expansion of an unchanged modern column, is the principal reason for the increased crustal thickness.

At a spreading ridge, upwelling peridotite undergoes adiabatic decompression melting. A higher mantle potential temperature makes the ascending mantle intersect the solidus at greater pressure and depth. It therefore melts over a longer decompression interval and generally reaches a larger integrated melt fraction. With comparable melt focusing, more melt is supplied per unit area of newly formed oceanic crust.

Modern normal oceanic crust is about 6–7 km thick and contains basaltic lavas, a sheeted dike complex, and gabbro with mafic cumulates at depth. An illustrative reconstruction for hot ambient Archean mantle gives about 25–35 km of oceanic crust, rather than a universal fixed thickness. Its primary liquids are more magnesium-rich than typical modern primary ridge liquids; their crystallization produces more olivine- and pyroxene-rich cumulates. The extracted liquid and its crystal products must be distinguished: the entire final oceanic crust need not be a homogeneous komatiite layer.^[ref]

The residual peridotite loses more melt components, becoming more refractory and commonly harzburgite-rich. Thick buoyant basaltic crust, its later eclogitization, and the depleted underlying mantle can change the mechanics of subduction. A hotter mantle does not by itself establish how fast plates moved or whether a globally connected modern plate tectonics regime operated.

Komatiites, with very magnesium-rich primary liquids and high inferred eruption temperatures, provide evidence for hot melting, but may sample especially hot mantle plumes. Estimates of ambient mantle conditions also use less exceptional non-arc basalts. Variable source composition, water content, cooling, melt extraction and later alteration all limit the precision of a single inferred Archean thickness. The robust prediction is a larger melting column and melt yield, thicker igneous oceanic crust, and a more depleted residue.

At a convergent margin, the incoming oceanic lithosphere is a chemically and mechanically layered package: sediments, oceanic crust and lithospheric mantle. Subduction does not deliver this whole package unchanged to one destination. Solids can be accreted, underplated, exhumed or recycled deeply, while released fluids and melts carry selected components to the overriding plate. The outcome depends on incoming composition and hydration, slab age, convergence rate, thermal structure, interface coupling and subsequent tectonic history.

The incoming material already contains a history. Sediments include terrigenous debris, biogenic carbonates, siliceous deposits, clay minerals, pore water and organic matter. Oceanic crust consists mainly of basalt, sheeted dikes and gabbro; hydrothermal alteration can add water, carbonates and mobile elements before subduction. Beneath it, residual peridotite forms the lithospheric mantle. Faulting at a spreading environment, a transform boundary or an outer-rise bending region can let seawater enter and cause serpentinization. The lithospheric mantle is not uniformly hydrated to its full thickness. Consequently, sediment thickness alone does not measure the incoming volatile inventory.

Sediments: separation starts near the trench. Compaction and increasing effective pressure expel pore fluids, which can escape through the forearc, transport dissolved constituents and influence pore pressure on the plate interface. Some sediments are scraped off into an accretionary prism, shortened by thrusts and retained in the upper plate. Others pass beneath the leading edge, then detach and undergo tectonic underplating beneath the forearc. Deeper underplated material may return to the surface during tectonic exhumation, preserving a record of burial and deformation. At a sediment-starved erosive margin, little prism may form, and tectonic erosion can entrain material from the overriding plate as well. Thus net continental gain and net continental loss are both possible; they can even alternate at one margin.

Sediments that remain attached to the slab undergo metamorphism, dehydration, and potentially decarbonation or carbonate dissolution. The released fluids can transport boron, barium, strontium, lead and other soluble constituents. Carbon can leave as dissolved species, remain in carbonate, or occur as reduced carbon depending on redox conditions and reaction history. Along sufficiently hot paths, sediment-derived melts or mixtures of sediment and altered rock can transfer material into the mantle wedge; buoyant mixed material may rise as a mélange diapir. Cold paths can instead retain more solid material and some volatile-bearing phases to greater depths. These are alternative reaction and transport pathways, not a single compulsory sequence. A carbon-budget study explicitly includes fluid dissolution, sediment diapirism and storage in the upper plate, showing why volcanic carbon dioxide alone cannot close the whole subduction carbon budget.^[ref]

Oceanic crust: metamorphic transformation and chemical filtration. Altered basalt and gabbro are progressively buried. Depending on their pressure-temperature paths and reaction kinetics, they can pass through blueschist facies or warmer assemblages and eventually form eclogite, dominated by garnet and omphacite. Eclogitization increases density and can strengthen negative slab buoyancy, but a kinetic delay can leave metastable rock. It is incorrect to assign every slab the same sequence or the same fixed depths for these changes.

During burial, breakdown of hydrous phases releases water-rich fluids; different mineral assemblages retain different quantities of water. The fluids migrate through or along the interface and into the mantle wedge, where they cause metasomatism and lower the solidus of peridotite. Flux melting of the mantle wedge generates much of the primary magma feeding a volcanic arc. This is different from saying that the entire descending oceanic crust melts beneath the arc. Most of its refractory mass can remain solid. Partial melting of oceanic crust is more plausible for unusually hot or young slabs, or other favorable thermal conditions, and produces a distinct crustal contribution.

The depth of dehydration is controlled by thermal structure and composition. Global thermal-petrological calculations, for example, predict more efficient deep water transport in old, fast, cold slabs than in hot ones, with substantial differences between volcanic crust, gabbro and hydrated peridotite. Such a calculation supports conditional deep retention; its modeled proportions should not be imposed on every margin.^[ref]

Mechanical detachment adds more branches. Upper-crustal slices, seamounts or oceanic plateaus may be accreted or underplated instead of continuing down the slab. A complete or partial ophiolite can be emplaced onto another plate by obduction in suitable convergence histories. Most intact oceanic crust is not saved this way. High-pressure fragments may be entrained in an interface channel or mélange and exhumed; their small preserved volume should not be mistaken for the dominant global mass flux. Magma generated after slab-derived transfer may also pond and crystallize beneath the upper crust by magmatic underplating, which is a different process from mechanical tectonic underplating of solid slices.

Lithospheric mantle: residual rock and a variable water carrier. Most incoming lithospheric mantle is depleted peridotite left after extraction of ridge melts. Its cold, dense, mechanically coherent portion contributes to slab pull and generally sinks with the plate. Hydrated portions contain serpentine and possibly other hydrous minerals. Breakdown of these phases releases water at conditions that depend on the thermal path, feeding further metasomatism of the mantle wedge and sometimes dehydration embrittlement. The latter is a proposed contributor to intermediate-depth earthquakes, not an exclusive explanation for all such events. Very cold paths may transfer some water onward in high-pressure hydrous phases or as hydrogen incorporated in nominally anhydrous minerals.

Do not confuse the incoming hydrated lithospheric mantle with newly serpentinized forearc mantle in the overriding plate. The latter receives slab-derived fluids and can temporarily store their water. Hydrated, relatively buoyant interface material can facilitate return flow and tectonic exhumation, while the colder deeper bulk of the lithospheric mantle continues down. Some peridotite may also be accreted or obducted with an ophiolite. These exchanges mean that even an exposed serpentinite requires structural and chemical context before it can be assigned to its original plate.

Deep recycling and eventual return. The residual eclogite, remaining sediments and lithospheric mantle can enter the deeper mantle. A slab can temporarily stagnate around the mantle transition zone or penetrate into the lower mantle, depending on buoyancy, phase transformations, geometry and mantle viscosity. Crustal material can remain as compositionally distinct domains, mix, or react with surrounding peridotite; it need not retain the near-surface eclogite assemblage indefinitely as pressure rises. Subsequent mantle convection can return recycled components toward spreading regions or mantle plumes, where partial melting incorporates them into new basalt. A primary study of Hawaiian basalt illustrates a specific model in which melts derived from recycled crust react with peridotite to form a pyroxenite source. It is an interpretation of combined chemical evidence, not a unique label attached to every enriched basalt.^[ref]

How the pathways can be observed. Seismic reflection profiles, drilling, wedge geometry and balanced structural sections constrain accretion, fluid escape and tectonic erosion. High-pressure metamorphism, mineral zoning, radiometric dating and pressure-temperature-time paths in exhumed rocks reveal burial, reaction and return. Their reconstructable histories are especially useful when deformation and fluid inclusions are linked to the dated mineral growth. Earthquake distributions, seismic tomography, seismic attenuation and electrical conductivity constrain the geometry and physical state of the slab and mantle wedge, although fluids, melts and temperature can produce nonunique signals.

Arc magmas carry complementary chemical evidence: relatively mobile barium, boron, lead or strontium can trace fluid-mediated transfer, whereas sediment-sensitive thorium and isotope ratios can help identify other contributions. Element mobility changes with fluid composition and pressure-temperature conditions, so no one element is an infallible tracer. Boron isotope ratios, radiogenic strontium isotope ratios and lead isotope ratios should be interpreted with measured incoming sediments and altered oceanic crust. Arc gases and forearc fluid discharge constrain volatile return; ocean-island basalts, mantle xenoliths and diamond inclusions provide possible evidence of deeper recycling. Each observation samples only part of the system.

For any element $X$, a useful control-volume statement isHere permanent accretion is counted outside the chosen reactive control volume, return includes fluid, melt and exhumed-solid exports, and transient underplating or reservoir growth contributes to storage until exported. For a different boundary choice, the terms must be regrouped consistently. The three incoming layers can have very different transfer efficiencies for bulk rock, water, carbon and a particular trace element. A small mass of escaping fluid can dominate an element's return while most of the solid slab continues into the deep mantle.