# Quantifying cosmic variance

## Abstract

We determine an expression for the cosmic variance of any “normal” galaxy survey based on examination of mag galaxies in the SDSS DR7 data cube. We find that cosmic variance will depend on a number of factors principally: total survey volume, survey aspect ratio, and whether the area surveyed is contiguous or comprised of independent sight-lines. As a rule of thumb cosmic variance falls below 10% once a volume of Mpc is surveyed for a single contiguous region with a 1:1 aspect ratio. Cosmic variance will be lower for higher aspect ratios and/or non-contiguous surveys. Extrapolating outside our test region we infer that cosmic variance in the entire SDSS DR7 main survey region is % to

The equation obtained from the SDSS DR7 region can be generalised to estimate the cosmic variance for any density measurement determined from normal galaxies (e.g., luminosity densities, stellar mass densities and cosmic star-formation rates) within the volume range to Mpc.

We apply our equation to show that 2 sightlines are required to ensure cosmic variance is % in any ASKAP galaxy survey (divided into intervals, i.e., Gyr intervals for ). Likewise 10 MeerKAT sightlines will be required to meet the same conditions. GAMA, VVDS, and zCOSMOS all suffer less than 10% cosmic variance ( 3%-8%) in intervals of 0.1, 0.25, and 0.5 respectively. Finally we show that cosmic variance is potentially at the 50-70% level, or greater, in the HST Ultra Deep Field depending on assumptions as to the evolution of clustering. 100 or 10 independent sightlines will be required to reduce cosmic variance to a manageable level (%) for HST ACS or HST WFC3 surveys respectively (in intervals). Cosmic variance is therefore a significant factor in the HST studies currently underway.

###### keywords:

galaxies: general — galaxies: luminosity functions, mass functions — galaxies: statistics — cosmology: large-scale structure of the Universe^{1}

^{2}

^{3}

## 1 Introduction

The Universe is not homogeneous except on the largest scales (Gpc,
Davis et al. 1985). As a consequence number and density measurements
derived from within modest volumes will show greater than Poisson
variation (see Szapudi & Colombi 1996 for example). This cosmic
variance^{4}

In this paper we aim to use the largest volume survey to date, the Sloan Digital Sky Survey (SDSS), to empirically determine some generically useful formulae for estimating the cosmic variance as a function of survey volume and survey shape. These formulae should also assist in the design of future surveys where trade offs between area and depth need to be made. Throughout we adopt a standard cosmology with the following parameter set: kmsMpc although as our analysis is based on very local volumes only the value of the adopted Hubble constant is significant.

## 2 The calibration sample

The Sloan Digital Sky Survey Data Release 7 (SDSS DR7; Abazajian et
al. 2009) represents the final release of the original main galaxy
survey programme. As our starting point we downloaded the main galaxy
survey spectroscopic compendium
catalogs^{5}

In our selected sub-region we find 463k galaxies with mag of which 415k (90 per cent) have reliable redshifts (zconffinal). To assess the uniformity of the data over the entire region we construct the integrated number-counts to mag in 0.1 sq deg cells (as shown in Fig. 1 lower) and perform a bilinear fit to a simple flat projection of this region (i.e., we interpret Right Ascension and Declination as Cartesian but keep cell sizes as strictly to generate uniform 0.1 sq deg cells). The fit reveals a weak gradient in both Right Ascension and Declination, caused by the presence of a series of superclusters in the lower right corner of Fig. 1 (lower). This illustrates that even the SDSS is not immune to cosmic variance exhibiting a detectable large scale structure across the entire target region. Note that the density gradient extends over a scale of 100 sq deg which at the median redshift of equates to a 1Gpc scale structure. We do not explore the nature of this trend any further but do note some orthogonal indication of extremely large Gpc-scale foreground structures arising from WMAP studies (see for example Hansen, Banday & Górski 2004 and Tangen 2010).

### 2.1 Managing spectroscopic incompleteness

