Title: The architecture of intrusions in magmatic mush
Magmatic reservoirs located in the upper crust have been shown to result from the repeated intrusions of new magmas, and spend much of the time as a crystal-rich mush. The geometry of the intrusion of new magmas may greatly affect the thermal and compositional evolution of the reservoir. Despite advances in our understanding of the physical processes that may occur in a magmatic reservoir, the resulting architecture of the composite system remains poorly constrained. Here we performed numerical simulations coupling a computational fluid dynamics and a discrete element method in order to illuminate the geometry and emplacement dynamics of a new intrusion into mush and the relevant physical parameters controlling it. Our results show that the geometry of the intrusion is to first order controlled by the density contrast that exists between the melt phases of the intrusion and resident mush rather than the bulk density contrast as is usually assumed. When the intruded melt is denser than the host melt, the intrusion pounds at the base of the mush and emplaced as a horizontal layer. The occurrence of Rayleigh-Taylor instability leading to the rapid ascent of the intruded material through the mush was observed when the intruded melt was lighter than the host one and was also unrelated to the bulk density contrast. In the absence of density contrasts between the two melt phases, the intrusion may fluidize the host crystal network and slowly ascend through the mush. The effect of the viscosity contrast between the intruded and host materials was found to have a lesser importance on the architecture of intrusions in a mush. Analyzing the eruptive sequence of well documented eruptions involving an intrusion as the trigger shows a good agreement with our modeling results, highlighting the importance of specifically considering granular dynamics when evaluating magmas and mush physical processes. more »« less
Carrara, Alexandre; Burgisser, Alain; Bergantz, George W
(, Volcanica)
none
(Ed.)
Currently, our ability to interpret the mechanics of magma mingling and mixing is limited by an incomplete understanding of the modes of mixing across all melt fractions and compositions. Here, we present numerical simulations of the emplace- ment of crystal-free magma in crystal-rich reservoirs employing a computational fluid dynamics and discrete element method (CFD–DEM). We performed two runs corresponding to the emplacement of basalt into two end-member types of magmas mush (basaltic and dacitic). We found that the intruded volumes have similar shapes and are surrounded by a halo where the crys- tal volume fraction of the mush is lower. The dynamics of intruded melt are, however, different. Importantly, the mingling of the intruded and host materials starts after emplacement and consists in the incorporation of mush material into the intruded magma. Our findings imply that purely thermo-mechanical processes controlled by grain-scale dynamics are sufficient to explain fundamental aspects of recharge.
Carbotte, Suzanne M.; Arnulf, Adrien; Spiegelman, Marc; Lee, Michelle; Harding, Alistair; Kent, Graham; Canales, Juan Pablo; Nedimović, Mladen
(, Geology)
Abstract Magmatic systems are composed of melt accumulations and crystal mush that evolve with melt transport, contributing to igneous processes, volcano dynamics, and eruption triggering. Geophysical studies of active volcanoes have revealed details of shallow-level melt reservoirs, but little is known about fine-scale melt distribution at deeper levels dominated by crystal mush. Here, we present new seismic reflection images from Axial Seamount, northeastern Pacific Ocean, revealing a 3–5-km-wide conduit of vertically stacked melt lenses, with near-regular spacing of 300–450 m extending into the inferred mush zone of the mid-to-lower crust. This column of lenses underlies the shallowest melt-rich portion of the upper-crustal magma reservoir, where three dike intrusion and eruption events initiated. The pipe-like zone is similar in geometry and depth extent to the volcano inflation source modeled from geodetic records, and we infer that melt ascent by porous flow focused within the melt lens conduit led to the inflation-triggered eruptions. The multiple near-horizontal lenses are interpreted as melt-rich layers formed via mush compaction, an interpretation supported by one-dimensional numerical models of porous flow in a viscoelastic matrix.
Dick, Henry J.; Robinson, Paul T.; MacLeod, C.J. MacLeod; Kinoshita, Hajimu
(, Progress in earth and planetary science)
This paper presents the first detailed geologic map of in situ lower ocean crust; the product of six surveys of Atlantis Bank on the SW Indian Ridge. This combined with major and trace element compositions of primary magmatic phases in 99 seafloor gabbros shows there are both significant vertical and ridge-parallel variations in crustal composition and thickness, but a continuity of the basic stratigraphy parallel to spreading. This stratigraphy is not that of magmatic sedimentation in a large crustal magma chamber. Instead, it is the product of dynamic accretion where the lower crust formed by episodic intrusion, large-scale upward migration of interstitial melt due to crystal mush compaction, and continuous tectonic extension accompanied by hyper- and sub-solidus, crystal-plastic deformation. Five crossings of the gabbro-peridotite contact along the transform wall show that massive mantle peridotite is intruded by cumulate residues of moderately to highly evolved magmas, few of which could be even close to equilibrium with a primary mantle magma. This contact then does not represent the crust-mantle boundary as envisaged in the ophiolite analog for ocean crust. The residues of the magmas parental to the shallow crust must also lie beneath the center of the complex. This, and the nearly complete absence of dunites in peridotites from the transform wall, shows that melt transport through the shallow lithosphere was largely restricted to the central region of the paleo-ridge segment. There is almost no evidence for a melt lens or high-level storage of primitive melt in the upper 1500 m of Atlantis Bank. Thus, the composition of associated mid-ocean ridge basalt appears largely controlled by fractional crystallization of primitive cumulates at depth, near or at the base of the crust, modified somewhat by melt-rock reaction during transport through the overlying cumulate pile to the seafloor. Inliers of the dike-gabbro transition show that the uppermost gabbros crystallized at depth and were then emplaced upward, as they cooled, into the zone of diking. ODP and IODP drilling along the center of the gabbro massif also found few primitive gabbros that could have been in equilibrium with the original overlying lavas. Evidence of large-scale upward, permeable transport of interstitial melt through the gabbros is ubiquitous. Thus, post-cumulus processes, including extensive reaction, dissolution, and re-precipitation within the cumulate pile have obscured nearly all evidence of earlier primitive origins. We suggest that many of the gabbros may have started as primitive cumulates but were hybridized and transformed by later, migrating melts to evolved compositions, even as they ascended to higher levels, while new primitive cumulates were emplaced near the base of the crust. Mass balance for a likely parental melt intruded from the mantle to form the crust, however, requires that such primitive cumulates must exist at depth beneath Atlantis Bank at the center of the magmatic complex. The Atlantis Bank Gabbro Massif accreted by direct magma intrusion into the lower crust, followed by upward diapiric flow, first as a crystal mush, then by solid-state, crystal-plastic deformation, and finally by detachment faulting to the sea floor. The strongly asymmetric spreading to the south, parallel to the transform, was due to fault capture, with the bounding faults on the northern rift valley wall cut off by the detachment fault, which extended across the zone of intrusion causing rapid migration of the plate boundary to the north.
Abstract The Duluth Complex is a large mafic intrusive system located in northeastern Minnesota emplaced as part of the 1.1-Ga Midcontinent Rift. Several Fe–Ti oxide-bearing ultramafic intrusions are hosted along the Western Margin of the Duluth Complex, and are discordant bodies present in a variety of geometries, hosted in multiple rock types, and dominated by peridotite, pyroxenite, and semi-massive to massive Fe–Ti oxide rock types. Their origin has been debated, and here we present geochemical evidence and modeling that supports a purely magmatic origin for the Titac and Longnose Fe–Ti oxide-bearing ultramafic intrusions. Ilmenite and titanomagnetite textures indicate a protracted cooling process, and δ34S values of sulfides reveal little assimilation of the footwall Virginia Formation, a fine-grained pelitic unit that contains sulfide-rich bands. We model the crystallization of a hypothetical parental magma composition to the host intrusion of Longnose using Rhyolite-MELTS and demonstrate that the accumulation of Fe–Ti oxides in the discordant intrusions cannot be explained by density-driven segregation of crystallized Fe–Ti oxides. Instead, we show that the development of silicate liquid immiscibility, occurring by the unmixing of the silicate melt into conjugate Si- and Fe-rich melts, can result in the effective segregation and transportation of the Fe-rich melt. The Fe-rich melt is ~2 orders of magnitude less viscous than the Si-rich melt, allowing the Fe-rich melt to be more effectively segregated and transported in the mush regime (crystallinities >50%). This suggests that viscosity, in addition to density, plays a significant role in forming the discordant Fe–Ti oxide-bearing ultramafic intrusions. We propose a genetic model that could also be responsible for the Fe–Ti oxide-rich layers or bands that are hosted within the igneous stratigraphy of mafic intrusions of the Duluth Complex.
Sas, May; Shane, Phil; Kuritani, Takeshi; Zellmer, Georg F; Kent, Adam J; Nakagawa, Mitsuhiro
(, Journal of Petrology)
null
(Ed.)
Abstract The Okataina Volcanic Centre (OVC), located in the Taupo Volcanic Zone, New Zealand, is a dominantly rhyolitic magmatic system in an arc setting, where eruptions are thought to be driven by mafic recharge. Here, Sr–Pb isotopes, and compositional and textural variations in plagioclase phenocrysts from 10 rhyolitic deposits (two caldera, one immediately post-caldera, four intra-caldera, and three extra-caldera) are used to investigate the OVC magmatic system and identify the sources and assimilants within this diverse mush zone. Plagioclase interiors exhibit normal and reverse zoning, and are commonly in disequilibrium with their accompanying glass, melt inclusions, and whole-rock compositions. This indicates that the crystals nucleated in melts that differed from their carrier magma. In contrast, the outermost rims of crystals exhibit normal zoning that is compositionally consistent with growth in cooling and fractionating melts just prior to eruption. At the intra-crystal scale, the total suite of 87Sr/86Sr ratios are highly variable (0·7042–0·7065 ± 0·0004 average 2SE); however, the majority (95 %) of the crystals are internally homogeneous within error. At whole-crystal scale (where better precision is obtained), 87Sr/86Sr ratios are much more homogeneous (0·70512–0·70543 ± 0·00001 average 2SE) and overlap with their host whole-rock Sr isotopic ratios. Whole-crystal Pb isotopic ratios also largely overlap with whole-rock Pb ratios. The plagioclase and whole-rock isotopic compositions indicate significant crustal assimilation (≥20 %) of Torlesse-like metasediments (local basement rock) by a depleted mid-ocean ridge mantle magma source, and Pb isotopes require variable fluid-dominant subduction flux. The new data support previous petrogenetic models for OVC magmas that require crystal growth in compositionally and thermally distinct magmas within a complex of disconnected melt-and-mush reservoirs. These reservoirs were rejuvenated by underplating basaltic magmas that serve as an eruption trigger. However, the outermost rims of the plagioclase imply that interaction between silicic melts and eruption-triggering mafic influx is largely limited to heat and volatile transfer, and results in rapid mobilization and syn-eruption mixing of rhyolitic melts. Finally, relatively uniform isotopic compositions of plagioclase indicate balanced contributions from the crust and mantle over the lifespan of the OVC magmatic system.
Carrara, Alexandre, Burgisser, Alain, and Bergantz, George W. The architecture of intrusions in magmatic mush. Retrieved from https://par.nsf.gov/biblio/10555686. Earth and Planetary Science Letters 549.C Web. doi:10.1016/j.epsl.2020.116539.
Carrara, Alexandre, Burgisser, Alain, & Bergantz, George W. The architecture of intrusions in magmatic mush. Earth and Planetary Science Letters, 549 (C). Retrieved from https://par.nsf.gov/biblio/10555686. https://doi.org/10.1016/j.epsl.2020.116539
@article{osti_10555686,
place = {Country unknown/Code not available},
title = {The architecture of intrusions in magmatic mush},
url = {https://par.nsf.gov/biblio/10555686},
DOI = {10.1016/j.epsl.2020.116539},
abstractNote = {Magmatic reservoirs located in the upper crust have been shown to result from the repeated intrusions of new magmas, and spend much of the time as a crystal-rich mush. The geometry of the intrusion of new magmas may greatly affect the thermal and compositional evolution of the reservoir. Despite advances in our understanding of the physical processes that may occur in a magmatic reservoir, the resulting architecture of the composite system remains poorly constrained. Here we performed numerical simulations coupling a computational fluid dynamics and a discrete element method in order to illuminate the geometry and emplacement dynamics of a new intrusion into mush and the relevant physical parameters controlling it. Our results show that the geometry of the intrusion is to first order controlled by the density contrast that exists between the melt phases of the intrusion and resident mush rather than the bulk density contrast as is usually assumed. When the intruded melt is denser than the host melt, the intrusion pounds at the base of the mush and emplaced as a horizontal layer. The occurrence of Rayleigh-Taylor instability leading to the rapid ascent of the intruded material through the mush was observed when the intruded melt was lighter than the host one and was also unrelated to the bulk density contrast. In the absence of density contrasts between the two melt phases, the intrusion may fluidize the host crystal network and slowly ascend through the mush. The effect of the viscosity contrast between the intruded and host materials was found to have a lesser importance on the architecture of intrusions in a mush. Analyzing the eruptive sequence of well documented eruptions involving an intrusion as the trigger shows a good agreement with our modeling results, highlighting the importance of specifically considering granular dynamics when evaluating magmas and mush physical processes.},
journal = {Earth and Planetary Science Letters},
volume = {549},
number = {C},
publisher = {none},
author = {Carrara, Alexandre and Burgisser, Alain and Bergantz, George W},
editor = {none}
}
Warning: Leaving National Science Foundation Website
You are now leaving the National Science Foundation website to go to a non-government website.
Website:
NSF takes no responsibility for and exercises no control over the views expressed or the accuracy of
the information contained on this site. Also be aware that NSF's privacy policy does not apply to this site.
An official website of the United States government
