Conditional ground‐motion models (CGMMs) for estimating Arias intensity () for earthquakes in subduction zones are developed. The estimate of is conditioned in these models on the estimated peak ground acceleration (PGA), the spectral acceleration at (SA1), time‐averaged shear‐wave velocity in the top 30 m (), and magnitude (). Random‐effects regressions are used to develop CGMMs for Japan, Taiwan, South America, and New Zealand. By combining the conditional models of with the ground‐motion models (GMMs) for PGA and SA1, the conditional models are converted to scenario‐based GMMs that can be used to estimate the median and its standard deviation directly for a given earthquake scenario and site conditions. The conditional scaling approach ensures the estimated values are consistent with a design spectrum that may correspond to above‐average spectral values for the controlling scenario. In addition, this approach captures the complex ground‐motion scaling effects found in GMMs for spectral acceleration, such as sediment‐depth effects, soil nonlinearity effects, and regionalization effects, in the developed scenario‐based models for . Estimates from the new scenario‐based models are compared to those from traditional GMMs for in subduction zones.