arXiv:astro-ph/0601599v1 26 Jan 2006 hcg,54 ot li vne hcg,I 60637 IL Chicago, Avenue, Ellis South 5640 Chicago, ocoayh2 004Cao,Poland Cracow, 30-084 2, Podchorazych 12Bnign Switzerland Binningen, 4102 aai,I 60510 IL Batavia, rne al a rcsN 88003-8001 NM Cruces Las Mall, Frenger eH´tr5,108Gaaa Spain Granada, 18008 Hu´etor 50, de abig,C30A ntdKingdom; United 0HA, CB3 Cambridge, Germany Heidelberg, which models creating by observation and theory function. tween mass satellite initial the pop- of less region a sampling ulous later and observationally therefore 2004), were are Klypin we that & that subhalos Gnedin, larger, (Kravtsov, massive orig- ob- stripped or much more tidally 2002), the in in al. that formed et embedded possible inally (Stoehr fact subhalos is dark in massive It more are in 2002a). satellites al. formation served Somerville et (e.g., Benson inhibited be 2002; qualitatively. could subsystems least mass low at problem, the observations. to obstacles theory serious al. CDM most “miss- the et matching of the for Moore one discrepancy, is This 1999; problem, abundance satellite” 2002b). al. subha- ing al. observed et dark et Klypin the Benson (e.g., 1 low-mass than 1999; least epoch more dwarf at present magnitude of the predict mod- of at the scales, Such orders los large from 2 at systems. form – successful smaller M31 while of and els, accretion Way and Milky merger the like galaxies oy73,Cp on ot Africa South Town, Cape 7935, tory ailB Zucker B. Daniel Barentine 10 9 11 6 8 7 5 2 3 4 rpittpstuigL using 2021 typeset 25, Preprint July version Draft 1 be- discrepancy the resolve to aim solutions these of All address to proposed been have solutions of number A large models, (CDM) matter dark cold hierarchical In srnmshsIsiu,Universit¨at Venusstras Institut, Basel, Astronomisches eateto srnm,NwMxc tt nvriy 1320 University, State Mexico New Astronomy, of Department C/ikOsraoy 16Hg t,SnaCu,C 95064 CA Cruz, Santa St., High 1156 Observatory, UCO/Lick ot fia srnmclOsraoy OBx9 Observa- 9, Box PO Observatory, Astronomical African South nttt eAto´sc eAdlcı CI) /aioB C/Camino Andaluc´ıa Astrof´ısica (CSIC), de de Instituto pcePitOsraoy ..Bx5,Snpt M88349 NM Sunspot, 59, Box P.O. Observatory, Point Apache a-lnkIsiu ¨rAtooi,Koisul1,D-6 K¨onigstuhl 17, f¨ur Astronomie, Max-Planck-Institut urn drs:Isiueo srnm,MdnlyRoad, Madingley Astronomy, of Institute address: Current eateto srnm n srpyis nvriyof University Astrophysics, and Astronomy of Department t uoaOsraoy rcwPdggclUiest,ul University, Pedagogical Cracow Observatory, Suhora Mt. em ainlAclrtrLbrtr,PO o 500, Box P.O. Laboratory, Accelerator National Fermi Grebel siae t itneadohrpyia rpris efidthat find follo We Using s (SDSS). properties. spheroidal of Survey physical brightness dwarf other Sky surface and new central Digital distance a Sloan its X, the estimated Andromeda from of photometry discovery lar the report We hc ttedrvddsac ouu,( modulus, distance derived the at which M xrml an aeltsmyb lniu nteLclGroup. Local the Andr in and plentiful IX headings: be Andromeda Subject may of discoveries satellites faint The extremely satellite. M31 luminosity low V NRMD ,ANWDAFSHRIA AELT FM1 PHOTOM M31: OF SATELLITE SPHEROIDAL DWARF NEW A X, ANDROMEDA 9 − ∼ 5 oadBrewington Howard , o .Holtzman A. Jon , 8 .INTRODUCTION 1. . 1,2 1 a Long Dan ± lxiY Kniazev Y. Alexei , A T 0 E . ;teevle r ut oprbet hs fAdoeaI,ap a IX, Andromeda of those to comparable quite are values these 5; tl mltajv 10/10/03 v. emulateapj style X aais wr aais niiul(nrmd X nrmd X) Andromeda IX, (Andromeda individual galaxies: — dwarf galaxies: vlto oa Group Local — evolution 9 rcH ele,Jr. Neilsen, H. Eric , 6 eeA .Walterbos M. A. Rene , [email protected] 9 µ .Brinkmann J. , V, 0 1,3 ∼ ai Mart David , 26 . a arcsec mag 7 rf eso uy2,2021 25, July version Draft e7,CH- se m − 9117 ABSTRACT 8 11 ajo ihe Harvanek Michael , M tuoNitta Atsuko , . ´ inez-Delgado ) 0 6 − ∼ osac .Rockosi M. Constance , 2 utn20;Tce ta.20)adt eettar- select to 2004; and al. Ivezi´c et 2005) 2003; al. et al. Tucker 2002; et Pier al. 2005; properties et al. 2002; Stoughton Lupton et astrometric Hogg al. 1999; et and Smith Szalay 1998; & pipelines photometric Gunn, through (Lupton, processed al. measure are et 2005) Gunn to al. et 1996; Gunn 2001; al. et Fukugita mgn aai h v DSbnpse ( bandpasses SDSS five the in data Imaging cpcsre hthsmapped has that survey scopic 1998). Garnavich & Stanek 1998; Holland e.g., ue itnet 3 f73kc(( kpc 783 as- of have M31 we to work, distance paper this For a companion sumed 2005b). a al. in et presented information (Guhathakurta abundance are and therefrom, new this kinematic derived of the field to the and luminosity in satellite, in of comparable spectroscopy Keck appears IX. companion And which spheroidal this one dwarf In M31, new followup to a and of SDSS companions. using data, discovery, M31 photometric the additional or- report in we for letter, telescopes, search other select to on to observations der surroundings deeper its for and fields M31 of photometry objects SDSS point. two surface that these lowest to luminous, up of least found each the discovery, al. brightness be et its to (Willman determined of was respectively; Major Way, time Milky the the Ursa and at M31 and find of satellites to 2004b) 2005), IX; used (And al. been IX et Andromeda have the Zucker members, in SDSS Group and Local from new 2002), data two years al. et two Willman past possi- surface e.g., low it extremely (see, made as with brightnesses structures has such stellar (SDSS) surveys detect Survey to CCD ble Sky wide-field Digital that of Sloan in advent the by galaxies The hampered of are expected regime. so brightnesses do surface lumi- low to galaxy the attempts the but of function, end more faint nosity comple- the place on to A constraints be stringent would satellites. approach observational observable mentary directly fewer predict 24 n oa paetmgiueof magnitude apparent total a and DS(oke l 00 sa mgn n spectro- and imaging an is 2000) al. et (York SDSS ic h icvr fAdI,w aebe using been have we IX, And of discovery the Since . 12 .OSRAIN N AAANALYSIS DATA AND OBSERVATIONS 2. − 9 24 1,4 tpai .Snedden A. Stephanie , . rcF Bell F. Eric , 4 ilsa bouemgiueof magnitude absolute an yields 34, 9 .J Kleinman J. S. , nrmd a dereddened a has X Andromeda tlieo 3,bsdo stel- on based M31, of atellite -piaigdt ehave we data imaging w-up md ugs htsuch that suggest X omeda 7 oadG York G. Donald , 1 asWle Rix Hans-Walter , reviously-discovered 9 ue Krzesinski Jurek , 9 ∼ V galaxies: — tot m ETRY 1 / − ∼ ftesky. the of 4 M 16 8 .C. J. , ) . ,g ,i z i, r, g, u, 1, 0 1 v K. Eva , 24 = 9,10 . 47; , ; 2 Zucker et al. gets for spectroscopic followup. The SDSS scan used roughly NE – SW) at the center. The superior depth in this work was carried out on 5 October 2002 (see and resolution of the NOT data are evident in the re- Zucker et al. 2004a), processed with the same pipeline sulting CMDs of the central 2′-radius region (Fig. 3a and as Data Release 1 (Abazajian et al. 2003), and spans d), in which And X’s narrow, blue RGB clearly stands ∼ 18◦ × 2.5◦ along the major axis of M31. For deredden- out. The background-subtracted dereddened Hess dia- ing and conversion from SDSS magnitudes to V, I magni- grams (Figs.3c and f) are overplotted with the dered- tudes, we use Schlegel, Finkbeiner, & Davis (1998) and dened SDSS CMD of NGC 2419, a Galactic globular Smith et al. (2002), respectively. cluster ([Fe/H]= −2.12), shifted to a distance modulus By selecting stars in the SDSS M31 scan with colors of 24.24 (see §3); given the uncertainties in both And X’s and magnitudes characteristic of red giant branch (RGB) distance and reddening, the agreement between And X stars at the distance of M31, we discovered And IX, a new and NGC 2419 is remarkable, even hinting at the possible dwarf spheroidal satellite of M31 (Zucker et al. 2004b). presence of blue horizontal branch stars in the (g,g − r) Owing to their low metallicity (Harbeck et al. 2005), the CMD, and provides further evidence for And X’s low RGB stars in And IX are relatively blue; limiting our metallicity ([Fe/H] ∼−2). photometric criteria to select blue RGB stars at the dis- tance of M31, a large number of additional, less promi- nent overdensities were revealed. It is important to note that these stellar overdensities are not directly visible in SDSS images (see, e.g., Fig.1a), and are only detected as enhancements of the spatial density of stars within a specific color-magnitude bin (Fig.1b); in SDSS data, such enhancements may consist of a few stars spread across several square arcminutes. Deeper followup observations are therefore required in order to determine the nature of these overdensities, in particular to ascertain if they represent new M31 satellites. In the autumn of 2004, we used the 4.2m William Her- schel Telescope (WHT) and 2.5m Nordic Optical Tele- scope (NOT) on La Palma to image targets selected from these SDSS stellar overdensities. We observed at the WHT the nights of 14 and 15 September 2004, using the Prime Focus Imaging Camera (PFIP) with Harris V I filters. Unfortunately, we were only able to obtain short integrations of one of the satellite candidates un- der suitable observing conditions (seeing . 1.5′′). The Fig. 1.— as seen by SDSS and WHT: a) Sum of SDSS g,r,i images centered on And X. b) The spatial density WHT images were reduced with standard IRAF tasks, of presumed SDSS “blue” red giant stars in the same area, binned and DAOPHOTII/ALLSTAR (Stetson 1994) was used to ob- 1′ × 1′ and smoothed with a 1.′5 FWHM Gaussian. The grayscale tain V and I photometry, which we bootstrap-calibrated bar to the right indicates the significance of the spatial overden- sity, in units of the background standard deviation. c) Dered- with SDSS data. Figure 1c shows the resulting dered- ′ dened WHT CMD of stars within 2 of the center of And X, with dened (I, V − I) color-magnitude diagram (CMD), for ′ fiducial sequences overplotted for Galactic globular clusters with stars in a circular 2 -radius region centered on the loca- metallicities of (left to right) [Fe/H] = −2.2 (M15), −1.6 (M2), tion of the satellite candidate. Despite significant field −0.7 (47 Tuc), and −0.3 (NGC 6553) (Da Costa & Armandroff 1990; Sagar et al. 1999), shifted to a distance modulus of 24.47. contamination, an RGB is visible. Subtracting the area- ′ d) Dereddened WHT Hess diagram of stars within 2 of the center scaled CMD of an outer control region yields the Hess of And X, minus the area-scaled Hess diagram of a control field. diagram in Figure 1d. The narrowness, shape, and color The data are binned by 0m. 15 in (V − I) and I. The grayscale bar of the RGB indicate the presence of a spatially distinct, indicates the number of stars in each bin. Overplotted fiducials are metal-poor ([Fe/H] . −1.6) stellar population at roughly as for c). 3. PROPERTIES OF AND X the distance of M31 (the tip of the RGB (TRGB) ap- pears to be at I ∼ 20.2 ± 0.2), similar to that of And IX Using background-subtracted photometric data from (Zucker et al. 2004b; Harbeck et al. 2005). The addi- the NOT observations, we constructed a spatial den- tional similarities between the observed candidate and sity map of probable And X members in radial bins of And IX, with regard to low surface brightness (i.e., no 15′′, to which we fit two-dimensional King (1966) and obvious unresolved luminous component) and size (∼ few S´ersic (1968) profiles. We then applied the techniques of arcminutes), are further indications that the object is a Kniazev et al. (2004) to the NOT g,r,i images, masking new dwarf spheroidal (dSph) satellite of M31; we have foreground stars and subtracting a fitted sky level, bin- therefore given it the name Andromeda X (And X). ning the resulting data in circular annuli to measure the During our subsequent run at the NOT (6 – 11 Octo- central surface brightness and total magnitude in each ber 2004), we used the Mosaic Camera (MOSCA) and filter. No unresolved luminous component was detected SDSS gri filters to obtain deeper imaging of this candi- in And X, so we used the method of Zucker et al. (2004a) ′′ ′′ date with excellent seeing (typically 0. 6 − 0. 9). We used to estimate the total surface brightness and luminosity the same reduction and analysis tasks to obtain SDSS- by comparison with SDSS observations of the Pegasus calibrated gri photometry. Figure 2 shows the coaddi- dIrr; comparison to the luminosity function of the Draco tion of the NOT g,r,i imaging data, in which And X dSph (Odenkirchen et al. 2001) confirms that these es- is revealed as an elliptical stellar overdensity (oriented timates are reasonable for a range of stellar populations Andromeda X 3

Fig. 2.— And X as seen by the NOT: summed g,r,i images centered on And X. Fig. 3.— Dereddened CMDs of And X from NOT data: a) (g,g−r) CMD of stars within 2′ of the center of And X. b) (g,g−r) CMD of stars from an outer control region. c) (g,g − r) Hess (Zucker et al. 2004a). The undereddened central surface ′ brighnesses derived for And X in this way are 27m. 6±0.4, diagram of stars within 2 of the center of And X, minus the area- m m scaled Hess diagram of the control region. The data are binned by 26. 7 ± 0.3, and 26. 3 ± 0.3 in g, r, and i, respectively 0m. 15 in g and (g −r). The SDSS CMD of an inner annulus of NGC m (27. 1±0.5 in V ); the measured total magnitudes in g,r,i 2419 (a Galactic globular cluster; [Fe/H]= −2.12), dereddened and m m m m shifted to a distance modulus of 24.24 (see §3), is overplotted as are 17. 1 ± 0.4, 16. 1 ± 0.3, and 15. 8 ± 0.3 (16. 5 ± 0.5 in ′ V ). Table 1 lists the derived structural parameters and dots. d) (i,g − i) CMD of stars within 2 of the center of And X. e) (i,g − i) CMD of stars from the outer control region. f) (i,g − i) dereddened V -band magnitudes for And IX and And X. Hess diagram of the central region, minus the scaled control region, Determining a more precise distance to And X with binned by 0m. 15 in i and 0m. 10 in (g − i). The CMD of NGC the data in hand is problematic, due to the sparseness 2419 is overplotted, as in c). The dashed line in d) - f) shows the of And X’s CMD in the region of the TRGB. We closely approximate detection limit of our data; the grayscale bar indicates the number of stars per bin in c) and f). examined the field-star-subtracted I-band stellar lumi- nosity functions derived from both WHT and NOT data within the central 2′ radius of And X; however, even 10 convolution with an edge detection kernel did not yield And VII NGC 147 an unambiguous TRGB magnitude, with the measured And V ITRGB ranging from approximately 20.14 to 20.36, af- NGC 185 5 ter correction for a foreground reddening of AI = 0.24. And X