Spectroscopic incompleteness across the survey is also unlikely to be uniform but dependent on observing conditions (i.e., extinction, zenith angle etc) as well as the flux of the object in question. These cannot be decoupled and must be treated simultaneously. Here we divide the region into cells of regular sky coverage of 1 sq.deg (containing on average 75 galaxies). Within each cell we then determine the completeness as a function of apparent magnitude in 0.25 mag bins. Each galaxy, , within this cell with a known redshift, , is then assigned a weight, . This is the number of galaxies in some apparent magnitude bin, , divided by the number with known redshift in the same bin, i.e., . The global, or mean, weight for a cell (i.e., integrated over magnitude) is not particularly meaningful but can be obtained from the average assigned weight within each cell. Fig. 2 shows the mean weight in 1 sq.deg cells indicating that the incompleteness is not strongly clustered but distributed fairly randomly across the region except for a central swathe of exceptionally high completeness introduced by virtue of DR7 data.

### 2.2 Creating suitable test particles

Having established a sub-region and a mechanism of correcting for the spectroscopic incompleteness as a function of apparent magnitude and spatial location we now require a set of test particles to start exploring the cosmic variance within the test region. We adopt as test particles the most common galaxy in a locally observed sample, i.e., mag galaxies. These are both numerous and readily detectable to large distances. Selecting brighter systems potentially introduces additional variance as brighter sources are known to cluster more strongly (i.e., atypical; c.f., Norberg et al. 2002b), whereas low luminosity sources would restrict the depth of our volume and their intrinsic clustering properties are less well known. We adopt mag taken from the recent LF estimates of Hill et al. (2010). The r-band is used as this is the filter in which the SDSS main spectroscopic sample is selected. Fig. 3 shows the distribution of absolute magnitude versus redshift with the systems within our defined volume shown in green. In deriving the absolute magnitudes we adopted universal and corrections ( and ). The choice of and corrections are not critical nor is the assumption of universal rather than individual corrections as firstly the corrections are small at low-z (, see Hill et al. 2010), and secondly we will always compare cells constructed over identical redshift ranges. From Fig. 3 we can see that the maximum redshift we’re capable of sampling, due to the SDSS spectroscopic limit, is . Fig. 4 shows the variance along the line of sight by counting the space-density of test particles (i.e., galaxies) as a function of redshift, both cumulative (mauve curve) and differential (green data points with errorbars) distributions are shown. The differential counts show that the SDSS appears to suffer from an extreme underdensity locally () followed by a strong overdensity at () with the cumulative density well behaved from to (within 10% of the mean density in this range, see red lines on Fig. 4) after which incompleteness (i.e., traditional Malmquist bias) starts to affect the sample (r.f. Fig. 3). We therefore adopt a z range of 0.03 to 0.10 (minimum/maximum transverse co-moving scale of 2.2/7.3 Mpc/deg) which contains 100117 test galaxies distributed over an area of 5150 sq.deg with a co-moving radial length of Mpc, a lookback interval of 0.9Gyr, and a volume of Mpc. In the results that follow we will therefore only be sensitive to cosmic variance on scales below Mpc, however as cosmic variance is generally seen to be decreasing the loss of sensitivity to variance on scales greater than this volume is not expected to be a significant issue. In what follows the radial co-moving distance is typically comparable or larger (Mpc) than the transverse lengths (Mpc) and is therefore not the dominating dimension contributing to the cosmic variance. This would not be the case if the depth was quantised more discreetly, hence for what follows we have two caveats:

(1) We are not sensitive to the cosmic variance on volumes larger than Mpc

(2) We require that the shortest co-moving lengths defining the volume are tangential to the line-of-sight. In general this will be the case if the depth interval is Gyr in lookback time.

Finally Fig. 6 shows the histogram of the counts in 1 sq. deg cells for our test particles indicating a well behaved distribution somewhere between an ideal Normal and logNormal distribution — typically in previous studies a Normal distribution is assumed and we therefore follow convention. We elect to quantify cosmic variance using the simplest measure of standard deviation defined as: where , is the mean of the counts in cells for that particular cell size, and the counts in the cell. Note that this statistic will by definition include the intrinsic Poisson component which in all cases is the minor component (see blue line on Fig. 7).

## 3 Results

