Basic fluid equations are the main ingredient in the development of theories of Rayleigh–
Taylor buoyancy-induced instability. Turbulence arises in the late stage of the instability …

This hands-on introduction to numerical geodynamic modelling provides a solid grounding
in the necessary mathematical theory and techniques, including continuum mechanics and …

High‐pressure metamorphic rocks provide evidence that in subduction zones material can
return from depths of more than 100 km to the surface. The pressure‐temperature paths …

Calculating surface topography in geodynamic models is a common numerical problem.
Besides other approaches, the so-called 'sticky air'approach has gained interest as a free …

We have designed a 2D thermal–mechanical code, incorporating both a characteristics
based marker-in-cell method and conservative finite-difference (FD) schemes. In this paper …

Initiation of diapirs is one of the least understood aspects of salt tectonics. Differential
sedimentary loading or erosion are both effective but not universal. A survey of 18 major salt …

We have extended our previous 2D method [Gerya, TV, Yuen, DA, 2003. Characteristics-
based marker-in-cell method with conservative finite-differences schemes for modeling …

Numerically modelling the dynamics of a self-consistently subducting lithosphere is a
challenging task because of the decoupling problems of the slab from the free surface. We …

This paper presents a model for the genesis of migmatites and granites during orogenic
evolution based on the analysis of several Phanerozoic crustal segments and a review of …

There has never been a convincing explanation of the way in which diapirs of molten granite
can effectively rise through mantle and crust. We argue here that this is mainly because the …