We present a three dimensional velocity analysis of Milky Way disk kinematics using LAMOST K giant stars and the GPS1 proper motion catalogue . We find that Galactic disk stars near the anticenter direction ( in the range of Galactocentric distance between R = 8 and 13 kpc and vertical position between Z = - 2 and 2 kpc ) exhibit asymmetrical motions in the Galactocentric radial , azimuthal , and vertical components . Radial motions are not zero , thus departing from circularity in the orbits ; they increase outwards within R \lesssim 12 kpc , show some oscillation in the northern ( 0 < Z < 2 kpc ) stars , and have north-south asymmetry in the region corresponding to a well-known nearby northern structure in the velocity field . There is a clear vertical gradient in azimuthal velocity , and also an asymmetry that shifts from a larger azimuthal velocity above the plane near the solar radius to faster rotation below the plane at radii of 11-12 kpc . Stars both above and below the plane at R \gtrsim 9 kpc exhibit net upward vertical motions . We discuss some possible mechanisms that might create the asymmetrical motions , such as external perturbations due to dwarf galaxy minor mergers or dark matter sub-halos , warp dynamics , internal processes due to spiral arms or the Galactic bar , and ( most likely ) a combination of some or all of these components .