The NGA‐Sub (subduction zone earthquake) database developed by the Pacific Earthquake Engineering Research center is used to derive new correlation coefficients for a number of ground‐motion intensity measure (IM) parameters from ground motions in subduction zone earthquakes, considering both interface and intraslab tectonic settings. The IMs include peak ground acceleration, pseudospectral accelerations with periods from 0.01 to 10 s, Arias intensity, cumulative absolute velocity, peak ground velocity, and significant duration. Comparisons of the estimated correlation coefficients for ground motions from the interface and intraslab tectonic settings generally show a good agreement. Our estimations are also in good agreement with correlation coefficients estimated in previous studies that used ground motions from shallow crustal earthquakes, supporting the concept that any variation in correlation coefficients comes from spectral shape (i.e., the distribution of peaks and troughs) rather than tectonic region. This study also explores the influence of parameters such as magnitude, distance, and site conditions on the estimated correlation coefficients. We did not find apparent trends of the correlation coefficients with respect to these parameters. Finally, we propose analytical models to estimate correlation coefficients between the IMs explored in this study, considering both subduction interface and intraslab tectonic settings.