The Landau damping of geodesic acoustic mode (GAM) in a torodial rotating tokamak plasma is analytically investigated by taking into account the finite-orbit-width (FOW) resonance effect to the 3rd order. The analytical result is shown to agree well with the numerical solution. The dependence of the damping rate on the toroidal Mach number M relies on krρi. For sufficiently small krρi, the damping rate monotonically decreases with M . For relatively large krρi, the damping rate increases with M until approaching the maximum and then decreases with M .In the kinetic framework, keeping terms to the 1st order finite-Larmor-radius (FLR) effect, which represents the leading order polarization, and the 1st order finiteorbit-width (FOW) effect of passing particles δ i ∼ qρ i , where q is the safety factor of tokamaks and ρ i is the ion gyroradius, the dispersion relation of geodesic acoustic mode (GAM) [1] is derived. The classical Landau damping rate of GAM is found to be ∝ e −q 2 Ω 2 and independent of δ i , where Ω is the GAM frequency normalized by R/v T i with major radius R and ions thermal velocity v T i . Later, some theoretical analysis[2-4], numerical evaluation [5], and simulation [6,7] all indicate that the high-order FOW effect plays a key role in the collisionless damping of GAM, specifically in the large q region. The resonant damping rate is sensitive to and significantly enhanced by k r ρ i , where k r is the radial wave number. When only the 2-nd resonance is taken into account, a recent NEMORB simulation [8] performed a good agreement with the analytical result to the 2-nd FOW effects for about q < 3.5 [2]. It was shown that for q > 3.5, the discrepancy between the theoretical result with 2nd harmonics and the TEMPEST simulation data becomes remarkable [6]. That means higher-order resonance should be considered. The numerical evaluation performed by Gao et al. [5] shows that the damping rate with 3-rd resonance and the one with 4-th resonance have only slight discrepancy when q is about greater than 7. Meanwhile, Xu et al. numerically found that the damping rate with 4-th order resonance is almost the same as the rate with 10-th resonance [6].In a recent work[9], Guo and co-authors investigated the collisionless damping rate of GAM by taking into account the toroidal rotation. They numerically evaluated the influence of toroidal Mach number on the Landau damping rate by considering from 2nd to 5th FOW resonance effect, respectively. Similar to the case without toroidal rotation [5,6,10], they found that the damping rate was significantly enhanced by 3rd resonance and the damping rate with 3rd resonance and the one with 4th * Electronic address: hjren@ustc.edu.cn or 5th resonance are almost the same with each other. In this Letter, theoretical investigation on the collisionless damping of GAM in a toroidally rotating tokamak plasma is performed. We derive the analytical expression for the Landau damping rate by considering the FOW resonance effect to the 3rd order. Good agreement is found between our...