Electromagnetic effects of the kinetic geodesic acoustic modes (KGAMs) are numerically studied in low β(= plasma pressure/magnetic pressure) tokamak plasmas. The parallel component of the perturbed vector potential is considered along with the electrostatic potential perturbation. The finite Larmor radius and finite orbit width of the ions as well as electron parallel dynamics are all taken into account. Systematic harmonic and ordering analysis is performed for collisionless damping of the KGAMs, assuming β~(κρi)2, where κand ρiare the radial component of the KGAM wave vector and the Larmor radius of the ions, respectively. It is found that the electron parallel dynamics enhances the damping of the electrostatic KGAM modes when the safety factor q is high. In addition, the electromagnetic (finite β effect is revealed to enhance and weaken the damping of the modes in plasmas of low and high safety factor q~2.0 and 5.5, respectively. The harmonic features of the KGAMs are discussed as well.