This article explores in depth a watershed concept to partition the universe, introduced in \cite{2019MNRAS.489L...1D} and applied to the {\it Cosmicflows-3} observational dataset. We present a series of tests conducted with cosmological dark matter simulations. In particular we are interested in quantifying the evolution with redshift of large scale structures when defined as segmented basins of attraction. This new dynamical definition in the field of measuring standard rulers demonstrates robustness since all basins show a density contrast $\delta$ above one (mean universe density) independently of the simulation spatial resolution or the redshift. Another major finding is that density profiles of the basins show universality in slope. Consequently, there is a unique definition of what is a gravitational watershed at large scale, that can be further used as a probe for cosmology studies.

Type Ia Supernovae have yet again the opportunity to revolutionize the field of cosmology as the new generation of surveys are acquiring thousands of nearby SNeIa opening a new era in cosmology: the direct measurement of the growth of structure parametrized by $fD$. This method is based on the SNeIa peculiar velocities derived from the residual to the Hubble law as direct tracers of the full gravitational potential caused by large scale structure. With this technique, we could probe not only the properties of dark energy, but also the laws of gravity. In this paper we present the analytical framework and forecasts. We show that ZTF and LSST will be able to reach 5\% precision on $fD$ by 2027. Our analysis is not significantly sensitive to photo-typing, but known selection functions and spectroscopic redshifts are mandatory. We finally introduce an idea of a dedicated spectrograph that would get all the required information in addition to boost the efficiency to each SNeIa so that we could reach the 5\% precision within the first two years of LSST operation and the few percent level by the end of the survey.

This letter presents a new approach using the cosmic peculiar velocity field to characterize the morphology and size of large scale structures in the local Universe. The algorithm developed uses the three-dimensional peculiar velocity field to compute flow lines, or streamlines. The local Universe is then partitioned into volumes corresponding to gravitational basins, also called watersheds, among the different end-points of the velocity flow lines. This new methodology is first tested on numerical cosmological simulations, used as benchmark for the method, and then applied to the {\it Cosmic-Flows} project observational data in order to to pay particular attention to the nearby superclusters including ours. More extensive tests on both simulated and observational data will be discussed in an accompanying paper.