The linear eigenmode characterizations and the nonlinear turbulence energy spreading of the drift waves in a tokamak plasma with strong pedestal gradient are numerically investigated based on an electromagnetic Landau fluid model. By the linear eigenmode analysis, it is found that the dominant instability in the low β regime is the ion-temperature-gradient (ITGc) mode and the electron drift wave instability (eDWIp) in the core and edge region with strong density gradient, respectively. Multiple eigenstates of the eDWIp with different peak locations in the poloidal direction can be obtained by the eigenvalue problem solver. The dominant one is the high order eDWIp corresponding to the unconventional ballooning mode structure with multiple peaks in the poloidal position, in contrast to the conventional modes that peak at the outboard mid-plane, and has been verified through initial value simulation. In the high β regime, the dominant eigenmodes in the core and edge region are the conventional and unconventional kinetic ballooning modes respectively. In the nonlinear simulation, an inward turbulence spreading phenomenon during the quasi-saturation phase of the edge turbulence is clearly observed. The inward speed of the turbulence energy front in the high β regime is much faster than that in the low β regime. It is interestingly found that the speed of the turbulence energy front increases with the increase of the plasma β in the low β regime, while it is almost unchanged in the high β regime. It is identified that the turbulence spreading in the low and high β regimes are determined by the nonlinear dynamics and the linear toroidal coupling respectively.