This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

We describe the design and performance of the near-infrared (1.51--1.70 micron), fiber-fed, multi-object (300 fibers), high resolution (R = lambda/delta lambda ~ 22,500) spectrograph built for the Apache Point Observatory Galactic Evolution Experiment (APOGEE). APOGEE is a survey of ~ 10^5 red giant stars that systematically sampled all Milky Way populations (bulge, disk, and halo) to study the Galaxy's chemical and kinematical history. It was part of the Sloan Digital Sky Survey III (SDSS-III) from 2011 -- 2014 using the 2.5 m Sloan Foundation Telescope at Apache Point Observatory, New Mexico. The APOGEE-2 survey is now using the spectrograph as part of SDSS-IV, as well as a second spectrograph, a close copy of the first, operating at the 2.5 m du Pont Telescope at Las Campanas Observatory in Chile. Although several fiber-fed, multi-object, high resolution spectrographs have been built for visual wavelength spectroscopy, the APOGEE spectrograph is one of the first such instruments built for observations in the near-infrared. The instrument's successful development was enabled by several key innovations, including a "gang connector" to allow simultaneous connections of 300 fibers; hermetically sealed feedthroughs to allow fibers to pass through the cryostat wall continuously; the first cryogenically deployed mosaic volume phase holographic grating; and a large refractive camera that includes mono-crystalline silicon and fused silica elements with diameters as large as ~ 400 mm. This paper contains a comprehensive description of all aspects of the instrument including the fiber system, optics and opto-mechanics, detector arrays, mechanics and cryogenics, instrument control, calibration system, optical performance and stability, lessons learned, and design changes for the second instrument.

We present a search for new members of the 300 km s$^{-1}$ stream (300S) near the dwarf galaxy Segue 1 using wide-field survey data. We identify 11 previously unknown bright stream members in the APOGEE-2 and SEGUE-1 and 2 spectroscopic surveys. Based on the spatial distribution of the high-velocity stars, we confirm for the first time that this kinematic structure is associated with a 24$^{\circ}$-long stream seen in SDSS and Pan-STARRS imaging data. The 300S stars display a metallicity range of $-2.17 < {\rm [Fe/H]} < -1.24$, with an intrinsic dispersion of 0.21$_{-0.09}^{+0.12}$ dex. They also have chemical abundance patterns similar to those of Local Group dwarf galaxies, as well as that of the Milky Way halo. Using the open-source code galpy to model the orbit of the stream, we find that the progenitor of the stream passed perigalacticon about 70 Myr ago, with a closest approach to the Galactic Center of about 4.1 kpc. Using Pan-STARRS DR1 data, we obtain an integrated stream luminosity of $4 \times 10^3$ L$_{\odot}$. We conclude that the progenitor of the stream was a dwarf galaxy that is probably similar to the satellites that were accreted to build the present-day Milky Way halo.

We report the detection of diffuse starlight in an extragalactic HI cloud in the nearby Leo~I galaxy group. We detect the source, dubbed BST1047+1156, in both broadband optical and GALEX ultraviolet light. Spanning ~ 2 kpc in radius, it has a peak surface brightness of $\mu_B$=28.8 mag/arcsec$^2$, making it the lowest surface brightness object ever detected via integrated light. Although the object is extremely gas-rich, with a gas fraction of $f_g$=0.99, its peak HI column density is well below levels where star formation is typically observed in galaxies. Nonetheless, BST1047+1156 shows evidence for young stellar populations: along with the detected UV emission, the object is extremely blue, with B-V=0.14$\pm$0.09. The object sports two tidal tails and is found embedded within diffuse gas connecting the spiral galaxy M96 to the group's extended HI Leo Ring. The nature of BST1047+1156 is unclear. It could be a disrupting tidal dwarf, recently spawned from star formation triggered in the Leo I Group's tidal debris. Alternatively, the object may have been a pre-existing galaxy --- the most extreme example of a gas-rich field LSB known to date --- which had a recent burst of star formation triggered by encounters in the group environment.

We use deep Hubble Space Telescope imaging in the outskirts of the nearby spiral M101 to study stellar populations in the galaxy's outer disk and halo. Our ACS field lies 17.6 arcmin (36 kpc) from the center of M101 and targets the blue "NE Plume" of M101's outer disk, while the parallel WFC3 field lies at a distance of 23.3 arcmin (47 kpc) to sample the galaxy's stellar halo. The WFC3 halo field shows a well-defined red giant branch characterized by low metallicity ([M/H]=-1.7 $\pm$ 0.2), with no evidence of young stellar populations. In contrast, the ACS disk field shows multiple stellar populations, including a young main sequence, blue and red helium burning stars, and old RGB and AGB populations. The mean metallicity of these disk stars is quite low: [M/H]=-1.3 $\pm$ 0.2 for the RGB population, and -1.15 $\pm$ 0.2 for the younger helium burning sequences. Of particular interest is a bunching of stars along the BHeB sequence, indicative of an evolving cohort of massive young stars. We show that the young stellar populations in this field are well-described by a decaying burst of star formation that peaked ~ 300-400 Myr ago, along with a more extended star formation history to produce the older RGB and AGB populations. These results confirm and extend the results from our previous deep surface photometry of M101's outer disk, providing an important cross-check on stellar population studies using resolved stellar populations versus integrated light photometry. We discuss our results in the context of halo formation models and the interaction history of M101 and its companions.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present the discovery of a vast cloud of ionized gas 13$^{\prime}$ (32 kpc) north of the interacting system M51. We detected this cloud via deep narrow-band imaging with the Burrell Schmidt Telescope, where it appears as an extended, diffuse H$\alpha$-emitting feature with no embedded compact regions. The Cloud spans $\sim$10$^{\prime}\times$3$^{\prime}$ (25$\times$7.5 kpc) in size and has no stellar counterpart; comparisons with our previous deep broadband imaging show no detected continuum light to a limit of $\mu_{\rm lim, B} \sim$30 mag arcsec$^{-2}$. WIYN SparsePak observations confirm the cloud's kinematic association with M51, and the high NII/H$\alpha$, SII/H$\alpha$, and OI/H$\alpha$ line ratios we measure imply a hard ionization source such as AGN photoionization or shock heating rather than photoionization due to young stars. Given the strong NII emission, we infer roughly solar metallicity for the cloud, ruling out an origin due to infall of primordial gas. Instead we favor models where the gas has been expelled from the inner regions of the M51 system due to tidal stripping or starburst/AGN winds and has been subsequently ionized either by shocks or a fading AGN. This latter scenario raises the intriguing possibility that M51 may be the nearest example of an AGN fossil nebula or light echo, akin to the famous "Hanny's Voorwerp" in the IC 2497 system.

We present a multiwavelength study of star formation within the nearby M101 Group, including new deep H$\alpha$ imaging of M101 and its two companions. We perform a statistical analysis of the H$\alpha$ to FUV flux ratios in HII regions located in three different environments: M101's inner disk, M101's outer disk, and M101's lower mass companion galaxy NGC~5474. We find that, once bulk radial trends in extinction are taken into account, both the median and scatter in F$_{\rm H\alpha}/$F$_{\rm FUV}$ in HII regions are invariant across all of these environments. Also, using Starburst99 models, we are able to qualitatively reproduce the distributions of F$_{\rm H\alpha}/$F$_{\rm FUV}$ throughout these different environments using a standard Kroupa IMF, hence we find no need to invoke truncations in the upper mass end of the IMF to explain the young star-forming regions in the M101 Group even at extremely low surface density. This implies that star formation in low density environments differs from star formation in high density environments only by intensity and not by cloud-to-cloud physics.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

APOGEE-2 is a high-resolution, near-infrared spectroscopic survey observing roughly 300,000 stars across the entire sky. It is the successor to APOGEE and is part of the Sloan Digital Sky Survey IV (SDSS-IV). APOGEE-2 is expanding upon APOGEE's goals of addressing critical questions of stellar astrophysics, stellar populations, and Galactic chemodynamical evolution using (1) an enhanced set of target types and (2) a second spectrograph at Las Campanas Observatory in Chile. APOGEE-2 is targeting red giant branch (RGB) and red clump (RC) stars, RR Lyrae, low-mass dwarf stars, young stellar objects, and numerous other Milky Way and Local Group sources across the entire sky from both hemispheres. In this paper, we describe the APOGEE-2 observational design, target selection catalogs and algorithms, and the targeting-related documentation included in the SDSS data releases.

We report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We present a deep imaging survey of the Virgo Cluster, designed to study the connection between cluster galaxies and Virgo's diffuse intracluster light (ICL). Our observations span roughly 16 square degrees and reach a 3-sigma depth of mu(B)=29.5 and mu(V)=28.5 mag/arcsec^2. At these depths, the limiting systematic uncertainties are astrophysical: scattered starlight from foreground Galactic dust, and variations in faint background sources. The dust-scattered starlight is well-traced by deep far-infrared imaging, making it distinguishable from true Virgo diffuse light. Our imaging maps the Virgo core around M87 and the adjacent M86/M84 region, in subcluster B around M49, and in the more distant W' cloud around NGC 4365. Most of the detected ICL is found in the Virgo core and within the W' cloud, with little evidence for extensive ICL in subcluster B. The large amount of diffuse light seen in the infalling W' cloud likely illustrates the importance of the group environment for generating ICL. The bulk of the detected ICL is fairly red (B-V=0.7-0.9), indicative of old stellar populations. We estimate a total Virgo ICL fraction of 7-15%, somewhat smaller than expected for massive, evolved clusters, suggesting that Virgo is still growing its ICL component. We trace M87's extremely boxy halo out to ~ 150 kpc, and show that the current stripping rate of low luminosity galaxies is insufficient to have built M87's outer halo over a Hubble time. Finally, we identify another large ultra-diffuse galaxy in Virgo, likely in the process of being shredded by the cluster tidal field.

With the advent of the space missions CoRoT and Kepler, it has become feasible to determine precise asteroseismic masses and ages for large samples of red-giant stars. In this paper, we present the CoRoGEE dataset -- obtained from CoRoT lightcurves for 606 red giant stars in two fields of the Galactic disc which have been co-observed for an ancillary project of APOGEE. We have used the Bayesian parameter estimation code PARAM to calculate distances, extinctions, masses, and ages for these stars in a homogeneous analysis, resulting in relative statistical uncertainties of $\sim2\%$ in distance, $\sim4\%$ in radius, $\sim9\%$ in mass and $\sim25\%$ in age. We also assess systematic age uncertainties due to different input physics and mass loss. We discuss the correlation between ages and chemical abundance patterns of field stars over a large radial range of the Milky Way's disc (5 kpc $<R_{\rm Gal}<$ 14 kpc), focussing on the [$\alpha$/Fe]-[Fe/H]-age plane in five radial bins of the Galactic disc. We find an overall agreement with the expectations of chemical-evolution models computed before the present data were available, especially for the outer regions. However, our data also indicate that a significant fraction of stars now observed near and beyond the Solar Neighbourhood migrated from inner regions. Mock CoRoGEE observations of a chemo-dynamical Milky Way disc model show that the number of high-metallicity stars in the outer disc is too high to be accounted for even by the strong radial mixing present in the model. The mock observations also reveal that the age distribution of the [$\alpha$/Fe]-enhanced sequence in the CoRoGEE inner-disc field is much broader than expected from a combination of radial mixing and observational errors. We suggest that a thick disc/bulge component that formed stars for more than 3 Gyr may account for these discrepancies.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

We present deep, wide-field imaging of the M96 galaxy group (also known as the Leo I Group). Down to surface brightness limits of $\mu_{B}=30.1$ and $\mu_{V}=29.5$, we find no diffuse, large-scale optical counterpart to the "Leo Ring", an extended HI ring surrounding the central elliptical M105 (NGC 3379). However, we do find a number of extremely low surface-brightness ($\mu_{B} \gtrsim 29$) small-scale streamlike features, possibly tidal in origin, two of which may be associated with the Ring. In addition we present detailed surface photometry of each of the group's most massive members -- M105, NGC 3384, M96 (NGC 3368), and M95 (NGC 3351) -- out to large radius and low surface brightness, where we search for signatures of interaction and accretion events. We find that the outer isophotes of both M105 and M95 appear almost completely undisturbed, in contrast to NGC 3384 which shows a system of diffuse shells indicative of a recent minor merger. We also find photometric evidence that M96 is accreting gas from the HI ring, in agreement with HI data. In general, however, interaction signatures in the M96 Group are extremely subtle for a group environment, and provide some tension with interaction scenarios for the formation of the Leo HI Ring. The lack of a significant component of diffuse intragroup starlight in the M96 Group is consistent with its status as a loose galaxy group in which encounters are relatively mild and infrequent.

