We discuss a new formalism for light propagation which can be used within the regime of validity of geometric optics, but with no limitation on the scales of interest: from inside the Galaxy to the ultra-large scales of cosmology. One of our main results is that within this framework it is possible to calculate all relevant observables (image magnification, parallax, position drift or proper motion) by simply differentiating the photon trajectory with respect to the initial data. We then focus on a new observable, which we name the distance slip: it is defined as the relative difference between the angular diameter distance and the parallax distance. Its peculiarity lies in the fact that its value is independent of the momentary motions of both the source and the observer and that for short distances it shows a tomographic property, being proportional to the amount of matter along the line of sight. After describing further its properties and methods of measurement, we specialize our study of the distance slip to cosmology. We show that it does not depend on the Hubble constant H_0 and that its dependence on the other cosmological parameters is different from other distance indicators. This suggests that the distance slip may contain new information.

We consider the standard problem of observational astronomy, i.e. the observations of light emission from a distant region of spacetime in general relativity. The goal is to describe the changes between the measurements of the light performed by a sample of observers slightly displaced with respect to each other and moving with different 4-velocities and 4-accelerations. In our approach, all results of observations can be expressed as functions of the kinematic variables, describing the motions of the observers and the emitting bodies with respect to their local inertial frames, and four linear bilocal geodesic operators, describing the influence of the spacetime geometry on light propagation. The operators are functionals of the curvature tensor along the line of sight. The results are based on the assumption that the regions of emissions and observations are sufficiently small so that the spacetime curvature effects are negligible within each of them, although they are significant for the light propagation between them. The new formulation provides a uniform approach to optical phenomena in curved spacetimes and, as an application, we discuss the problem of a fully relativistic definition of the parallax and position drifts (or proper motions). We then use the results to construct combinations of observables which are completely insensitive to the motion of both the observer and the emitter. These combinations by construction probe the spacetime geometry between the observation and emission regions and in our formalism we may express them as functionals of the Riemann tensor along the line of sight. For short distances one of these combinations depends only on the matter content along the line of sight. This opens up the possibility to measure the matter content of a spacetime in a tomography-like manner irrespective of the motions of the emitter and the observer.

We consider the question of determining the optical drift effects in general relativity, i.e. the rate of change of the apparent position, redshift, Jacobi matrix, angular distance and luminosity distance of a distant object as registered by an observer in an arbitrary spacetime. We present a fully relativistic and covariant approach, in which the problem is reduced to a hierarchy of ODE's solved along the line of sight. The 4-velocities and 4-accelerations of the observer and the emitter and the geometry of the spacetime along the line of sight constitute the input data. We build on the standard relativistic geometric optics formalism and extend it to include the time derivatives of the observables. In the process we obtain two general, non-perturbative relations: the first one between the gravitational lensing, represented by the Jacobi matrix, and the apparent position drift, also called the cosmic parallax, and the second one between the apparent position drift and the redshift drift. The applications of the results include the theoretical study of the drift effects of cosmological origin (so-called real-time cosmology) in numerical or exact Universe models.

The apparent properties of distant objects encode information about the way the light they emit propagates to an observer, and therefore about the curvature of the underlying spacetime. Measuring the relationship between the redshift $z$ and the luminosity distance $D_{\rm L}$ of a standard candle, for example, yields information on the Universe's matter content. In practice, however, in order to decode this information the observer needs to make an assumption about the functional form of the $D_{\rm L}(z)$ relation; in other words, a cosmological model needs to be assumed. In this work, we use numerical-relativity simulations, equipped with a new ray-tracing module, to numerically obtain this relation for a few black-hole--lattice cosmologies and compare it to the well-known Friedmann-Lema\^itre-Robertson-Walker case, as well as to other relevant cosmologies and to the Empty-Beam Approximation. We find that the latter provides the best estimate of the luminosity distance and formulate a simple argument to account for this agreement. We also find that a Friedmann-Lema\^itre-Robertson-Walker model can reproduce this observable exactly, as long as a time-dependent cosmological constant is included in the fit. Finally, the dependence of these results on the lattice mass-to-spacing ratio $\mu$ is discussed: we discover that, unlike the expansion rate, the $D_{\rm L}(z)$ relation in a black-hole lattice does not tend to that measured in the corresponding continuum spacetime as $\mu \to 0$.

Astronomical observations reveal hierarchical structures in the Universe, from galaxies, groups of galaxies, clusters and superclusters, to filaments and voids. On the largest scales it seems that some kind of statistical homogeneity can be observed. As a result, modern cosmological models are based on spatially homogeneous and isotropic solutions of the Einstein equations, and the evolution of the universe is approximated by the Friedmann equations. In parallel to standard homogeneous cosmology, the field of inhomogeneous cosmology and backreaction is being developed. This field investigates whether small scale inhomogeneities via non-linear effects can backreact and alter the properties of the Universe on its largest scales, leading to a non-Friedmannian evolution. This paper presents the current status of inhomogeneous cosmology and backreaction. It also discusses future prospects of the field of inhomogeneous cosmology, which is based on a survey of 50 academics working in the field of inhomogeneous cosmology.

When do the nonlinear effects of general relativity matter in astrophysical situations? They are obviously relevant for very compact sources of the gravitational field, such as neutron stars or black holes. In this paper I discuss another, less obvious situation, in which large relativistic effects may arise due to a complicated, multiscale structure of the matter distribution. I present an exact solution with an inhomogeneous energy density distribution in the form of a hierarchy of nested voids and overdensities of various sizes, extending from the homogeneity scale down to arbitrary small scales. I show that although each of the voids and overdensities seems to be very weakly relativistic, and thus easy to describe using the linearized general relativity, the solution taken as a whole lies in fact in the nonlinear regime. Its nonlinear properties are most easily seen when we compare the ADM mass of the solution and the integral of the local mass density: the difference between them, i.e. the relativistic mass deficit, can be significant provided that the inhomogeneities extend to sufficiently small scales. The non-additivity of masses implies a large backreaction effect, i.e. significant discrepancy between the averaged, large-scale effective stress-energy tensor and the naive average of the local energy density. I show that this is a general relativistic effect arising due to the inhomogeneous, multiscale structure. I also discuss the relevance of the results in cosmology and relativistic astrophysics.