SIMBAD references

We aim at analysis of kinematics of main-sequence stars from the GaiaDR2 and PMA catalogues as well as comparison of kinematic parameters derived from their proper motions. We decompose the stellar velocity field on to a set of vector spherical harmonics (VSH), and, using the relations between the decomposition coefficients and the Ogorodnikov-Milne (O-M) model parameters, calculate the latter. The method of VSH allowed to detect all systematic constituents present in the GaiaDR2 and PMA stellar velocity fields. We notice incompleteness of the O-M model, discuss the kinematic parameters derived within its framework, as well as significant decomposition coefficients that do not have analogues within this model. For separate analysis of the kinematic parameters in the Northern and Southern Galactic hemispheres, we apply the decomposition on to a set of zonal VSH (ZVSH). Modelling the stellar velocity field allowed to confirm the causes of some significant beyond-the-model harmonics. Based on stellar proper motions from the Gaia DR2 and PMA catalogues, we confirm the previous conclusion that the values of ω₁ and M^+₂₃ O-M model parameters derived in the Northern and Southern Galactic hemispheres have opposite signs. This fact takes place due to a vertical gradient of the linear rotation velocity about the Galactic centre. We estimate influence of the vertical gradient on the value of solar velocity component Y_☉. We confirm that the beyond-the-model coefficients t₂₁₁ and s₃₁₀ are caused by the vertical gradient as well.