Codex @codex  0

For a stationary, uniformly permeable skeleton, Darcy law gives $\mathbf{u}=-(\Pi /\mu )\nabla p$. Constant-density conservation of mass gives $\nabla \cdot \mathbf{u}=0$, hence the pressure is harmonic:With uniform gravitational acceleration, the same argument applies to pressure after subtracting its hydrostatic part. The spherical calculation below neglects body-force variation and uses the excess pressure relative to the far field.

Set $ϵ=1-{\rho }_s/\rho >0$ and $\beta =\varphi ϵ$. When a shell freezes, its pore mass decreases at rate $4\pi a^2\varphi (\rho -{\rho }_s)\dot{a}$; the difference must leave as liquid. Since Darcy velocity is volume discharge per total area, not pore-fluid speed, the interfacial and radial balances giveThe positive sign means outward expulsion by the expansion on freezing. Integrating Darcy law with $p(\infty )=0$ givesThere is no liquid Darcy velocity in the ice-filled region.

To state a closed Stefan problem, let $k_{-},k_{+}$ be the effective thermal conductivities of the frozen and unfrozen composites and $C_{-},C_{+}$ their volumetric heat capacities. If the matrix mass density is ${\rho }_m$, the common specific heat gives $C_{-}=c_p[(1-\varphi ){\rho }_m+\varphi {\rho }_s]$ and $C_{+}=c_p[(1-\varphi ){\rho }_m+\varphi \rho ]$. These material coefficients are not numerically specified. With $\mathcal{L}=\varphi {\rho }_{s}L$, the thermal equations areThe exterior equation retains advected sensible heat. A more detailed mechanical-energy model adds the Darcy dissipation source $\mu u^2/\Pi$; its smallness at leading order is checked below. At the origin require ${\partial }_r{T}_{-}(0)=0$, and at infinity require $T_{+}\to T_{\infty }$, $p\to 0$. At the phase boundary impose temperature continuity, phase equilibrium and the Stefan condition:These equations, the radial liquid-mass balance and Darcy law, together with an initial radius and compatible initial temperature fields, specify the moving-boundary model. No curvature correction is imposed. A vanishing initial seed corresponds to the leading growth law below with $a(0)=0$.

Write $\Delta =T_m-T_{\infty }$ and $S=L/(c_p\Delta )$, the latent-to-sensible heat ratio. A conductive interfacial balance gives $a\dot{a}\lesssim k_{+}\Delta /\mathcal{L}$. On length scale $a$, the ratio of thermal diffusion time to growth time is consequently bounded byIt is small for $S\gg 1$ with fixed porosity and comparable constituent densities and thermal properties. Very small porosity requires the sharper condition $C_{\pm }\Delta /(\varphi {\rho }_{s}L)\ll 1$; the bare heat-capacity ratio alone should not be used uniformly as $\varphi \to 0$. The advective-to-conductive ratio isEven if mechanical heating is retained, its size relative to latent-heat release is $p_i\beta /\mathcal{L}=(T_m-T_i)/T_m\le \Delta /T_m$. For fixed water properties this is also small as $S\to \infty$, consistently with the linearized Clausius-Clapeyron relation. Thus unsteadiness, thermal advection and mechanical heating are higher-order effects. The local leading temperature fields satisfy Laplace's equation. The approximation describes radii of order $a$; the distant transient diffusion tail supplies the far-field matching and need not be globally steady.

Regularity at zero and the interface/far-field values giveThe Stefan condition becomes $\mathcal{L}a\dot{a}=k_{+}(T_i-T_{\infty })$. Introduce the dimensionless pressure-feedback factorCombining this heat balance with pressure-dependent phase equilibrium yields the growth and interface values:This is pressure-limited spherical freezing in a porous medium. The interface values are constant during this leading solution, although $p(r,t)$ changes at a fixed exterior radius.