We present results from deep, wide-field surface photometry of three nearby (D=4--7 Mpc) spiral galaxies: M94 (NGC 4736), M64 (NGC 4826), and M106 (NGC 4258). Our imaging reaches limiting surface brightnesses of $\mu_{B} \sim$ 28 -- 30 mag arcsec$^{-2}$ and probes colors down to $\mu_{B} \sim$ 27.5 mag arcsec$^{-2}$. We compare our broadband optical data to available ultraviolet and high column-density HI data to better constrain the star forming history and stellar populations of the outermost parts of each galaxy's disk. Each galaxy has a well-defined radius beyond which little star formation occurs and the disk light appears both azimuthally smooth and red in color, suggestive of old, well-mixed stellar populations. Given the lack of ongoing star formation or blue stellar populations in these galaxies' outer disks, the most likely mechanisms for their formation are dynamical processes such as disk heating or radial migration, rather than inside-out growth of the disks. This is also implied by the similarity in outer disk properties despite each galaxy showing distinct levels of environmental influence, from a purely isolated galaxy (M94) to one experiencing weak tidal perturbations from its satellite galaxies (M106) to a galaxy recovering from a recent merger (M64), suggesting that a variety of evolutionary histories can yield similar outer disk structure. While this suggests a common secular mechanism for outer disk formation, the large extent of these smooth, red stellar populations---which reach several disk scalelengths beyond the galaxies' spiral structure---may challenge models of radial migration given the lack of any non-axisymmetric forcing at such large radii.

We statistically quantify the amount of substructure in the Milky Way stellar halo using a sample of 4568 halo K giant stars at Galactocentric distances ranging over 5-125 kpc. These stars have been selected photometrically and confirmed spectroscopically as K giants from the Sloan Digital Sky Survey's SEGUE project. Using a position-velocity clustering estimator (the 4distance) and a model of a smooth stellar halo, we quantify the amount of substructure in the halo, divided by distance and metallicity. Overall, we find that the halo as a whole is highly structured. We also confirm earlier work using BHB stars which showed that there is an increasing amount of substructure with increasing Galactocentric radius, and additionally find that the amount of substructure in the halo increases with increasing metallicity. Comparing to resampled BHB stars, we find that K giants and BHBs have similar amounts of substructure over equivalent ranges of Galactocentric radius. Using a friends-of-friends algorithm to identify members of individual groups, we find that a large fraction (~33%) of grouped stars are associated with Sgr, and identify stars belonging to other halo star streams: the Orphan Stream, the Cetus Polar Stream, and others, including previously unknown substructures. A large fraction of sample K giants (more than 50%) are not grouped into any substructure. We find also that the Sgr stream strongly dominates groups in the outer halo for all except the most metal-poor stars, and suggest that this is the source of the increase of substructure with Galactocentric radius and metallicity.

We present griz observations for the clusters M92, M13 and NGC 6791 and gr photometry for M71, Be 29 and NGC 7789. In addition we present new membership identifications for all these clusters, which have been observed spectroscopically as calibrators for the SDSS/SEGUE survey; this paper focuses in particular on the red giant branch stars in the clusters. In a number of cases, these giants were too bright to be observed in the normal SDSS survey operations, and we describe the procedure used to obtain spectra for these stars. For M71, also present a new variable reddening map and a new fiducial for the gr giant branch. For NGC 7789, we derived a transformation from Teff to g-r for giants of near solar abundance, using IRFM Teff measures of stars with good ugriz and 2MASS photometry and SEGUE spectra. The result of our analysis is a robust list of known cluster members with correctly dereddened and (if needed) transformed gr photometry for crucial calibration efforts for SDSS and SEGUE.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), one of three core programs in the Sloan Digital Sky Survey-IV (SDSS-IV), is an integral-field spectroscopic (IFS) survey of roughly 10,000 nearby galaxies. It employs dithered observations using 17 hexagonal bundles of 2 arcsec fibers to obtain resolved spectroscopy over a wide wavelength range of 3,600-10,300A. To map the internal variations within each galaxy, we need to perform accurate {\it spectral surface photometry}, which is to calibrate the specific intensity at every spatial location sampled by each individual aperture element of the integral field unit. The calibration must correct only for the flux loss due to atmospheric throughput and the instrument response, but not for losses due to the finite geometry of the fiber aperture. This requires the use of standard star measurements to strictly separate these two flux loss factors (throughput versus geometry), a difficult challenge with standard single-fiber spectroscopy techniques due to various practical limitations. Therefore, we developed a technique for spectral surface photometry using multiple small fiber-bundles targeting standard stars simultaneously with galaxy observations. We discuss the principles of our approach and how they compare to previous efforts, and we demonstrate the precision and accuracy achieved. MaNGA's relative calibration between the wavelengths of H$\alpha$ and H$\beta$ has a root-mean-square (RMS) of 1.7%, while that between [NII] $\lambda$6583A and [OII] $\lambda$3727A has an RMS of 4.7%. Using extinction-corrected star formation rates and gas-phase metallicities as an illustration, this level of precision guarantees that flux calibration errors will be sub-dominant when estimating these quantities. The absolute calibration is better than 5% for more than 89% of MaNGA's wavelength range.

The Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey (SDSS) III project, has provided the largest survey of galaxy redshifts available to date, in terms of both the number of galaxy redshifts measured by a single survey, and the effective cosmological volume covered. Key to analysing the clustering of these data to provide cosmological measurements is understanding the detailed properties of this sample. Potential issues include variations in the target catalogue caused by changes either in the targeting algorithm or properties of the data used, the pattern of spectroscopic observations, the spatial distribution of targets for which redshifts were not obtained, and variations in the target sky density due to observational systematics. We document here the target selection algorithms used to create the galaxy samples that comprise BOSS. We also present the algorithms used to create large scale structure catalogues for the final Data Release (DR12) samples and the associated random catalogues that quantify the survey mask. The algorithms are an evolution of those used by the BOSS team to construct catalogues from earlier data, and have been designed to accurately quantify the galaxy sample. The code used, designated MKSAMPLE, is released with this paper.

Astrophysical sources are now observed by many different instruments at different wavelengths, from radio to high-energy gamma-rays, with an unprecedented quality. Putting all these data together to form a coherent view, however, is a very difficult task. Each instrument has its own data format, software and analysis procedure, which are difficult to combine. It is for example very challenging to perform a broadband fit of the energy spectrum of the source. The Multi-Mission Maximum Likelihood framework (3ML) aims to solve this issue, providing a common framework which allows for a coherent modeling of sources using all the available data, independent of their origin. At the same time, thanks to its architecture based on plug-ins, 3ML uses the existing official software of each instrument for the corresponding data in a way which is transparent to the user. 3ML is based on the likelihood formalism, in which a model summarizing our knowledge about a particular region of the sky is convolved with the instrument response and compared to the corresponding data. The user can choose between a frequentist analysis, and a Bayesian analysis. In the former, parameters of the model are optimized in order to obtain the best match to the data (i.e., the maximum of the likelihood). In the latter, the priors specified by the user are used to build the posterior distribution, which is then sampled with Markov Chain Monte Carlo or Multinest. Our implementation of this idea is very flexible, allowing the study of point sources as well as extended sources with arbitrary spectra. We will review the problem we aim to solve, the 3ML concepts and its innovative potential.

We report the discovery of three large (R29 >~ 1 arcminute) extremely low surface brightness (mu_(V,0) ~ 27.0) galaxies identified using our deep, wide-field imaging of the Virgo Cluster from the Burrell Schmidt telescope. Complementary data from the Next Generation Virgo Cluster Survey do not resolve red giant branch stars in these objects down to i=24, yielding a lower distance limit of 2.5 Mpc. At the Virgo distance, these objects have half-light radii 3-10 kpc and luminosities L_V=2-9x10^7 Lsun. These galaxies are comparable in size but lower in surface brightness than the large ultradiffuse LSB galaxies recently identified in the Coma cluster, and are located well within Virgo's virial radius; two are projected directly on the cluster core. One object appears to be a nucleated LSB in the process of being tidally stripped to form a new Virgo ultracompact dwarf galaxy. The others show no sign of tidal disruption, despite the fact that such objects should be most vulnerable to tidal destruction in the cluster environment. The relative proximity of Virgo makes these objects amenable to detailed studies of their structural properties and stellar populations. They thus provide an important new window onto the connection between cluster environment and galaxy evolution at the extremes.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We report the discovery of a group of apparently young CoRoT red-giant stars exhibiting enhanced [alpha/Fe] abundance ratios (as determined from APOGEE spectra) with respect to Solar values. Their existence is not explained by standard chemical evolution models of the Milky Way, and shows that the chemical-enrichment history of the Galactic disc is more complex. We find similar stars in previously published samples for which isochrone-ages could be robustly obtained, although in smaller relative numbers, which could explain why these stars have not received prior attention. The young [alpha/Fe]-rich stars are much more numerous in the CoRoT-APOGEE (CoRoGEE) inner-field sample than in any other high-resolution sample available at present, as only CoRoGEE can explore the inner-disc regions and provide ages for its field stars. The kinematic properties of the young [$\alpha$/Fe]-rich stars are not clearly thick-disc like, despite their rather large distances from the Galactic mid-plane. Our tentative interpretation of these and previous intriguing observations in the Milky Way is that these stars were formed close to the end of the Galactic bar, near corotation -- a region where gas can be kept inert for longer times, compared to other regions shocked more frequently by the passage of spiral arms. Moreover, that is where the mass return from older inner-disc stellar generations should be maximal (according to an inside-out disc-formation scenario), further diluting the in-situ gas. Other possibilities to explain these observations (e.g., a recent gas-accretion event) are also discussed.

Using a sample of 69,919 red giants from the SDSS-III/APOGEE Data Release 12, we measure the distribution of stars in the [$\alpha$/Fe] vs. [Fe/H] plane and the metallicity distribution functions (MDF) across an unprecedented volume of the Milky Way disk, with radius $3<R<15$ kpc and height $|z|<2$ kpc. Stars in the inner disk ($R<5$ kpc) lie along a single track in [$\alpha$/Fe] vs. [Fe/H], starting with $\alpha$-enhanced, metal-poor stars and ending at [$\alpha$/Fe]$\sim0$ and [Fe/H]$\sim+0.4$. At larger radii we find two distinct sequences in [$\alpha$/Fe] vs. [Fe/H] space, with a roughly solar-$\alpha$ sequence that spans a decade in metallicity and a high-$\alpha$ sequence that merges with the low-$\alpha$ sequence at super-solar [Fe/H]. The location of the high-$\alpha$ sequence is nearly constant across the disk, however there are very few high-$\alpha$ stars at $R>11$ kpc. The peak of the midplane MDF shifts to lower metallicity at larger $R$, reflecting the Galactic metallicity gradient. Most strikingly, the shape of the midplane MDF changes systematically with radius, with a negatively skewed distribution at $3<R<7$ kpc, to a roughly Gaussian distribution at the solar annulus, to a positively skewed shape in the outer Galaxy. For stars with $|z|>1$ kpc or [$\alpha$/Fe]$>0.18$, the MDF shows little dependence on $R$. The positive skewness of the outer disk MDF may be a signature of radial migration; we show that blurring of stellar populations by orbital eccentricities is not enough to explain the reversal of MDF shape but a simple model of radial migration can do so.

