We investigate number count statistics as measures for transition to homogeneity of the matter distribution in the Universe and analyse how such statistics might be `dressed' by the assumed survey selection function. Since the estimated survey selection function -- which ideally accounts for selection bias in the observed distribution -- is degenerate with the estimated underlying distribution of galaxies, the ability to identify the correct survey selection function is of importance for obtaining reliable estimates for clustering statistics. Selection functions of existing galaxy catalogues are modelled from data to resemble the redshift distribution and mean density of the observed galaxies. Proposed estimates of the selection function for upcoming surveys in addition resemble the angular distribution of the observed galaxies. We argue that such modelling of the selection function could potentially underestimate the deviance from homogeneity at scales probed by existing catalogues. We investigate the impact of conventionally applied methods for estimation of the survey selection function on number count in sphere statistics in a toy model setting. The example density distribution is asymptotically homogeneous, while non-linear density fluctuations are present regionally. We find that density oscillations with period comparable to characteristic scales of the survey are suppressed when conventional estimates of the survey selection function are invoked, resulting in number count statistics which are biased towards homogeneity. For our concrete toy model with maximum density contrasts of 1 and period of the density oscillation comparable in size to the survey radius, we find that the homogeneity scale is underestimated by ~40%, however this quantitative result is dependent on the model setup including the form of the simplistic density field considered and the survey geometry.

Recently it has been noted by Di Valentino, Melchiorri and Silk (2019) that the enhanced lensing signal relative to that expected in the spatially flat $\Lambda$CDM model poses a possible crisis for the Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) class of models usually used to interpret cosmological data. The 'crisis' amounts to inconsistencies between cosmological datasets arising when the FLRW curvature parameter $\Omega_k$ is determined from the data rather than constrained to be zero a priori. Moreover, the already substantial discrepancy between the Hubble parameter as determined by Planck and local observations increases to the level of $5\sigma$. While such inconsistencies might arise from systematic effects of astrophysical origin affecting the Planck Cosmic Microwave Background (CMB) power spectra at small angular scales, it is an option that the inconsistencies are due to the failure of the FLRW assumption. In this paper we recall how the FLRW curvature ansatz is expected to be violated for generic relativistic spacetimes. We explain how the FLRW conservation equation for volume-averaged spatial curvature is modified through structure formation, and we illustrate in a simple framework how the curvature tension in a FLRW spacetime can be resolved - and is even expected to occur - from the point of view of general relativity. Requiring early-time convergence towards a Friedmannian model with a spatial curvature parameter $\Omega_{k0}$ equal to that preferred from the Planck power spectra resolves the Hubble tension within our dark energy-free model.

We investigate - in a generic setting - the regime of applicability of the Alcock-Paczynski (AP) scaling conventionally applied to test different cosmological models, given a fiducial measurement of the baryon acoustic oscillation (BAO) characteristic scale in the galaxy 2-point correlation function. We quantify the error in conventional AP scaling methods, for which our ignorance about the true cosmology is parameterised in terms of two constant AP scaling parameters. We propose a new, and as it turns out, improved version of the constant AP scaling, also consisting of two scaling parameters. The two constant AP scaling methods are almost indistinguishable when the fiducial model used in data reduction and the "true" underlying cosmology are not differing substantially in terms of metric gradients, but are otherwise expected to differ. Our new methods can be applied to existing analyses through a reinterpretation of the results of the conventional AP scaling. This reinterpretation might be important in model universes where curvature gradients above the scale of galaxies are significant. We test our theoretical findings on $\Lambda$CDM mock catalogues. The conventional constant AP scaling methods are surprisingly successful for pairs of large-scale metrics, but eventually break down when toy models allowing for large metric gradients are tested. The new constant AP scaling methods proposed in this paper are efficient for all test models examined. We find systematic errors of ~1% in the recovery of the BAO scale when the true model is distant from the fiducial, which cannot be attributed to any constant AP approximation. The level of systematic uncertainty is robust to the exact fitting method employed. This indicates that caution must be taken with the error budget when extrapolating the BAO acoustic scale measurements obtained in the standard literature.

Few statements in cosmology can be made without assuming a cosmological model within which to interpret data. Statements about cosmic acceleration are no exception to this rule, and the inferred positive volume acceleration of our Universe often quoted in the literature is valid in the context of the standard Friedmann-Lemaitre-Robertson-Walker (FLRW) class of space-times. Using the Joint Light-curve Analysis (JLA) catalogue of type Ia supernovae (SNIa), we examine the fit of a class of exact scaling solutions with dynamical spatial curvature formulated in the framework of a scalar averaging scheme for relativistic inhomogeneous space-times. In these models, global volume acceleration may emerge as a result of the non-local variance between expansion rates of clusters and voids, the latter gaining volume dominance in the late-epoch Universe. We find best-fit parameters for a scaling model of backreaction that are reasonably consistent with previously found constraints from SNIa, CMB, and baryon acoustic oscillations data. The quality of fit of the scaling solutions is indistinguishable from that of the LCDM model and the timescape cosmology from an Akaike Information Criterion (AIC) perspective. This indicates that a broad class of models can account for the $z\lesssim 1$ expansion history.

We develop methods for investigating baryon acoustic oscillation (BAO) features in cosmological models with non-trivial (but slowly varying) averaged spatial curvature: models that are not necessarily flat, close to flat, nor with constant spatial curvature. The class of models to which our methods apply include Lemaitre-Tolman-Bondi models, modified gravity cosmologies, and inhomogeneous cosmologies with backreaction - in which we do not have a prediction of the shape of the spatial 2-point correlation function, but where we nevertheless expect to see a BAO feature in the present-day galaxy distribution, in form of an excess in the galaxy 2-point correlation function. We apply our methods to the Baryon Oscillation Spectroscopic Survey (BOSS) dataset, investigating both the Lambda Cold Dark Matter ($\Lambda$CDM) and timescape cosmological models as case studies. The correlation functions measured in the two fiducial models contain a similarly-pronounced BAO feature. We use the relative tangential and radial BAO scales to measure the anisotropic Alcock-Paczy\'nski distortion parameter, $\epsilon$, which is independent of the underlying BAO preferred scale. We find that $\epsilon$ is consistent with zero in both fiducial cosmologies, indicating that models with a different spatial curvature behaviour can account for the relative positions of the tangential and radial BAO scale. We validate our methods using $\Lambda$CDM mocks.

We introduce a generalization of the 4-dimensional averaging window function of Gasperini, Marozzi and Veneziano (2010) that may prove useful for a number of applications. The covariant nature of spatial scalar averaging schemes to address the averaging problem in relativistic cosmology is an important property that is implied by construction, but usually remains implicit. We employ here the approach of Gasperini et al. for two reasons. First, the formalism and its generalization presented here are manifestly covariant. Second, the formalism is convenient for disentangling the dependencies on foliation, volume measure, and boundaries in the averaged expressions entering in scalar averaging schemes. These properties will prove handy for simplifying expressions, but also for investigating extremal foliations and for comparing averaged properties of different foliations directly. The proposed generalization of the window function allows for choosing the most appropriate averaging scheme for the physical problem at hand, and for distinguishing between the role of the foliation itself and the role of the volume measure in averaged dynamic equations. We also show that one particular window function obtained from this generalized class results in an averaging scheme corresponding to that of a recent investigation by Buchert, Mourier and Roy (2018) and, as a byproduct, we explicitly show that the general equations for backreaction derived therein are covariant.

Parameters that quantify the acceleration of cosmic expansion are conventionally determined within the standard Friedmann-Lemaitre-Robertson-Walker (FLRW) model, which fixes spatial curvature to be homogeneous. Generic averages of Einstein's equations in inhomogeneous cosmology lead to models with non-rigidly evolving average spatial curvature, and different parametrizations of apparent cosmic acceleration. The timescape cosmology is a viable example of such a model without dark energy. Using the largest available supernova data set, the JLA catalogue, we find that the timescape model fits the luminosity distance-redshift data with a likelihood that is statistically indistinguishable from the standard spatially flat $\Lambda$ cold dark matter cosmology by Bayesian comparison. In the timescape case cosmic acceleration is non-zero but has a marginal amplitude, with best-fitting apparent deceleration parameter, $q_0=-0.043^{+0.004}_{-0.000}$. Systematic issues regarding standardization of supernova light curves are analysed. Cuts of data at the statistical homogeneity scale affect light curve parameter fits independent of cosmology. A cosmological model dependence of empirical changes to the mean colour parameter is also found. Irrespective of which model ultimately fits better, we argue that as a competitive model with a non-FLRW expansion history, the timescape model may prove a useful diagnostic tool for disentangling selection effects and astrophysical systematics from the underlying expansion history.