Long memory processes are widely used in many scientific fields, such as economics, physics, and engineering. Change point detection problems have received considerable attention in the literature because of their wide range of possible applications. Here we describe a wavelet-based Bayesian procedure for the estimation and location of multiple change points in the long memory parameter of Gaussian autoregressive fractionally integrated moving average models (ARFIMA(p, d, q)), with unknown autoregressive and moving average parameters. Our methodology allows the number of change points to be unknown. The reversible jump Markov chain Monte Carlo algorithm is used for posterior inference. The method also produces estimates of all model parameters. Performances are evaluated on simulated data and on the benchmark Nile river dataset.