Number counts observations available with new surveys such as the Euclid mission will be an important source of information about the metric of the Universe. We compute the low red-shift expansion for the energy density and the density contrast using an exact spherically symmetric solution in presence of a cosmological constant. At low red-shift the expansion is more precise than linear perturbation theory prediction. We then use the local expansion to reconstruct the metric from the monopole of the density contrast. We test the inversion method using numerical calculations and find a good agreement within the regime of validity of the red-shift expansion. The method could be applied to observational data to reconstruct the metric of the local Universe with a level of precision higher than the one achievable using perturbation theory.

Recent measurements of the cosmic microwave background (CMB) radiation have shown an apparent tension with the present value of the Hubble parameter inferred from local observations of supernovae, which look closer, i.e. brighter, than what is expected in a homogeneous model with a value of $H_0$ equal to the one estimated from CMB observations. We examine the possibility that such a discrepancy is the consequence of the presence of a local inhomogeneity seeded by primordial curvature perturbations, finding that a negative peak of the order of less than two standard deviations could allow to fit low red-shift supernovae observations without the need of using a value of the Hubble parameter different from $H_0^{CMB}$. The type of inhomogeneity we consider does not modify the distance to the last scattering, making it compatible with the constraints of the PLANCK mission data. Our calculation is fully relativistic and non perturbative, and for this reason shows important effects which were missed in the previous investigations using relativistic perturbations or Newtonian approximations, because the structures seeded by primordial curvature perturbations can be today highly non linear, and relativist Doppler terms cannot be neglected. Comparison with studies of local structure based on galaxy surveys and luminosity density analysis reveals that the density profile we find could in fact be compatible with the one obtained for the same region of sky where is located most of the Cepheids used to calibrate the luminosity distance of the supernovae employed for the local $H_0$ estimation, suggesting a possible directional dependence and calibration bias which could be partially attributed to the presence of the Sloan Great Wall.

In order to estimate the effects of local structure on the Hubble parameter we calculate the low-redshift expansion for $H(z)$ and $\frac{\delta H}{H}$ for an observer at the center of a spherically symmetric matter distribution in presence of a cosmological constant. We then test the accuracy of the formulae comparing them with fully relativistic non perturbative numerical calculations for different cases for the density profile. The low red-shift expansion we obtain gives results more precise than perturbation theory since is based on the use of an exact solution of Einstein's field equations. For larger density contrasts the low red-shift formulae accuracy improves respect to the perturbation theory accuracy because the latter is based on the assumption of a small density contrast, while the former does not rely on such assumption. The formulae can be used to take into account the effects on the Hubble expansion parameter due to the monopole component of the local structure. If the $H(z)$ observations will show deviations from the $\Lambda CDM$ prediction compatible with the formulae we have derived, this could be considered an independent evidence of the existence of a local inhomogeneity, and the formulae could be used to determine the characteristics of this local structure.