The damage to structures during an earthquake can be varied depending on the frequency characteristics of seismic waves and the geological properties of the ground. Therefore, considering such attributes in the design ground motions is crucial. The Korean seismic design standard (KDS 17 10 00) provides design response spectra for various ground classifications. If required for time-domain analysis, ground motion time series can be either selected and adjusted from motions recorded at rock sites in intraplate regions or artificially synthesized. Ground motion time series at soil sites should be obtained from site response analysis. However, in practice, selecting suitable ground motion records is challenging due to the overall lack of large earthquakes in intraplate regions, and artificially synthesized time series often leads to unrealistic responses of structures. As an alternative approach, this study provides a case study of generating ground motion time series based on the hybrid broadband ground motion simulation of selected scenario earthquakes at sites in the Nakdonggang delta region. This research is significant as it provides a novel method for generating ground motion time series that can be used in seismic design and response analysis. For large-magnitude earthquake scenarios close to the epicenter, the simulated response spectra surpassed the 1000-year design response spectra in some specific frequency ranges. Subsequently, the acceleration time series at each location were used as input motions to perform nonlinear 1D site response analysis through the PySeismoSoil Package to account for the site response characteristics at each location. The results of the study revealed a tendency to amplify ground motion in the mid to long-period range in most places within the study area. Additionally, significant amplification in the short-period range was observed in some locations characterized by a thin soil layer and relatively high shear wave velocity soil near the upper bedrock.