The interaction of a subducting lithospheric slab with a chemical or phase boundary

We have used a finite element model of time‐dependent convection to determine the conditions for penetration of the subducted plate into the lower mantle. A temperature‐dependent and non‐Newtonian rheology is applied to achieve platelike behavior of the upper and sinking thermal boundary layer of convection. The 650‐km discontinuity is taken as either a chemical or phase boundary or as a combination of both. It is represented by a marker chain which effects additional buoyancy when distorted out of its equilibrium position. When the compositional density contrast is greater than about 5%, the descending slab is deflected sidewards at the boundary and two‐layer convection prevails. A resulting depression of the boundary in the range of 50–200 km should be detectable with seismic methods. Below 5% density difference the slab plunges several hundred kilometers into the lower mantle, and below 2% it will probably not stop before reaching the core‐mantle boundary and extensive mixing would be expected. With a pure phase change a negative Clapeyron slope of about −6 MPa/K (−60 bar/K) is required to establish a type of “leaky” double‐layer convection. A more moderate slope can aid a small compositional density difference to prevent slab penetration into the lower mantle. With the present uncertainties about the physical nature of the 650‐km discontinuity, a variety of convective styles appears possible on dynamical grounds.