Two-dimensional electromagnetic particle-in-cell simulations in a magnetized, homogeneous, collisionless electron-proton plasma demonstrate the forward cascade of whistler turbulence. The simulations represent decaying turbulence, in which an initial, narrowband spectrum of fluctuations at wavenumbers kc∕ωe≃0.1 cascades toward increased damping at kc∕ωe≃1.0, where c∕ωe is the electron inertial length. The turbulence displays magnetic energy spectra that are relatively steep functions of wavenumber and are anisotropic with more energy in directions relatively perpendicular to the background magnetic field Bo=x̂Bo than at the same wavenumbers parallel to Bo. In the weak turbulence regime, the primary new results of the simulations are as follows: (1) Magnetic spectra of the cascading fluctuations become more anisotropic with increasing fluctuation energy; (2) the wavevector dependence of the three magnetic energy ratios, ∣δBj∣2∕∣δB∣2 with j=x,y,z, show good agreement with linear dispersion theory for whistler fluctuations; (3) the magnetic compressibility summed over the cascading modes satisfies 0.3≲∣δBx∣2∕∣δB∣2≲0.6; and (4) the turbulence heats electrons in directions both parallel and perpendicular to Bo, with stronger heating in the parallel direction.