We present results from two-dimensional (2-D) numerical experiments on the thermal and dynamical evolution of the subducting slab and of the overlying mantle wedge for a range in subduction parameters. These include subduction rate and the age and rheology of both subducting and overriding plates. Experiments also consider the influence of slab forcing conditions (from purely kinematic to purely dynamic) on the evolution of both the slab and mantle wedge. One goal is to determine how different parameters control thermal evolution of the slab-wedge interface, from just after subduction initiation up through roughly 500--600 km of subduction, where temperatures are approaching steady state. An additional goal is to define optimal conditions for the melting of slab sediments and crust. Results show slab surface temperatures (SSTs) depend strongly on subduction velocity, plate thermal structure, and upper mantle (or wedge) viscosity structure. Fast subduction beneath a thick (>70 km) overriding plate results in the coolest SSTs. Maximum SSTs are recorded as an early transient event for cases of slow subduction (100 km) which deflects the zone of maximum shear away from slab-wedge interface.¿ 1997 American Geophysical Union