**Abstract** : Full sky surveys of peculiar velocity are arguably the best way to map the large scale structure out to distances of a few times 100 Mpc/h. Using the largest and most accurate ever catalog of galaxy peculiar velocities "Cosmicflows-2", the large scale structure has been reconstructed by means of the Wiener filter and constrained realizations assuming as a Bayesian prior model the LCDM model with the WMAP inferred cosmological parameters. The present paper focuses on studying the bulk flow of the local flow field, defined as the mean velocity of top-hat spheres with radii ranging out to R=500 Mpc/h. The estimated large scale structures, in general, and the bulk flow, in particular, are determined by the tension between the observational data and the assumed prior model. A prerequisite for such an analysis is the requirement that the estimated bulk flow is consistent with the prior model. Such a consistency is found here. At R=50(150) Mpc/h the estimated bulk velocity is 250+/-21 (239+/-38) km/s. The corresponding cosmic variance at these radii is 126(60)km/s, which implies that these estimated bulk flows are dominated by the data and not by the assumed prior model. The estimated bulk velocity is dominated by the data out to R~200 Mpc/h, where the cosmic variance on the individual Supergalactic Cartesian components (of the r.m.s. values) exceeds the variance of the Constrained Realizations by at least a factor of 2. The supergalactic SGX and SGY components of the CMB dipole velocity are recovered by the Wiener filter velocity field down to a very few km/s. The SGZ component of the estimated velocity, the one that is most affected by the Zone of Avoidance, is off by 126 km/s (an almost 2 sigma discrepancy).