We present deep, wide-field imaging of the M51 system using CWRU's Burrell Schmidt telescope at KPNO to study the faint tidal features that constrain its interaction history. Our images trace M51's tidal morphology down to a limiting surface brightness of $\mu_{B,lim}\sim $30 mag arcsec$^{-2}$, and provide accurate colors ($\sigma_{B-V} < 0.1$) down to $\mu_B\sim 28$. We identify two new tidal streams in the system (the South and Northeast Plumes) with surface brightnesses of $\mu_B =29$ and luminosities of $\sim 10^6 L_{\odot,B}$. While the Northeast Plume may be a faint outer extension of the tidal "crown" north of NGC 5195 (M51b), the South Plume has no analogue in any existing M51 simulation and may represent a distinct tidal stream or disrupted dwarf galaxy. We also trace the extremely diffuse Northwest Plume out to a total extent of 20' (43 kpc) from NGC 5194(M51a), and show it to be physically distinct from the overlapping bright tidal streams from M51b. The Northwest Plume's morphology and red color ($B-V=0.8$) instead argue that it originated from tidal stripping of M51a's extreme outer disk. Finally, we confirm the strong segregation of gas and stars in the Southeast Tail, and do not detect any diffuse stellar component in the HI portion of the tail. Extant simulations of M51 have difficulty matching both the wealth of tidal structure in the system and the lack of stars in the HI tail, motivating new modeling campaigns to study the dynamical evolution of this classic interacting system.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

We describe the design, manufacture, and performance of bare-fiber integral field units (IFUs) for the SDSS-IV survey MaNGA (Mapping Nearby Galaxies at APO) on the the Sloan 2.5 m telescope at Apache Point Observatory (APO). MaNGA is a luminosity-selected integral-field spectroscopic survey of 10,000 local galaxies covering 360-1030 nm at R ~ 2200. The IFUs have hexagonal dense packing of fibers with packing regularity of 3 um (RMS), and throughput of 96+/-0.5% from 350 nm to 1 um in the lab. Their sizes range from 19 to 127 fibers (3-7 hexagonal layers) using Polymicro FBP 120:132:150 um core:clad:buffer fibers to reach a fill fraction of 56%. High throughput (and low focal-ratio degradation) is achieved by maintaining the fiber cladding and buffer intact, ensuring excellent surface polish, and applying a multi-layer AR coating of the input and output surfaces. In operations on-sky, the IFUs show only an additional 2.3% FRD-related variability in throughput despite repeated mechanical stressing during plate plugging (however other losses are present). The IFUs achieve on-sky throughput 5% above the single-fiber feeds used in SDSS-III/BOSS, attributable to equivalent performance compared to single fibers and additional gains from the AR coating. The manufacturing process is geared toward mass-production of high-multiplex systems. The low-stress process involves a precision ferrule with hexagonal inner shape designed to lead inserted fibers to settle in a dense hexagonal pattern. The ferrule inner diameter is tapered at progressively shallower angles toward its tip and the final 2 mm are straight and only a few um larger than necessary to hold the desired number of fibers. This process scales to accommodate other fiber sizes and to IFUs with substantially larger fiber count. (Abridged)

We present the first APOKASC catalog of spectroscopic and asteroseismic properties of 1916 red giants observed in the Kepler fields. The spectroscopic parameters provided from the Apache Point Observatory Galactic Evolution Experiment project are complemented with asteroseismic surface gravities, masses, radii, and mean densities determined by members of the Kepler Asteroseismology Science Consortium. We assess both random and systematic sources of error and include a discussion of sample selection for giants in the Kepler fields. Total uncertainties in the main catalog properties are of order 80 K in Teff , 0.06 dex in [M/H], 0.014 dex in log g, and 12% and 5% in mass and radius, respectively; these reflect a combination of systematic and random errors. Asteroseismic surface gravities are substantially more precise and accurate than spectroscopic ones, and we find good agreement between their mean values and the calibrated spectroscopic surface gravities. There are, however, systematic underlying trends with Teff and log g. Our effective temperature scale is between 0-200 K cooler than that expected from the Infrared Flux Method, depending on the adopted extinction map, which provides evidence for a lower value on average than that inferred for the Kepler Input Catalog (KIC). We find a reasonable correspondence between the photometric KIC and spectroscopic APOKASC metallicity scales, with increased dispersion in KIC metallicities as the absolute metal abundance decreases, and offsets in Teff and log g consistent with those derived in the literature. We present mean fitting relations between APOKASC and KIC observables and discuss future prospects, strengths, and limitations of the catalog data.

Using G dwarfs from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey, we have determined a vertical metallicity gradient over a large volume of the Milky Way's disk, and examined how this gradient varies for different [a/Fe] subsamples. This sample contains over 40,000 stars with low-resolution spectroscopy over 144 lines of sight. We employ the SEGUE Stellar Parameter Pipeline (SSPP) to obtain estimates of effective temperature, surface gravity, [Fe/H], and [a/Fe] for each star and extract multiple volume-complete subsamples of approximately 1000 stars each. Based on the survey's consistent target-selection algorithm, we adjust each subsample to determine an unbiased picture of the disk in [Fe/H] and [a/Fe]; consequently, each individual star represents the properties of many. The SEGUE sample allows us to constrain the vertical metallicity gradient for a large number of stars over a significant volume of the disk, between ~0.3 and 1.6 kpc from the Galactic plane, and examine the in situ structure, in contrast to previous analyses which are more limited in scope. This work does not pre-suppose a disk structure, whether composed of a single complex population or a distinct thin and thick disk component. The metallicity gradient is -0.243 +0.039 -0.053 dex/kpc for the sample as a whole, which we compare to various literature results. Each [a/Fe] subsample dominates at a different range of heights above the plane of the Galaxy, which is exhibited in the gradient found in the sample as a whole. Stars over a limited range in [a/Fe] show little change in median [Fe/H] with height. If we associate [a/Fe] with age, our consistent vertical metallicity gradients with [a/Fe] suggest that stars formed in different epochs exhibit comparable vertical structure, implying similar star-formation processes and evolution.

Fundamental stellar properties, such as mass, radius, and age, can be inferred using asteroseismology. Cool stars with convective envelopes have turbulent motions that can stochastically drive and damp pulsations. The properties of the oscillation frequency power spectrum can be tied to mass and radius through solar-scaled asteroseismic relations. Stellar properties derived using these scaling relations need verification over a range of metallicities. Because the age and mass of halo stars are well-constrained by astrophysical priors, they provide an independent, empirical check on asteroseismic mass estimates in the low-metallicity regime. We identify nine metal-poor red giants (including six stars that are kinematically associated with the halo) from a sample observed by both the Kepler space telescope and the Sloan Digital Sky Survey-III APOGEE spectroscopic survey. We compare masses inferred using asteroseismology to those expected for halo and thick-disk stars. Although our sample is small, standard scaling relations, combined with asteroseismic parameters from the APOKASC Catalog, produce masses that are systematically higher (<{\Delta}M>=0.17+/-0.05 Msun) than astrophysical expectations. The magnitude of the mass discrepancy is reduced by known theoretical corrections to the measured large frequency separation scaling relationship. Using alternative methods for measuring asteroseismic parameters induces systematic shifts at the 0.04 Msun level. We also compare published asteroseismic analyses with scaling relationship masses to examine the impact of using the frequency of maximum power as a constraint. Upcoming APOKASC observations will provide a larger sample of ~100 metal-poor stars, important for detailed asteroseismic characterization of Galactic stellar populations.

We present an online catalog of distance determinations for $\rm 6036$ K giants, most of which are members of the Milky Way's stellar halo. Their medium-resolution spectra from SDSS/SEGUE are used to derive metallicities and rough gravity estimates, along with radial velocities. Distance moduli are derived from a comparison of each star's apparent magnitude with the absolute magnitude of empirically calibrated color-luminosity fiducials, at the observed $(g-r)_0$ color and spectroscopic [Fe/H]. We employ a probabilistic approach that makes it straightforward to properly propagate the errors in metallicities, magnitudes, and colors into distance uncertainties. We also fold in ${\it prior}$ information about the giant-branch luminosity function and the different metallicity distributions of the SEGUE K-giant targeting sub-categories. We show that the metallicity prior plays a small role in the distance estimates, but that neglecting the luminosity prior could lead to a systematic distance modulus bias of up to 0.25 mag, compared to the case of using the luminosity prior. We find a median distance precision of $16\%$, with distance estimates most precise for the least metal-poor stars near the tip of the red-giant branch. The precision and accuracy of our distance estimates are validated with observations of globular and open clusters. The stars in our catalog are up to 125 kpc distant from the Galactic center, with 283 stars beyond 50 kpc, forming the largest available spectroscopic sample of distant tracers in the Galactic halo.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We present radio, optical/NIR, and X-ray observations of the afterglow of the short-duration 130603B, and uncover a break in the radio and optical bands at 0.5 d after the burst, best explained as a jet break with an inferred jet opening angle of 4-8 deg. GRB 130603B is only the third short GRB with a radio afterglow detection to date, and the first time that a jet break is evident in the radio band. We model the temporal evolution of the spectral energy distribution to determine the burst explosion properties and find an isotropic-equivalent kinetic energy of (0.6-1.7) x 10^51 erg and a circumburst density of 5 x 10^-3-30 cm^-3. From the inferred opening angle of GRB 130603B, we calculate beaming-corrected energies of Egamma (0.5-2) x 10^49 erg and EK (0.1-1.6) x 10^49 erg. Along with previous measurements and lower limits we find a median short GRB opening angle of 10 deg. Using the all-sky observed rate of 10 Gpc^-3 yr^-1, this implies a true short GRB rate of 20 yr^-1 within 200 Mpc, the Advanced LIGO/VIRGO sensitivity range for neutron star binary mergers. Finally, we uncover evidence for significant excess emission in the X-ray afterglow of GRB 130603B at >1 d and conclude that the additional energy component could be due to fall-back accretion or spin-down energy from a magnetar formed following the merger.

In the past decade, the notion that globular clusters (GCs) are composed of coeval stars with homogeneous initial chemical compositions has been challenged by growing evidence that they host an intricate stellar population mix, likely indicative of a complex history of star formation and chemical enrichment. Several models have been proposed to explain the existence of multiple stellar populations in GCs, but no single model provides a fully satisfactory match to existing data. Correlations between chemistry and global parameters such as cluster mass or luminosity are fundamental clues to the physics of GC formation. In this Letter, we present an analysis of the mean abundances of Fe, Mg, C, N, and Ca for 72 old GCs from the Andromeda galaxy. We show for the first time that there is a correlation between the masses of GCs and the mean stellar abundances of nitrogen, spanning almost two decades in mass. This result sheds new light on the formation of GCs, providing important constraints on their internal chemical evolution and mass loss history.

We present the design and performance of the multi-object fiber spectrographs for the Sloan Digital Sky Survey (SDSS) and their upgrade for the Baryon Oscillation Spectroscopic Survey (BOSS). Originally commissioned in Fall 1999 on the 2.5-m aperture Sloan Telescope at Apache Point Observatory, the spectrographs produced more than 1.5 million spectra for the SDSS and SDSS-II surveys, enabling a wide variety of Galactic and extra-galactic science including the first observation of baryon acoustic oscillations in 2005. The spectrographs were upgraded in 2009 and are currently in use for BOSS, the flagship survey of the third-generation SDSS-III project. BOSS will measure redshifts of 1.35 million massive galaxies to redshift 0.7 and Lyman-alpha absorption of 160,000 high redshift quasars over 10,000 square degrees of sky, making percent level measurements of the absolute cosmic distance scale of the Universe and placing tight constraints on the equation of state of dark energy. The twin multi-object fiber spectrographs utilize a simple optical layout with reflective collimators, gratings, all-refractive cameras, and state-of-the-art CCD detectors to produce hundreds of spectra simultaneously in two channels over a bandpass covering the near ultraviolet to the near infrared, with a resolving power R = \lambda/FWHM ~ 2000. Building on proven heritage, the spectrographs were upgraded for BOSS with volume-phase holographic gratings and modern CCD detectors, improving the peak throughput by nearly a factor of two, extending the bandpass to cover 360 < \lambda < 1000 nm, and increasing the number of fibers from 640 to 1000 per exposure. In this paper we describe the original SDSS spectrograph design and the upgrades implemented for BOSS, and document the predicted and measured performances.

