We present a new iterative model order reduction method for large-scale linear time-invariant dynamical systems, based on a combined singular value decomposition–adaptive-order rational Arnoldi (SVD-AORA) approach. This method is an extension of the SVD-rational Krylov method. It is based on two-sided projections: the SVD side depends on the observability Gramian by the resolution of the Lyapunov equation, and the Krylov side is generated by the adaptive-order rational Arnoldi based on moment matching. The use of the SVD provides stability for the reduced system, and the use of the AORA method provides numerical efficiency and a relative lower computation complexity. The reduced model obtained is asymptotically stable and minimizes the error ($H_{2}$ and $H_{\infty }$) between the original and the reduced system. Two examples are given to study the performance of the proposed approach.