The (low frequency) zonal flow-ion temperature gradient (ITG) wave system, constructed on Braginskii's fluid model in tokamak, is shown to be a reaction-diffusion-advection system; it is derived by making use of a multiple spatiotemporal scale technique and two-dimensional (2D) ballooning theory. For real regular group velocities of ITG waves, two distinct temporal processes, sharing a very similar meso-scale radial structure, are identified in the nonlinear self-organized stage. The stationary and quasi-stationary structures reflect a particular feature of the poloidal group velocity. The equation set posed to be an initial value problem is numerically solved for JET low mode parameters; the results are presented in several figures and two movies that show the spatiotemporal evolutions as well as the spectrum analysis-frequency-wave number spectrum, auto power spectrum, and Lissajous diagram. This approach reveals that the zonal flow in tokamak is a local traveling wave. For the quasi-stationary process, the cycle of ITG wave energy is composed of two consecutive phases in distinct spatiotemporal structures: a pair of Cavitons growing and breathing slowly without long range propagation, followed by a sudden decay into many Instantons that carry negative wave energy rapidly into infinity. A spotlight onto the motion of Instantons for a given radial position reproduces a Blob-Hole temporal structure; the occurrence as well as the rapid decay of Caviton into Instantons is triggered by zero-crossing of radial group velocity. A sample of the radial profile of zonal flow contributed from 31 nonlinearly coupled rational surfaces near plasma edge is found to be very similar to that observed in the JET Ohmic phase [J. C. Hillesheim et al., Phys. Rev. Lett. 116, 165002 (2016)]. The theory predicts an interior asymmetric dipole structure associated with the zonal flow that is driven by the gradients of ITG turbulence intensity. Published by AIP Publishing.