Confusions about the melt migration and the ‘porosity’ compositional field

Hi all,

I’m trying to study the influence of lithospheric properties (such as lithospheric thickness, fracture zones) on plume material melting and melt migration.

I want to model melt migration with the Melt_global material model. In this subsection only melt and background rock seem to be distinguished. However, I need to further define the properties of the background solid.

I wonder if it is possible to represent different background solid by setting varying porosity in the initial compositional field. For example, setting high porosity in a local area of the lithosphere to add a magma channel. I have tried a simple setup-- set the lithospheric velocity to 0 and the initial composition model is defined as shown below:
compositional fields

However, the self-defined high porosity decreases rapidly in the early stage of simulation and approaches 0. :sob: :sob:

I don’t know if I misunderstood the porosity here. Do my Settings fail to meet the principles on which the program is based?

Best wishes,


Hi Judy,

Sorry for the slow reply (and I guess this also answers my question from the other thread; you do need magma transport for your problem). For what you’re trying to do, I think you need to pick a material model that is closest to what you want to do and then modify it yourself. For example, if you wanted to have a channel that makes it easy for melt to flow though, you could add a compositional field to one of the melt material models and give that a higher permeability, making it easier for melt to flow though a specific place.

The porosity is simply the volume fraction of melt at a given location, so if, for example, your temperature was below the solidus then all your melt would start to crystallize.

But studying melt migration in the setting of a plume interacting with the overlying lithosphere that has a visco-elasto-plastic rheology is a challenging problem; you would expect interactions between the melt and the fractures in the lithosphere, for example, and this is something that we haven’t implemented in ASPECT yet. So can you be a bit more specific about what features you would actually need for your problem, and then we can see what’s currently in ASPECT and what you would have to add and who to talk to. It could also be good to come to one of the ASPECT user meetings on Wednesday.

Best wishes,

Hi Juliane,

Thank you for your attention to this matter.

I read your paper published in 2016 and found that my previous understanding of porosity was wrong.

I want to calculate the magmatic yield (in other words, accretive crustal thickness). Actually, at present my main concern is how does melt get extracted through the lithosphere, rather than the complex interactions between the melt and the lithosphere.
I would like to know if and how Melt_Global implements the melt extraction process. Is it based on extraction model by Schmeling (2006)? How is the critical melt fraction defined?

I’m interested in what you say about the compositional field with higher permeability.

But I’m not quite sure what you mean. To define different permeabilities for different compositional fields, can it be achieved directly by defining the reference permeability parameter in the Melt_Global material model?

Thank you again for your patience and cooperation.

Best wishes,