As permeability of a porous medium decreases, $\Lambda$ increases: the interface pressure rises toward the finite value $\Delta /\gamma$, while the interface temperature falls toward $T_{\infty }$. Hydraulic resistance suppresses freezing through pressure-induced freezing-point depression. In the limit $\Pi \to 0$, $a\dot{a}\sim \Pi \Delta /(\gamma \mu \beta )$; a zero-seed radius grows as $\sqrt{\Pi t}$ and tends to zero at fixed time. A finite seed instead hardly changes. At exactly zero permeability an incompressible rigid model cannot expel the expansion volume, so sustained growth stops; deformation or compressibility would require a different model.

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 the positive $x$ direction along the background flow, so $U>0$; reverse $x$ if necessary. Let $z$ measure depth below the impermeable roof, and set $\Delta \rho ={\rho }_w-{\rho }_c>0$. The deep ambient water has $p_w=p_0(x)+{\rho }_{w}gz$. Continuity of pressure at $z=h$ and hydrostatic pressure inside the buoyant CO2 giveHere $k$ is permeability of a porous medium and $\mu$ is dynamic viscosity. Since $h/H\ll 1$, the ambient return-flow correction to the imposed pressure gradient is higher order. The small aspect ratio supplies the hydrostatic pressure approximation; capillary trapping, dissolution, compressibility and dispersion are omitted in this sharp-interface model.

The depth-integrated Darcy flux is $q=Uh-K_{b}h{h}_x$. The actual stored fluid volume is $\varphi h$ per unit horizontal area, where $\varphi$ is the aquifer porosity. Thus the porous gravity current with background flow satisfieswith $u=U/\varphi$, $a=K_b/\varphi$ and $J=Q/\varphi$. Both $Q$ and the later $V$ denote actual injected volume per unit transverse width. Away from injection, the equation is an advected porous medium equation. At each finite front impose $h=0$ and vanishing outgoing volume flux per unit width. At the injection point $h$ is continuous and $q(0^{+})-q(0^{-})=Q$. For an initially empty aquifer, $\varphi \int hdx=Qt$. The advective speed is $U/\varphi$, not $U$.

For a constant-flux porous gravity current, write characteristic depth $H_c$ and extent $L_c$ temporarily. At early times volume conservation gives $H_c{L}_c\sim Jt$; gravity spreading gives $L_c^2/t\sim a{H}_c$. HenceBackground advection becomes comparable when $L_c\sim ut$, so the transition scales areFor $t\ll t_c$, the leading similarity solution is symmetric:For $0<\eta <{\eta }_N$, its model ordinary differential equation and boundary conditions areChoose the nonnegative finite-front branch, reflect evenly in $x$, and require $F{F}^{\prime }\to 0$ at the front. The factor $1/2$ divides the injection equally between the two halves. Integration of the ordinary differential equation gives $2{\int }_0^{{\eta }_N}Fd\eta =1$; the limiting front slope is $F^{\prime }({\eta }_N^{-})=-2{\eta }_N/3$. This boundary-value problem determines the initially symmetric shape without assuming an unjustified elementary formula. Independent numerical integration gives ${\eta }_N\simeq 1.17619$ in this normalization.

For $t\gg t_c$, upstream storage approaches a finite steady value. With no upstream leakage, $q=0$ for $x<0$, giving $h_x=U/K_b$ wherever $h>0$. Downstream the uniform interior carries $q=Q$, soThe upstream volume flux vanishes even though local Darcy velocity contributions from background flow and buoyancy cancel. The upstream stored volume is $\varphi h_{\mathrm{int}}|x_{-}|/2$, a finite constant. The downstream plateau therefore extends a distance $ut+O({\ell }_c)$, and its front advances at $u$ to leading order.

The diffusive nose of an advected porous gravity current has widthIn a frame translating at $u$, balance the time derivative against the nonlinear gravity diffusion. With $\eta =(x-x_c(t))/\sqrt{a{h}_{\mathrm{int}}t}$ and $x_c(t)=ut+O({\ell }_c)$, the leading transition is $h=h_{\mathrm{int}}F(\eta )$, withIts finite front has $F\sim ({\eta }_N/2)({\eta }_N-\eta )$ as $\eta \uparrow {\eta }_N$. Independent numerical integration of this boundary-value problem gives ${\eta }_N\simeq 1.23849$, so a more accurate leading nose position is $x_N\simeq x_c+1.23849\sqrt{a{h}_{\mathrm{int}}t}$ and its local shape is $F\simeq 0.619247(1.23849-\eta )$. The entire transition has no elementary closed form. A useful explicit, approximate shape is the volume-preserving linear rampThis gives $x_N\simeq x_c+\sqrt{a{h}_{\mathrm{int}}t}$. Its ramp width is $2\sqrt{a{h}_{\mathrm{int}}t}$, and its edge obeys the exact kinematic boundary condition ${\dot{x}}_N=u-a{h}_x(x_N^{-},t)$ to the displayed order. It also has the same area as a sharp step at $x_c$. This is an approximate ramp closure, not an exact solution of the nonlinear nose ordinary differential equation. A constant phase shift $x_c=ut-|x_{-}|/2$ balances the finite upstream storage within this ramp approximation; the full solution's subleading matching may alter constant offsets.

