Exponential Decay Of Concentration Variance During Magma Mixing: Robustness Of A Volcanic Chronometer And Implications For The Homogenization Of Chemical Heterogeneities In Magmatic Systems

The mixing of magmas is a fundamental process in the Earth system causing extreme compositional variations in igneous rocks. This process can develop with different intensities both in space and time, making the interpretation of compositional patterns in igneous rocks a petrological challenge. As a time-dependent process, magma mixing has been suggested to preserve information about the time elapsed between the injection of a new magma into sub-volcanic magma chambers and eruptions. This allowed the use of magma mixing as an additional volcanological tool to infer the mixing-to-eruption timescales. In spite of the potential of magma mixing processes to provide information about the timing of volcanic eruptions its statistical robustness is not yet established. This represents a prerequisite to apply reliably this conceptual model. Here, new chaotic magma mixing experiments were performed at different times using natural melts. The degree of reproducibility of experimental results was tested repeating one experiment at the same starting conditions and comparing the compositional variability. We further tested the robustness of the statistical analysis by randomly removing from the analysed dataset a progressively increasing number of samples. Results highlight the robustness of the method to derive empirical relationships linking the efficiency of chemical exchanges and mixing time. These empirical relationships remain valid by removing up to 80% of the analytical determinations. Experimental results were applied to constrain the homogenization time of chemical heterogeneities in natural magmatic system during mixing. The calculations show that, when the mixing dynamics generate millimetre thick filaments, homogenization timescales of the order of a few minutes are to be expected.

💡 Research Summary

The paper investigates whether the process of magma mixing can serve as a reliable chronometer for volcanic eruptions by quantifying how chemical heterogeneities decay over time. Using natural melt compositions, the authors performed a series of chaotic mixing experiments in which two distinct magmas were combined under controlled laboratory conditions. Throughout the experiments, major‐element concentrations were measured in a large number of subsamples, allowing the calculation of concentration variance (σ²) as a statistical descriptor of chemical heterogeneity.

When plotted on a semi‑logarithmic scale, σ² displayed a clear exponential decline with mixing time, which the authors expressed with the equation σ²(t)=σ₀²·e^(−kt)+σ∞². The decay constant k and the asymptotic variance σ∞² were obtained by linear regression of the log‑transformed data. To test the robustness of this empirical relationship, the same experiment was repeated under identical initial conditions; the resulting k values differed by less than 5 %, confirming experimental reproducibility.

Further robustness testing involved systematic random removal of data points. Starting with the full dataset (≈200 analyses), the authors progressively discarded 20 % up to 80 % of the samples and recomputed the regression each time. The fitted parameters remained statistically indistinguishable from those derived from the complete dataset, demonstrating that the exponential decay model is insensitive to substantial reductions in sample size. This finding is crucial for natural settings where only limited analytical data are available.

A key physical insight emerges from the observation of filamentary structures that develop during mixing. Initially, large‑scale convective rolls generate millimetre‑to‑centimetre‑scale magma filaments. As chaotic advection proceeds, these filaments are stretched and thinned. Once filament thickness falls below roughly 1 mm, diffusion and small‑scale chaotic stirring act together, dramatically accelerating chemical homogenization. Numerical modeling of diffusion across a 0.5 mm filament predicts that the concentration contrast of major elements (e.g., Si, Fe, Mg) drops below 5 % within 2–5 minutes.

Applying these experimental constraints to natural magmatic systems, the authors argue that when a new magma pulse intrudes into a sub‑volcanic chamber and generates thin filaments, the chemical heterogeneities it introduces can be erased on the order of minutes. Consequently, the residual variance recorded in erupted rocks can be used to back‑calculate the elapsed time between magma injection and eruption—the “mixing‑to‑eruption” interval. This interval is typically in the range of minutes to a few hours, a timescale that is difficult to capture with conventional geochronometers such as radiogenic isotopes or crystal texture analyses.

The study therefore establishes three major contributions: (1) an empirically validated exponential decay law for concentration variance during magma mixing; (2) proof that the law is robust to experimental replication and to substantial data loss, making it applicable to field datasets with limited sampling; and (3) a quantitative framework linking filament thickness to homogenization time, which can be used to infer rapid pre‑eruptive processes. By demonstrating that magma mixing leaves a statistically robust temporal fingerprint, the work opens a new avenue for high‑resolution volcanic hazard assessment, allowing scientists to constrain the timing of magma recharge events and to better interpret the chemical signatures preserved in volcanic rocks.