Assuming a metallicity of [Fe/H]= −2.2 and a TRGB And IX color of (V − I)TRGB =1.2 (as for And IX; Zucker et al. NGC 205 2004b), from the calibration of Da Costa & Armandroff 0 (1990) we obtain the relation (m − M) = I +3.98, M 32 0 TRGB Y (degrees) And I yielding a distance modulus range of (m − M)0 ∼ 24.12 And III to 24.34, with uncertainties on the order of 0.1 mag (667 −5 kpc to 738 kpc, ±30 to 35 kpc). At an M31 distance of 783 kpc, the angular separation of And X from the And II ◦ center of M31, ∼ 5.5, gives a projected separation of M 33 75 kpc; considering the above range of And X distance −10 50 kpc estimates puts And X at 87 kpc to 138 kpc from the center of M31, well within M31’s estimated virial ra- 10 5 0 −5 −10 dius (e.g., Benson et al. 2002a) and the radius at which X (degrees) M31 halo stars have been found (Guhathakurta et al. Fig. 4.— Projected distribution of M31’s nearest satellites (after 2005a). Thus And X is in all likelihood a bound satel- Ferguson et al. 2002), with the addition of And IX and And X. The lite of M31, with a dereddened integrated magnitude of inner ellipse indicates the rough optical size of M31’s disk, and the Vtot ∼ 16.1, translating to an of outer ellipse the approximate extent of the INT survey (see, e.g., m Ferguson et al. 2002). The coverage of the SDSS M31 scan is shown MV ∼−8.1 ± 0. 5. by the long stripe along the major axis of M31. 4. DISCUSSION As shown in Table 1, And X is comparable in size, nosity dSph companion of M31, a result confirmed by surface brightness, and to And IX, kinematic data (Guhathakurta et al. 2005b). The ear- and, despite the uncertainty in its distance, And X ap- lier discovery of And IX raised the question of whether pears to be somewhat lower in luminosity (MV ∼−8.1). such low luminosity galaxies were a rarity in the Local We conclude that And X is a new extremely low lumi- Group, or if And IX represented the tip of an iceberg, 4 Zucker et al.