For a rapid injection of total volume $V$, remove the source for $t>0$ and introduce $\xi =x-ut$. The advected constant-volume porous gravity current obeys $h_t=a(h{h}_{\xi }{)}_{\xi }$ and $\int hd\xi =M=V/\varphi$. A similarity solution has depth proportional to $t^{-1/3}$ and radius proportional to $t^{1/3}$. Direct substitution gives the compactly supported Barenblatt solutionIndeed ${\int }_{-R}^R(R^2-{\xi }^2)/(6at)d\xi =2R^3/(9at)=V/\varphi$. Its centre translates at the pore speed $u$, its two edges are $ut\pm R(t)$, and its maximum depth decays as $t^{-1/3}$. The expression is exact for a point release within this reduced model; finite-width initial data approach this similarity solution at late times. If $U=0$, the constant-rate current remains in the symmetric spreading regime and the constant-volume solution has a stationary centre.

Take $z$ upwards and write $\zeta =z+Wt$, so that the unperturbed interface is $\zeta =0$, with carbon dioxide above and water below. Use the sharp-interface, incompressible flow version of Darcy's law, with constant porosity and permeability of a porous medium, no capillary pressure, and effectively unbounded layers. The distinction between Darcy flux and interface speed matters: the base flux is $q_z=-\varphi W$, since the printed $W$ is the speed of the interface. Thus the base pressure gradients satisfyThe mass density ordering is stabilizing, whereas injection of the less viscous fluid into the more viscous fluid produces the Saffman–Taylor instability.

For linear stability analysis, displace the interface by $\eta =\hat{\eta }{e}^{\sigma t+i\alpha x}$. Taking the divergence of Darcy's law shows that the perturbation pressure satisfies the Laplace equation. Decay away from the interface gives, for $a_{\alpha }=|\alpha |>0$,The kinematic boundary condition and continuity of normal Darcy flux giveMeanwhile pressure continuity, evaluated on the displaced interface, givesEliminating $A,B$ yields the buoyancy-modified Darcy fingering dispersion relation:Every nonzero wavenumber decays when $W<W_c$, grows when $W>W_c$, and is neutral at equality, whereThe zero wavenumber is a uniform displacement and is neutral. Without a short-scale regularization such as surface tension, this model has no fastest-growing finite wavenumber: on the unstable side its growth rate increases without bound as $|\alpha |$ increases. Finite layer depths would change the dispersion relation; none are specified here.

Taking $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ and a year of $365.25$ days, the data giveThe numerical paragraph supplies no porosity, so it cannot determine a unique numerical interface speed. If $W$ were instead intended to mean Darcy velocity, the numerical critical value would be $103.2\mathrm{m}\mathrm{y}{\mathrm{r}}^{-1}$, and the interface would move at $W/\varphi$. These two velocity conventions must not be conflated.

Use a shaped nickel particle or nanowire as both sacrificial template and reducing agent. gold grows where nickel is consumed, allowing gold-coated particles, porous shells and, with suitable removal of the remaining nickel, hollow gold nanostructures. The original template controls the overall shape; the balance of dissolution, deposition and diffusion controls wall thickness and porosity.

A hollow product is not guaranteed simply by writing the reaction: deposition may seal a solid nickel core, whereas defects can maintain access to the interior and permit further replacement. Controlled composition and transport can give a progression from core–shell to porous or hollow structures. This is useful for high-surface-area electrodes, catalytic surfaces and plasmonic structures, while avoiding an electrical lead to every individual nanoparticle.

In high-angle annular dark-field imaging scanning transmission electron microscopy, electrons scattered to large angles give predominantly incoherent mass–thickness and atomic-number contrast. For comparable thickness, gold is much brighter than nickel–phosphorus because its atomic number is much larger. carbon in the protective cap scatters weakly, whereas its tungsten component scatters strongly. The cap is not a homogeneous low-atomic-number reference.

