The standard cosmological model predicts the existence of a Cosmic Neutrino Background, which has not yet been observed directly. Some experiments aiming at its detection are currently under development, despite the tiny kinetic energy of the cosmological relic neutrinos, which makes this task incredibly challenging. Since massive neutrinos are attracted by the gravitational potential of our Galaxy, they can cluster locally. Neutrinos should be more abundant at the Earth position than at an average point in the Universe. This fact may enhance the expected event rate in any future experiment. Past calculations of the local neutrino clustering factor only considered a spherical distribution of matter in the Milky Way and neglected the influence of other nearby objects like the Virgo cluster, although recent $N$-body simulations suggest that the latter may actually be important. In this paper, we adopt a back-tracking technique, well established in the calculation of cosmic rays fluxes, to perform the first three-dimensional calculation of the number density of relic neutrinos at the Solar System, taking into account not only the matter composition of the Milky Way, but also the contribution of the Andromeda galaxy and the Virgo cluster. The effect of Virgo is indeed found to be relevant and to depend non-trivially on the value of the neutrino mass. Our results show that the local neutrino density is enhanced by 0.53% for a neutrino mass of 10 meV, 12% for 50 meV, 50% for 100 meV or 500% for 300 meV.

Massive neutrinos, due to their free streaming, produce a suppression in the matter power spectrum at intermediate and small scales which could be probed by galaxy clustering and/or weak lensing observables. This effect happens at scales that are also influenced by baryon feedback, i.e. galactic winds or Active Galactic Nuclei (AGN) feedback, which in realistic hydrodynamic simulations has also been shown to produce a suppression of power. Leaving aside, for the moment, the complex issue of galaxy bias, we focus here on matter clustering and tomographic weak lensing, we investigate the possible degeneracy between baryon feedback and neutrinos showing that it is not likely to degrade significantly the measurement of neutrino mass in future surveys. To do so, we generate mock data sets and fit them using the Markov Chain Monte Carlo (MCMC) technique and explore degeneracies between feedback parameters and neutrino mass. We model baryon feedback through fitting functions, while massive neutrinos are accounted for, also in the non-linear regime, using Halofit calibrated against accurate N-body neutrino simulations. In the error budget, we include the uncertainty in the modelling of non-linearities. For both matter clustering and weak lensing, we always recover the input neutrino mass within $\sim 0.25\sigma$ confidence level. Finally, we also take into account the intrinsic alignment effect in the weak lensing mock data. Even in this case, we are able to recover the right parameters: in particular, we find a significant degeneracy pattern between $M_\nu$ and the intrinsic alignment parameter $A_\mathrm{IA}$ .