An official website of the United States government
Abstract The rheology of crustal mushes is a crucial parameter controlling melt segregation and magma flow. However, the relations between mush dynamics and crystal size and shape distribution remain poorly understood because of the complexity of melt‐crystal and crystal‐crystal interactions. We performed analog experiments to characterize the mechanisms that control pore space reduction associated with repacking. Three suspensions of monodisperse particles with different geometries and aspect ratios (1:1, 2:1, 4:1) in a viscous fluid were tested. Our results show that particle aspect ratios strongly control the melt extraction processes. We identify two competing mechanisms that enable melt extraction at grain scale. The first mechanism leads to continuous deformation and melt extraction and is associated with “diffuse” frictional dissipation between neighboring particles. The second is stochastic, localized, and nearly instantaneous and is associated with the development and destruction of force chains percolating through the granular assembly.
more »
« less
Repacking in Compacting Mushes at Intermediate Melt Fractions: Constraints From Numerical Modeling and Phase Separation Experiments on Granular Media
Abstract Before large volumes of crystal poor rhyolites are mobilized as melt, they are extracted through the reduction of pore space within their corresponding crystal matrix (compaction). Petrological and mechanical models suggest that a significant fraction of this process occurs at intermediate melt fractions (ca. 0.3–0.6). The timescales associated with such extraction processes have important ramifications for volcanic hazards. However, it remains unclear how melt is redistributed at the grain‐scale and whether using continuum scale models for compaction is suitable to estimate extraction timescales at these melt fractions. To explore these issues, we develop and apply a two‐phase continuum model of compaction to two suites of analog phase separation experiments—one conducted at low and the other at high temperatures, T, and pressures, P. We characterize the ability of the crystal matrix to resist porosity change using parameterizations of granular phenomena and find that repacking explains both data sets well. A transition between compaction by repacking to melt‐enhanced grain boundary diffusion‐controlled creep near the maximum packing fraction of the mush may explain the difference in compaction rates inferred from high T + P experiments and measured in previous deformation experiments. When upscaling results to magmatic systems at intermediate melt fractions, repacking may provide an efficient mechanism to redistribute melt. Finally, outside nearly instantaneous force chain disruption events occasionally recorded in the low T + P experiments, melt loss is continuous, and two‐phase dynamics can be solved at the continuum scale with an effective matrix viscosity.
more »
« less
- Award ID(s):
- 2054414
- PAR ID:
- 10517809
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 7
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We performed a series of extrusion experiments on partially molten samples of forsterite plus 10 vol% of an anorthite‐rich melt to investigate melt segregation in a pipe‐extrusion geometry and test the predictions of two‐phase flow theory with viscous anisotropy. The employed flow geometry has not been experimentally investigated for partially molten rocks; however, numerical solutions for a similar, pipe‐Poiseuille geometry are available. Samples were extruded from a 6‐mm diameter reservoir into a 2‐mm diameter channel under a fixed normal stress at 1350°C and 0.1 MPa. The melt distribution in the channel was subsequently mapped with optical and backscattered electron microscopy and analyzed via quantitative image analysis. Melt segregated from the center toward the outer radius of the channel. The melt fraction at the wall increased with increasing extrusion duration and with increasing shear stress. The melt fraction profiles are parabolic with the melt fraction at the wall reaching 0.17–0.66, values 2 to 16 times higher than at the channel center. Segregation of melt toward the wall of the channel is consistent with base‐state melt segregation as predicted by two‐phase flow theory with viscous anisotropy. However, melt‐rich sheets inclined at a low angle to the wall, which are anticipated from two‐phase flow theory, were not observed, indicating that the compaction length is larger than the channel diameter. The results of our experiments are a test of two‐phase flow theory that includes viscous anisotropy, an essential theoretical frame work needed for modeling large‐scale melt migration and segregation in the upper mantle.more » « less
-
Abstract We conducted experiments to study melt migration in crystal‐rich mushes, with application to magma ascent within transcrustal magma reservoirs. Mushes with crystal volume fractions of 0.59–0.83 were prepared by hot‐pressing crushed borosilicate glass mixed with different proportions of quartz sand particles. Each experimental sample comprises stacked disks of mush and soda‐lime glass, a proxy for crystal‐free magma. Samples were subjected to confining pressures of 100–300 MPa and a temperature of 900°C (above the glass transition temperatures of the borosilicate and soda‐lime glasses) for up to 6 h. The bottom and circumference of the mush and soda lime disks experience the confining pressure, but the top of the mush disks is at room pressure, resulting in a pore‐pressure gradient across the mush layer. Following cooling and decompression, we determined the area fraction and morphology of soda‐lime melt that migrated into the mush layer during experiments. Melt fraction is more strongly correlated to crystal fraction than pore‐pressure gradient, increasing with crystal fraction before sharply decreasing as crystal fractions exceed 0.8. This change at 0.8 coincides with the transition from crystals in the mush moving during soda‐lime migration to crystals forming a continuous rigid network. In our experiments, melt migration occurred by viscous fingering, but near the mobile‐to‐rigid transition, melt migration is enhanced by additional capillary action. Our results indicate that magma migration may peak when rigid mushes “unlock” to become mobile. This transition may mark an increase in magma migration, a potential precursor to volcanic unrest and eruption.more » « less
-
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.more » « less
-
Abstract. Fluid and melt transport in the solid mantle can be modeled as a two-phase flow in which the liquid flow is resisted by the compaction of the viscously deforming solid mantle. Given the wide impact of liquid transport on the geodynamical and geochemical evolution of the Earth, the so-called “compaction equations” are increasingly being incorporated into geodynamical modeling studies. When implementing these equations, it is common to use a regularization technique to handle the porosity singularity in the dry mantle. Moreover, it is also common to enforce a positive porosity (liquid fraction) to avoid unphysical negative values of porosity. However, the effects of this “capped” porosity on the liquid flow and mass conservation have not been quantitatively evaluated. Here, we investigate these effects using a series of 1- and 2-dimensional numerical models implemented using the commercial finite-element package COMSOL Multiphysics®. The results of benchmarking experiments against a semi-analytical solution for 1- and 2-D solitary waves illustrate the successful implementation of the compaction equations. We show that the solutions are accurate when the element size is smaller than half of the compaction length. Furthermore, in time-evolving experiments where the solid is stationary (immobile), we show that the mass balance errors are similarly low for both the capped and uncapped (i.e., allowing negative porosity) experiments. When Couette flow, convective flow, or subduction corner flow of the solid mantle is assumed, the capped porosity leads to overestimations of the mass of liquid in the model domain and the mass flux of liquid across the model boundaries, resulting in intrinsic errors in mass conservation even if a high mesh resolution is used. Despite the errors in mass balance, however, the distributions of the positive porosity and peaks (largest positive liquid fractions) in both the uncapped and capped experiments are similar. Hence, the capping of porosity in the compaction equations can be reasonably used to assess the main pathways and first-order distribution of fluids and melts in the mantle.more » « less
