The proposed work integrates the Arnoldi algorithm within the Krylov Subspace Method (KSM) framework to construct the reduced-order model (ROM) of a commensurate-order fractional-order system. First, the original commensurate fractional-order model (FOM) is transformed into an equivalent integer-order system (IOS) through a suitable domain transformation. This conversion enables the use of well-established linear system reduction techniques. On the transformed IOS, a Krylov subspace is generated to approximate the system’s dominant dynamics. The Arnoldi algorithm is then employed to construct an orthonormal basis for this Krylov subspace. By iteratively building a Hessenberg matrix representation of the large system matrix, Arnoldi efficiently captures the most influential spectral properties while preserving numerical stability. Thus, Arnoldi complements KSM by ensuring stable and computationally efficient basis construction, while KSM provides the theoretical foundation for projection-based model reduction. After reduction, the IOS-based ROM is transformed back into fractional-order form using the inverse transformation, yielding a compact fractional-order ROM. The main contribution of the work is to propose an Arnoldi-enhanced Krylov subspace framework tailored for commensurate-order fractional-order systems via domain transformation. The method achieves improved approximation accuracy in both time and frequency domains, as validated through performance indices such as Integral Absolute Error (IAE), Integral of Time Weighted Absolute Error (ITAE), Integral Time Square Error (ITSE), Integral Square Error (ISE), and \(\mathcal {H}_{\infty }\) norm. Comparative studies demonstrate that the proposed Arnoldi–KSM approach produces reduced models that closely replicate the dynamics of high-order fractional systems, outperforming existing ROM techniques in terms of accuracy and compactness. To test the effectiveness of the proposed method, two examples of higher order FOM are simulated in MATLAB.