We present a numerical model of the main phase (370-335 Ma) of the Variscan orogeny in the central part of the Bohemian Massif. The crustal deformation in our model is driven by radiogenic heating in the felsic lower crust, the lateral contraction of the Moldanubian domain due to convergence with the Saxothuringian plate (in the early stage of orogeny), and the indentation of the Brunovistulian basement into the weakened orogenic root (in the late stage).
Our model explains the main geological events inferred from the geological record in the Moldanubian domain: formation of the orogenic plateau and onset of sedimentation at about 345 Ma, rapid exhumation of the orogenic lower crust at about 340 Ma and subsurface flow of crustal material (similar to 335 Ma and later). The results of our modeling suggest that delamination of the lithosphere, often invoked to explain the high temperature metamorphism in the orogenic lower crust of the Bohemian Massif, is not the only physical mechanism which can transfer a sufficient amount of heat to the crust to trigger its overturn.