in reconciling the predictions of current CDM galaxy for- Table 1. Properties of And IX and And X mation models with the observable local .

Parametera And IX And X R.A.(J2000.0) 00 52 52.8 01 06 33.7 DBZ acknowledges support from a National Science DEC.(J2000.0) +43 12 00 +44 48 15.8 ′ b ′ Foundation International Postdoctoral Fellowship. EFB King rc 1.33 1.33 ′ b ′ acknowledges the financial support provided through the King rt 5.9 7.2 S´ersic n 1.6b 1.2 European Community’s Human Potential Program un- m m AV 0. 26 0. 41 der contract HPRN-CT-2002-00316, SISCO. We thank m m m µ0,V 26. 77 ± 0. 09 26.7 ± 0. 5 m m m M. Schirmer and T. Pursimo for assistance in reducing Vtot 16. 17 ± 0. 06 16.1 ± 0. 5 m m m m m WHT and NOT data, respectively. This work utilized (m−M)0 24. 48 ± 0. 20 24. 12 − 24. 34 ± 0. 10 m m Mtot,V −8. 3 −8. 1 observations made with the Nordic Optical Telescope, aSurface brightnesses and integrated magnitudes are corrected for operated on the island of La Palma jointly by Denmark, the mean Galactic foreground reddenings, AV, shown. Unless oth- Finland, Iceland, Norway, and Sweden, in the Spanish erwise noted, values for And IX are from Zucker et al. (2004b). Observatorio del Roque de los Muchachos of the Insti- b From Harbeck et al. (2005). tuto de Astrof´isica de Canarias; these observations were funded by the Optical Infrared Coordination Network (OPTICON), a major international collaboration sup- one of a large population of faint dwarf satellites that ported by the Research Infrastructures Programme of the have remained undetected owing to their extremely low European Commission’s Sixth Framework Programme. surface brightnesses and low luminosities. And X’s dis- Funding for the creation and distribution of the SDSS covery, in conjunction with the recent discovery of the Archive has been provided by the Alfred P. Sloan Foun- Ursa Major dwarf (Willman et al. 2005), strongly sug- dation, the Participating Institutions, the National Aero- gests that the latter scenario is closer to the truth – that nautics and Space Administration, the National Sci- M31 and the Milky Way may in fact have a large number ence Foundation, the U.S. Department of Energy, the of low luminosity satellites. Japanese Monbukagakusho, and the Max Planck Soci- Further analysis of SDSS data is therefore likely to ety. The SDSS Web site is http://www.sdss.org/. yield more such dwarfs. To test this scenario, we are The SDSS is managed by the Astrophysical Research following up other stellar overdensities in the SDSS M31 Consortium (ARC) for the Participating Institutions. scan with deeper observations on other telescopes. In- The Participating Institutions are The University of triguingly, in projection And IX and And X lie relatively Chicago, Fermilab, the Institute for Advanced Study, close to each other and to the major axis of M31, north- the Japan Participation Group, The Johns Hopkins Uni- east of M31’s center (Fig. 4); as more new M31 satellites versity, the Korean Scientist Group, Los Alamos Na- are found (as appears quite likely), we will be able to tional Laboratory, the Max-Planck-Institute for Astron- place constraints on the spatial distribution of M31 satel- omy (MPIA), the Max-Planck-Institute for Astrophysics lites and their (an)isotropy (e.g., Koch & Grebel 2005), (MPA), New Mexico State University, University of as well as on the satellite luminosity function. Thus the Pittsburgh, University of Portsmouth, Princeton Uni- discovery and characterization of extremely low luminos- versity, the United States Naval Observatory, and the ity galaxies like And IX and And X play a critical role University of Washington.

