It is generally accepted that tokamak edge plasma turbulence causes anomalous diffusion. Theoretical studies addressing the anomalous diffusion in these potential structures are usually based on the test-particle drift approximation and on the electrostatic field resulting from the Hasegawa-Mima model (see, e.g. [1]) or the Hasegawa-Wakatani model [2].In our preceding papers (e.g.[3]), we discussed the effect of the anomalous ion diffusion in question on radial electric field generation. There we used a very simplified model of the turbulent potential structures, namely a spatially periodic and time-independent potential. Using a Hamiltonian approach (which takes into account the cyclotron motion), we have found, for impurity ions C + and usual potential amplitudes, a substantial increase in the diffusion of these ions (both of the Gaussian and Lévy-walk forms) [4], resulting in the generation of a radial electric field [3]. Using a drift approximation for this case, no diffusion and no electric field is observed. Radial electric fields play an important role in the formation of transport barriers in tokamaks. It seems that our mechanism, resulting from anomalous diffusion of impurities can play an interesting role.Recently, a more realistic form of this potential, namely the Hasegawa-Wakatani potential (HW), was considered. The HW model [5] is based on the numerical solution of mode-coupling equations for the resistive drift-wave instability, modelling in this way the turbulent processes in the edge plasma. We have followed their approach, using the parameters of the Czech CASTOR tokamak [6] for the numerical simulation. A typical instantaneous potential relief inside a space sample with dimensions 4 cm radially 8 cm poloidally is presented in Fig. 1. We use the 2-dimensional Particle-in-Cell (PIC) code BIT2 to follow C + ions. These particles are initially distributed uniformly over the whole sample. We assume periodic boundary conditions in the poloidal direction. In the radial direction, the boundary conditions are periodic for the calculation of the variance and diffusion coefficient and non-periodic for the self-consistent calculation of the radial electric field generation. The electrons and plasma ions are treated as a fixed background and the HW potential is taken from a fluid simulation.We calculate the time histories of the average variance X 2 (t) and the running diffusion coefficient(and similarly for the y direction), where x j , y j are the radial and poloidal co-ordinates of the jth particle, respectively, t is time, and N is the number of particles. In Fig. 3 we see that the diffusion is larger in the poloidal direction than in the radial. Contrary to the results presented in [4], the diffusion in the radial direction is rather Gaussian, and only the diffusion in the poloidal direction has a parabolic-type variance time trace (Fig. 2) indicating Lévy-walk dynamics. As already proposed in our recent papers ([3], [7]), this anomalous diffusion must necessarily results in a positive-charge deficit....