In the early stages of aerospace design, reduced-order models (ROMs) are crucial for minimizing computational costs associated with using physics-rich field information in many-query scenarios requiring multiple evaluations. The intricacy of aerospace design demands the use of high-dimensional design spaces to capture detailed features and design variability accurately. However, these spaces introduce significant challenges, including the curse of dimensionality, which stems from both high-dimensional inputs and outputs necessitating substantial training data and computational effort. To address these complexities, this study introduces a novel multifidelity, parametric, and nonintrusive ROM framework designed for high-dimensional contexts. It integrates machine learning techniques for manifold alignment and dimension reduction—employing proper orthogonal decomposition and model-based active subspace—with multifidelity regression for ROM construction. Our approach is validated through two test cases: the 2D RAE 2822 airfoil and the 3D NASA CRM wing, assessing various fidelity levels, training data ratios, and sample sizes. Compared to the single-fidelity principal component–active subspace (PCAS) method, our multifidelity solution offers improved cost-accuracy benefits and achieves better predictive accuracy with reduced computational demands. Moreover, our methodology outperforms the manifold-aligned ROM method by 50% in handling scenarios with large input dimensions, underscoring its efficacy in addressing the complex challenges of aerospace design.