Earthquake source time functions deconvolved from teleseismic broadband P wave recordings are used to examine rupture variations for 417 underthrusting earthquakes located on the interplate interface in circum-Pacific subduction zones. Moment-scaled duration of significant moment release varies with depth, with longer-duration events occurring in the shallowest 20 km of the megathrusts. The source time functions are also used to estimate radiated seismic energy. Two estimates are obtained: a simple scaled triangle substitution for the moment release history provides a minimum estimate, while integration of the time function shape provides an estimate limited only by the bandwidth of the teleseismic deconvolutions. While these energy calculations underestimate total energy, they enable systematic comparisons of rupture process within and between subduction zones. We do not find significant depth dependence for radiated energy overall, but some regions do show mild trends of increasing energy/seismic moment ratios (E/Mo) with increasing source depth that correspond to rupture duration variations in those regions. The observations of longer rupture duration at shallow depth with moderate E/Mo may be due to heterogeneous friction and structural features on the shallow plate interface.