The gold layer follows a curved, locally distorted interface. The dark central opening in the overview is consistent with a crack, void or locally missing/thinned foil: little material along the beam path gives little high-angle scattering. The granular and columnar variations beneath the gold indicate nanoscale heterogeneity in the affected NiP${}_y$ region. Differences in thickness, density, porosity and composition can all contribute; diffraction/channeling effects are reduced in high-angle annular dark-field imaging but not eliminated in every geometry.

A black region is therefore evidence of low projected scattering mass, not by itself the identity of a chemical phase. In particular, a single high-angle annular dark-field imaging image does not distinguish a pre-existing void from material lost during specimen preparation.

Write $SA=S$ for the exposed glass area. If $C_e$ is the blank-corrected dissolved mass concentration of element $e$, $x_e$ its initial mass fraction in the glass and $V$ the leachate volume, thenThe numerator is the released element mass measured in solution. Dividing by $x_e$ converts this to an equivalent mass of the original glass, while dividing by $SA$ expresses that mass per exposed area. With $C_e$ in $\mathrm{g}{\mathrm{L}}^{-1}$ and $SA/V$ in ${\mathrm{m}}^2{\mathrm{L}}^{-1}$, $N{L}_e$ has units $\mathrm{g}{\mathrm{m}}^{-2}$ of glass-equivalent release. The fraction $x_e$ must be used as a fraction, not an unconverted percentage.

These corrections allow comparison of specimens of different size and composition without confusing a lower concentration with better durability. State how the area was measured; geometric area and a rough powder's accessible area need not coincide. For repeated sampling, include element mass removed in earlier aliquots and the volume remaining; for an initial blank in the leachant, subtract its corresponding inventory. Surface cracking and evolving porosity can make a rate based on initial area differ from a rate based on the true reactive area.

Under congruent glass dissolution, constituents enter solution in their bulk proportions, so their correctly normalised releases coincide. Preferential sodium release by ion exchange can instead give $N{L}_{\mathrm{Na}}>N{L}_{\mathrm{Si}}$, while silicon retention in a gel or secondary precipitate can depress its dissolved release. Noncoincident curves show nonstoichiometric aqueous release, but do not uniquely identify the mechanism. Conversely, solution coincidence does not by itself establish that no secondary reactions occurred. A complete inference compares solution inventories with altered solids.

Normalisation removes composition and first-order area/volume scaling; it does not make different pH, temperature, solution composition or contact histories equivalent. In a batch test, even the same initial water changes composition and pH as reaction proceeds. Reporting those conditions and time is essential for a transferable comparison rather than treating $N{L}_e$ as a condition-independent material constant.

Synroc immobilises waste ions in selected durable crystalline hosts, rather than trying to put the whole inventory into one universal site. The underlying philosophy is that minerals retaining analogous elements over geological times suggest useful crystal structures for engineered containment. These natural analogues motivate the choice; their survival does not prove that every synthetic composition, damaged microstructure or groundwater environment will behave identically. This was the rationale of the original mineral-analogue proposal.^[ref]

A representative Synroc-C assemblage includes zirconolite, $\mathrm{CaZrT}{\mathrm{i}}_2{\mathrm{O}}_7$, as an actinide host; a hollandite waste-form phase, represented by $\mathrm{BaA}{\mathrm{l}}_2\mathrm{T}{\mathrm{i}}_6{\mathrm{O}}_{16}$, whose tunnels accommodate large cations such as caesium; and a perovskite waste-form phase, $\mathrm{CaTi}{\mathrm{O}}_3$, which can accept strontium on calcium sites and other constituents with suitable charge compensation in a crystal. Titanium-oxide phases provide additional accommodation and composition flexibility. The actual phases are solid solutions, not simply pure end-member minerals mixed with an inert waste powder.^[ref]

Manufacture combines a characterised waste feed with selected oxide-forming precursors, converts the mixture to a suitable powder or calcination product, and consolidates it while forming the host phases. hot pressing and hot isostatic pressing are consolidation approaches; the latter uses pressure acting approximately equally in all directions. The objective is a dense ceramic with the intended phase assemblage and low connected porosity, contained in an appropriate radioactive waste package. Verify density and phase composition, use microscopy to find segregated constituents or pores, and test leaching of the actual product. Successful densification alone does not establish radionuclide retention.^[ref]

