3D-1D broadband simulations of ground motion are normally performed by using a hybrid approach where full wave numerical methods are used for the low frequencies and a stochastic contribution has added a posteriori to achieve high frequencies. The cross-over for hybridization is established empirically in a band between about 0.5 Hz and 2 Hz depending on the earthquake magnitude size, distance, local structure and shape of Fourier spectra of low frequency synthetic seismograms.