Averaging and evolving inhomogeneities are non-commuting operations. This implies the existence of deviations of an averaged model from the standard Friedmann-Lemaitre cosmologies. We quantify these deviations, encoded in a backreaction parameter, in the framework of Newtonian cosmology. We employ the linear theory of gravitational instability in the Eulerian and Lagrangian approaches, as well as the spherically- and plane-symmetric solutions as standards of reference. We propose a model for the evolution of the average characteristics of a spatial domain for generic initial conditions that contains the spherical top-hat model and the planar collapse model as exact sub cases. A central result is that the backreaction term itself, calculated on sufficiently large domains, is small but, still, its presence can drive the cosmological parameters on the averaging domain far away from their global values of the standard model. We quantify the variations of these parameters in terms of the fluctuations in the initial data as derived from the power spectrum of initial cold dark matter density fluctuations. E.g. in a domain with a radius of 100Mpc today and initially one-sigma fluctuations, the density parameters deviate from their homogeneous values by 15%; three-sigma fluctuations lead to deviations larger than 100%.

Most cosmological models studied today are based on the assumption of homogeneity and isotropy. Observationally one can find evidence that supports these assumptions on very large scales, the strongest being the almost isotropy of the Cosmic Microwave Background radiation after assigning the whole dipole to our proper motion relative to this background. However, on small and on intermediate scales up to several hundreds of Mpcs, there are strong deviations from homogeneity and isotropy. Here the problem arises how to relate the observations with the homogeneous and isotropic models. The usual proposal for solving this problem is to assume that Friedmann-Lemaitre models describe the mean observables. Such mean values may be identified with spatial averages. For Newtonian fluid dynamics the averaging procedure has been discussed in detail in Buchert and Ehlers (1997), leading to an additional backreaction term in the Friedmann equation. We use the Eulerian linear approximation and the `Zel'dovich approximation' to estimate the effect of the backreaction term on the expansion. Our results indicate that even for domains matching the background density in the mean, the evolution of the scale factor strongly deviates from the Friedmann solution, critically depending on the velocity field inside.