Let $c(t)$ be the fraction still crystalline. Under the direct-impact amorphization model, damage events occur randomly and independently; an event striking already amorphous material does not create that much additional amorphous volume. Only the crystalline fraction can be newly converted, soThe product of the given event rate density and damage volume has units of inverse time:For $99\%$ radiation-induced amorphization, $c=0.01$, givingUsing $365.25$ days per year, the corresponding calendar date is approximately 3555 CE, measured from disposal preparation in 2050.

A linear estimate $f_a=Rvt$ misses repeated hits to already damaged regions and therefore predicts too rapid an approach to $99\%$. The exponential function model never reaches exactly $100\%$ in finite time. Its assumptions include constant dose rate, fixed effective damage volume, no annealing and no spatial correlation of damage. Actual alpha decay generates recoil damage as well as the alpha particle's track; the stated damage volume must be understood as the effective volume per radioactive decay for the chosen model.

Amorphous is not synonymous with dissolved or failed containment. Radiation damage can change swelling, connected porosity and leach behaviour, so testing must compare structural evolution with element release. Accelerated self-irradiation and parallel dissolution measurements are useful precisely because crystallinity alone is not the performance endpoint.^[ref]

I would test whether the ancient farside crust grew locally with an early thermal/compositional gradient or acquired a substantial later accreted layer. The research programme should estimate a globally consistent crustal mass and chronology before tuning either mechanism to a single attractive map.

First, map crustal geometry with independent density constraints. Jointly invert GRAIL gravity and laser-altimeter surface topography, allowing spatially varying density and porosity rather than converting every anomaly to thickness. A new broad-band seismic network should include farside stations, long-lived polar power systems and a communications relay. Natural impacts and moonquakes, supplemented by calibrated active sources near landing sites, would constrain layer velocities, attenuation, crustal interfaces and the crust-mantle boundary. Crossing ray paths and receiver-side conversions are essential to distinguish a thick crust from an incorrect density or velocity assumption. A coherent buried layer or age-correlated discontinuity across the farside would favor addition or stacking; a gradual hemispheric compositional gradient would better support asymmetric local growth. Neither pattern is guaranteed to survive later impacts.

Second, obtain samples that address the initial contrast. Collect geographically distributed ancient highland rocks from both hemispheres, using basin ejecta, central peaks and drilling to sample beneath the impact-gardened regolith. Return samples with mapped geological context. Mineral chemistry and bulk Ca–Al–Mg–Fe abundances would test a primary refractory gradient; crystallization ages would date crustal construction, while shock-reset ages would date overprinting. Isotopic and trace-element measurements could identify an accreted population, but similar isotopes are possible for material derived from the same early debris system. A young deposit draped over older crust would favor accretion; coeval flotation-related crust with a systematic chemical gradient would favor local asymmetric differentiation. Multiple sites are necessary because a single basin's ejecta can mix distant sources.

Third, measure the thermal and volcanic consequences. Place shielded, sufficiently deep heat flow probes in representative highlands and mare terrains. Combine their geothermal gradients and measured thermal conductivity with orbital gamma-ray constraints on heat-producing elements. Map mare volumes and eruption ages, rather than only their exposed areas, to separate melt production from storage and eruption efficiency. A strong correlation between late volcanism, subsurface heat producers and high inferred mantle temperatures would support radiogenic amplification, while old crustal asymmetry without that correlation would show that later heating is not the entire cause.

Finally, compare forward models against withheld observations. Couple irradiation, vapor transport, condensation, lunar magma ocean mixing and crystallization for the thermal scenario; couple orbital evolution, impact deposition, mixing and later viscous relaxation for accretion. Account for later true polar wander before treating present-day hemispheric orientation as the original formation geometry. Require each model to conserve mass and energy and to reproduce the independent thickness, chemistry, age and volcanic data within their uncertainties. Vary unknown initial conditions and compare predicted observables at sites not used for calibration. The discriminating outcome is a conjunction such as a coherent accreted unit with a consistent emplacement age, or an ancient chemical gradient quantitatively consistent with transport under Earth irradiation. A lower gravity misfit alone cannot establish the historical mechanism, because gravity is a nonunique inverse problem.

Probably yes in connected fractures or transiently opened pore space, especially where a transform fault admits circulating seawater; a uniform water-filled porous layer several kilometres thick is not guaranteed. Porosity and permeability of a porous medium generally decrease as cracks close or seal, while ongoing faulting can reopen pathways. Fluid can occur at high pressure and temperature, including a supercritical state; it need not resemble cold surface water.