We use deep surface photometry of the giant elliptical M49 (NGC 4472), obtained as part of our survey for diffuse light in the Virgo Cluster, to study the stellar populations in its outer halo. Our data trace M49's stellar halo out to ~ 100 kpc (7 Re), where we find that the shallow color gradient seen in the inner regions becomes dramatically steeper. The outer regions of the galaxy are quite blue (B-V ~ 0.7); if this is purely a metallicity effect, it argues for extremely metal poor stellar populations with [Fe/H] < -1. We also find that the extended accretion shells around M49 are distinctly redder than the galaxy's surrounding halo, suggesting that we are likely witnessing the buildup of both the stellar mass and metallicity in M49's outer halo due to late time accretion. While such growth of galaxy halos is predicted by models of hierarchical accretion, this growth is thought to be driven by more massive accretion events which have correspondingly higher mean metallicity than inferred for M49's halo. Thus the extremely metal-poor nature of M49's extended halo provides some tension against current models for elliptical galaxy formation.

We have used deep, wide-field optical imaging to study the faint outskirts of the luminous spiral galaxy M101 (NGC 5457), as well as its surrounding environment. Over six square degrees, our imaging has a limiting surface brightness of mu_B ~ 29.5 mag/arcsec^2, and has revealed the stellar structure of M101's disk out to nearly 25 arcminutes (50 kpc), three times our measured R25 isophotal size of the optical disk. At these radii, the well-known asymmetry of the inner disk slews 180 degrees, resulting in an asymmetric plume of light at large radius which follows the very extended HI disk to the northeast of M101. This plume has very blue colors (B-V ~ 0.2), suggesting it is the somewhat more evolved (few hundred Myr to ~ 1 Gyr) counterpart of the young far ultraviolet emitting population traced by GALEX imaging. We also detect another, redder spur of extended light to the east of the disk, and both structures are reminiscent of features produced during fly-by galaxy interactions. However, we see no evidence of very extended tidal tails around M101 or any of its companions which might be expected from a recent encounter with a massive companion. We consider the properties of M101's outer disk in light of possible past interactions with the nearby companion galaxies NGC 5477 and NGC 5474. The detection of optical starlight at such large radii gives us the ability to study star formation histories and stellar populations in outer disks over a longer timescales than those traced by the UV or Halpha emitting populations. Our data suggest ongoing buildup of the M101's outer disk due to encounters in the group environment triggering extended star formation and tidal heating of existing disk populations.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

We present the metallicity distribution function (MDF) for 24,270 G and 16,847 K dwarfs at distances from 0.2 to 2.3 kpc from the Galactic plane, based on spectroscopy from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. This stellar sample is significantly larger in both number and volume than previous spectroscopic analyses, which were limited to the solar vicinity, making it ideal for comparison with local volume-limited samples and Galactic models. For the first time, we have corrected the MDF for the various observational biases introduced by the SEGUE target selection strategy. The SEGUE sample is particularly notable for K dwarfs, which are too faint to examine spectroscopically far from the solar neighborhood. The MDF of both spectral types becomes more metal-poor with increasing |Z|, which reflects the transition from a sample with small [alpha/Fe] values at small heights to one with enhanced [alpha/Fe] above 1 kpc. Comparison of our SEGUE distributions to those of two different Milky Way models reveals that both are more metal-rich than our observed distributions at all heights above the plane. Our unbiased observations of G and K dwarfs provide valuable constraints over the |Z|-height range of the Milky Way disk for chemical and dynamical Galaxy evolution models, previously only calibrated to the solar neighborhood, with particular utility for thin- and thick-disk formation models.

The oxygen abundance gradients among nebular emission line regions in spiral galaxies have been used as important constraints for models of chemical evolution. We present the largest ever full-wavelength optical spectroscopic sample of emission line nebulae in a spiral galaxy (M31). We have collected spectra of 253 HII regions and 407 planetary nebulae with the Hectospec multi-fiber spectrograph of the MMT. We measure the line-of-sight extinction for 199 HII regions and 333 PNe; we derive oxygen abundance directly, based on the electron temperature, for 51 PNe; and we use strong line methods to estimate oxygen abundance for 192 HII regions and nitrogen abundance for 52 HII regions. The relatively shallow oxygen abundance gradient of the more extended HII regions in our sample is generally in agreement with the result of Zaritsky et al. (1994), based on only 19 M31 HII regions, but varies with the strong-line diagnostic employed. Our large sample size demonstrates that there is significant intrinsic scatter around this abundance gradient, as much as 3 times the systematic uncertainty in the strong line diagnostics. The intrinsic scatter is similar in the nitrogen abundances, although the gradient is significantly steeper. On small scales (deprojected distance < 0.5 kpc), HII regions exhibit local variations in oxygen abundance that are larger than 0.3 dex in 33% of neighboring pairs. We do not identify a significant oxygen abundance gradient among PNe, but we do find a significant gradient in the [N II] ratio that varies systematically with surface brightness. Our results underscore the complex and inhomogeneous nature of the ISM of M31, and our dataset illustrates systematic effects relevant to future studies of the metallicity gradients in nearby spiral galaxies.

We study two nearby, early-type galaxies, NGC4342 and NGC4291, that host unusually massive black holes relative to their low stellar mass. The observed black hole-to-bulge mass ratios of NGC4342 and NGC4291 are ~6.9% and ~1.9%, respectively, which significantly exceed the typical observed ratio of ~0.2%. As a consequence of the exceedingly large black hole-to-bulge mass ratios, NGC4342 and NGC4291 are ~5.1 sigma and ~3.4 sigma outliers from the M_BH - M_bulge scaling relation, respectively. In this paper, we explore the origin of the unusually high black hole-to-bulge mass ratio. Based on Chandra X-ray observations of the hot gas content of NGC4342 and NGC4291, we compute gravitating mass profiles, and conclude that both galaxies reside in massive dark matter halos, which extend well beyond the stellar light. The presence of dark matter halos around NGC4342 and NGC4291 and a deep optical image of the environment of NGC4342 indicate that tidal stripping, in which >90% of the stellar mass was lost, cannot explain the observed high black hole-to-bulge mass ratios. Therefore, we conclude that these galaxies formed with low stellar masses, implying that the bulge and black hole did not grow in tandem. We also find that the black hole mass correlates well with the properties of the dark matter halo, suggesting that dark matter halos may play a major role in regulating the growth of the supermassive black holes.

We examine the \alpha-element abundance ratio, [\alpha/Fe], of 5620 stars, observed by the Sloan Extension for Galactic Understanding and Exploration survey in the region 6 kpc < R < 16 kpc, 0.15 kpc < |Z| < 1.5 kpc, as a function of Galactocentric radius R and distance from the Galactic plane |Z|. Our results show that the high-\alpha\ thick disk population has a short scale length (L_thick ~ 1.8 kpc) compared to the low-\alpha population, which is typically associated with the thin disk. We find that the fraction of high-\alpha\ stars in the inner disk increases at large |Z|, and that high-\alpha\ stars lag in rotation compared to low-\alpha\ stars. In contrast, the fraction of high-\alpha\ stars in the outer disk is low at all |Z|, and high- and low-\alpha\ stars have similar rotational velocities up to 1.5 kpc from the plane. We interpret these results to indicate that different processes were responsible for the high-\alpha\ populations in the inner and outer disk. The high-\alpha\ population in the inner disk has a short scale length and large scale height, consistent with a scenario in which the thick disk forms during an early gas-rich accretion phase. Stars far from the plane in the outer disk may have reached their current locations through heating by minor mergers. The lack of high-\alpha\ stars at large R and |Z| also places strict constraints on the strength of radial migration via transient spiral structure.

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5\sigma in both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7\sigma. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

We present absorption line indices measured in the integrated spectra of globular clusters both from the Galaxy and from M 31. Our samples include 41 Galactic globular clusters, and more than 300 clusters in M 31. The conversion of instrumental equivalent widths into the Lick system is described, and zero-point uncertainties are provided. Comparison of line indices of old M 31 clusters and Galactic globular clusters suggests an absence of important differences in chemical composition between the two cluster systems. In particular, CN indices in the spectra of M 31 and Galactic clusters are essentially consistent with each other, in disagreement with several previous works. We reanalyze some of the previous data, and conclude that reported CN differences between M 31 and Galactic clusters were mostly due to data calibration uncertainties. Our data support the conclusion that the chemical compositions of Milky Way and M 31 globular clusters are not substantially different, and that there is no need to resort to enhanced nitrogen abundances to account for the optical spectra of M 31 globular clusters.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We employ measurements of the [alpha/Fe] ratio derived from low-resolution (R~2000) spectra of 17,277 G-type dwarfs from the SEGUE survey to separate them into likely thin- and thick-disk subsamples. Both subsamples exhibit strong gradients of orbital rotational velocity with metallicity, of opposite signs, -20 to -30 km/s/dex for the thin-disk and +40 to +50 km/s/dex for the thick-disk population. The rotational velocity is uncorrelated with Galactocentric distance for the thin-disk subsample, and exhibits a small trend for the thick-disk subsample. The rotational velocity decreases with distance from the plane for both disk components, with similar slopes (-9.0 {\pm} 1.0 km/s/kpc). Thick-disk stars exhibit a strong trend of orbital eccentricity with metallicity (about -0.2/dex), while the eccentricity does not change with metallicity for the thin-disk subsample. The eccentricity is almost independent of Galactocentric radius for the thin-disk population, while a marginal gradient of the eccentricity with radius exists for the thick-disk population. Both subsamples possess similar positive gradients of eccentricity with distance from the Galactic plane. The shapes of the eccentricity distributions for the thin- and thick-disk populations are independent of distance from the plane, and include no significant numbers of stars with eccentricity above 0.6. Among several contemporary models of disk evolution we consider, radial migration appears to have played an important role in the evolution of the thin-disk population, but possibly less so for the thick disk, relative to the gas-rich merger or disk heating scenarios. We emphasize that more physically realistic models and simulations need to be constructed in order to carry out the detailed quantitative comparisons that our new data enable.

We present the initial results of a 40 night contiguous ground-based campaign of time series photometric observations of a 1.39 sq. deg field located within the NASA Kepler mission field of view. The goal of this pre-launch survey was to search for transiting extrasolar planets and to provide independent variability information of stellar sources. We have gathered a data set containing light curves of 54,687 stars from which we have created a statistical sub-sample of 13,786 stars between 14< r <18.5 and have statistically examined each light curve to test for variability. We present a summary of our preliminary photometric findings including the overall level and content of stellar variability in this portion of the Kepler field and give some examples of unusual variable stars found within. We present a preliminary catalog of 2,457 candidate variable stars, of which 776 show signs of periodicity. We also present three potential exoplanet candidates, all of which should be observable in detail by the Kepler mission.

We continue our deep optical imaging survey of the Virgo cluster using the CWRU Burrell Schmidt telescope by presenting B-band surface photometry of the core of the Virgo cluster in order to study the cluster's intracluster light (ICL). We find ICL features down to mu_b ~ 29 mag sq. arcsec, confirming the results of Mihos et al. (2005), who saw a vast web of low-surface brightness streams, arcs, plumes, and diffuse light in the Virgo cluster core using V-band imaging. By combining these two data sets, we are able to measure the optical colors of many of the cluster's low-surface brightness features. While much of our imaging area is contaminated by galactic cirrus, the cluster core near the cD galaxy, M87, is unobscured. We trace the color profile of M87 out to over 2000 arcsec, and find a blueing trend with radius, continuing out to the largest radii. Moreover, we have measured the colors of several ICL features which extend beyond M87's outermost reaches and find that they have similar colors to the M87's halo itself, B-V ~ 0.8. The common colors of these features suggests that the extended outer envelopes of cD galaxies, such as M87, may be formed from similar streams, created by tidal interactions within the cluster, that have since dissolved into a smooth background in the cluster potential.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

We present new high signal-to-noise spectroscopic data on the M31 globular cluster system, obtained with the Hectospec multifiber spectrograph on the 6.5m MMT. More than 300 clusters have been observed at a resolution of 5A. The primary focus of this paper is the determination of mean cluster metallicities, ages and reddenings. Metallicities were estimated using a calibration of Lick indices with [Fe/H] provided by Galactic GCs. The metallicity distribution is not generally bimodal, in strong distinction with the bimodal Galactic globular distribution. Rather, the M31 distribution shows a broad peak, centered at [Fe/H]=-1, suggesting that the cluster systems of M31 and the Milky Way had different formation histories. Ages for clusters with [Fe/H] > -1 were determined using the automatic stellar population analysis program EZ_Ages. We find no evidence for massive clusters in M31 with intermediate ages. Moreover, we find that the mean ages of the old GCs are remarkably constant over about a decade in metallicity (-0.95 < [Fe/H] < 0.0). (abridged)

