This study aims to systematically investigate the effects of spectral scaling and matching procedures of earthquake records on 1D, 2D, and 3D Ground Response Analyses (GRA). In this study, analyses were performed on the same soil profile using finite element software such as Deepsoil, PLAXIS, and FLAC. Selected earthquake records were scaled using a simple scaling method based on a structural period range of 0.5-1s and were also spectrally matched; the effects of these procedures on ground behavior were evaluated in terms of surface spectral accelerations. In light of the findings, the influence of the scaling factor and model dimensionality (1D, 2D, 3D) on the obtained spectral acceleration values was investigated. While it was determined that model dimensionality produces different performance outputs even under same scaling conditions, the application of different scaling ratios to the same model dimension also led to significant changes in the analysis results. These findings indicate that both model dimensionality and scaling factor have independent and cumulative effects on spectral acceleration, revealing that the optimization of these parameters is of critical importance for analysis accuracy.