Abstract-In this work, a new method has been proposed for the finite-difference time-domain (FDTD) analysis of the transient grounding resistance (TGR) of large grounding systems. To calculate the TGR, a coarse grid is occupied to model the earthing conductor, the convolution PML (CPML) is chose to truncate the computational domain, and the parallel implementation is involved to overcome the memory limit of the serial FDTD. With this model, the effect of the earthing conductor number and topology structure, the buried depth, and the ground permittivity and conductivity on the TGR is tested to find an optimized program to decrease the TGR of the lightning protection grounding systems.