We analyze our accurate kinematical data for the old clusters in the inner regions of M31. These velocities are based on high S/N Hectospec data (Caldwell et al 2010). The data are well suited for analysis of M31's inner regions because we took particular care to correct for contamination by unresolved field stars from the disk and bulge in the fibers. The metal poor clusters show kinematics which are compatible with a pressure-supported spheroid. The kinematics of metal-rich clusters, however, argue for a disk population. In particular the innermost region (inside 2 kpc) shows the kinematics of the x2 family of bar periodic orbits, arguing for the existence of an inner Lindblad resonance in M31.

We present new limits on the ejection of metal-rich old-population hypervelocity stars from the Galactic center (GC) as probed by the SEGUE-2 survey. Our limits are a factor of 3-10 more stringent than previously reported, depending on stellar type. Compared to the known population of B-star ejectees, there can be no more than 30 times more metal-rich old-population F/G stars ejected from the GC. Because B stars comprise a tiny fraction of a normal stellar population, this places significant limits on a combination of the GC mass function and the ejection mechanism for hypervelocity stars. In the presence of a normal GC mass function, our results require an ejection mechanism that is about 5.5 times more efficient at ejecting B-stars compared to low-mass F/G stars.

We use deep V-band surface photometry of five of the brightest elliptical galaxies in the Virgo cluster to search for diffuse tidal streams, shells, and plumes in their outer halos (r > 50 kpc). We fit and subtract elliptical isophotal models from the galaxy images to reveal a variety of substructure, with surface brightnesses in the range mu_V= 26-29 mag/arcsec^2. M49 possesses an extended, interleaved shell system reminiscent of the radial accretion of a satellite companion, while M89's complex system of shells and plumes suggests a more complicated accretion history involving either multiple events or a major merger. M87 has a set of long streamers as might be expected from stripping of low luminosity dwarfs on radial orbits in Virgo. M86 also displays a number of small streams indicative of stripping of dwarf companions, but these comprise much less luminosity than those of M87. Only M84 lacks significant tidal features. We quantify the photometric properties of these structures, and discuss their origins in the context of each galaxy's environment and kinematics within the Virgo cluster.

We report on a single star, B030D, observed as part of a large survey of objects in M31, which has the unusual radial velocity of -780 km/s. Based on details of its spectrum, we find that the star is an F supergiant, with a circumstellar shell. The evolutionary status of the star could be one of a post-mainsequence close binary, a symbiotic nova, or less likely, a post-AGB star, which additional observations could help sort out. Membership of the star in the Andromeda Giant Stream can explain its highly negative velocity.

We present a means of characterizing and removing internal reflections between the CCD and other optical surfaces in an astronomical camera. The stellar reflections appear as out-of-focus images and are not necessarily axisymmetric about the star. Using long exposures of very bright stars as calibration images we are able to measure the position, size, and intensity of reflections as a function of their position on the field. We also measure the extended stellar point-spread function out to one degree. Together this information can be used to create an empirical model of the excess light from bright stars and reduce systematic artifacts in deep surface photometry. We then reduce a set of deep observations of the Virgo cluster with our method to demonstrate its efficacy and to provide a comparison with other strategies for removing scattered light.

We compare the distribution of diffuse intracluster light detected in the Virgo Cluster via broadband imaging with that inferred from searches for intracluster planetary nebulae (IPNe). We find a rough correspondence on large scales (~ 100 kpc) between the two, but with very large scatter (~ 1.3 mag/arcsec^2). On smaller scales (1 -- 10 kpc), the presence or absence of correlation is clearly dependent on the underlying surface brightness. On these scales, we find a correlation in regions of higher surface brightness (mu_V < ~27) which are dominated by the halos of large galaxies such as M87, M86, and M84. In those cases, we are likely tracing PNe associated with galaxies rather than true IPNe. In true intracluster fields, at lower surface brightness, the correlation between luminosity and IPN candidates is much weaker. While a correlation between broadband light and IPNe is expected based on stellar populations, a variety of statistical, physical, and methodological effects can act to wash out this correlation and explain the lack of a strong correlation at lower surface brightness found here. [abridged]

We have measured the amount of kinematic substructure in the Galactic halo using the final data set from the Spaghetti project, a pencil-beam high latitude sky survey. Our sample contains 101 photometrically selected and spectroscopically confirmed giants with accurate distance, radial velocity and metallicity information. We have developed a new clustering estimator: the "4distance" measure, which when applied to our data set leads to the identification of 1 group and 7 pairs of clumped stars. The group, with 6 members, can confidently be matched to tidal debris of the Sagittarius dwarf galaxy. Two pairs match the properties of known Virgo structures. Using models of the disruption of Sagittarius in Galactic potentials with different degrees of dark halo flattening, we show that this favors a spherical or prolate halo shape, as demonstrated by Newberg et al. (2007) using SDSS data. One additional pair can be linked to older Sagittarius debris. We find that 20% of the stars in the Spaghetti data set are in substructures. From comparison with random data sets we derive a very conservative lower limit of 10% to the amount of substructure in the halo. However, comparison to numerical simulations shows that our results are also consistent with a halo entirely built up from disrupted satellites, provided the dominating features are relatively broad due to early merging or relatively heavy progenitor satellites.

We show that the Sagittarius dwarf tidal stream can be traced with very red K/M-giant stars selected from SDSS photometry. A subset of these stars are spectroscopically confirmed with SEGUE and SDSS spectra, and the distance scale of 2MASS and SDSS M giants is calibrated to the RR Lyrae distance scale. The absolute g band magnitude of the K/M-giant stars at the tip of the giant branch is M_g=-1.0. The line-of-sight velocities of the M giant and BHB stars that are spatially coincident with the Sgr dwarf tidal stream are consistent with those of previous authors, reinforcing the need for new models that can explain all of the Sgr tidal debris stream observations. We estimate stellar densities along the tidal tails that can be used to help constrain future models. The K/M-giant, BHB, and F-turnoff stars in the lower surface brightness tidal stream that is adjacent to the main leading Sgr dwarf tidal tail have velocities and metallicities that are similar to those of the stars in the leading tidal tail. The ratio of K/M giants to BHBs and BHBs to F-turnoff stars are also similar for both branches of the leading tidal tail. We show that there is an additional low-metallicity tidal stream near the Sgr trailing tidal tail.

Using a sample of 248 metal-poor stars (RR Lyraes, red giants and RHB stars) which is remarkable for the accuracy of its 6-D kinematical data, we find a new component for the local halo which has an axial ratio c/a ~ 0.2, a similar flattening to the thick disk. It has a small prograde rotation but is supported by velocity anisotropy, and contains more intermediate-metallicity stars (with -1.5 < [Fe/H] < -1.0) than the rest of our sample. We suggest that this component was formed quite late, during or after the formation of the disk. It formed either from the gas that was accreted by the last major mergers experienced by the Galaxy, or by dynamical friction of massive infalling satellite(s) with the halo and possibly the stellar disk or thick disk. The remainder of the stars in our sample exhibit a clumpy distribution in energy and angular momentum, suggesting that the early, chaotic conditions under which the inner halo formed were not violent enough to erase the record of their origins. The clumpy structure suggests that a relatively small number of progenitors were responsible for building up the inner halo, in line with theoretical expectations. We find a difference in mean binding energy between the RR Lyrae variables and the red giants in our sample, suggesting that more of the RR Lyraes in the sample belong to the outer halo, and that the outer halo may be somewhat younger, as first suggested by Searle and Zinn (1978). We also find that the RR Lyrae mean rotation is more negative than the red giants, which is consistent with the recent result of Carollo et al.(2007) that the outer halo has a retrograde rotation and with the difference in kinematics seen between RR Lyraes and BHB stars by Kinman et al.(2007).

We present an updated catalog of 1300 objects in the field of M31, including 670 likely star clusters of various types. Archival images were inspected to confirm cluster classifications where possible, but most of the classifications were based on spectra taken of ~1000 objects with the Hectospec fiber positioner and spectrograph on the 6.5m MMT. The spectra and images of young clusters are analyzed in detail in this paper. Ages, reddenings and masses of 140 young clusters are derived by comparing the observed spectra and photometry with model spectra. We find these to have masses as great as 10^5 with a median of 10^4 M_sun, and a median age of 0.25 Gyr. Thus these clusters are similar in mass to the LMC young massive clusters, and are in between Milky Way open clusters and the globulars of M31 and the Milky Way. Most of the M31 young clusters have the low concentration typical of Milky Way open clusters, and we expect that most of these will be disrupted in the next Gyr, but a few have high concentrations and will likely survive longer. The spatial distribution of the young clusters is well correlated with the star-forming regions as mapped out by mid-IR emission. A kinematic analysis likewise confirms the spatial association of the young clusters with the young disk in M31.

We present photometry for globular and open cluster stars observed with the Sloan Digital Sky Survey (SDSS). In order to exploit over 100 million stellar objects with r < 22.5 mag observed by SDSS, we need to understand the characteristics of stars in the SDSS ugriz filters. While star clusters provide important calibration samples for stellar colors, the regions close to globular clusters, where the fraction of field stars is smallest, are too crowded for the standard SDSS photometric pipeline to process. To complement the SDSS imaging survey, we reduce the SDSS imaging data for crowded cluster fields using the DAOPHOT/ALLFRAME suite of programs and present photometry for 17 globular clusters and 3 open clusters in a SDSS value-added catalog. Our photometry and cluster fiducial sequences are on the native SDSS 2.5-meter ugriz photometric system, and the fiducial sequences can be directly applied to the SDSS photometry without relying upon any transformations. Model photometry for red giant branch and main-sequence stars obtained by Girardi et al. cannot be matched simultaneously to fiducial sequences; their colors differ by ~0.02-0.05 mag. Good agreement (< ~0.02 mag in colors) is found with Clem et al. empirical fiducial sequences in u'g'r'i'z' when using the transformation equations in Tucker et al.

We use N-body simulations to explore the origin and a plausible orbit for the Orphan Stream, one of the faintest substructures discovered so far in the outer halo of our Galaxy. We are able to reproduce its position, velocity and distance measurements by appealing to a single wrap of a double-component satellite galaxy. We find that the progenitor of the Orphan Stream could have been an object similar to today's Milky Way dwarfs, such as Carina, Draco, Leo II or Sculptor; and unlikely to be connected to Complex A or Ursa Major II. Our models suggest that such progenitors, if accreted on orbits with apocenters smaller than ~35 kpc, are likely to give rise to very low surface brightness streams, which may be hiding in the outer halo and remain largely undetected with current techniques. The systematic discovery of these ghostly substructures may well require wide field spectroscopic surveys of the Milky Way's outer stellar halo.

In our study of M31's globular cluster system with MMT/Hectospec, we have obtained high-quality spectra of 85 clusters with ages less than 1 Gyr. With the exception of Hubble V, the young cluster in NGC 205, we find that these young clusters have kinematics and spatial distribution consistent with membership in M31's young disk. Preliminary estimates of the cluster masses and structural parameters, using spectroscopically derived ages and HST imaging, confirms earlier suggestions that M31 has clusters similar to the LMC's young populous clusters.

We have assembled a sample of halo stars in the solar neighborhood to look for halo substructure in velocity and angular momentum space. Our sample includes red giants, RR Lyrae, and red horizontal branch stars within 2.5 kpc of the Sun with [Fe/H] less than -1.0. It was chosen to include stars with accurate distances, space velocities, and metallicities as well as well-quantified errors. We confirm the existence of the streams found by Helmi and coworkers, which we refer to as the H99 streams. These streams have a double-peaked velocity distribution in the z direction. We use the results of modeling of the H99 streams by Helmi and collaborators to test how one might use v_z velocity information and radial velocity information to detect kinematic substructure in the halo. We find that detecting the H99 streams with radial velocities alone would require a large sample. We use the velocity distribution of the H99 streams to estimate their age. From our model of the progenitor of the H99 streams, we determine that it was accreted between 6 and 9 Gyr ago. The H99 streams have [alpha/Fe] abundances similar to other halo stars in the solar neighborhood, suggesting that the gas that formed these stars were enriched mostly by Type II SNe. We have also discovered in angular momentum space two other possible substructures, which we refer to as the retrograde and prograde outliers. The retrograde outliers are likely to be halo substructure, but the prograde outliers are most likely part of the smooth halo. The retrograde outliers have significant structure in the v_phi direction and show a range of [alpha/Fe]. The methods presented in this paper can be used to exploit the kinematic information present in future large databases like RAVE, SDSSII/SEGUE, and Gaia.

