No abstract

Pegmatites are shallow, coarse-grained magmatic intrusions with crystals occasionally approaching meters in length. Compared to their plutonic hosts, pegmatites are thought to have cooled rapidly, suggesting that these large crystals must have grown fast. Growth rates and conditions, however, remain poorly constrained. Here we investigate quartz crystals and their trace element compositions from miarolitic cavities in the Stewart pegmatite in southern California, USA, to quantify crystal growth rates. Trace element concentrations deviate considerably from equilibrium and are best explained by kinetic effects associated with rapid crystal growth. Kinetic crystal growth theory is used to show that crystals accelerated from an initial growth rate of 10−6–10−7 m s−1 to 10−5–10−4 m s−1 (10-100 mm day−1 to 1–10 m day−1), indicating meter sized crystals could have formed within days, if these rates are sustained throughout pegmatite formation. The rapid growth rates require that quartz crystals grew from thin (micron scale) chemical boundary layers at the fluid-crystal interfaces. A strong advective component is required to sustain such thin boundary layers. Turbulent conditions (high Reynolds number) in these miarolitic cavities are shown to exist during crystallization, suggesting that volatile exsolution, crystallization, and cavity generation occur together.

No abstract

Trace element zoning is often used to unravel the crystallization history of phenocrysts in magmatic systems, but interpretation requires quantifying the relative importance of equilibrium versus disequilibrium. Published partition coefficients for phosphorous (P) in olivine vary by more than a factor of ten. After considering kinetic effects, a new equilibrium partition coefficient was extrapolated from a re-examination of natural and experimental systems, indicating that P partition coefficients in olivine are significantly over-estimated. These new partitioning constraints allow us to establish a theoretical P Equilibrium Fractionation Array (PEFA) for mid-ocean ridge basalts (MORBs), revealing that most olivines from MORBs have excess P (2–15 times PEFA) and are thus in disequilibrium. Using an independent case study of natural dendritic olivines, we show that such P enrichments can be explained by diffusion-limited incorporation of P during rapid crystal growth. If growth rate can be related to cooling, the rapid growth rates of olivines have implications for magma system dynamics, such as the size of magma bodies or where crystallization occurs within the body.

Summary The 2011-12 eruption at Cordón Caulle in Chile produced crystal-poor rhyolitic magma with crystal-rich mafic enclaves whose interstitial glass is of identical composition to the host rhyolite. Eruptible rhyolites are thought to be genetically associated with crystal-rich magma mushes, and the enclaves within the Cordón Caulle rhyolite support the existence of a magma mush from which the erupted magma was derived. Moreover, towards the end of the 2011-12 eruption, subsidence gave way to inflation that has on average been continuous through at least 2020. We hypothesize that magma segregation from a crystal mush could be the source of the observed inflation. Conceptually, magma withdrawal from a crystal-poor rhyolite reservoir caused its depressurization, which could have led to upward flow of interstitial melt within an underlying crystal mush, causing a new batch of magma to segregate and partially recharge the crystal-poor rhyolite body. Because the compressibility of the crystalline matrix of the mush is expected to be lower than that of the interstitial melt, which likely contains some fraction of volatile bubbles, this redistribution of melt would result in a net increase in volume of the system and in the observed inflation. We use numerical modeling of subsurface magma flow and storage to show under which conditions such a scenario is supported by geodetic and petrologic observations.