The temperature-dependent edge magnetic susceptibility [Formula: see text] and the uniform magnetic susceptibility χ in zigzag graphene nanoribbons is studied within the Hubbard model on a honeycomb lattice. By using the determinant quantum Monte Carlo (DQMC) method, it is found that the ferromagnetic fluctuations at the zigzag edge dominate around half-filling, and that the fluctuations are strengthened markedly by the on-site Coulomb interaction U, which may lead to a possible high-temperature edge ferromagnetic behaviour in low-doped zigzag graphene nanoribbons.