The permeability of magma mush assembled from anisotropic tabular crystals

<p>The extraction of melt from a mush in a magma reservoir is of wide interest. All models for melt extraction from a mush require knowledge of mush permeability, and yet this remains poorly constrained. This permeability is typically calculated using the Kozeny-Carman model or variants thereof, which require a priori knowledge of the microstructural geometry. Such models are not calibrated or tested for packs of crystals of a range of shapes found in natural mush piles, leading to the potential for oversimplification of complex natural systems.</p><p>Essentially, a magma mush with minimal crystal-crystal intergrowth is composed of packed crystals where the pore space is filled with interstitial melt. Therefore, this can be studied as a granular medium. We use numerical methods to create domains of closely packed, randomly oriented cuboids in which we keep the short and intermediate axes lengths equal (i.e. square cross section) and vary the long axis magnitude. Our synthetic ‘crystals’ therefore cover the range from oblate to prolate, passing through a cubic shape. We supplement these with 3D numerical packs of spherical particles in cubic lattice arrangements or random arrangements. For the sphere packs we use various polydispersivity of sphere sizes. The permeability of all of these pack types is calculated using a numerical simulation (both LBflow and Avizo-based algorithms) with imposed periodic boundary conditions. The preliminary results suggest that the permeability of a granular medium scales with the specific surface area exclusively, without requiring prior knowledge of the geometry and size distribution of the particles.</p><p>We suggest that the model toward which we are working will allow magma mush permeability to be modelled more accurately. If our approach is embedded in existing continuum models for mush compaction and melt extraction, then more accurate estimates of melt accumulation rates prior to very large eruptions could be found.</p><p>Keywords: melt segregation, compaction, granular media, fluid flow, numerical simulation</p>

Hydrological Modeling of Highly Glacierized Basins (Andes, Alps, and Central Asia)

Soil and Water Assessment Tool (SWAT) was used to simulate five glacierized river basins that are global in coverage and vary in climate. The river basins included the Narayani (Nepal), Vakhsh (Central Asia), Rhone (Switzerland), Mendoza (Central Andes, Argentina), and Central Dry Andes (Chile) with a total area of 85,000 km2. A modified SWAT snow algorithm was applied in order to consider spatial variation of associated snow melt /accumulation by elevation band across each subbasin. In the previous studies, melt rates varied as a function of elevation resulting from an air temperature gradient while the snow parameters were constant throughout the entire basin. A major improvement of the new snow algorithm is separating the glaciers from seasonal snow based on their characteristics. Two SWAT snow algorithms were evaluated in simulation of monthly runoff from glaciered watershed: 1) the snow parameters are lumped (i.e. constant throughout the entire basin) and 2) the snow parameters are spatially variable based on elevation band-subbasin (i.e. modified snow algorithm). Applying the distributed SWAT snow algorithm improved the model performance in simulation of monthly runoff with snow-glacial regime, so that mean RSR decreased to 0.49 from 0.55 and NSE increased to 0.75 from 0.69. Improvement of model performance was negligible in simulation of monthly runoff from the basins with monsoon runoff regime.

Granite formation: Stepwise accumulation of melt or connected networks?

ABSTRACTSeveral authors have proposed that granitic melt accumulation and transport from the source region occurs in networks of connected melt-filled veins and dykes. These models envisage the smallest leucosomes as ‘rivulets’ that connect to feed larger dykes that form the ‘rivers’ through which magma ascends through the sub-solidus crust. This paper critically reviews this ‘rivulets-feeding-rivers’ model. It is argued that such melt-filled networks are unlikely to develop in nature, because melt flows and accumulates well before a fully connected network can be established. In the alternative stepwise accumulation model, flow and accumulation is transient in both space and time. Observations on migmatites at Port Navalo, France, that were used to support the existence of melt-filled networks are discussed and reinterpreted. In this interpretation, the structures in these migmatites are consistent with the collapse and draining of individual melt batches, supporting the stepwise accumulation model.