Unexpected Deep Migration of Sediment in Coupled ASPECT and Fastscape Simulations

Hello,everyone
I am using the coupled software ASPECT and Fastscape. The model’s surface is initially set as continental crust , and the surface boundary condition is defined as sediment. After Fastscape returns a velocity, ASPECT deforms the mesh. However, if the mesh deforms only upwards, the sediment should thicken upwards while the points beneath the mesh should remain as continental crust. Nevertheless, in the actual model, I observe that the sediment thickens and extends beyond one grid thickness.


Are there any ideas for this?
Thanks!
Xiao

This is sediment composition.

@Xiao The general approach towards debugging issues such as this is to make the problem simpler – in fact, as simple as possible to still show the issue. What have you already investigated to figure out what is going on?
Best
W.

@Xiao - In addition to @bangerth suggestion regarding simplifying the problem, can you provide us a bit more information about the model that produced the image you posted? A short summary of the overall model design and prm file would be very helpful.

My initial reaction to the images you posted is that:

  1. There is some numerical diffusion occurring associated with compositional fields
  2. There are localized deformation structures advecting sediment downwards
  3. A higher resolution near the surface may be prudent

But as @bangerth pointed out, it would be ideal if you can reproduce the issue with a very simple model.

Cheers,
John