Скачать или смотреть Case0: Base Case with Magmatic Crust Formation, High Melt Damage and Subaerial Lava Flow

This video presents a 2D high-resolution visco-elasto-plastic geodynamic model of lithospheric extension and magmatic crust formation with subaerial lava flows simulated using a cellular-automata approach. This is the base-case model in the magmatic crust formation play list. The melt damage factor is equal to 10 (High Melt Damage), the initial thermal anomaly in the asthenosphere (ΔTₘ) is 100 °C, the eruption interval (Δtₑᵣᵤₚₜ) is 50 thousand years and the fraction of melt extruded at the surface (𝑓ₑₓₜ,ₘₐₓ) ranges from 6% to 50% as the characteristic magmatic crust increases from 6km to 7.5km.

The model is initialized with a thermal anomaly in the asthenosphere (ΔTₘ), representing a rapidly spreading plume head beneath the lithosphere. This elevated thermal structure enhances partial melting and leads to the formation of thick magmatic crust during extension.

This model incorporates probabilistic frictional–plastic melt damage above melt-focusing regions to capture the weakening effects of melt networks and dike intrusions. This mechanism is essential for producing buried, seaward-dipping lava flows that accumulate toward a central zone of localized deformation. The localized cloud of transparent red colors in the thick volcanic flows above the shallow partially molten mantle are areas with relatively high plastic-strain rate associated with melt damage.

Model setup:

Domain size: 500 km × 160 km
Sticky-air layer: 10 km thick
Local mesh refinement: 200 m resolution from 150–350 km horizontally and to 50 km depth
Displacement limit for markers as a fraction of cell dimensions: 0.3
Convergence Criteria (L2 Norm of Velocity Solution): 1e-4
Rock properties: temperature dependent thermal conductivity and heat capacity, melt depletion density
Viscous Rheology: Composite, partial melt
Frictional-plastic Model: Strain weakening, melt damage, fluid pressure, time-dependent randomization using a randomization factor of 10

Visualization: This video shows a zoomed-in view centered on the zone where deformation localizes and magmatic crust builds.

The simulation was produced using the EarthBox geodynamic modeling framework.
Model input files for this case are available here (case_name=case0):
🔗 https://github.com/eakneller/EarthBox.jl