We quantify the variability of faint unresolved optical sources using a catalog based on multiple SDSS imaging observations. The catalog covers SDSS Stripe 82, and contains 58 million photometric observations in the SDSS ugriz system for 1.4 million unresolved sources. In each photometric bandpass we compute various low-order lightcurve statistics and use them to select and study variable sources. We find that 2% of unresolved optical sources brighter than g=20.5 appear variable at the 0.05 mag level (rms) simultaneously in the g and r bands. The majority (2/3) of these variable sources are low-redshift (<2) quasars, although they represent only 2% of all sources in the adopted flux-limited sample. We find that at least 90% of quasars are variable at the 0.03 mag level (rms) and confirm that variability is as good a method for finding low-redshift quasars as is the UV excess color selection (at high Galactic latitudes). We analyze the distribution of lightcurve skewness for quasars and find that is centered on zero. We find that about 1/4 of the variable stars are RR Lyrae stars, and that only 0.5% of stars from the main stellar locus are variable at the 0.05 mag level. The distribution of lightcurve skewness in the g-r vs. u-g color-color diagram on the main stellar locus is found to be bimodal (with one mode consistent with Algol-like behavior). Using over six hundred RR Lyrae stars, we demonstrate rich halo substructure out to distances of 100 kpc. We extrapolate these results to expected performance by the Large Synoptic Survey Telescope and estimate that it will obtain well-sampled 2% accurate, multi-color lightcurves for ~2 million low-redshift quasars, and will discover at least 50 million variable stars.

We present a cautionary study exploring the reliability of the H-delta line in the integrated spectra of galaxies for determining galaxy ages. Our database consists of the observed integrated spectra of ~120 early-type galaxies, of 7 metal-rich globular clusters in M31 and the Galactic globular cluster 47 Tuc, and of the open cluster M67. We have measured H-delta using index definitions designed to assess contamination from the CN molecule in and around H-delta by choosing combinations of bandpasses that both avoid and include a region of CN molecular lines redward of H-delta. We find systematic differences in the ages derived from H-delta measurements among the various definitions when extracting ages from H-delta in old stellar populations with enhanced CN bands due to non-solar abundance ratios. We propose that neighboring CN lines have a strong effect on pseudocontinuum and central bandpass levels. For stellar populations which have non-solar abundance ratios in C and/or N, population synthesis models that do not account for abundance ratio variations cannot reproduce accurately the CN 4216 \AA band, which leads to a corresponding inaccuracy in reproducing the various H-delta indices. Hence, caution must be used when extracting galaxy ages from the H-delta line in old stellar populations with significant non-solar abundance ratios.

We compare Stetson's photometric standards with measurements listed in a standard star catalog constructed using repeated SDSS imaging observations. The SDSS catalog includes over 700,000 candidate standard stars from the equatorial stripe 82 (|Dec|<1.266 deg) in the RA range 20h 34' to 4h 00', and with the $r$ band magnitudes in the range 14--21. The distributions of measurements for individual sources demonstrate that the SDSS photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in ugriz, respectively (about twice as good as for individual SDSS runs). We derive mean photometric transformations between the SDSS gri and the BVRI system using 1165 Stetson stars found in the equatorial stripe 82, and then study the spatial variation of the difference in zeropoints between the two catalogs. Using third order polynomials to describe the color terms, we find that photometric measurements for main-sequence stars can be transformed between the two systems with systematic errors smaller than a few millimagnitudes. The spatial variation of photometric zeropoints in the two catalogs typically does not exceed 0.01 magnitude. Consequently, the SDSS Standard Star Catalog for Stripe 82 can be used to calibrate new data in both the SDSS ugriz and the BVRI systems with a similar accuracy.

We present a detailed analysis of the star formation history (SFH) of three fields in M33 located ~ 4 - 6 visual scale lengths from its nucleus. These fields were imaged with the Advanced Camera for Surveys on the Hubble Space Telescope and reach ~ 2.5 magnitudes below the red clump of core helium burning stars. The observed color-magnitude diagrams are modeled as linear combinations of individual synthetic populations with different ages and metallicities. To gain a better understanding of the systematic errors we have conducted the analysis with two different sets of stellar evolutionary tracks which we designate as Padova (Girardi et al. 2000) and Teramo (Pietrinferni et al. 2004). The precise details of the results depend on which tracks are used but we can make several conclusions that are fairly robust despite the differences. Both sets of tracks predict the mean age to increase and the mean metallicity to decrease with radius. Allowing age and metallicity to be free parameters and assuming star formation began ~ 14 Gyr ago, we find that the mean age of all stars and stellar remnants increases from ~ 6 Gyr to ~ 8 Gyr and the mean global metallicity decreases from ~ -0.7 to ~ -0.9. The fraction of stars formed by 4.5 Gyr ago increases from ~ 65% to ~ 80%. The mean star formation rate 80 - 800 Myr ago decreases from ~ 30% of the lifetime average to just ~ 5%. The random errors on these estimates are ~ 10%, 1.0 Gyr, and 0.1 dex. By comparing the results of the two sets of stellar tracks for the real data and for test populations with known SFH we have estimated the systematic errors to be 15%, 1.0 Gyr, and 0.2 dex. These do not include uncertainties in the bolometric corrections or variations in alpha-element abundance which deserve future study.

Studying the stellar populations in the outskirts of spiral galaxies can provide important constraints on their structure, formation, and evolution. To that end, we present VI photometry obtained with the Advanced Camera for Surveys for three fields located ~ 20' - 30' in projected distance southeast of M33's nucleus (corresponding to ~ 4 - 6 visual scale lengths or ~ 9 - 13 kpc in deprojected radius). The color-magnitude diagrams reveal a mixed stellar population whose youngest constituents have ages no greater than ~ 100 Myr and whose oldest members have ages of at least several Gyr. The presence of stars as massive as 3 - 5 Msun is consistent with global star formation thresholds in disk galaxies but could argue for a threshold in M33 that is on the low end of observational and theoretical expectations. The metallicity gradient as inferred by comparing the observed red giant branch (RGB) to the Galactic globular clusters is consistent with M33's inner disk gradient traced by several other studies. The surface density of RGB stars drops off exponentially with a radial scale length of 4.7' +/- 0.1'. The scale length increases with age in a manner similar to the vertical scale height of several nearby late-type spirals. Based on the metallicity gradient, density gradient, and mixed nature of the stellar populations, we conclude these fields are dominated by a disk population although we cannot rule out the presence of a small halo component.

We present integrated photometry and color-magnitude diagrams (CMDs) for 24 star clusters in M33, of which 12 were previously uncataloged. These results are based on Advanced Camera for Surveys observations from the Hubble Space Telescope of two fields in M33. Our integrated V magnitudes and V-I colors for the previously identified objects are in good agreement with published photometry. We are able to estimate ages for 21 of these clusters using features in the CMDs, including isochrone fitting to the main sequence turnoffs for 17 of the clusters. Comparisons of these ages with the clusters' integrated colors and magnitudes suggest that simple stellar population models perform reasonably well in predicting these properties.

We present observations of RR Lyrae variables in the Local Group late-type spiral galaxy M33. Using the Advanced Camera for Surveys on the Hubble Space Telescope, we have identified 64 ab-type RR Lyraes in M33. We have estimated reddenings for these stars based on their minimum light V-I colors and metallicities based on their periods. From the distributions of these properties, we conclude that the RR Lyraes belong to two populations - one associated with the halo of M33 and the other with its disk. Given that RR Lyraes are produced by populations older than ~10 Gyr, this suggests that not only does the field halo of M33 contain an old component, but so does its disk. This is one of the best pieces of evidence for the existence of a halo field component in M33. Using a relation between RR Lyrae absolute magnitude and metallicity (Mv(RR) = 0.23[Fe/H] + 0.93), we estimate a mean distance modulus of <(m-M)o> = 24.67 +/- 0.08 for M33. This places M33 approximately 70 kpc beyond M31 in line-of-sight distance.

We present deep optical imaging of the inner 1.5 x 1.5 degrees of the Virgo cluster to search for diffuse intracluster light (ICL). Our image reaches a 1 sigma depth of mu_v=28.5 mag/arcsec^2 -- 1.5 mag/arcsec^2 deeper than previous surveys -- and reveals an intricate web of diffuse intracluster light. We see several long (>100 kpc) tidal streamers, as well as a myriad of smaller-scale tidal tails and bridges between galaxies. The diffuse halo of M87 is traced out to nearly 200 kpc, appearing very irregular on these scales, while significant diffuse light is also detected around the M84/M86 pair. Several galaxies in the core are embedded in common envelopes, suggesting they are true physical subgroups. The complex substructure of Virgo's diffuse ICL reflects the hierarchical nature of cluster assembly, rather than being the product of smooth accretion around a central galaxy.

We have compiled a new sample of 240 halo objects with accurate distance and radial velocity measurements, including globular clusters, satellite galaxies, field blue horizontal branch stars and red giant stars from the Spaghetti survey. The new data lead to a significant increase in the number of known objects for Galactocentric radii beyond 50 kpc, which allows a reliable determination of the radial velocity dispersion profile out to very large distances. The radial velocity dispersion shows an almost constant value of 120 km/s out to 30 kpc and then continuously declines down to 50 km/s at about 120 kpc. This fall-off puts important constraints on the density profile and total mass of the dark matter halo of the Milky Way. For a constant velocity anisotropy, the isothermal profile is ruled out, while both a dark halo following a truncated flat model of mass $1.2^{+1.8}_{-0.5}\times 10^{12}$ M_sun and an NFW profile of mass $0.8^{+1.2}_{-0.2}\times 10^{12}$ M_sun and c=18 are consistent with the data. The significant increase in the number of tracers combined with the large extent of the region probed by these has allowed a more precise determination of the Milky Way mass in comparison to previous works. We also show how different assumptions for the velocity anisotropy affect the performance of the mass models.

We present kinematics of 135 planetary nebulae in M31 from a survey covering 3.9 square degrees and extending out to 15 kpc from the southwest major axis and more than 20 kpc along the minor axis. The majority of our sample, even well outside the disk, shows significant rotational support (mean line-of-sight velocity 116 km/s). We argue that these PN belong to the outer part of M31's large de Vaucouleurs bulge. Only five PN have velocities clearly inconsistent with this fast rotating bulge. All five may belong to tidal streams in M31's outer halo. One is projected on the Northern Spur, and is counter-rotating with respect to the disk there. Two are projected along the major axis at X=-10 kpc and have M32-like velocities; they could be debris from that galaxy. The remaining two halo PN are located near the center of the galaxy and their velocities follow the gradient found by Ibata et al. (2004), implying that these PN could belong to the Southern Stream. If M31 has a non-rotating, pressure-supported halo, we have yet to find it, and it must be a very minor component of the galaxy.

We report the search for intracluster light in four Abell Type II/III (non-cD) galaxy clusters: Abell 801, 1234, 1553, & 1914. We find on average that these clusters contain $\sim$ 10% of their detected stellar luminosity in a diffuse component. We show that for two of the clusters the intracluster light closely follows the galaxy distribution, but in the other two cases, there are noticeable differences between the spatial distribution of the galaxies and the intracluster light. We report the results of a search for intracluster tidal debris in each cluster, and note that Abell 1914 in particular has a number of strong tidal features likely due to its status as a recent cluster merger. One of the Abell 1914 features appears to be spatially coincident with an extension seen in weak lensing maps, implying the feature traces a large amount of mass. We compare these results to numerical simulations of hierarchically-formed galaxy clusters, and find good general agreement between the observed and simulated images, although we also find that our observations sample only the brightest features of the intracluster light. Together, these results suggest that intracluster light can be a valuable tool in determining the evolutionary state of galaxy clusters.