REFERENCES Abazajian, K., et al. 2003, AJ, 126, 2081 Lupton, R., Gunn, J., & Szalay, A. 1999, AJ, 118, 1406 Benson, A. J., Lacey, C. G., Baugh, C. M., Cole, S., & Frenk, C. S. Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Stadel, 2002a, MNRAS, 333, 156 J., & Tozzi, P. 1999, ApJ, 524, 19 Benson, A. J., Frenk, C. S., Lacey, C. G., Baugh, C. M., & Cole, Odenkirchen, M. et al. 2001, AJ, 122, 2538 S. 2002b, MNRAS, 333, 177 Pier, J.R., Munn, J.A., Hindsley, R.B, Hennessy, G.S., Kent, S.M., Da Costa, G. S., & Armandroff, T. E. 1990, AJ, 100, 162 Lupton, R.H., & Ivezic, Z. 2003, AJ, 125, 1559 Ferguson, A. M. N., Irwin, M. J., Ibata, R. A., Lewis, G. F., & Sagar, R., Subramaniam, A., Richtler, T., & Grebel, E. K. 1999, Tanvir, N. R. 2002, AJ, 124, 1452 A&AS, 135, 391 Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 & Schneider, D. P. 1996, AJ, 111, 1748 S´ersic, J. 1968, Atlas de galaxias australes, C´ordova: Observatorio Guhathakurta, P., et al. 2005a, Nature, submitted Astron´omico Guhathakurta, P. et al. 2005b, in preparation Smith, J.A. et al. 2002, AJ, 123, 2121 Gunn, J.E. et al. 1998, AJ, 116, 3040 Somerville, R. S. 2002, ApJ, 572, L23 Gunn, J.E. et al. 2005, AJ, submitted Stanek, K., & Garnavich, P. 1998, ApJ, 503, L131 Harbeck, D., Gallagher, J., Grebel, E., Koch, A., & Zucker, D. Stetson, P. B. 1994, PASP, 106, 250 2005, ApJ, 623, 159 Stoehr, F., White, S. D. M., Tormen, G., & Springel, V. 2002, Hogg, D.W., Finkbeiner, D.P., Schlegel, D.J., and Gunn, J.E. 2001, MNRAS, 335, L84 AJ, 122, 2129 Stoughton, C. et al. 2002, AJ, 123, 485 Holland, S. 1998, AJ, 115, 1916 Tucker, D., et al. 2005, AJ, submitted Ivezi´c, Z.ˇ et al., AN, 325, 583 van den Bergh, S. 1999, A&A Rev., 9, 273 King, I. R. 1966, AJ, 71, 64 Willman, B., Dalcanton, J., Ivezi´c, Z.,ˇ Jackson, T., Lupton, R., Kniazev, A. Y., Grebel, E. K., Pustilnik, S. A., Pramskij, A. G., Brinkmann, J., Hennessy, G., & Hindsley, R. 2002, AJ, 123, 848 Kniazeva, T. F., Prada, F., & Harbeck, D. 2004, AJ, 127, 704 Willman, B., et al. 2005, ApJ, 626, L85 Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, York, D.G. et al. 2000, AJ, 120, 1579 ApJ, 522, 82 Zucker et al. 2004a, ApJ, 612, L117 Koch, A., & Grebel, E. 2005, AJ, in press (astro-ph/0509258) Zucker et al. 2004b, ApJ, 612, L121 Kravtsov, A., Gnedin, O., & Klypin, A. 2004, ApJ, 609, 482 Lupton, R. H. 2005, AJ, submitted