A modified lattice Boltzmann model based on the two-dimensional, nine-velocity lattice-Bhatnagar-Gross-Krook fluid is presented for axisymmetric flows. A spatially and temporally varying source term is incorporated into the evolution equation for the momentum distribution function on a two-dimensional Cartesian lattice. The precise form of the source term is derived through a Chapman-Enskog analysis, so that the additional axisymmetric contributions in the Navier-Stokes equations are furnished when written in the cylindrical polar coordinate system.