Having defined an extensive test region (see Fig. 5) containing test particles which reflect the underlying local structures contained in the volume, we can now sample the variance by repetitively extracting counts in fixed sized cells (truncated square-based pyramids) at random locations and repeat for increasing cell sizes. In the next two subsections we firstly explore the basic variance in square cells from 1 sq.deg to 2048 sq.deg, and secondly the variance in rectangular cells where the aspect ratio is varied from 1:1 to 1:128. In the sample extractions which follow we only sample individual test particles once, resulting in a fully uncorrelated measure of the cosmic variance but a reduction in the statistical accuracy for the larger samples and aspect ratios where fewer independent samples can be constructed.

### 3.1 The variance in square regions

Fig. 7 shows the (cosmic) variance (data points) derived from our test region. This is shown as a percentage versus volume sampled. The range of individual values measured is shown by the grey bars and the accuracy to which the mean variance is measured is shown as errorbars. We can see that the variance decreases steadily from 60 per cent for volumes of Mpc to 10 per cent at our sampling limit of Mpc. We fit a simple second order polynomial to the data and find a good ( percent) fit given by:

(1) |

where represents the cosmic variance for a volume, , as a percentage (i.e., ). This expression can be used to provide a robust estimate of the cosmic variance for any given square shaped survey given the total volume sampled.

### 3.2 Contiguous versus sparse sampling

Given the scale of cosmic variance it is worth asking whether it can be more easily overcome by multiple independent sight-lines rather than a single contiguous survey. To explore this we build up a larger survey by combining independent regions each of sq.deg to obtain a survey of area . Fig. 8 shows the results as red lines originating from the building block area (). Unsurprisingly the cosmic variance of a larger survey is significantly reduced if comprised of multiple smaller blocks rather than a single contiguous survey. For example for a survey of 32 sq.deg constructed from multiple 1 sq deg. blocks the cosmic variance is reduced from 31 per cent to 11 per cent. The cosmic variance empirically decreases with multiple sight-lines by , as one would expect if the sightlines are indeed decoupled. This holds regardless of the base survey area (i.e., independent of ). Hence for multiple sight-lines Eqn. 1 can now be modified to:

(2) |

where is the number of independent sight lines each of volume . Similarly the cosmic variance can be determined for independent regions of differing area and combined assuming Poisson statistics as long as the blocks are fully independent and non-contiguous.

### 3.3 Dependence on aspect ratio

Survey shape and in particular the aspect ratio, or window function, is likely to impact on the cosmic variance. Not all surveys will be square/circular and may have significantly different dimensions in length and width. For example the Millennium Galaxy Catalogue (Liske et al., 2003; Driver et al 2005) has a width of 0.5 deg and a length of 75 sq deg giving an aspect ratio of 1:150. Similarly the GAMA survey (Driver et al. 2009; Baldry et al. 2010) consists of 3 chunks of 4 deg by 12 deg regions for an aspect ratio of 1:3. This aspect ratio, when taken to the extreme, can help significantly to reduce cosmic variance. Observationally long thin strips are often easier to observe because of the Earth’s rotation and are more robust to cosmic variance however for the same reasons are less suitable for studies of large scale structure. Fig. 9 shows the outcome of modifying the aspect ratio by reproducing Fig. 7 for various aspect ratios ranging from 1:1 to 1:128. The cosmic variance follows an almost identical trend as for Fig. 7 but offset in amplitude. Note that the data becomes noisier for higher aspect ratios because of the limited number of independent samplings possible. For very large volumes the aspect range no longer appears to provide a gain in cosmic variance, this is because the longer tangential length approaches and exceeds the depth of the test volume resulting in a cap to the gain in cosmic variance as any volumes cosmic variance is always dominated by the two shortest lengths. In these cases the data points are shown as open symbols and the connecting lines as dotted. Finally we can incorporate the aspect ratio into Eqn. 2 by simply allowing the amplitude to vary, thus:

(3) | |||||

where is the aspect ratio, e.g., 128 for 1:128. Eqn. 3 now provides a robust estimate of the cosmic variance for the interval for almost any survey in terms of the sampling volume, , the aspect ratio of the survey window, , and the number of independent volumes, .

