Cooling of a Magmatic System Under Thermal Chaotic Mixing

The cooling of a basaltic melt undergoing chaotic advection is studied numerically for a magma with a temperature-dependent viscosity in a two-dimensional (2D) cavity with moving boundary. Different statistical mixing and energy indicators are used to characterize the efficiency of cooling by thermal chaotic mixing. We show that different cooling rates can be obtained during the thermal mixing of a single basaltic magmatic batch undergoing chaotic advection. This process can induce complex temperature patterns inside the magma chamber. The emergence of chaotic dynamics strongly modulates the temperature fields over time and greatly increases the cooling rates. This mechanism has implications for the thermal lifetime of the magmatic body and may favor the appearance of chemical heterogeneities in the igneous system as a result of different crystallization rates. Results from this study also highlight that even a single magma batch can develop, under chaotic thermal advection, complex thermal and therefore compositional patterns resulting from different cooling rates, which can account for some natural features that, to date, have received unsatisfactory explanations, including the production of magmatic enclaves showing completely different cooling histories compared with the host magma, compositional zoning in mineral phases, and the generation of large-scale compositional zoning observed in many plutons worldwide.