We exploit the 1+1+2 formalism to covariantly describe the inhomogeneous and anisotropic Szekeres models. It is shown that an \emph{average scale length} can be defined \emph{covariantly} which satisfies a 2d equation of motion driven from the \emph{effective gravitational mass} (EGM) contained in the dust cloud. The contributions to the EGM are encoded to the energy density of the dust fluid and the free gravitational field $E_{ab}$. We show that the quasi-symmetric property of the Szekeres models is justified through the existence of 3 independent \emph{Intrinsic Killing Vector Fields (IKVFs)}. In addition the notions of the Apparent and Absolute Apparent Horizons are briefly discussed and we give an alternative gauge-invariant form to define them in terms of the kinematical variables of the spacelike congruences. We argue that the proposed program can be used in order to express Sachs' optical equations in a covariant form and analyze the confrontation of a spatially inhomogeneous irrotational overdense fluid model with the observational data.

"Diagonal" spatially inhomogeneous (SI) models are introduced under the assumption of the existence of (proper) intrinsic symmetries and can be seen, in some sense, complementary to the Szekeres models. The structure of this class of spacetimes can be regarded as a generalization of the (twist-free) Locally Rotationally Symmetric (LRS) geometries without any global isometry containing, however, these models as special cases. We consider geometries where a six-dimensional algebra $\mathcal{IC}$ of Intrinsic Conformal Vector Fields (ICVFs) exists acting on a $2-$dimensional (pseudo)-Riemannian manifold. Its members $\mathbf{X}_{\alpha }$, constituted of 3 Intrinsic Killing Vector Fields (IKVFs) and 3 \emph{proper} and \emph{gradient} ICVFs, as well as the specific form of the gravitational field are given explicitly. An interesting consequence, in contrast with the Szekeres models, is the immediate existence of \emph{conserved quantities along null geodesics}. We check computationally that the magnetic part $H_{ab}$ of the Weyl tensor vanishes whereas the shear $\sigma_{ab}$ and the electric part $E_{ab}$ share a common eigenframe irrespective of the fluid interpretation of the models. A side result is the fact that the spacetimes are foliated by a set of \emph{conformally flat } $3-$dimensional \emph{timelike} slices when the anisotropy of the \emph{flux-free} fluid is described only in terms of the 3 \emph{principal inhomogeneous "pressures"} $p_{\alpha}$ or equivalently when the Ricci tensor shares the same basis of eigenvectors with $\sigma_{ab}$ and $E_{ab}$. The conformal flatness also indicates that a \emph{10-dimensional algebra} of ICVFs $\mathbf{\Xi}$ acting on the $3-$dimensional timelike slices is highly possible to exist enriching in that way the set of conserved quantities admitted by the SI models found in the present paper.

We expand the study of generalized brane cosmologies by allowing for an $f(\tilde{\cal R})$ gravity term on the brane, with $\tilde{\cal R}$ the curvature scalar derived from the induced metric. We also include arbitrary matter components on the brane and in the five-dimensional bulk. At low energies, the effect of the bulk on the brane evolution can be described through a mirage component, termed generalized dark radiation, in the effective four-dimensional field equations. Using the covariant formalism, we derive the exact form of these equations. We also derive an effective conservation equation involving the brane matter and the generalized dark radiation. At low energies the coupled brane-bulk system has a purely four-dimensional description. The applications of the formalism include generalizations of the Starobinsky model and the Dvali-Gabadadze-Porrati cosmology.

We extend the covariant analysis of the brane cosmological evolution in order to take into account, apart from a general matter content and an induced-gravity term on the brane, a Gauss-Bonnet term in the bulk. The gravitational effect of the bulk matter on the brane evolution can be described in terms of the total bulk mass as measured by a bulk observer at the location of the brane. This mass appears in the effective Friedmann equation through a term characterized as generalized dark radiation that induces mirage effects in the evolution. We discuss the normal and self-accelerating branches of the combined system. We also derive the Raychaudhuri equation that can be used in order to determine if the cosmological evolution is accelerating.

The consistency of the constraint with the evolution equations for spatially inhomogeneous and irrotational silent (SIIS) models of Petrov type I, demands that the former are preserved along the timelike congruence represented by the velocity of the dust fluid, leading to \emph{new} non-trivial constraints. This fact has been used to conjecture that the resulting models correspond to the spatially homogeneous (SH) models of Bianchi type I, at least for the case where the cosmological constant vanish. By exploiting the full set of the constraint equations as expressed in the 1+3 covariant formalism and using elements from the theory of the spacelike congruences, we provide a direct and simple proof of this conjecture for vacuum and dust fluid models, which shows that the Szekeres family of solutions represents the most general class of SIIS models. The suggested procedure also shows that, the uniqueness of the SIIS of the Petrov type D is not, in general, affected by the presence of a non-zero pressure fluid. Therefore, in order to allow a broader class of Petrov type I solutions apart from the SH models of Bianchi type I, one should consider more general ``silent'' configurations by relaxing the vanishing of the vorticity and the magnetic part of the Weyl tensor but maintaining their ``silence'' properties i.e. the vanishing of the curls of $E_{ab},H_{ab}$ and the pressure $p$.

We study the cosmological evolution on a brane with induced gravity within a bulk with arbitrary matter content. We consider a Friedmann-Robertson-Walker brane, invariantly characterized by a six-dimensional group of isometries. We derive the effective Friedmann and Raychaudhuri equations. We show that the Hubble expansion rate on the brane depends on the covariantly defined integrated mass in the bulk, which determines the energy density of the generalized dark radiation. The Friedmann equation has two branches, distinguished by the two possible values of the parameter $\ex=\pm 1$. The branch with $\ex=1$ is characterized by an effective cosmological constant and accelerated expansion for low energy densities. Another remarkable feature is that the contribution from the generalized dark radiation appears with a negative sign. As a result, the presence of the bulk corresponds to an effective negative energy density on the brane, without violation of the weak energy condition. The transition from a period of domination of the matter energy density by non-relativistic brane matter to domination by the generalized dark radiation corresponds to a crossing of the phantom divide $w=-1$.

The acceleration parameter defined through the local volume expansion is negative for a pressureless, irrotational fluid with positive energy density. In the presence of inhomogeneities or anisotropies the volume expansion rate results from averaging over various directions. On the other hand, the observation of light from a certain source in the sky provides information on the expansion along the direction to that source. If there are preferred directions in the underlying geometry one can define several expansion parameters. We provide such definitions for the case of the Tolman-Bondi metric. We then examine the effect of a localized inhomogeneity on the surrounding cosmological fluid. Our framework is similar in spirit to the model of spherical collapse. For an observer in the vicinity of a central overdensity, the perceived local evolution is consistent with acceleration in the direction towards the center of the overdensity, and deceleration perpendicularly to it. A negative mass leads to deceleration along the radial direction, and acceleration perpendicularly to it. If the observer is located at the center of an overdensity the null geodesics are radial. The form of the luminosity distance as a function of the redshift is consistent with acceleration for a certain range of redshifts.

We derive exact solutions of the Einstein equations in the context of the Randall-Sundrum model with matter both on the brane and in the bulk. The bulk metric is a generalization of the static metrics describing the interior of stellar objects. We study the cosmological evolution on the brane. Our solutions describe energy ouflow from the brane, and the appearance of ``mirage'' contributions to the Hubble expansion that alter the standard cosmological evolution.