  We calculate the D-wave charmonium production in $e^+e^-$ annihilation at
BELLE and BABAR at $\sqrt{s}=10.6$ GeV using the nonrelativistic QCD
factorization formalism, including color singlet and color octet contributions.
We analyze the contributions of various processes for $\delta_J(J=1,2)$
production. We find that both the color singlet and color octet channels may
give substantial contributions, and the production rates are estimated to be
$\sigma (\delta_1)\simeq 0.043-0.16pb$ and $\sigma (\delta_2)\simeq
0.094-0.29pb$, which, however, are very sensitive to the choice of the color
octet matrix elements of the D-wave charmonium states. The measurement of
D-wave charmonium production at BABAR and BELLE in the future will be very
helpful to test the color octet mechanism and to determine the color octet
matrix elements.
