Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe role of volatiles in reactive melt transport in the asthenosphere
67
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Experimental studies of mantle petrology find that small concentrations of water and carbon dioxide have a large effect on the solidus temperature and distribution of melting in the upper mantle. However, it has remained unclear what effect small fractions of deep, volatile-rich melts have on melt transport and reactive melting in the shallow asthenosphere. Here we present theory and computations indicating that low-degree, reactive, volatile-rich melts cause channelisation of magmatic flow at depths approximately corresponding to the anhydrous solidus temperature. These results are obtained with a novel method to simulate the thermochemical evolution of the upper mantle in the presence of volatiles. The method uses a thermodynamically consistent framework for reactive, disequilibrium, multi-component melting. It is coupled with a system of equations representing conservation of mass, momentum, and energy for a partially molten grain aggregate. Application of this method in two-phase, three-component upwelling-column models demonstrates that it reproduces leading-order features of hydrated and carbonated peridotite melting; in particular, it captures the production of low-degree, volatile-rich melt at depths far below the volatile-free solidus. The models predict that segregation of volatile-rich, deep melts promotes a reactive channeling instability that creates fast and chemically isolated pathways of melt extraction. Reactive channeling occurs where volatile-rich melts flux the base of the silicate melting region, enhancing dissolution of fusible components from the ambient mantle. We find this effect to be similarly expressed for models of both hydrated and carbonated mantle melting. These findings indicate that despite their small concentrations, water and carbon dioxide have an important control on the extent and style of magma genesis, as well as on the dynamics of melt transport.
rate research
Read More
Multi-phase reactive transport processes are ubiquitous in igneous systems. A challenging aspect of modelling igneous phenomena is that they range from solid-dominated porous to liquid-dominated suspension flows and therefore entail a wide spectrum of rheological conditions, flow speeds, and length scales. Most previous models have been restricted to the two-phase limits of porous melt transport in deforming, partially molten rock and crystal settling in convecting magma bodies. The goal of this paper is to develop a framework that can capture igneous system from source to surface at all phase proportions including not only rock and melt but also an exsolved volatile phase. Here, we derive an n-phase reactive transport model building on the concepts of Mixture Theory, along with principles of Rational Thermodynamics and procedures of Non-equilibrium Thermodynamics. Our model operates at the macroscopic system scale and requires constitutive relations for fluxes within and transfers between phases, which are the processes that together give rise to reactive transport phenomena. We introduce a phase- and process-wise symmetrical formulation for fluxes and transfers of entropy, mass, momentum, and volume, and propose phenomenological coefficient closures that determine how fluxes and transfers respond to mechanical and thermodynamic forces. Finally, we demonstrate that the known limits of two-phase porous and suspension flow emerge as special cases of our general model and discuss some ramifications for modelling pertinent two- and three-phase flow problems in igneous systems.
Deep-Earth volatile cycles couple the mantle with near-surface reservoirs. Volatiles are emitted by volcanism and, in particular, from mid-ocean ridges, which are the most prolific source of basaltic volcanism. Estimates of volatile extraction from the asthenosphere beneath ridges typically rely on measurements of undegassed lavas combined with simple petrogenetic models of the mean degree of melting. Estimated volatile fluxes have large uncertainties; this is partly due to a poor understanding of how volatiles are transported by magma in the asthenosphere. Here, we assess the fate of mantle volatiles through numerical simulations of melting and melt transport at mid-ocean ridges. Our simulations are based on two-phase, magma/mantle dynamics theory coupled to idealised thermodynamic model of mantle melting in the presence of water and carbon dioxide. We combine simulation results with catalogued observations of all ridge segments to estimate a range of likely volatile output from the global mid-ocean ridge system. We thus predict global MOR crust production of 66-73 Gt/yr (22-24 km3/yr) and global volatile output of 52-110 Mt/yr, corresponding to mantle volatile contents of 100--200~ppm. We find that volatile extraction is limited: up to half of deep, volatile-rich melt is not focused to the axis but is rather deposited along the LAB. As these distal melts crystallise and fractionate, they metasomatise the base of the lithosphere, creating rheological heterogeneity that could contribute to the seismic signature of the LAB.
The reactive-infiltration instability, which develops when a porous matrix is dissolved by a flowing fluid, contains two important length scales. Here we outline a linear stability analysis that simultaneously incorporates both scales. We show that the commonly used thin-front model is a limiting case of a more general theory, which also includes convection-dominated dissolution as another special case. The wavelength of the instability is bounded from below, and lies in the range 1mm to 1km for physically reasonable flow rates and reaction rates. We obtain a closed form for the growth rate when the change in porosity is small.
Reactive infiltration instability (RII) drives the development of many natural and engineered flow systems. These are encountered e.g. in hydraulic fracturing, geologic carbon storage and well stimulation in enhanced oil recovery. The surface area of the rocks changes as the pore structure evolves. We combined a reactor network model with grey scale tomography to seek the morphological interpretation for differences among geometric, reactive and apparent surface areas of dissolving natural porous materials. The approach allowed us to delineate the experimentally convoluted variables and study independently the effects of initial geometry and macroscopic flowrate. Simulations based on North Sea chalk microstructure showed that geometric surface not only serves as the interface for water-rock interactions but also represents the regional transport heterogeneities that can be amplified indefinitely by dissolutive percolation. Hence, RII leads to channelization of the solid matrix, which results in fluid focusing and an increase in geometric surface area. Fluid focusing reduces the reactive surface area and the residence time of reactants, both of which amplify the differences in question, i.e. they are self-supporting. Our results also suggested that the growing and merging of microchannels near the fluid entrance leads to the macroscopic fast initial dissolution of chemically homogeneous materials.
The tendency of irreversible processes to generate entropy is the ultimate driving force for the evolution of nature. In engineering, entropy production is often used as a measure of usable energy losses. In this study we show that the analysis of the entropy production patterns can help understand the vastly diversified experimental observations of water-rock interactions in natural porous media. We first present a numerical scheme for the analysis of entropy production in dissolving porous media. Our scheme uses a greyscale digital model of natural chalk obtained by X-ray nanotomography. Greyscale models preserve structural heterogeneities with very high fidelity, which is essential for simulating a system dominated by infiltration instability. We focus on the coupling between two types of entropy production: the percolative entropy generated by dissipating the kinetic energy of fluid flow and the reactive entropy that originates from the consumption of chemical free energy. Their temporal patterns pinpoint three stages of microstructural evolution. We then show that the regional mixing deteriorates infiltration instability by reducing local variations in reactant distribution. In addition, we show that the microstructural evolution can be particularly sensitive to the initially present transport heterogeneities when the global flowrate is small. This dependence on flowrate indicates that the need to resolve the structural features of a porous system is greater when the residence time of the fluid is long.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
