AIMS. We calculate the stellar density using star counts obtained from Gaia DR2 up to a Galactocentric distance R=20 kpc with a deconvolution technique for the parallax errors. Then we analyse the density in order to study the structure of the outer Galactic disc, mainly the warp. METHODS. In order to carry out the deconvolution, we used the Lucy inversion technique for recovering the corrected star counts. We also used the Gaia luminosity function of stars with $M_G<10$ to extract the stellar density from the star counts. RESULTS. The stellar density maps can be fitted by an exponential disc in the radial direction $h_r=2.07\pm0.07$ kpc, with a weak dependence on the azimuth, extended up to 20 kpc without any cut-off. The flare and warp are clearly visible. The best fit of a symmetrical S-shaped warp gives $z_w= z_\odot+(37\pm 4.2(stat.)-0.91(syst.))$ pc $(R/R_\odot )^{2.42\pm 0.76(stat.) + 0.129 (syst.)} sin(\phi+9.3\pm 7.37 (stat.) +4.48 (syst.))$ for the whole population. When we analyse the northern and southern warps separately, we obtain an asymmetry of an $\sim25\%$ larger amplitude in the north. This result may be influenced by extinction because the Gaia G band is quite prone to extinction biases. However, we tested the accuracy of the extinction map we used, which shows that the extinction is determined very well in the outer disc. Nevertheless, we recall that we do not know the full extinction error, and neither do we know the systematic error of the map, which may influence the final result. The analysis was also carried out for very luminous stars alone ($M_G<-2$), which on average represents a younger population. We obtain similar scale-length values, while the maximum amplitude of the warp is $20-30\%$ larger than with the whole population. The north-south asymmetry is maintained.

Context: In our Paper I, by using statistical deconvolution methods, extended kinematics maps of Gaia-DR2 data have been produced in a range of heliocentric distances that are a factor of two to three larger than those analyzed previously by the Gaia Collaboration with the same data. It added the range of Galactocentric distances between 13 kpc and 20 kpc to the previous maps. Aims: Here, we investigate the dynamical effects produced by different mechanisms that can explain the radial and vertical components of these extended kinematic maps, including a decomposition of bending and breathing of the vertical components. This paper as a whole tries to be a compendium of different dynamical mechanisms whose predictions can be compared to the kinematic maps. Methods: Using analytical methods or simulations, we are able to predict the main dynamical factors and compare them to the predictions of the extended kinematic maps of Gaia-DR2. Results: The gravitational influence of Galactic components that are different from the disk, such as the long bar or bulge, the spiral arms, or a tidal interaction with Sagittarius dwarf galaxy, may explain some features of the velocity maps, especially in the inner parts of the disk. However, they are not sufficient in explaining the most conspicuous gradients in the outer disk. Vertical motions might be dominated by external perturbations or mergers, although a minor component may be due to a warp whose amplitude evolves with time. Here, we show with two different methods, which analyze the dispersion of velocities, that the mass distribution of the disk is flared. Despite these partial explanations, the main observed features can only be explained in terms of out-of-equilibrium models, which are either due to external perturbers or to the fact that the disk has not had time to reach equilibrium since its formation.

Physics of the gravitational effect of the galactic bar and spiral structure is presented. Physical equations differ from the conventionally used equations. Application to the motion of the Sun is treated. The speed of the Sun is taken to be consistent with the Oort constants. Galactic radial migration of the Sun is less than +- 0.4 kpc for the four-armed spiral structure. The Sun remains about 75 % of its existence within galactocentric distances (7.8 - 8.2) kpc and the results are practically independent on the spiral structure strength. Thus, the radial distance changes only within 5 % from the value of 8 kpc. Galactic radial migration of the Sun is less than +- (0.3 - 1.2) kpc, for the two-armed spiral structure. The Sun remains (29 - 95) % of its existence within galactocentric distances (7.8 - 8.2) kpc and the results strongly depend on the spiral structure strength and the angular speed of the spiral arms. The radial distance changes within (3.8 - 15.0) % from the value of 8 kpc. If observational arguments prefer relevant radial migration of the Sun, then the Milky Way is characterized by the two-arm spiral structure.