This is the second of a series of papers extending the 1+3 covariant and gauge invariant treatment of kinetic theory to an examination of Cosmic Microwave Background temperature anisotropies arising from inhomogeneities in the early universe. The first paper dealt with algebraic issues. Here we derive the mode form of the integrated Boltzmann equations, first, giving a covariant version of the standard derivation using the mode recursion relations, second, demonstrating the link to the multipole divergence equations and finally various analytic ways of solving the resulting equations are discussed. A general integral form of solution is obtained for the equations with Thomson scattering. The covariant Friedmann-Lemaitre multipole form of the transport equations are found using the covariant and gauge-invariant generalization of the Peebles and Yu expansion in Thompson scattering time. The dispersion relations and damping scale are then obtained from the covariant approach. The equations are integrated to give the covariant and gauge-invariant equivalent of the canonical scalar sourced anisotropies. We carry out a simple treatment of the matter dominated free-streaming projection, slow decoupling, and tight-coupling cases, with the aim both giving a unified transparent derivation of this range of results and clarifying the connection between the more usual approaches (for example that of Hu and Sugiyama) and the treatment for scalar perturbations (for example the treatment of Challinor and Lasenby).

I outline some nonperturbative relativistic effects that arise from gravitational corrections to the Boltzmann equations. These may be important for the study of CMB temperature anisotropies, particularly their interpretation. These terms are not included in the canonical treatment as they arise from the exact equations. Here a weakly nonlinear investigation of these effects is defined and investigated with an emphasis on a Rees-Sciama sourced effect -- the imprint of structure evolution on the CMB. It is shown that gravitational nonlinearity in the weakly nonlinear extension of almost-FLRW temperature anisotropies leads to cancellation on small-scales when threading in the Newtonian frame. In the general frame this cancellation does not occur. In the context of a flat almost-FLRW CDM model I provide a heuristic argument for a nonperturbative small scale correction, due to the Rees-Sciama effect, of not more than $\Delta T/T \sim 10^{-6}-10^{-5}$ near $\ell \sim 100 - 300$. The effect of mild gravitational nonlinearity could be more sophisticated than previously expected.

This is the first of a series of papers extending a covariant and gauge invariant (CGI) treatment of kinetic theory in curved space-times to a treatment of Cosmic Background Radiation (CBR) temperature anisotropies arising from inhomogeneities in the early universe. This paper deals with algebraic issues, both generically and in the context of models linearised about RW geometries. The approach represents radiation anisotropies by PSTF tensors. The Angular correlation functions for the mode coefficients are found in terms of these quantities, following the Wilson-Silk approach, but derived and dealt with in CGI form. The covariant multipole and mode-expanded angular correlation functions are related to the usual treatments in the literature. The CGI mode expansion is related to the coordinate approach by linking the Legendre functions to the PSTF representation, using a covariant addition theorem for the tensors to generate the Legendre Polynomial recursion relation. This paper lays the foundation for further papers in the series, which use this formalism in a CGI approach to develop solutions of the Boltzmann and Liouville equations for the CBR before and after decoupling, thus providing a unified CGI derivation of the variety of approaches to CBR anisotropies in the current literature.

Assuming that the cosmological principle holds, Maartens, Ellis and Stoeger (MES) recently constructed a detailed scheme linking anisotropies in the cosmic background radiation (CMB) with anisotropies and inhomogeneities in the large scale structure of the universe and showed how to place limits on those anisotropies and inhomogeneities simply by using CMB quadrupole and octupole limits. First we indicate and discuss the connection between the covariant multipole moments of the temperature anisotropy used in the MES scheme and the quadrupole and octupole results from COBE. Then we introduce those results into the MES limit equations to obtain definite quantitative limits on the complete set of cosmological measures of anisotropy and inhomogeneity. We find that all the anisotropy measures are less than 10^{-4} in the case of those not affected by the expansion rate H, and less than 10^{-6} Mpc^{-1} in the case of those which are. These results quantitatively demonstrate that the observable universe is indeed close to Friedmann-Lemaitre-Robertson-Walker (FLRW) on the largest scales, and can be adequately modelled by an almost-FLRW model -- that is, the anisotropies and inhomogeneities characterizing the observable universe on the largest scales are not too large to be considered perturbations to FLRW.

We develop a new approach to local nonlinear effects in cosmic microwave background anisotropies, and discuss the qualitative features of these effects. New couplings of the baryonic velocity to radiation multipoles are found, arising from nonlinear Thomson scattering effects. We also find a new nonlinear shear effect on small angular scales. The full set of evolution and constraint equations is derived, including the nonlinear generalizations of the radiation multipole hierarchy, and of the dynamics of multi-fluids. These equations govern radiation anisotropies in any inhomogeneous spacetime, but their main application is to second-order effects in a universe that is close to the Friedmann models. Qualitative analysis is given here, and quantitative calculations are taken up in further papers.