We study the expected variance of measurements of the Hubble constant, $H_0$, as calculated in either linear perturbation theory or using non-linear velocity power spectra derived from $N$-body simulations. We compare the variance with that obtained by carrying out mock observations in the N-body simulations, and show that the estimator typically used for the local Hubble constant in studies based on perturbation theory is different from the one used in studies based on N-body simulations. The latter gives larger weight to distant sources, which explains why studies based on N-body simulations tend to obtain a smaller variance than that found from studies based on the power spectrum. Although both approaches result in a variance too small to explain the discrepancy between the value of $H_0$ from CMB measurements and the value measured in the local universe, these considerations are important in light of the percent determination of the Hubble constant in the local universe.

With the upcoming sky survey with the Large Synoptic Survey Telescope a great sample of type Ia supernovae will be observed, allowing for a precise mapping of the velocity structure of the universe. Since the source of peculiar velocities is variations in the density field, cosmological parameters related to the matter distribution can subsequently be extracted from the velocity power spectrum. One way to quantify this is through the angular power spectrum of radial peculiar velocities on spheres at different redshifts. We investigate how well this observable can be measured, despite the problems caused by areas with no information. To obtain a realistic distribution of supernovae, we create mock supernova catalogs by using a semi-analytical code for galaxy formation on the merger trees extracted from N-body simulations. We measure the cosmic variance in the velocity power spectrum by repeating the procedure many times for differently located observers, and vary different aspects of the analysis, such as the observer environment, to see how this affects the measurements. Our results confirm the findings from earlier studies regarding the precision with which the angular velocity power spectrum can be determined in the near future. This level of precision has been found to imply, that the angular velocity power spectrum from type Ia supernovae is competitive in its potential to measure parameters such as $\sigma_8$. This makes the peculiar velocity power spectrum from type Ia supernovae a promising new observable, which deserves further attention.

In the current state of cosmology, where cosmological parameters are being measured to percent accuracy, it is essential to understand all sources of error to high precision. In this paper we present the results of a study of the local variations in the Hubble constant measured at the distance scale of the Coma Cluster, and test the validity of correcting for the peculiar velocities predicted by gravitational instability theory. The study is based on N-body simulations, and includes models featuring a coupling between dark energy and dark matter, as well as two $\Lambda$CDM simulations with different values of $\sigma_8$. It is found that the variance in the local flows is significantly larger in the coupled models, which increases the uncertainty in the local measurements of the Hubble constant in these scenarios. By comparing the results from the different simulations, it is found that most of the effect is caused by the higher value of $\sigma_8$ in the coupled cosmologies, though this cannot account for all of the additional variance. Given the discrepancy between different estimates of the Hubble constant in the universe today, cosmological models causing a greater cosmic variance is something that we should be aware of.

The effects of local inhomogeneities on low redshift $H_0$ determinations are studied by estimating the redshift-distance relation of mock sources in N-body simulations. The results are compared to those obtained using the standard approach based on Hubble's law. The comparison shows a clear tendency for the standard approach to yield lower values of $H_0$ than the approach based on the scheme using light rays. The difference is, however, small. More precisely, it is found that the overall effect of inhomogeneities on the determination of $H_0$ is a small increase in the local estimates of about $0.3\%$ compared to the results obtained with Hubble's law, when based on a typical distribution of supernovae in the redshift range $0.01 < z < 0.1$. The overall conclusion of the study is a verification of the results that have earlier been obtained by using Hubble's law: The effects of inhomogeneities on local $H_0$ estimates are not significant enough to make it plausible that differences in high- and low-redshift estimates of $H_0$ are due to small inhomogeneities within the setting of standard cosmology.

We have carefully studied how local measurements of the Hubble constant, $H_0$, can be influenced by a variety of different parameters related to survey geometry, depth, and size, as well as observer position in space. Our study is based on N-body simulations of structure in the standard $\Lambda$CDM model and our conclusion is that the expected variance in measurements of $H_0$ is far too small to explain the current discrepancy between the low value of $H_0$ inferred from measurements of the cosmic microwave background (CMB) by the Planck collaboration and the value measured directly in the local universe by use of Type Ia supernovae. This conclusion is very robust and does not change with different assumptions about survey geometry and effective sky coverage or observer position in space.