### 3.4 Generalising over all redshift for any survey

Extrapolating the current method beyond becomes non-trivial for a number of reasons: firstly, and foremost, the clustering signature of the population is evolving, with the galaxy population expected to be less clustered towards higher redshift, secondly one needs to consider the three dimensional volume shape, (i.e., in the earlier section we kept the redshift baseline constant at which equates to a physical co-moving distance of Mpc). The first of these issues cannot be addressed using the SDSS and, as no suitable dataset exists at higher-z, it is currently empirically intractable. However one can adopt the values from Eqn. 3 as a robust upper limit. The issue of survey shape is also itself problematic in two ways: firstly as one moves to higher-z, for a fixed survey window, the volume becomes less conic/pyramidic and more cylindrical/cuboid due to the tendency towards a nearly constant angular-diameter-transverse length relation; secondly the freedom allowed by modifying all three dimensions of the sampling volume makes the derivation of a direct empirical expression valid beyond our maximum length impossible. In a future paper (Robotham & Driver, in prep.) we will address the first of these by providing an online Cosmic Variance calculator which allows the user to specify the precise dimensions of their survey cuboid upto approximately Mpc volumes. However this may still not cover the very deep pencil beam surveys where one might wish to bin results over radial co-moving lengths greater than Mpc. However we can make two relatively simple assumptions to provide an analytical workaround.

(1) Over scales greater than Mpc one expects no correlation of structure (although note our apparent detection of weak structure over a 1Gpc linear scale across the entire SDSS). Certainly if one axis is significantly greater than the other two one expects the cosmic variance across the longer dimension to contribute least to the total variance. Hence the cosmic variance in long cuboids should scale according to Poisson statistics if the long radial length is increased/decreased over a range Mpc. i.e., a volume of Mpc should have less variance than that defined by Mpc.

(2) We equate any survey volume to a cuboid where we preserve the total volume, range, and aspect ratio and derive the appropriate transverse lengths. Thus our Mpc volume sampled by 1 sq.deg to can be equated to a cuboid of dimensions Mpc.

The first of these assumptions is reasonable assuming the radial length is always the greater of the three cuboid sides. The second can be tested by extracting equal area and equal radial length square based pyramids and cones. Adopting a constant radial length of 291Mpc we reproduce the initial results shown on Fig. 7 by now sampling our volume with cuboids (red line) rather than the original square-based pyramids (data points). To derive our cuboids we fix the total volume, the range and derive the required transverse co-moving lengths. The red curve closely follows the original data points indicating no correction is required for the change in geometric shape.

Given these two caveats one can now trivially determine an approximate cosmic variance for any survey volume. We achieve this by replacing the survey volume, , with the product of the median redshift transverse lengths (for the survey bin in question), & , and a radial depth all expressed in Mpc, then, assuming that Mpc, we find from Eqn. 3 and the caveats above that:

(4) | |||||

Note that the ratio of the transverse lengths () at the median redshift has replaced the Aspect ratio () in Eqn. 3. For conic/cylindrical surveys one can replace and with where is the transverse radius of the survey at the median redshift. This equation is strictly valid over the range to Mpc

### 3.5 Comparison to predictions from simulations

In a recent study Moster et al. (2010) derived cosmic variance values from a purely numerical/analytical route using simulations coupled with the adoption of a Halo Mass Function. We reproduce on Table. 1 columns from their Table. 5 for galaxies of mass M, the approximate turn-over point of the stellar mass function (see Baldry, Glazebrook & Driver 2008), and compare to our predictions using Eqn. 4 for the GOODs, GEMs, and COSMOS surveys. In general the values shown in Table. 3 paint a remarkably consistent picture of the cosmic variance with estimates based on the two entirely distinct methods producing remarkably consistent values for . One interesting offset is the tendency for the empirical method to produce lower cosmic variance values at higher redshift. While this initially appears counter-intuitive it can be explained in the context that while general clustering decreases as a function of redshift it presumably increases for fixed stellar mass. In other words the systems at high-z are destined to become the highly clustered superluminous population at low redshift. This is borne out if we compare our cosmic variance values to the lower mass values of Moster et al. at higher redshift. For example for the cosmic variance in a GOODS field in redshift interval 3.58 to 4.00 for stellar masses of Moster et al find consistent with our value of 30.6%. We therefore conclude that our method continues to provide a reasonable estimate of the cosmic variance at any redshift and that the discrepancy between Moster et al. and our work at high redshift is fully explained by the change in stellar mass as a function of redshift. One should therefore use our formulae if the point is sampled or the Moster et al. formulae for know stellar mass ranges.

