Subduction zones are one of the most important tectonic environments on our planet, with much volcanism occuring above the subducting slab. In this project we study thermo-chemical structures and mantle flow above the subducting slab, with the goal of understanding the pattern of surface volcanism and crustal growth and seismological observations of these regions. This investigation is performed using numerical models. A key feature of our models is three-dimensionality (3-D), which is important because the structures that occur in nature are three-dimensional. This research builds on our 3-D petrological-thermomechanical models funded by a previous SNF grant, which showed that 3-D thermal-chemical convection patterns (and hence melt productivity distribution) in intra-oceanic subduction zones can be strongly variable and depend on both the magnitude and spatial distribution of chemical buoyancy above slabs. The objective of this new project is to develop a robust 3-D model of oceanic-continental subduction with magma generation and crustal growth, and to achieve a better understanding the fluid/melt generation and propagation atop the subducting slab in the southern Alaskan subduction zone, for which good seismic and petrological information exists. Our 3-D simulations will therefore account for: (i) mantle flow associated with the specific plate convergence rate, (ii) slab water release, (iii) slab fluid propagation that will trigger partial melting at the slab surface, (iv) melt ascent, (vi) variations in density and viscosity of both fluids and rocks as a function of local conditions and (v) relations between melting and melt transport dynamics in 3-D and crustal growth (volcanic activity) at the surface. By comparing with seismic structure and petrological observations in the southern Alaskan subduction zone, we will get a robust oceanic-continental subduction model with well- constrained controlling parameters particularly applicable for this zone. This project funds a postdoctoral scholar for 1.5 years.