Ongoing developments in earthquake source inversions incorporate nonplanar fault geometries as inputs to the inversion process, improving previous approaches that relied solely on planar fault surfaces. This evolution motivates advancing the existing framework for constraining fault geometry, particularly in subduction zones where plate boundary surfaces that host highly hazardous earthquakes are clearly nonplanar. Here, we improve upon the existing framework for the constraint of the seismic rupture plane of subduction interfaces by incorporating active seismic and seafloor sediment thickness data with existing independent data sets and inverting for the most probable nonplanar subduction geometry. Constraining the rupture interface a priori with independent geological and seismological information reduces the uncertainty in the derived earthquake source inversion parameters over models that rely on simpler assumptions, such as the moment tensor inferred fault plane. Examples are shown for a number of wellconstrained global locations. We expand the coverage of previous analyses to a more uniform global data set and show that even in areas of sparse data this approach is able to accurately constrain the approximate subduction geometry, particularly when aided with the addition of data from local active seismic surveys. In addition, we show an example of the integration of many two-dimensional profiles into a threedimensional surface for the Sunda subduction zone and introduce the development of a new global threedimensional subduction interface model: Slab1.0. ?? 2009 by the American Geophysical Union.