Impurity accumulation in the core plasma leads to fuel dilution and higher radiative losses that can lead to loss of the H-mode, to thermal collapse of the plasma, and eventually even to a disruption in tokamaks. In present experiments, it has been shown that Edge Localized Modes (ELMs) at sufficiently high frequency are required to prevent W accumulation in the core, by expelling impurities from the edge plasma region, thus preventing their penetration into the plasma core. We present a full-orbit particle extension of the MHD code JOREK suitable for simulating impurity transport during ELMs. This model has been applied to the simulation of an ELM crash in ASDEX Upgrade, where we have quantified the displacement of W particles across flux surfaces. The transport mechanism is shown to be the particle E × B-drifts due to the electric field created by the MHD instability underlying the ELM. In- and outward transport is observed as a series of interchange motions, leading to a superdiffusive behavior. This causes not only the particles near the plasma pedestal to move outwards but also the particles outside of the pedestal to move inwards. This has important consequences for operation with W in ITER, where it is expected to be screened in the pedestal, and ELMs are shown here to increase the core W density. A comparison with existing diffusive modeling is made, showing a qualitative agreement and the limitations of this simplified modeling approach.