We present a survey of macroscopic excitations of harmonically confined Bose-Einstein condensates (BEC), described by Gross-Pitaevskii (GP) equation, in search of routes to develop quantum turbulence. These excitations can all be created by phase imprinting techniques on an otherwise equilibrium BoseEinstein condensate. We analyze two crossed vortices, two parallel anti-vortices, a vortex ring, a vortex with topological charge Q = 2, and a tangle of 4 vortices. Since GP equation is time-reversal invariant, we are careful to distinguish time intervals in which this symmetry is preserved and those in which rounding errors play a role. We find that the system tends to reach stationary states that may be widely classified as having either an array of vortices with collective excitations at different length scales or an agitated state composed mainly of Bogoliubov phonons.