The motion of the Adriatic microplate is thought to be highly sensitive to the surrounding subduction zones and the convergence of Africa and Eurasia. However, our understanding of the mantle dynamics in the Mediterranean region and its effect on plate motion remains incomplete. Here, we present a large set of 3D thermomechanical models of the entire Mediterranean region over the last 35 Myr to understand what controls the motion of the Adriatic microplate. The simulations take the convergence of the African and Arabian plates with the Eurasian plate into account, along with the dynamics of the subduction systems in the western (Apennines-Calabria), central (Dinarides-Hellenides) Mediterranean and in the Alpine-Carpathian region. Our results demonstrate that the subduction systems around Adria are highly coupled, which gives rise to complex asthenospheric flow in the central Mediterranean. We find that the plate motion of the Adriatic microplate over the last 35 Myr is controlled by interplay of three main factors: 1) the convergence between the African and Eurasian plates, 2) the retreat of the Alpine subduction zone to the north of Adria, and 3) the distance between the Calabrian and Hellenic subduction zones around Adria. Furthermore, in a system characterized by active convergence between Africa and Eurasia, the slab pull exerted by nearby subduction zones can only notably influence the motion of the Adriatic microplate if these subduction zones are located within a few hundred kilometers of Adria.

The Yellowstone magmatic system is one of the largest magmatic systems on Earth, and thus an ideal location to study magmatic processes. Whereas previous seismic tomography results could only image a shallow magma reservoir, a recent study using more seismometers showed that a second and massive partially molten mush reservoir exists above the Moho \citep{huang2015yellowstone}. To understand the measurable surface response of this system to visco-elasto-plastic deformation, it is thus important to take the whole system from the mantle plume up to the shallow magma reservoirs into account. Here, we employ lithospheric-scale 3D visco-elasto-plastic geodynamic models to test the influence of parameters such as the connectivity of the reservoirs and rheology of the lithosphere on the dynamics of the system. A gravity inversion is used to constrain the effective density of the magma reservoirs, and an adjoint modelling approach reveals the key model parameters affecting the surface velocity. Model results show that a combination of connected reservoirs with plastic rheology can explain the recorded slow vertical surface uplift rates of around 1.2 cm/yr, as representing a long term background signal. A geodynamic inversion to fit the model to observed GPS surface velocities reveals that the magnitude of surface uplift varies strongly with the viscosity difference between the reservoirs and the crust. Even though stress directions have not been used as inversion parameters, modelled stress orientations are consistent with observations. However, phases of larger uplift velocities can also result from magma reservoir inflation which is a short term effect. We consider two approaches: 1) overpressure in the magma reservoir in the asthenosphere and 2) inflation of the uppermost reservoir prescribed by an internal kinematic boundary condition. We demonstrate that the asthenosphere inflation has a smaller effect on the surface velocities in comparison with the uppermost reservoir inflation. We show that the pure buoyant uplift of magma bodies in combination with magma reservoir inflation can explain (varying) observed uplift rates at the example of the Yellowstone volcanic system.