Survey | z | Moster et al. (2010) | This work |

range | for | (%) | |

GOODS ( | 0.00—1.12 | 0.126 | 15.5% |

1.12—1.58 | 0.194 | 23.0% | |

1.58—1.99 | 0.241 | 25.1% | |

1.99—2.39 | 0.295 | 26.5% | |

2.39—2.78 | 0.365 | 28.0% | |

2.78—3.17 | 0.446 | 29.2% | |

3.17—3.58 | 0.534 | 29.7% | |

3.58—4.00 | 0.647 | 30.6% | |

GEMS ( | 0.00—1.12 | 0.098 | 11.3 |

1.12—1.58 | 0.140 | 16.0% | |

1.58—1.99 | 0.169 | 17.3% | |

1.99—2.39 | 0.203 | 18.1% | |

2.39—2.78 | 0.247 | 19.0% | |

2.78—3.17 | 0.299 | 19.7% | |

3.17—3.58 | 0.354 | 20.0% | |

3.58—4.00 | 0.425 | 20.5% | |

COSMOS ( | 0.00—1.12 | 0.057 | 6.5% |

1.12—1.58 | 0.069 | 8.5% | |

1.58—1.99 | 0.080 | 9.0% | |

1.99—2.39 | 0.093 | 9.3% | |

2.39—2.78 | 0.110 | 9.7% | |

2.78—3.17 | 0.131 | 10.0% | |

3.17—3.58 | 0.153 | 10.1% | |

3.58—4.00 | 0.181 | 10.4% |

## 4 Cosmic variance values for specific surveys.

The main purpose of this paper is to derive a simple credible path to empirically based cosmic variance estimates for recent, ongoing, and upcoming surveys. Table. 2 shows a variety of cosmic variance estimates based on Eqn. 4. The main conclusion is the need for multiple independent sightlines to reduce the impact of cosmic variance. In particular any deep ASKAP survey should have a minimum of 2 ultra-deep fields and any deep MeerKAT survey should have a minimum of 10 ultra-deep fields to probe the HI universe in intervals and keep cosmic variance below 10%. Likewise approximately 100 HST ACS or 10 HST WFC3 ultra-deep fields are required to keep cosmic variance below 10% in intervals of 1. Finally we note that the GAMA, VVDS and zCOSMOS surveys all have cosmic variance below the 10% level in intervals of of 0.1, 0.25 and 0.5 respectively.

Survey | Redshift | Area | Aspect | Cos. Var. per | Number of | Final |

Name | range | sq.deg. | ratio | pointing | pointings | Cos. Var |

MGC | 0.0 — 0.1 | 30 | 1:128 | 19% | 1 | 19% |

MGC | 0.0 — 0.2 | 30 | 1:128 | 10% | 1 | 10% |

GAMA I | 0.0 — 0.1 | 48 | 1:3 | 24% | 3 | 14% |

GAMA I | 0.1 — 0.2 | 48 | 1:3 | 15% | 3 | 8% |

GAMA I | 0.2 — 0.3 | 48 | 1:3 | 11% | 3 | 6% |

GAMA I | 0.3 — 0.4 | 48 | 1:3 | 9% | 3 | 5% |

GAMA I | 0.4 — 0.5 | 48 | 1:3 | 8% | 3 | 5% |

GAMA I | 0.0 — 0.5 | 48 | 1:3 | 5% | 3 | 3% |

ASKAP | 0.0 — 0.1 | 36 | 1:1 | 27% | 2 | 19% |

ASKAP | 0.1 — 0.2 | 36 | 1:1 | 17% | 2 | 12% |

ASKAP | 0.2 — 0.3 | 36 | 1:1 | 13% | 2 | 9% |

ASKAP | 0.3 — 0.4 | 36 | 1:1 | 11% | 2 | 8% |

ASKAP | 0.4 — 0.5 | 36 | 1:1 | 10% | 2 | 7% |

MeerKAT | 0.0 — 0.1 | 1 | 1:1 | 58% | 10 | 18% |

MeerKAT | 0.1 — 0.2 | 1 | 1:1 | 41% | 10 | 13% |

MeerKAT | 0.2 — 0.3 | 1 | 1:1 | 34% | 10 | 11% |

MeerKAT | 0.3 — 0.4 | 1 | 1:1 | 30% | 10 | 10% |

MeerKAT | 0.4 — 0.5 | 1 | 1:1 | 28% | 10 | 9% |

MeerKAT | 0.5 — 0.6 | 1 | 1:1 | 26% | 10 | 8% |

MeerKAT | 0.6 — 0.7 | 1 | 1:1 | 25% | 10 | 8% |

MeerKAT | 0.7 — 0.8 | 1 | 1:1 | 24% | 10 | 8% |

MeerKAT | 0.8 — 0.9 | 1 | 1:1 | 24% | 10 | 7% |

MeerKAT | 0.9 — 1.0 | 1 | 1:1 | 23% | 10 | 7% |

VVDS | 0.00 — 0.50 | 4 | 1:1 | 10% | 4 | 5% |

VVDS | 0.50 — 0.75 | 4 | 1:1 | 10% | 4 | 5% |

VVDS | 0.75 — 1.00 | 4 | 1:1 | 10% | 4 | 5% |

VVDS | 1.00 — 1.25 | 4 | 1:1 | 9% | 4 | 5% |

VVDS | 1.25 — 1.50 | 4 | 1:1 | 9% | 4 | 5% |

VVDS | 1.50 — 1.75 | 4 | 1:1 | 9% | 4 | 5% |

VVDS | 1.75 — 2.00 | 4 | 1:1 | 9% | 4 | 5% |

zCOSMOS | 0.0 — 0.5 | 2 | 1:1 | 12% | 1 | 12% |

zCOSMOS | 0.5 — 1.0 | 2 | 1:1 | 9% | 1 | 9% |

zCOSMOS | 1.0 — 1.5 | 2 | 1:1 | 8% | 1 | 8% |

zCOSMOS | 1.5 — 2.0 | 2 | 1:1 | 8% | 1 | 8% |

zCOSMOS | 2.0 — 2.5 | 2 | 1:1 | 8% | 1 | 8% |

HST ACS | 1.0 — 1.5 | 4.8E-05 | 1:1 | 71% | 10 | 22% |

HST ACS | 1.5 — 2.0 | 4.8E-05 | 1:1 | 77% | 10 | 24% |

HST ACS | 2.0 — 3.0 | 4.8E-05 | 1:1 | 61% | 10 | 19% |

HST ACS | 3.0 — 4.0 | 4.8E-05 | 1:1 | 70% | 10 | 22% |

HST ACS | 4.0 — 5.0 | 4.8E-05 | 1:1 | 79% | 10 | 25% |

HST ACS | 5.0 — 6.0 | 4.8E-05 | 1:1 | 87% | 10 | 28% |

HST ACS | 6.0 — 7.0 | 4.8E-05 | 1:1 | 96% | 10 | 30% |

HST ACS | 7.0 — 20.0 | 4.8E-05 | 1:1 | 39% | 10 | 12% |

HST WFC3 | 1.0 — 1.5 | 1.3E-03 | 1:1.1 | 44% | 10 | 14% |

HST WFC3 | 1.5 — 2.0 | 1.3E-03 | 1:1.1 | 46% | 10 | 15% |

HST WFC3 | 2.0 — 3.0 | 1.3E-03 | 1:1.1 | 36% | 10 | 11% |

HST WFC3 | 3.0 — 4.0 | 1.3E-03 | 1:1.1 | 41% | 10 | 13% |

HST WFC3 | 4.0 — 5.0 | 1.3E-03 | 1:1.1 | 46% | 10 | 15% |

HST WFC3 | 5.0 — 6.0 | 1.3E-03 | 1:1.1 | 51% | 10 | 16% |

HST WFC3 | 6.0 — 7.0 | 1.3E-03 | 1:1.1 | 55% | 10 | 17% |

HST WFC3 | 7.0 — 20.0 | 1.3E-03 | 1:1.1 | 22% | 10 | 7% |

## 5 Conclusions

We have derived a simple empirical expression for calculating cosmic variance for almost any extragalactic survey. The results are entirely empirical and based on resampling the SDSS DR7. The resulting equations agree extremely well with the recent numerical results by Moster et al. (2010). The two resulting equations provide corrections for robustly and for under the following caveats:

(1) Te derived cosmic variance is for mag population only and assumed not to evolve with lookback time – this is clearly incompatible with our understanding of the evolution of structure and hence beyond the derived values should be taken as indicative only.

(2) That above Mpc cosmic variance scales with radial co-moving length according to Poisson statistics.

The two equations are then used to determine cosmic variance values for a number of recent, ongoing and planned surveys.

## Acknowledgments

SPD thanks the University of Western Australia and the International Centre for Radio Astronomy Research (ICRAR) for financial support and to Profs Peter Quinn and Lister Staveley Smith, and Dr Martin Meyer for stimulating discussions on this topic during his Sabbatical stay.

Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/.

The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.

## References

Abazajian K.N., et al., 2009, ApJS, 182, 543

Baldry I. K., et al., 2010, MNRAS, 404, 86

Baldry I.K., Glazebrook K., Driver S.P., 2008, MNRAS, 388, 945

Baugh C., 2006, RPPh, 69, 3101

Berlind A.A., Weinberg D.H., 2002, ApJ, 575, 587

Cole S., Lacey C. G. Baugh C. M., Frenk C.S., 2000, MNRAS, 319, 168

Davies M, Efstathiou G., Frenk C.S., White S.D.M., 1985, ApJ, 292, 371

Driver S.P., Odewahn S.C., Echevarria L., Cohen S.H., Windhorst R.A., Phillipps S., 2003, AJ, 126, 2662

Driver S.P., Liske J., Cross N.J.G., De Propris R. Allen P.D., 2005, MNRAS, 360, 81

Driver S.P., et al., 2009, A&G, 50, 12

Hill D., Driver S.P., Cameron E., Cross N.J.G., Liske J., Robotham A., 2010, MNRAS, 404, 1215

Hansen F.K., Banday A.J., Górski K.M., 2004, 354, 641

Hopkins A.M., Beacom J.F., 2006, ApJ, 651, 142

Liske J., Lemon D.J., Driver S.P., Cross N.J.G., Couch W.J., 2003, MNRAS, 344, 307

Moster B.P., Somerville R.S., Newman J.A., Rix, H-W, 2010, ApJ, submitted (astro-ph/1001.1737)

Newman J.A., Davis M., 2002, ApJ, 564, 567

Norberg P., et al., 2002a, MNRAS, 336, 907

Norberg P., et al., 2002b, MNRAS, 332, 827

Somerville R.S., Lee K., Ferguson H.C., Gardner J.P., Moustakas LA., Giavalisco M., 2004, ApJ, 600, 171

Springel V et al., 2005, Nature, 435, 629

Szapudi I., Columbi S., 1996, ApJ, 470, 131

Tangen K., 2010, submitted, (astro-ph/0910.4164)

Trenti M., Stiavelli M., 2008, ApJ, 676, 767

Wilkins S.M. Trentham N., Hopkins A.M., MNRAS, 385, 687

Zwaan M., Meyer M., Staveley-Smith L., Webster R., 2005, MNRAS, 359, 30

### Footnotes

- pubyear: 2006
- volume: 000
- pagerange: Quantifying cosmic variance–References
- Technically the term sample variance is more correct but here we adhere to the current convention of using the term cosmic variance to describe perturbations in measurements within our Universe due to sampling size.
- http://www.sdss.org/dr7/products/spectra/getspectra.html main galaxy files.