Spatiotemporal quenches are efficient at preparing ground states of critical Hamiltonians that have emergent low-energy descriptions with Lorentz invariance [Agarwal, Bhatt, and Sondhi, Phys. Rev. Lett. 120, 210604 (2018)]. The critical transverse field Ising model with nearest-neighbor interactions, for instance, maps to free fermions with a relativistic low energy dispersion. However, spin models realized in artificial quantum simulators based on neutral Rydberg atoms, or trapped ions, generically exhibit long-range power-law decay of interactions with for a wide range of . In this paper, we study the fate of spatiotemporal quenches in these models with a fixed velocity for the propagation of the quench front, using the numerical time-dependent variational principle. For , where the critical theory is suggested to have a dynamical critical exponent , our simulations show that optimal cooling is achieved when the front velocity approaches , the effective speed of excitations in the critical model. The energy density is inhomogeneously distributed in space, with prominent hot regions populated by excitations copropagating with the quench front, and cold regions populated by counterpropagating excitations. Lowering largely blurs the boundaries between these regions. For , we find that the Doppler cooling effect disappears, as expected from renormalization group results for the critical model, which suggest a dispersion with . Instead, we show that excitations are controlled by two relevant length scales whose ratio is related to that of the front velocity to a threshold velocity that ultimately determines the adiabaticity of the quench.