We give an update of our ongoing survey for intracluster light (ICL), in a sample of distant Abell clusters. We find that the amount of intracluster starlight is comparable to that seen in nearby clusters, and that tidal debris appears to be common.

We report the detection of a new satellite of M31, projected close to M32. And VIII is tidally distorted, with length ~10 kpc and width a few kpc. It contains 5-12 planetary nebulae (PNe) and 1-3 globular clusters, and has a velocity of -204 km/s with respect to M31, some 350 km/s away from M32's velocity. There are also about 4 x 10^5 solar masses of HI, well-separated from the disk, at the same position and velocity. The satellite has luminosity of 1.2-2.4 x 10^8 solar luminosities, and a central surface brightness of order \mu_V=24. Both these values are typical of Local Group dwarf galaxies. Its surface brightness is some 6 magnitudes brighter than any of the stellar streams found in the Milky Way or M31. The three associated globular clusters have reddening consistent with foreground reddening from the Milky Way only, making it likely that the satellite is in front of M31, unlike the giant tidal stream of Ibata et al.(2001), which is behind M31 in the SE quadrant. However, the major axis of And VIII is aligned with the western edge of this giant stream, and we suggest that its unusual fan shape is caused by superposition of two streams, the westernmost of which was tidally stripped from And VIII.

We show that there is a subsystem of the M31 globular clusters with THIN DISK kinematics. These clusters span the entire metallicity range of the M31 globular cluster system, in contrast to the (thick) disk globulars in the Milky Way which are predominantly metal-rich. Disk globular clusters are found across the entire disk of M31 and form around 40% of the clusters projected on its disk. The existence of such a disk system suggests that there was a relatively large thin disk in place very early in M31's history. Accurate measures of the ages of these clusters will constrain the epoch of disk formation in M31. There is currently no strong evidence for differences in age between Milky Way and M31 globulars. While age differences are subtle for old populations, it is unlikely that disk clusters with [Fe/H] around -2.0 were formed after significant star formation began in the galaxy, as the proto-cluster gas would be enriched by supernova ejecta. Thus it is likely that M31 had a rather large disk in place at early epochs. The very existence of such a cold disk means that M31 has suffered no mergers with an object of 10% or more of the disk mass since the clusters were formed. This makes Brown et al (2003)'s suggestion that M31 could have suffered an equal-mass merger 6-8 Gyr ago less viable.

We present early observations of the afterglow of the Gamma-Ray Burst (GRB) 030329 and the spectroscopic discovery of its associated supernova SN 2003dh. We obtained spectra of the afterglow of GRB 030329 each night from March 30.12 (0.6 days after the burst) to April 8.13 (UT) (9.6 days after the burst). The spectra cover a wavelength range of 350 nm to 850 nm. The early spectra consist of a power-law continuum (F_nu ~ nu^{-0.9}) with narrow emission lines originating from HII regions in the host galaxy, indicating a low redshift of z=0.1687. However, our spectra taken after 2003 Apr. 5 show broad peaks in flux characteristic of a supernova. Correcting for the afterglow emission, we find the spectrum of the supernova is remarkably similar to the type Ic `hypernova' SN 1998bw. While the presence of supernovae have been inferred from the light curves and colors of GRB afterglows in the past, this is the first direct, spectroscopic confirmation that a subset of classical gamma-ray bursts originate from supernovae.

We present some early results from our deep imaging survey of galaxy clusters intended to detect and study intracluster light (ICL). From our observations to date, we find that ICL is common in galaxy clusters, and that substructure in the ICL also appears to be common as well. We also discuss some initial comparisons of our imaging results to high-resolution numerical simulations of galaxy clusters, and give avenues for future research.

We report the initial results of a deep imaging survey of galaxy clusters. The primary goals of this survey are to quantify the amount of intracluster light as a function of cluster properties, and to quantify the frequency of tidal debris. We outline the techniques needed to perform such a survey, and we report findings for the first two galaxy clusters in the survey: Abell 1413, and MKW 7 . These clusters vary greatly in richness and structure. We show that our surface photometry reliably reaches to a surface brightness of \mu_v = 26.5 mags per arcsec. We find that both clusters show clear excesses over a best-fitting r^{1/4} profile: this was expected for Abell 1413, but not for MKW 7. Both clusters also show evidence of tidal debris in the form of plumes and arc-like structures, but no long tidal arcs were detected. We also find that the central cD galaxy in Abell 1413 is flattened at large radii, with an ellipticity of $\approx 0.8$, the largest measured ellipticity of any cD galaxy to date.

The last 10-20 years has seen a profound shift in views of how the Galaxy's halo formed. The idea of a monolithic early collapse of a single system (Eggen, Lynden-Bell and Sandage 1962) has been challenged by observations at high redshift and by cosmological models of structure formation. These findings imply that we should see clear evidence of hierarchical formation processes in nearby galaxies. Recent studies of our Galaxy, made possible by large-scale CCD surveys such as the Sloan Digital Sky Survey (SDSS), have begun to reveal tantalizing evidence of substructure in the outer halo. We review evidence for tidal streams associated with known Milky Way satellites and for star streams whose progenitors are still unknown. This includes results from the SDSS and our own ongoing pencil-beam halo survey, the Spaghetti survey.

As part of the Spaghetti Project Survey (SPS) we have detected a concentration of giant stars well above expectations for a smooth halo model. The position (l~350, b~50) and distance (~50 kpc) of this concentration match those of the Northern over-density detected by SDSS (Yanny et al. 2000, Ivezic et al. 2000). We find additional evidence for structure at ~80 kpc in the same direction. We present radial velocities for many of these stars, including the first published results from the 6.5m Magellan telescope. The radial velocities for stars in these structures are in excellent agreement with models of the dynamical evolution of the Sgr dwarf tidal debris, whose center is 60 degrees away. The metallicity of stars in these streams is lower than that of the main body of the Sgr dwarf, which may indicate a radial metallicity gradient prior to disruption.

We have simulated the evolution of tidal debris in the Galactic halo in order to guide our ongoing survey to determine the fraction of halo mass accreted via satellite infall. Contrary to naive expectations that the satellite debris will produce a single narrow velocity peak on a smooth distribution, there are many different signatures of substructure, including multiple peaks and broad but asymmetrical velocity distributions. Observations of the simulations show that there is a high probability of detecting the presence of tidal debris with a pencil beam survey of 100 square degrees. In the limiting case of a single 10^7 Msun satellite contributing 1% of the luminous halo mass the detection probability is a few percent using just the velocities of 100 halo stars in a single 1 square degree field. The detection probabilities scale with the accreted fraction of the halo and the number of fields surveyed. There is also surprisingly little dependence of the detection probabilities on the time since the satellite became tidally disrupted, or on the initial orbit of the satellite, except for the time spent in the survey volume.

We critically examine the use of the Washington photometric system (with the 51 filter) for identifying distant halo giants. While this is the most powerful photometric technique for isolating G and K giant stars, spectroscopic follow-up of giant candidates is vital. There are two situations in which interlopers outnumber genuine giants in the diagnostic M-51/M-T2 plot, and are indistinguishable photometrically from the giants. (1) In deep surveys covering tens of square degrees, very metal-poor halo dwarfs are a significant contaminant. An example is our survey of the outer halo (Morrison et al. 2000, Dohm-Palmer et al. 2000), where these metal-poor dwarfs dominate the number of photometric giant candidates at magnitudes fainter than V = 18 and cannot be isolated photometrically. (2) In deep surveys of smaller areas with low photometric precision, most objects in the giant region of the color-color plot are dwarfs whose photometric errors have moved them there. Color errors in M-51 and M-T2 need to be smaller than 0.03 mag to avoid this problem. An example of a survey whose photometric errors place the giant identifications under question is the survey for extra-tidal giants around the Carina dwarf spheroidal of Majewski et al. (2000a). Accurate photometry and spectroscopic follow-up of giant candidates are essential when using the Washington system to identify the rare outer-halo giants.

We present imaging results from a high Galactic latitude survey designed to examine the structure of the Galactic halo. The objective of the survey is to identify candidate halo stars which can be observed spectroscopically to obtain radial velocities and confirm halo membership. The Washington filter system is used for its ability to distinguish between dwarfs and giants, as well as provide a metallicity indicator. Our most successful imaging run used the BTC camera on the CTIO 4m telescope in April 1999. Photometric conditions during these observations provided superb photometry, with average errors for a star at $M=18.5$ of 0.009, 0.008, 0.011, and 0.009 for $C$, $M$, $DDO51$, and $T2$ respectively. These data are available with the electronic version of this paper, as well as through ADC (http://adc.gsfc.nasa.gov/). We use these data as a template to describe the details of our photometric reduction process. It is designed to perform CCD reductions and stellar photometry automatically during the observation run without the aid of external packages, such as IRAF and IDL. We describe necessary deviations from this procedure for other instruments used in the survey up to June 2000. Preliminary results from spectroscopic observations indicate a 97% efficiency in eliminating normal dwarfs from halo giant candidates for $M<18.5$. Unfortunately, low-metallicity subdwarfs cannot be photometrically distinguished from giants using the Washington filters. These major contaminates unavoidably reduced the overall giant identification efficiency to 66% for $M<18.5$. Our improved knowledge of these stars will increase this efficiency for future spectroscopic observations.

We describe a major survey of the Milky Way halo designed to test for kinematic substructure caused by destruction of accreted satellites. We use the Washington photometric system to identify halo stars efficiently for spectroscopic followup. Tracers include halo giants (detectable out to more than 100 kpc), blue horizontal branch stars, halo stars near the main sequence turnoff, and the ``blue metal-poor stars'' of Preston et al (1994). We demonstrate the success of our survey by showing spectra of stars we have identified in all these categories, including giants as distant as 75 kpc. We discuss the problem of identifying the most distant halo giants. In particular, extremely metal-poor halo K dwarfs are present in approximately equal numbers to the distant giants for V fainter than 18, and we show that our method will distinguish reliably between these two groups of metal-poor stars. We plan to survey 100 square degrees at high galactic latitude, and expect to increase the numbers of known halo giants, BHB stars and turnoff stars by more than an order of magnitude. In addition to the strong test that this large sample will provide for the question `was the Milky Way halo accreted from satellite galaxies?', we will improve the accuracy of mass measurements of the Milky Way beyond 50 kpc via the kinematics of the many distant giants and BHB stars we will find. We show that one of our first datasets constrains the halo density law over galactocentric radii of 5-20 kpc and z heights of 2-15 kpc. The data support a flattened power-law halo with b/a of 0.6 and exponent -3.0. More complex models with a varying axial ratio may be needed with a larger dataset.

We have identified twenty-seven candidate very low mass members of the relatively young Alpha Persei open cluster from a six square degree CCD imaging survey. Based on their I magnitudes and the nominal age and distance to the cluster, these objects should have masses less than 0.1 Msunif they are cluster members. We have subsequently obtained intermediate resolution spectra of seventeen of these objects using the Keck II telescope and LRIS spectrograph. We have also obtained near-IR photometry for many of the stars. Our primary goal was to determine the location of the "lithium depletion boundary" and hence to derive a precise age for the cluster. We detect lithium with equivalent widths greater than or equal to 0.4 \AA in five of the program objects. We have constructed a color-magnitude diagram for the faint end of the Alpha Persei main sequence. These data allow us to accurately determine the Alpha Persei single-star lithium depletion boundary at M(I$_C$) = 11.47, M(Bol) = 11.42, (R-I)$_{C0}$ = 2.12, spectral type M6.0. By reference to theoretical evolutionary models, this converts fairly directly into an age for the Alpha Persei cluster of 90 $\pm$ 10 Myr. At this age, the two faintest of our spectroscopically confirmed members should be sub-stellar (i.e., brown dwarfs) according to theoretical models.

