  We present a detailed numerical study of solutions to the Zakharov-Kuznetsov
equation in three spatial dimensions. The equation is a three-dimensional
generalization of the Korteweg-de Vries equation, though, not completely
integrable. This equation is $L^2$-subcritical, and thus, solutions exist
globally, for example, in the $H^1$ energy space.
  We first study stability of solitons with various perturbations in sizes and
symmetry, and show asymptotic stability and formation of radiation, confirming
the asymptotic stability result in \cite{FHRY2020} for a larger class of
initial data. We then investigate the solution behavior for different
localizations and rates of decay including exponential and algebraic decays,
and give positive confirmation toward the soliton resolution conjecture in this
equation. Finally, we investigate soliton interactions in various settings and
show that there is both a quasi-elastic interaction and a strong interaction
when two solitons merge into one, in all cases always emitting radiation in the
conic-type region of the negative $x$-direction.
