One of the goals of this study was to simulate earthquake shaking in the wider Dubrovnik area. To achieve this, an approximate 3D local crustal structure model, DF_3D, was first created. The DF_3D model comprises three main layers: the upper crust, lower crust, and mantle, and includes topography and bathymetry. Each element within the model contains information on P- and S-wave velocities and material densities. This newly defined model was then used and compared with a local 1D model and a regional 3D model in deterministic simulations of low-frequency ground motion recordings. Through visual inspection and analysis of the quality metrics—such as peak ground velocities, correlation between simulated and recorded ground motions, and the duration of strong shaking—the local 3D crustal model was determined to be the most suitable model currently available for calculating low-frequency simulations in the wider Dubrovnik area. Using a hybrid method, broadband ground motion recordings were subsequently calculated for ten moderate to strong earthquakes. These recordings were generated by merging separately computed low- and high-frequency contributions using deterministic and stochastic approaches. In calculating the high-frequency contributions, a base of Green’s functions was defined and calculated for a 1D model (DF_1D), derived from the DF_3D model. Through visual inspection, quality metric analysis of peak ground velocities and shaking duration, and time-frequency analysis of envelope and phase deviations in simulated and recorded recordings, the effectiveness of the hybrid approach for broadband ground motion simulation in the studied area remains unclear. Further detailed validation and parameterization are required to fully assess its applicability. Regardless of this outcome, low-frequency simulations of the Great Dubrovnik Earthquake of 1667 were conducted to enhance understanding of ground motion during strong earthquakes in this area. Six earthquake shaking scenarios were simulated for this earthquake, using six different kinematic finite-source models. For these scenarios, a database of simulated low-frequency seismograms and simulation videos was created, and their main characteristics were analyzed. Scenario validation was based on a comparison of simulated ground motion intensity metrics, including horizontal peak ground velocity, spectral acceleration for periods of T = 1.0 s, and Arias intensity, with observed macroseismic intensities. It was found that a combination of non-uniform slip distribution along the fault plane and bilateral temporal rupture symmetry yields the highest correlation coefficient, potentially providing the best description of the seismic source for the Great Dubrovnik Earthquake. Further modeling of the seismic source, crustal structure model, simulation methods, and validation is necessary to ensure that earthquake shaking scenarios accurately reflect expected ground motion characteristics, thereby contributing to an improved deterministic seismic hazard assessment in the wider Dubrovnik area.