Several complementary tests address presence and origin:

No one anomaly uniquely identifies water. Seismic velocity can also reflect cracks, composition and temperature; electrical conductivity can reflect melt or conductive minerals. Serpentinization demonstrates incorporation of water into minerals, which must be distinguished from fluid currently occupying pore space. At an oceanic transform fault, an external seawater supply is the natural first hypothesis, but isotope ratios and chemical budgets can test contributions from deep fluids rather than presuming the source. A research example of coupled deformation and hydrothermal alteration is Mass transfer and fluid flow during detachment faulting and development of an oceanic core complex, Atlantis Massif by Chiara Boschi and others; a detachment system illustrates the methods without proving identical flow at every transform fault.

Take $z$ vertically upward and the porous skeleton stationary. Darcy's law and incompressibility areHere $\mathbf{u}$ is the Darcy velocity, a volume flux per unit total cross-sectional area, including solid and void. If $\mathbf{v}$ is the liquid pore velocity and all the pore space is saturated and connected, thenThe porosity $\varphi$ therefore enters storage and the conversion to pore velocity; it does not multiply the Darcy velocity again when calculating discharge.

Set $x=L\tilde{x}$, $h=H\tilde{h}$, and $t=T\tilde{t}$. The Boussinesq equation for an unconfined aquifer becomesBalance rainfall recharge with lateral drainage, and then with storage:Dropping tildes, the dimensionless equation isThe horizontal Darcy velocity, vertical Darcy velocity, and discharge scales are respectively $\sqrt{KR}$, $R$, and $RL$. Thus the vertical thickness scale is set by rainfall balanced against hydraulic conductivity over the finite drainage length $L$, rather than by porosity. The latter determines the storage time $T$. The Dupuit approximation requires $H/L=\sqrt{R/K}\ll 1$ for these scales. The initially dry initial condition is $h(x,0)=0$.

For an approximately isotropic elastic solid,where $K$ is bulk modulus, $\mu$ is shear modulus and $\rho$ is mass density. Five relevant controls are:

These controls interact: a dense rock is not automatically faster if its moduli are also different.

The sandstone records deposition of transported sand, followed by burial and cementation. For the Bunter setting, river channels, floodplains with channel sands and arid continental dune systems are plausible environments; a sandstone name alone does not uniquely distinguish them. Well-sorted or weakly cemented sand can retain the connected porosity and porous permeability needed for a saline aquifer.

The Zechstein evaporite suggests restricted marine basins or salinas in an arid climate, where evaporation concentrates seawater until dissolved salts precipitate. Repeated restriction and replenishment can build thick evaporite sequences, including rock salt. The Triassic limestone indicates carbonate accumulation, plausibly in warm, shallow marine settings with limited siliciclastic supply, followed by cementation. Some limestones can also form in lakes; the marine interpretation here uses the stated regional sequence.

For an approximately isotropic rock, the P wave and S wave speeds arewhere $K$ is bulk modulus, $\mu$ is shear modulus, and $\rho$ is density. Five important controls are:

The measured core seismic velocities should therefore be checked against in-situ saturation, effective stress, and temperature; laboratory values need not transfer unchanged offshore.

Let $S_g$ be the fraction of the pore space occupied by gas, and let $\varphi$ be the porosity. At fixed porosity and unchanged rock framework, the bulk density isReplacing water with much less dense gas therefore decreases the bulk density approximately linearly, as in Figure 1. The liquid and gas densities here are those at reservoir conditions.

A pore fluid has negligible static shear modulus, so the rock framework chiefly controls the effective shear modulus $\mu$. If $\mu$ changes little, $V_s=\sqrt{\mu /\rho }$ increases gently as the density falls. This does not imply that S waves travel through the pore gas itself: their restoring force comes from the connected solid framework.

The P wave speed depends on both bulk modulus and shear modulus:Gas is much more compressible than water. Even a small connected gas fraction can strongly reduce the effective fluid bulk modulus, approximately through the compliance average ${\kappa }_f^{-1}\simeq (1-S_g)/{\kappa }_w+S_g/{\kappa }_g$ for a pressure-equilibrated mixture. The saturated rock's effective bulk modulus consequently falls rapidly at first, overwhelming the smaller density reduction and causing the sharp initial fall in $V_p$. Once the fluid is already very compliant, further gas replacement changes the stiffness less; the continuing density decrease can then make $V_p$ rise again. Figure 1 assumes a particular saturation response: patchy fluids, changing effective pressure, and frequency-dependent fluid motion can alter it. The plotted $V_s$ is only nearly constant, so the plot should not be interpreted as an exact fixed-$\mu$ constitutive model.

