We present a theoretical study of a fundamentally new wave–mean or wave–vortex interaction effect able to force persistent, cumulative change in mean flows in the absence of wave breaking or other kinds of wave dissipation. It is associated with the refraction of non-dissipating waves by inhomogeneous mean (vortical) flows. The effect is studied in detail in the simplest relevant model, the two-dimensional compressible flow equations with a generic polytropic equation of state. This includes the usual shallow-water equations as a special case. The refraction of a narrow, slowly varying wavetrain of small-amplitude gravity or sound waves obliquely incident on a single weak (low Froude or Mach number) vortex is studied in detail. It is shown that, concomitant with the changes in the waves' pseudomomentum due to the refraction, there is an equal and opposite recoil force that is felt, in effect, by the vortex core. This effective force is called a ‘remote recoil’ to stress that there is no need for the vortex core and wavetrain to overlap in physical space. There is an accompanying ‘far-field recoil’ that is still more remote, as in classical vortex-impulse problems. The remote-recoil effects are studied perturbatively using the wave amplitude and vortex weakness as small parameters. The nature of the remote recoil is demonstrated in various set-ups with wavetrains of finite or infinite length. The effective recoil force ${\bm R}_V$ on the vortex core is given by an expression resembling the classical Magnus force felt by moving cylinders with circulation. In the case of wavetrains of infinite length, an explicit formula for the scattering angle $\theta_*$ of waves passing a vortex at a distance is derived correct to second order in Froude or Mach number. To this order ${\bm R}_V\,{\propto}\,\theta_*$. The formula is cross-checked against numerical integrations of the ray-tracing equations. This work is part of an ongoing study of internal-gravity-wave dynamics in the atmosphere and may be important for the development of future gravity-wave parametrization schemes in numerical models of the global atmospheric circulation. At present, all such schemes neglect remote-recoil effects caused by horizontally inhomogeneous mean flows. Taking these effects into account should make the parametrization schemes significantly more accurate.