With the ongoing progress of computing power made available not only by large supercomputer facilities but also by relatively common workstations and desktops, physics-based source-to-site 3D numerical simulations of seismic ground motion will likely become the leading and most reliable tool to construct ground shaking scenarios from future earthquakes. This paper aims at providing an overview of recent progress on this subject, by taking advantage of the experience gained during a recent research contract between Politecnico di Milano, Italy, and Munich RE, Germany, with the objective to construct ground shaking scenarios from hypothetical earthquakes in large urban areas worldwide. Within this contract, the SPEED computer code was developed, based on a spectral element formulation enhanced by the Discontinuous Galerkin approach to treat non-conforming meshes. After illustrating the SPEED code, different case studies are overviewed, while the construction of shaking scenarios in the Po river Plain, Italy, is considered in more detail. Referring, in fact, to this case study, the comparison with strong motion records allows one to derive some interesting considerations on the pros and on the present limitations of such approach.