This paper reports on a numerical study of two-dimensional decaying turbulence in a square domain with no-slip walls. The generation of strong small-scale vortices near the no-slip walls have been observed in the lattice Boltzmann simulations just like in earlier pseudospectral calculations. Due to these vortices the enstrophy is not a monotone decaying function of time. Considering a number of simulations and taking their ensemble average, we have found that the decay of enstrophy and that of the kinetic energy can be described well by power-laws. The exponents of these laws depend on the Reynolds number in a similar manner than was observed before in pseudospectral simulations. Considering the ensemble averaged 1D Fourier energy spectra calculated along the walls, we could not find a simple power-law, which fits well to the simulation data. These spectra change in time and reveal an exponent close to -3 in the intermediate and an exponent -5/3 at low wavenumbers. On the other hand, the two-dimensional energy spectra, which remain almost steady in the intermediate decay stage, show clear power-law behavior with exponent larger than -3 depending on the initial Reynolds number. © 2010 World Scientific Publishing Company.