  (Abridged) By using an N-body simulation of a bulge that was formed via a bar
instability mechanism, we analyse the imprints of the initial (i.e. before bar
formation) location of stars on the bulge kinematics, in particular on the
heliocentric radial velocity distribution of bulge stars. Four different
latitudes were considered: $b=-4^\circ$, $-6^\circ$, $-8^\circ$, and
$-10^\circ$, along the bulge minor axis as well as outside it, at
$l=\pm5^\circ$ and $l=\pm10^\circ$. The bulge X-shaped structure comprises
stars that formed in the disk at different locations. Stars formed in the outer
disk, beyond the end of the bar, which are part of the boxy peanut-bulge
structure may show peaks in the velocity distributions at positive and negative
heliocentric radial velocities with high absolute values that can be larger
than 100 $\rm km$ $\rm s^{-1}$, depending on the observed direction. In some
cases the structure of the velocity field is more complex and several peaks are
observed. Stars formed in the inner disk, the most numerous, contribute
predominantly to the X-shaped structure and present different kinematic
characteristics. Our results may enable us to interpret the cold high-velocity
peak observed in the APOGEE commissioning data, as well as the excess of
high-velocity stars in the near and far arms of the X-shaped structure at
$l$=$0^\circ$ and $b$=$-6^\circ$. When compared with real data, the kinematic
picture becomes more complex due to the possible presence in the observed
samples of classical bulge and/or thick disk stars. Overall, our results point
to the existence of complex patterns and structures in the bulge velocity
fields, which are generated by the bar. This suggests that caution should be
used when interpreting the bulge kinematics: the presence of substructures,
peaks and clumps in the velocity fields is not necessarily a sign of past
accretion events.
