A central-moments-based lattice Boltzmann model for large-eddy simulation of neutrallystratified turbulent flows is described. Through comparative simulations of the airflow within and above a homogeneous plant canopy, the performance of the model is evaluated with respect to a conventional large-eddy-simulation model based on the incompressible Navier-Stokes equations. Simulated turbulence statistics, such as the mean velocity, velocity variances, velocity skewness, and power spectra, are shown to be almost identical between the two models. The spatial structure of coherent eddies and their maintenance processes are also confirmed to be properly represented by the lattice Boltzmann method through analysis of the turbulence kinetic energy budget and spatial two-point correlation functions. Using the simulated results, the energetics of the streamwise-elongated streaky structures commonly observed over vegetation and urban canopies are examined. While the short-wavelength components of the shear-generated streamwise kinetic energy are redirected rapidly by pressure to the lateral and vertical velocity components, long-wavelength energy tends to remain in the streamwise velocity component, which is dissipated in relatively slower processes. Consequently, the streaky structures persist in the streamwise velocity component.