We investigate the weakly non-linear evolution of cosmic gravitational clustering in phase space by looking at the Zeldovich solution in the discrete wavelet transform (DWT) representation. We show that if the initial perturbations are Gaussian, the relation between the evolved DWT mode and the initial perturbations in the weakly non-linear regime is quasi-local. That is, the evolved density perturbations are mainly determined by the initial perturbations localized in the same spatial range. Furthermore, we show that the evolved mode is monotonically related to the initial perturbed mode. Thus large (small) perturbed modes statistically correspond to the large (small) initial perturbed modes. We test this prediction by using QSO Ly$alpha$ absorption samples. The results show that the weakly non-linear features for both the transmitted flux and identified forest lines are quasi-localized. The locality and monotonic properties provide a solid basis for a DWT scale-by-scale Gaussianization reconstruction algorithm proposed by Feng & Fang (Feng & Fang, 2000) for data in the weakly non-linear regime. With the Zeldovich solution, we find also that the major non-Gaussianity caused by the weakly non-linear evolution is local scale-scale correlations. Therefore, to have a precise recovery of the initial Gaussian mass field, it is essential to remove the scale-scale correlations.