The current paper is concerned with the forced waves of Keller-Segel chemoattraction systems in shifting environments of the form, u t = u xx − χ(uv x) x + u(r(x − ct) − bu), x ∈ R 0 = v xx − νv + µu, x ∈ R, (0.1) where χ, b, ν, and µ are positive constants, c ∈ R, the resource function r(x) is globally Hölder continuous, bounded, r * = sup x∈R r(x) > 0, r(±∞) := lim x→±∞ r(x) exist, and either r(−∞) < 0 < r(∞), or r(±∞) < 0. Assume that b > 2χµ. In the case that r(−∞) < 0 < r(∞), it is shown that (0.1) has a forced wave solution connecting (r * b , µ ν r * b) and (0, 0) with speed c provided that c > χµr * 2 √ ν(b−χµ) − 2 r * (b−2χµ) b−χµ. In the case that r(±∞) < 0, it is shown that (0.1) has a forced wave solution connecting (0, 0) and (0, 0) with speed c provided that χ is sufficiently small and λ ∞ > 0, where λ ∞ is the generalized principal eigenvalue of the operator u(•) → u xx (•) + cu x (•) + r(•)u(•) on R in certain sense. Some numerical simulations are also carried out. The simulations indicate the existence of forced wave solutions in some parameter regions which are not covered in the theoretical results, induce several problems to be further studied, and also provide some illustration of the theoretical results.