Deciphering Igneous Rock Crystals: Unveiling Multifractal Patterns in Crystal Size Dynamics

Understanding magma plumbing systems hinges upon an intricate comprehension of crystal populations concerning size, chemistry, and origin. We introduce an innovative, yet elegantly simple approach—the ‘number–length of crystals (N-LoC) multifractal model’—to classify crystal sizes, unveiling compelling insights into their distribution dynamics. This model, a departure from conventional crystal size distribution (CSD) diagrams, reveals multifractal patterns indicative of distinct class sizes within igneous rock crystals. By synthesizing multiple samples from experimental studies, natural occurrences, and numerical models, we validate this method’s efficacy. Our bi-logarithmic N-LoC diagrams for cooling-driven crystallized samples transcend the confines of traditional CSD plots, identifying variable thresholds linked to cooling rates and quenching temperatures. These thresholds hint at pulsative nucleation and size-dependent growth events, offering glimpses into crystallization regimes and post-growth modifications like coalescence and coarsening. Examining multifractal log–log plots across time-series samples unravels crystallization histories during cooling or decompression. Notably, microlites within volcano conduits delineate thresholds influenced by decompression rate and style, mirroring nucleation and growth dynamics observed in experimental studies. Our fractal methodology, presenting a more direct approach with fewer assumptions than the classic CSD method, stands poised as a potent alternative or complementary tool. We delve into its potential, facilitating comparisons between eruptive styles in volcanoes while deliberating on inherent limitations. This work not only advances crystal size analysis methodologies but also holds promise for inferring nuanced volcanic processes and offers a streamlined avenue for crystal size evaluation in igneous rocks.