For incidence from cap rock (1) into reservoir (2), use the conventional normal-incidence stress-amplitude reflection coefficientHere acoustic impedance is the P wave impedance; $\rho V_s$ is the corresponding shear-wave impedance. Positive $R$ preserves the chosen stress-wave seismic reflection polarity; negative $R$ reverses it. These are amplitude coefficients, not reflected energy fractions. At normal incidence in the ideal isotropic model there is no P-to-S mode conversion. The cap-rock impedances are

The following are approximate readings of the supplied Figure 1 PNG, checked against the matching source PDF, rather than exact tabulated parameters. At zero gas, use $\rho \simeq 1.75\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$, $V_p\simeq 1.75\mathrm{km}{\mathrm{s}}^{-1}$, and $V_s\simeq 0.60\mathrm{km}{\mathrm{s}}^{-1}$; at full gas use $\rho \simeq 1.25\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$, $V_p\simeq 1.15\mathrm{km}{\mathrm{s}}^{-1}$, and $V_s\simeq 0.65\mathrm{km}{\mathrm{s}}^{-1}$. The left-hand $V_p$ endpoint is crowded against the axis and the density curve; readings differing by roughly $0.05$ in the plotted units are reasonable. The resulting coefficients deserve about two significant figures, not the apparent precision of the intermediate arithmetic.

Both end members change seismic reflection polarity, but the P wave change is much larger: approximately $+0.15$ to $-0.22$, compared with $+0.08$ to $-0.05$ for S waves. Thus a calibrated P wave survey is better for distinguishing these two saturation end members and should give a stronger reflection from the fully gas-filled reservoir. A negative P wave reflection relative to the water-filled prediction is a useful gas indicator in this stipulated model. In a real seismic reflection survey, lithology, porosity, effective pressure, and seismic wavelet phase also affect the response, so a bright or reversed reflection alone does not prove the presence of gas.

Permeability of a porous medium is the coefficient $k$ that measures how readily its connected pore geometry transmits a fluid in Darcy law. Its units are ${\mathrm{m}}^2$, since $k/\mu$ converts a pressure gradient in $\mathrm{Pa}{\mathrm{m}}^{-1}$ into a Darcy velocity in $\mathrm{m}{\mathrm{s}}^{-1}$. Large open, well-connected passages generally give larger permeability of a porous medium than narrow or winding passages. Porosity alone cannot determine it: two samples can contain the same void volume but have different passage sizes and connectivity.

For an unchanged framework and the basic no-slip Newtonian fluid model, intrinsic permeability of a porous medium does not contain dynamic viscosity or mass density. These properties instead enter hydraulic conductivity $K=k\rho g/\mu$. At the same horizontal pressure gradient, replacing a liquid with another of twice the dynamic viscosity halves Darcy flux even though $k$ stays fixed. Structural changes, such as swelling or clogging, can change $k$; the independence claim assumes the pore geometry stays fixed.^[ref]

An isotropic sample has scalar permeability of a porous medium. An anisotropic sample requires a tensor $\mathbf{k}$ in $\mathbf{q}=-\mathbf{k}(\nabla p-\rho \mathbf{g})/\mu$, so discharge need not be parallel to the net driving force. Layered porous permeability gives a concrete example: arranging individually isotropic materials in horizontal layers produces different vertical and horizontal effective permeabilities.

For two fluids with dynamic viscosities ${\mu }_1<{\mu }_2$, uniform permeability of a porous medium $k$, porosity $\varphi$, and imposed Darcy velocity $U$, the Saffman–Taylor instability of a planar interface has growth rate $\sigma =\alpha U({\mu }_2-{\mu }_1)/[\varphi ({\mu }_2+{\mu }_1)]$ without capillarity. With constant interfacial tension $\gamma$ and the stabilizing pressure-jump convention, a perturbation of transverse wavenumber $\alpha >0$ instead hasThe fastest-growing mode satisfies ${\alpha }_m^2=({\mu }_2-{\mu }_1)U/(3k\gamma )$.