Abstract. Using three-dimensional magnetohydrodynamic simulations, we investigate the propagation of low-entropy magnetic flux tubes ("bubbles") in the magnetotail. Our simulations address fundamental properties of the propagation and dynamics of such flux tubes rather than the actual formation process. We find that the early evolution, after a sudden reduction of pressure and entropy on a localized flux tube, is governed by re-establishing the balance of the total pressure in the dawn-dusk and north-south directions through compression on a time scale less than about 20 s for the typical magnetotail. The compression returns the equatorial pressure to its original unperturbed value, due to the fact that the magnetic field contributes only little to the total pressure, while farther away from the equatorial plane the magnetic field compression dominates. As a consequence the pressure is no longer constant along a flux tube. The subsequent evolution is characterized by earthward propagation at speeds of the order of 200-400 km/s, depending on the initial amount of depletion and the cross-tail extent of a bubble. Simple acceleration without depletion does not lead to significant earthward propagation. It hence seems that both the entropy reduction and the plasma acceleration play an important role in the generation of fast plasma flows and their propagation into the near tail. Earthward moving bubbles are found to be associated with field-aligned current systems, directed earthward on the dawnward edge and tailward on the duskward edge. This is consistent with current systems attributed to observed bursty bulk flows and their auroral effects.