The distance to NGC 3319 has been determined from Cepheid variable stars as part of the Hubble Space Telescope Key Project on the Extragalactic Distance Scale. Thirteen and four epochs of observations, using filters F555W (V) and F814W (I) respectively, were made with the Wide Field Planetary Camera 2. Thirty-three Cepheid variables between periods of 8 and 47 days were discovered. Adopting a Large Magellanic Cloud distance modulus of 18.50 +- 0.10 mag and extinction of E(V-I)=0.13 mag, a true reddening-corrected distance modulus (based on an analysis employing the ALLFRAME software package) of 30.78 +- 0.14 (random) +- 0.10 (systematic) mag and the extinction of E(V-I) = 0.06 mag were determined for NGC 3319. This galaxy is the last galaxy observed for the HST H0 Key Project.

We have obtained deep surface photometry of the edge-on spiral galaxy NGC 4244. Our data reliably reach 27.5 R magnitude arcsec^{-2}, a significant improvement on our earlier deep CCD surface photometry of other galaxies. NGC 4244 is a nearby Scd galaxy whose total luminosity is approximately one magnitude fainter than the peak of the Sc luminosity function. We find that it has a simple structure: a single exponential disk, with a scale height h_Z = 246 +/- 2 pc, a scale length h_R = 1.84 +/- 0.02 kpc and a disk cutoff at a radius R(max) = 10.0 kpc (5.4 scale lengths). We confirm a strong cutoff in the stellar disk at R(max), which happens over only 1 kpc. We do not see any statistically significant evidence for disk flaring with radius. Unlike the more luminous Sc galaxies NGC 5907 and M 33, NGC 4244 does not show any evidence for a second component, such as a thick disk or halo, at mu(R) < 27.5 magnitude arcsec^{-2}.

We report on the discovery of Cepheids in the Virgo spiral galaxy NGC 4535, based on observations made with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope. NGC 4535 is one of 18 galaxies observed as a part of The HST Key Project on the Extragalactic Distance Scale which aims to measure the Hubble constant to 10% accuracy. NGC 4535 was observed over 13 epochs using the F555W filter, and over 9 epochs using the F814W filter. The HST F555W and F814W data were transformed to the Johnson V and Kron-Cousins I magnitude systems, respectively. Photometry was performed using two independent programs, DoPHOT and DAOPHOT II/ALLFRAME. Period-luminosity relations in the V and I bands were constructed using 39 high-quality Cepheids present in our set of 50 variable candidates. We obtain a distance modulus of 31.02+/-0.26 mag, corresponding to a distance of 16.0+/-1.9 Mpc. Our distance estimate is based on values of mu = 18.50 +/- 0.10 mag and E(V-I) = 0.13 mag for the distance modulus and reddening of the LMC, respectively.

Using the Hubble Space Telescope (HST) 37 long-period Cepheid variables have been discovered in the Fornax Cluster spiral galaxy NGC 1365. The resulting V and I period-luminosity relations yield a true distance modulus of 31.35 +/- 0.07 mag, which corresponds to a distance of 18.6 +/- 0.6 Mpc. This measurement provides several routes for estimating the Hubble Constant. (1) Assuming this distance for the Fornax Cluster as a whole yields a local Hubble Constant of 70 +/-18_{random} [+/-7]_{systematic} km/s/Mpc. (2) Nine Cepheid-based distances to groups of galaxies out to and including the Fornax and Virgo clusters yield Ho = 73 (+/-16)_r [+/-7]_s km/s/Mpc. (3) Recalibrating the I-band Tully-Fisher relation using NGC 1365 and six nearby spiral galaxies, and applying it to 15 galaxy clusters out to 100 Mpc gives Ho = 76 (+/-3)_r [+/-8]_s km/s/Mpc. (4) Using a broad-based set of differential cluster distance moduli ranging from Fornax to Abell 2147 gives Ho = 72 (+/-)_r [+/-6]_s km/s/Mpc. And finally, (5) Assuming the NGC 1365 distance for the two additional Type Ia supernovae in Fornax and adding them to the SnIa calibration (correcting for light curve shape) gives Ho = 67 (+/-6)_r [+/-7]_s km/s/Mpc out to a distance in excess of 500 Mpc. All five of these Ho determinations agree to within their statistical errors. The resulting estimate of the Hubble Constant combining all these determinations is Ho = 72 (+/-5)_r [+/-12]_s km/s/Mpc.

The distance to NGC 4725 has been derived from Cepheid variables, as part of the Hubble Space Telescope Key Project on the Extragalactic Distance Scale. Thirteen F555W (V) and four F814W (I) epochs of cosmic-ray-split Wide Field and Planetary Camera 2 observations were obtained. Twenty Cepheids were discovered, with periods ranging from 12 to 49 days. Adopting a Large Magellanic Cloud distance modulus and extinction of 18.50+/-0.10 mag and E(V-I)=0.13 mag, respectively, a true reddening-corrected distance modulus (based on an analysis employing the ALLFRAME software package) of 30.50 +/- 0.16 (random) +/- 0.17 (systematic) mag was determined for NGC 4725. The corresponding of distance of 12.6 +/- 1.0 (random) +/- 1.0 (systematic) Mpc is in excellent agreement with that found with an independent analysis based upon the DoPHOT photometry package. With a foreground reddening of only E(V-I)=0.02, the inferred intrinsic reddening of this field in NGC 4725, E(V-I)=0.19, makes it one of the most highly-reddened, encountered by the HST Key Project, to date.

We have constructed color-magnitude diagrams (CMDs) for ten globular clusters in the halo of the nearby spiral galaxy M33 based on Hubble Space Telescope Wide-Field Planetary Camera 2 observations in the F555W (~V) and F814W (~I) filters. These data reveal the morphology of the HB and allow us to estimate the cluster metallicity using the shape and color of the red giant branch. The principal result we report herein is that 8 of the 10 clusters possess exclusively red HB morphologies yet their metallicities are as metal-poor as [Fe/H] ~ -1.6. Indeed, these 8 clusters present basically only giant branch clumps reminiscent of intermediate age star clusters in the Magellanic Clouds. In addition, two of the clusters form a second parameter pair which have similar metallicities but very disparate HB types. Under the assumption that cluster age is the global second parameter, the average age of halo globular clusters in M33 appears to be a few Gyr younger than halo clusters in the Milky Way. Using the observed properties of HB stars in M31 and M33 along with published main-sequence turnoff ages for the globular clusters in the Milky Way, LMC, SMC, and the Sagittarius dwarf spheroidal (Sgr), we attempt to sketch the early formation history of these galaxies. This indicates that the Milky Way, M31, M33, the LMC, and Sgr all experienced their first epoch of cluster formation soon after the Big Bang. Three to four Gyr later, the SMC began to form its first generation of clusters; the bulk of the M33 clusters formed later still. We note that the halo clusters in M33 formed over a much longer time period than those in the Milky Way and M31.

The Hubble Space Telescope is being used to measure accurate Cepheid distances to nearby galaxies with the ultimate aim of determining the Hubble constant, H_0. For the first time, it has become feasible to use Cepheid variables to derive a distance to a galaxy in the southern hemisphere cluster of Fornax. Based on the discovery of 37 Cepheids in the Fornax galaxy NGC 1365, a distance to this galaxy of 18.6 +/- 0.6 Mpc (statistical error only) is obtained. This distance leads to a value of H_0 = 70 +/- 7 (random) +/- 18 (systematic) km/sec/Mpc in good agreement with estimates of the Hubble constant further afield.

We report the detection of Cepheid variable stars in the barred spiral galaxy NGC 1365, located in the Fornax cluster, using the Hubble Space Telescope Wide Field and Planetary Camera 2. Twelve V (F555W) and four I (F814W) epochs of observation were obtained. The two photometry packages, ALLFRAME and DoPHOT, were separately used to obtain profile-fitting photometry of all the stars in the HST field. The search for Cepheid variable stars resulted in a sample of 52 variables, with periods between 14 and 60 days, in common with both datasets. ALLFRAME photometry and light curves of the Cepheids are presented. A subset of 34 Cepheids were selected on the basis of period, light curve shape, similar ALLFRAME and DoPHOT periods, color, and relative crowding, to fit the Cepheid period-luminosity relations in V and I for both ALLFRAME and DoPHOT. The measured distance modulus to NGC 1365 from the ALLFRAME photometry is 31.31 +/- 0.20 (random) +/- 0.18 (systematic) mag, corresponding to a distance of 18.3 +/- 1.7 (random) +/- 1.6 (systematic) Mpc. The reddening is measured to be E(V-I) = 0.16 +/- 0.08 mag. These values are in excellent agreement with those obtained using the DoPHOT photometry, namely a distance modulus of 31.26 +/- 0.10 mag, and a reddening of 0.15 +/- 0.10 mag (internal errors only).

The distance to NGC 7331 has been derived from Cepheid variables observed with HST/WFPC2, as part of the Extragalactic Distance Scale Key Project. Multi-epoch exposures in F555W (V) and F814W (I), with photometry derived independently from DoPHOT and DAOPHOT/ALLFRAME programs, were used to detect a total of 13 reliable Cepheids, with periods between 11 and 42 days. The relative distance moduli between NGC 7331 and the LMC, imply an extinction to NGC 7331 of A_V = 0.47+-0.15 mag, and an extinction-corrected distance modulus to NGC 7331 of 30.89+-0.14(random) mag, equivalent to a distance of 15.1 Mpc. There are additional systematic uncertainties in the distance modulus of +-0.12 mag due to the calibration of the Cepheid Period-Luminosity relation, and a systematic offset of +0.05+-0.04 mag if we applied the metallicity correction inferred from the M101 results of Kennicutt et al 1998.

Uncertainty in the metal abundance dependence of the Cepheid variable period- luminosity (PL) relation remains one of the outstanding sources of systematic error in the extragalactic distance scale and Hubble constant. To test for such a metallicity dependence, we have used the WFPC2 camera on the Hubble Space Telescope (HST) to observe Cepheids in two fields in the nearby spiral galaxy M101, which span a range in oxygen abundance of 0.7+-0.15 dex. A differential analysis of the PL relations in V and I in the two fields yields a marginally significant change in the inferred distance modulus on metal abundance, with d(m-M)/d[O/H] = -0.24+-0.16 mag/dex. The trend is in the theoretically predicted sense that metal-rich Cepheids appear brighter and closer than metal-poor stars. External comparisons of Cepheid distances with those derived from three other distance indicators, in particular the tip of the red giant branch method, further constrain the magnitude of any Z-dependence of the PL relation at V and I. The overall effects of any metallicity dependence on the distance scale derived with HST will be of the order of a few percent or less for most applications, though distances to individual galaxies at the extremes of the metal abundance range may be affected at the 10% level.

We report on the discovery of Cepheids in the field spiral galaxy NGC 3621, based on observations made with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope (HST). NGC 3621 is one of 18 galaxies observed as a part of The HST Key Project on the Extragalactic Distance Scale, which aims to measure the Hubble constant to 10% accuracy. Sixty-nine Cepheids with periods in the range 9--60 days were observed over 12 epochs using the F555W filter, and 4 epochs using the F814W filter. The HST F555W and F814W data were transformed to the Johnson V and Kron-Cousins I magnitude systems, respectively. Photometry was performed using two independent packages, DAOPHOT II/ALLFRAME and DoPHOT. Period-luminosity relations in the V and I bands were constructed using 36 fairly isolated Cepheids present in our set of 69 variables. Extinction-corrected distance moduli relative to the LMC of 10.63 +/- 0.07 mag and 10.56 +/- 0.10 mag were obtained using the ALLFRAME and DoPHOT data, respectively. True distance moduli of 29.13 +/- 0.18 mag and 29.06 +/- 0.18 mag, corresponding to distances of 6.3 Mpc and 6.1 Mpc, were obtained by assuming values of 18.50 +/- 0.10 mag for the distance modulus of the LMC and E(V-I) = 0.13 mag for the reddening of the LMC.

The presence of unseen halos of ``dark matter'' has long been inferred from the high rotation speeds of gas and stars in the outer parts of spiral galaxies$^{1}$. The volume density of this dark matter decreases less quickly from the galactic center than does that of the luminous mass (such as that in stars), meaning that the dark matter dominates the mass far from the center$^{1,2}$. While searching for faint starlight away from the plane of the edge-on disk galaxy \gal$^{3}$, we have found that the galaxy is surrounded by a faint luminous halo. The intensity of light from this halo falls less steeply than any known luminous component of spiral galaxies, but is consistent with the distribution of dark mass inferred from the galaxy's rotation curve.