The Universe on scales $10-100~h^{-1}$ Mpc is dominated by a cosmic web of voids, filaments, sheets and knots of galaxy clusters. These structures participate differently in the global expansion of the Universe: from non-expanding clusters to the above average expansion rate of voids. In this paper we characterize Hubble expansion anisotropies in the COMPOSITE sample of 4534 galaxies and clusters. We concentrate on the dipole and quadrupole in the rest frame of the Local Group. These both have statistically significant amplitudes. These anisotropies, and their redshift dependence, cannot be explained solely by a boost of the Local Group in the Friedmann-Lema\^itre-Robertson-Walker (FLRW) model which expands isotropically in the rest frame of the cosmic microwave background (CMB) radiation. We simulate the local expansion of the Universe with inhomogeneous Szekeres models, which match the standard FLRW model on $> 100~h^{-1}$ Mpc scales but exhibit nonkinematic relativistic differential expansion on small scales. We restrict models to be consistent with observed CMB temperature anisotropies, while simultaneously fitting the redshift variation of the Hubble expansion dipole. We include features to account for both the Local Void and the "Great Attractor". While this naturally accounts for the Hubble expansion and CMB dipoles, the simulated quadrupoles are smaller than observed. Further refinement to incorporate additional structures may improve this. This would enable a test of the hypothesis that some large angle CMB anomalies result from failing to treat the relativistic differential expansion of the background geometry, a natural feature of solutions to Einstein's equations not included in the current standard model of cosmology.

We analyze the spectrum of cosmic microwave background (CMB) anisotropies in the timescape cosmology: a potentially viable alternative to homogeneous isotropic cosmologies without dark energy. We exploit the fact that the timescape cosmology is extremely close to the standard cosmology at early epochs to adapt existing numerical codes to produce CMB anisotropy spectra, and to match these as closely as possible to the timescape expansion history. A variety of matching methods are studied and compared. We perform Markov Chain Monte Carlo analyses on the parameter space, and fit CMB multipoles $50 \le \ell \le 2500$ to the Planck satellite data. Parameter fits include a dressed Hubble constant, $H_0 = 61.0$ km/s/Mpc ($\pm$ 1.3% stat) ($\pm$ 8% sys), and a present void volume fraction $f_{v0} = 0.627$ ($\pm$ 2.3% stat) ($\pm$ 13% sys). We find best fit likelihoods which are comparable to that of the best fit Lambda CDM cosmology in the same multipole range. In contrast to earlier results, the parameter constraints afforded by this analysis no longer admit the possibility of a solution to the primordial lithium abundance anomaly. This issue is related to a strong constraint between the ratio of baryonic to nonbaryonic dark matter and the ratio of heights of the second and third acoustic peaks, which cannot be changed as long as the standard cosmology is assumed up to the surface of last scattering. These conclusions may change if backreaction terms are also included in the radiation-dominated primordial plasma.

The timescape cosmology represents a potentially viable alternative to the standard homogeneous cosmology, without the need for dark energy. Although average cosmic evolution in the timescape scenario only differs substantially from that of Friedmann-Lemaitre model at relatively late epochs when the contribution from the energy density of radiation is negligible, a full solution of the Buchert equations to incorporate radiation is necessary to smoothly match parameters to the epoch of photon decoupling and to obtain constraints from cosmic microwave background data. Here we extend the matter-dominated solution found in earlier work to include radiation, providing series solutions at early times and an efficient numerical integration strategy for generating the complete solution. The numerical solution is used to directly calculate the scale of the sound horizon at decoupling, and at the baryon drag epoch. The constraints on these scales from the Planck satellite data yield bounds on the timescape cosmological parameters, which are found to also agree with the best-fit values from a recent analysis of SDSS-II supernova data, while avoiding the problem of a primordial lithium-7 abundance anomaly.