The process of ground motion selection and scaling is an integral part of hazardand risk-consistent seismic demand analysis of structures. Due to the lack of ground motion records that naturally possess high amplitude and intensity, the research community generally relies on scaling the records to match a target hazard intensity level. The scaling factors used are frequently as high as 10. Due to the criticism received in previous research studies, the extent of amplitude scaling and its process has become a matter of debate, and various constraints on the scaling factors have been proposed. The primary argument against unrestricted amplitude scaling is the unrealistic nature of the scaled records and the possible biases caused in the engineering demand parameters (EDPs) of structures. This study presents a framework to utilize machine-learning and statistical techniques for the assessment of ground motion amplitude scaling for nonlinear time-history analysis (NTHA) of structures. The framework utilizes Bayesian non-parametric Gaussian process regressions (GPRs) as surrogate models to obtain statistical estimates of EDPs for scaled and unscaled ground motions. The GPR surrogate models are developed based on a large-scale analysis of five steel moment frames (SMFs) using 200 unscaled as-recorded ground motions for ten spectral acceleration levels, (𝑆 𝑎 (𝑇 1 )) (ranging from 0 g to maximum considered earthquake, MCE) and 2500 scaled ground motions representing 50 scale factors (𝑆𝐹), and the 10 𝑆 𝑎 (𝑇 1 ) levels for each SMF. For each building, two types of EDPs are considered: i) peak inter-story drift ratio (PIDR) and ii) peak floor acceleration (PFA). To provide a better interpretation of the GPR surrogate models, the concept of explainable artificial intelligence (i.e., Shapley additive explanation, SHAP) is used to obtain insights into the decision-making process of the GPR models with respect to the 𝑆𝐹 and 𝑆 𝑎 (𝑇 1 ). Then, for the 10 𝑆 𝑎 (𝑇 1 ) levels, the GPRbased EDP estimates under scaled ground motions corresponding to 50 different SFs are compared with the EDP estimates of unscaled ground motions. TheThis is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.