We consider the evolution of multi-pulse patterns in an extended Klausmeier equation -as generalisation of the wellstudied Gray-Scott system, a prototypical singularly perturbed reaction-diffusion equation -with parameters that change in time and/or space. As a first step we formally show -under certain conditions on the parameters -that the full PDE dynamics of a N-pulse configuration can be reduced to a N-dimensional dynamical system that describes the dynamics on a N-dimensional manifold M N . Next, we determine the local stability of M N via the quasi-steady spectrum associated to evolving N-pulse patterns. This analysis provides explicit information on the boundary ∂M N of M N . Following the dynamics on M N , a N-pulse pattern may move through ∂M N and 'fall off' M N . A direct nonlinear extrapolation of our linear analysis predicts the subsequent fast PDE dynamics as the pattern 'jumps' to another invariant manifold M M , and specifically predicts the number N − M of pulses that disappear during this jump. Combining the asymptotic analysis with numerical simulations of the dynamics on the various invariant manifolds yields a hybrid asymptotic-numerical method that describes the full pulse interaction process that starts with a N-pulse pattern and typically ends in the trivial homogeneous state without pulses. We extensively test this method against PDE simulations and deduce a number of general conjectures on the nature of pulse interactions with disappearing pulses. We especially consider the differences between the evolution of (sufficiently) irregular and regular patterns. In the former case, the disappearing process is gradual: irregular patterns loose their pulses one by one, jumping from manifold M k to M k−1 (k = N, . . . , 1). In contrast, regular, spatially periodic, patterns undergo catastrophic transitions in which either half or all pulses disappear (on a fast time scale) as the patterns pass through ∂M N . However, making a precise distinction between these two drastically different processes is quite subtle, since irregular N-pulse patterns that do not cross ∂M N typically evolve towards regularity. sometimes also called the generalised Klausmeier-Gray-Scott system [34,40]. This model is a generalization of the original ecological model by Klausmeier on the interplay between vegetation and water in semi-arid regions [23] -which was proposed to describe the appearance of vegetation patterns as crucial intermediate step in the desertification process that begins with a homogeneously vegetated terrain and ends with the non-vegetated bare soil state: the desert -see [8,28,32] and the references therein for observations of these patterns and their relevance for the desertification process. In (1.1), U(x, t) represents (the concentration of) water and V(x, t) vegetation; for simplicity -and as in [34,35,36,38, 40] -we consider the system in a 1-dimensional unbounded domain, i.e. x ∈ R; parameter a models the rainfall and m the mortality of the vegetation. Since the diffusion of water occurs on a mu...