High-resolution direct numerical simulations (DNSs) of incompressible turbulence based on an alias-free spectral method were performed on the Earth Simulator. Statistics of turbulence are studied by a DNS on 1024(3) grid points with a special emphasis on the spectra of moments fourth order in velocity. A brief review is given on some results of the preliminary analysis of the data of DNSs with up to 2048(3) grid points.