In this paper we present the results of the first low frequency all-sky search of continuous gravitational wave signals conducted on Virgo VSR2 and VSR4 data. The search covered the full sky, a frequency range between 20 Hz and 128 Hz with a range of spin-down between $-1.0 \times 10^{-10}$ Hz/s and $+1.5 \times 10^{-11}$ Hz/s, and was based on a hierarchical approach. The starting point was a set of short Fast Fourier Transforms (FFT), of length 8192 seconds, built from the calibrated strain data. Aggressive data cleaning, both in the time and frequency domains, has been done in order to remove, as much as possible, the effect of disturbances of instrumental origin. On each dataset a number of candidates has been selected, using the FrequencyHough transform in an incoherent step. Only coincident candidates among VSR2 and VSR4 have been examined in order to strongly reduce the false alarm probability, and the most significant candidates have been selected. Selected candidates have been subject to a follow-up by constructing a new set of longer FFTs followed by a further incoherent analysis, still based on the FrequencyHough transform. No evidence for continuous gravitational wave signals was found, therefore we have set a population-based joint VSR2-VSR4 90$\%$ confidence level upper limit on the dimensionless gravitational wave strain in the frequency range between 20 Hz and 128 Hz. This is the first all-sky search for continuous gravitational waves conducted, on data of ground-based interferometric detectors, at frequencies below 50 Hz. We set upper limits in the range between about $10^{-24}$ and $2\times 10^{-23}$ at most frequencies. Our upper limits on signal strain show an improvement of up to a factor of $\sim$2 with respect to the results of previous all-sky searches at frequencies below $80~\mathrm{Hz}$.

Open-access telescopes of all apertures are needed to operate a competitive and efficient national science program. While larger facilities contribute light-gathering power and angular resolution, smaller ones dominate for field of view, time-resolution, and especially, total available observing time, thereby enabling our entire, diversely-expert community. Smaller aperture telescopes therefore play a critical and indispensable role in advancing science. Thus, the divestment of NSF support for modest-aperture (1 - 4 m) public telescopes poses a serious threat to U.S. scientific leadership, which is compounded by the unknown consequences of the shift from observations driven by individual investigators to survey-driven science. Given the much higher cost efficiency and dramatic science returns for investments in modest aperture telescopes, it is hard to justify funding only the most expensive facilities. We therefore urge the Astro2020 panel to explicitly make the case for modest aperture facilities, and to recommend enhancing this funding stream to support and grow this critical component of the OIR System. Further study is urgently needed to prioritize the numerous exciting potential capabilities of smaller facilities,and to establish sustainable, long-term planning for the System.

We evaluate the optical/near-infrared (OIR) color variability of 3C 279 in both gamma-ray flaring and non-flaring states over 7-year timescales using the Small and Medium Aperture Research Telescope System (SMARTS) in Cerro Tololo, Chile and gamma-ray fluxes obtained from the Fermi Gamma-ray Space Telescope. This observing strategy differs from previous blazar color variability studies in two key ways: 1) the reported color variability is assessed across optical through near-infrared wavelengths, and 2) the color variability is assessed over timescales significantly longer than an individual flare or ground-based observing season. We highlight 3C 279 because of its complex color variability, which is difficult to reconcile with the simple "redder when brighter" behavior often associated with Flat Spectrum Radio Quasar (FSRQ) color variability. We suggest that the observed OIR color changes depend on a combination of the jet and disk emission. We parameterize this behavior in terms of a single variable, $\zeta^m_n$, representing a smooth transition from disk-dominated, to a mixed contribution, to a jet-dominated system, which provides an explanation of the long-term OIR color variability in the same blazar over time. This suggests a general scheme that could apply to OIR color variability in other blazars.

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

We report results of a wideband search for periodic gravitational waves from isolated neutron stars within the Orion spur towards both the inner and outer regions of our Galaxy. As gravitational waves interact very weakly with matter, the search is unimpeded by dust and concentrations of stars. One search disk (A) is $6.87^\circ$ in diameter and centered on $20^\textrm{h}10^\textrm{m}54.71^\textrm{s}+33^\circ33'25.29"$, and the other (B) is $7.45^\circ$ in diameter and centered on $8^\textrm{h}35^\textrm{m}20.61^\textrm{s}-46^\circ49'25.151"$. We explored the frequency range of 50-1500 Hz and frequency derivative from $0$ to $-5\times 10^{-9}$ Hz/s. A multi-stage, loosely coherent search program allowed probing more deeply than before in these two regions, while increasing coherence length with every stage. Rigorous followup parameters have winnowed initial coincidence set to only 70 candidates, to be examined manually. None of those 70 candidates proved to be consistent with an isolated gravitational wave emitter, and 95% confidence level upper limits were placed on continuous-wave strain amplitudes. Near $169$ Hz we achieve our lowest 95% CL upper limit on worst-case linearly polarized strain amplitude $h_0$ of $6.3\times 10^{-25}$, while at the high end of our frequency range we achieve a worst-case upper limit of $3.4\times 10^{-24}$ for all polarizations and sky locations.

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

We present multi-epoch optical spectroscopy of seven southern Fermi-monitored blazars from 2008 - 2013 using the Small and Medium Aperture Research Telescope System (SMARTS), with supplemental spectroscopy and polarization data from the Steward Observatory. We find that the emission lines are much less variable than the continuum; 4 of 7 blazars had no detectable emission line variability over the 5 years. This is consistent with photoionization primarily by an accretion disk, allowing us to use the lines as a probe of disk activity. Comparing optical emission line flux with Fermi $\gamma$-ray flux and optical polarized flux, we investigate whether relativistic jet variability is related to the accretion flow. In general, we see no such dependence, suggesting the jet variability is likely caused by internal processes like turbulence or shock acceleration rather than a variable accretion rate. However, three sources showed statistically significant emission line flares in close temporal proximity to very large Fermi $\gamma$-ray flares. While we do not have sufficient emission line data to quantitatively assess their correlation with the $\gamma$-ray flux, it appears that in some cases, the jet might provide additional photoionizing flux to the broad line region, which implies some gamma-rays are produced within the broad line region, at least for these large flares.

We describe directed searches for continuous gravitational waves in data from the sixth LIGO science data run. The targets were nine young supernova remnants not associated with pulsars; eight of the remnants are associated with non-pulsing suspected neutron stars. One target's parameters are uncertain enough to warrant two searches, for a total of ten. Each search covered a broad band of frequencies and first and second frequency derivatives for a fixed sky direction. The searches coherently integrated data from the two LIGO interferometers over time spans from 5.3-25.3 days using the matched-filtering F-statistic. We found no credible gravitational-wave signals. We set 95% confidence upper limits as strong (low) as $4\times10^{-25}$ on intrinsic strain, $2\times10^{-7}$ on fiducial ellipticity, and $4\times10^{-5}$ on r-mode amplitude. These beat the indirect limits from energy conservation and are within the range of theoretical predictions for neutron-star ellipticities and r-mode amplitudes.

We present results of a search for continuously-emitted gravitational radiation, directed at the brightest low-mass X-ray binary, Scorpius X-1. Our semi-coherent analysis covers 10 days of LIGO S5 data ranging from 50-550 Hz, and performs an incoherent sum of coherent $\mathcal{F}$-statistic power distributed amongst frequency-modulated orbital sidebands. All candidates not removed at the veto stage were found to be consistent with noise at a 1% false alarm rate. We present Bayesian 95% confidence upper limits on gravitational-wave strain amplitude using two different prior distributions: a standard one, with no a priori assumptions about the orientation of Scorpius X-1; and an angle-restricted one, using a prior derived from electromagnetic observations. Median strain upper limits of 1.3e-24 and 8e-25 are reported at 150 Hz for the standard and angle-restricted searches respectively. This proof of principle analysis was limited to a short observation time by unknown effects of accretion on the intrinsic spin frequency of the neutron star, but improves upon previous upper limits by factors of ~1.4 for the standard, and 2.3 for the angle-restricted search at the sensitive region of the detector.

In this paper we present the results of a coherent narrow-band search for continuous gravitational-wave signals from the Crab and Vela pulsars conducted on Virgo VSR4 data. In order to take into account a possible small mismatch between the gravitational wave frequency and two times the star rotation frequency, inferred from measurement of the electromagnetic pulse rate, a range of 0.02 Hz around two times the star rotational frequency has been searched for both the pulsars. No evidence for a signal has been found and 95$\%$ confidence level upper limits have been computed both assuming polarization parameters are completely unknown and that they are known with some uncertainty, as derived from X-ray observations of the pulsar wind torii. For Vela the upper limits are comparable to the spin-down limit, computed assuming that all the observed spin-down is due to the emission of gravitational waves. For Crab the upper limits are about a factor of two below the spin-down limit, and represent a significant improvement with respect to past analysis. This is the first time the spin-down limit is significantly overcome in a narrow-band search.

The remarkable scientific return and legacy of LSST, in the era that it will define, will not only be realized in the breakthrough science that will be achieved with catalog data. This Big Data survey will shape the way the entire astronomical community advances -- or fails to embrace -- new ways of approaching astronomical research and data. In this white paper, we address the NRC template questions 4,5,6,8 and 9, with a focus on the unique challenges for smaller, and often under-resourced, institutions, including institutions dedicated to underserved minority populations, in the efficient and effective use of LSST data products to maximize LSST's scientific return.

We perform time-dependent, spatially-resolved simulations of blazar emission to evaluate several flaring scenarios related to magnetic-field amplification and enhanced particle acceleration. The code explicitly accounts for light-travel-time effects and is applied to flares observed in the flat spectrum radio quasar (FSRQ) PKS 0208-512, which show optical/{\gamma}-ray correlation at some times, but orphan optical flares at other times. Changes in both the magnetic field and the particle acceleration efficiency are explored as causes of flares. Generally, external Compton emission appears to describe the available data better than a synchrotron self-Compton scenario, and in particular orphan optical flares are difficult to produce in the SSC framework. X-ray soft-excesses, {\gamma}-ray spectral hardening, and the detections at very high energies of certain FSRQs during flares find natural explanations in the EC scenario with particle acceleration change. Likewise, optical flares with/without {\gamma}-ray counterparts can be explained by different allocations of energy between the magnetization and particle acceleration, which may be related to the orientation of the magnetic field relative to the jet flow. We also calculate the degree of linear polarization and polarization angle as a function of time for a jet with helical magnetic field. Tightening of the magnetic helix immediately downstream of the jet perturbations, where flares occur, can be sufficient to explain the increases in the degree of polarization and a rotation by >= 180 degree of the observed polarization angle, if light-travel-time effects are properly considered.

We present multi-epoch optical observations of the blazar 3C 454.3 (z = 0.859) from 2008 August through 2011 December, using the SMARTS Consortium 1.5m+RCSpectrograph and 1.3m+ANDICAM in Cerro Tololo, Chile. The spectra reveal that the broad optical emission lines Mg II, H-beta and H-gamma are far less variable than the optical or gamma-ray continuum. Although, the gamma-rays varied by a factor of 100 above the EGRET era flux, the lines generally vary by a factor of 2 or less. Smaller variations in the gamma-ray flux did not produce significant variation in any of the observed emission lines. Therefore, to first order, the ionizing flux from the disk changes only slowly during large variations of the jet. However, two exceptions in the response of the broad emission lines are reported during the largest gamma-ray flares in 2009 December and 2010 November, when significant deviations from the mean line flux in H-gamma and Mg II were observed. H-gamma showed a maximum 3-sigma and 4-sigma deviation in each flare, respectively, corresponding to a factor of 1.7 and 2.5 increase in flux. Mg II showed a 2-sigma deviation in both flares; no variation was detected in H-beta during either flare. These significant deviations from the mean line flux also coincide with 7mm core ejections reported previously (Jorstad et al. 2012). The correlation of the increased emission line flux with mm core ejections, and gamma-ray, optical and UV flares suggests that the broad line region extends beyond the gamma-emitting region during the 2009 and 2010 flares.

We report multi-wavelength (IR-O, UV, X-ray, gamma-ray) data of the flat spectrum radio quasar (FSRQ) PKS B1424-418 (z=1.52) gathered during an active phase in 2013 April . As for a few other cases of FSRQs reported in literature, the detection by the Large Area Telescope onboard Fermi of gamma rays with energy above 10 GeV indicates that the emission likely occurs beyond the highly opaque (tau ~10) broad line region. This conclusion is strengthened by a model fit to the spectral energy distribution, which allows us to locate the emission region even beyond the distance generally assumed for the dusty torus. The consequent large size (~1 pc) inferred for the emission region cannot account for the observed daily-scale variability of the gamma-ray flux. We discuss the possibility that short-term variability results from fast magnetic reconnection events, as proposed in recent works.

The Yale/SMARTS optical-near-IR monitoring program has followed the variations in emission of the Fermi-LAT monitored blazars in the southern sky with closely spaced observations since 2008. We report the discovery of an optical-near-IR (OIR) outburst with no accompanying gamma-rays in the blazar PKS 0208-512, one of the targets of this program. While the source undergoes three outbursts of 1 mag or more at OIR wavelengths lasting for longer than 3 months during 2008-2011, only interval 1 and 3 have corresponding bright phases in GeV energies lasting longer than 1 month. The OIR outburst during interval 2 is comparable in brightness and temporal extent to the OIR flares during intervals 1 and 3 which do have gamma-ray counterparts. Gamma-ray and OIR variability are very well-correlated in most cases in the Fermi blazars and the lack of correlation in this case is anomalous. By analyzing the gamma-ray, OIR, and supporting multi-wavelength variability data in details, we speculate that the location of the outburst in the jet during interval 2 was closer to the black hole where the jet is more compact and the magnetic field strength is higher, and the bulk Lorentz factor of the material in the jet is smaller. These result in a much lower Compton dominance and no observable gamma-ray outburst during interval 2.

We report the discovery of an anomalous flare in a bright blazar, namely, PKS 0208-512, one of the targets of the Yale/SMARTS optical-near-infrared (OIR) monitoring program of Fermi blazars. We identify three intervals during which PKS 0208-512 undergoes outbursts at OIR wavelengths lasting for longer than 3 months. Its brightness increases and then decreases again by at least 1 magnitude in these intervals. In contrast, the source undergoes bright phases in GeV energies lasting for longer than 1 month during intervals 1 and 3 only. The OIR outburst during interval 2 is comparable in brightness and temporal extent to the OIR flares during intervals 1 and 3 which do have gamma-ray counterparts. By analyzing the gamma-ray, OIR, and supporting multi-wavelength variability data in details, we speculate that the OIR outburst during interval 2 was caused by a change in the magnetic field without any change in the total number of emitting electrons or Doppler factor of the emitting region. Alternatively, it is possible that the location of the outburst in the jet during interval 2 was closer to the black hole where the jet is more compact and the bulk Lorentz factor of the material in the jet is smaller. We also discuss the complex OIR spectral behavior during these three intervals.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We present the time variability properties of a sample of six blazars, AO 0235+164, 3C 273, 3C 279, PKS 1510-089, PKS 2155-304, and 3C 454.3, at optical-IR as well as gamma-ray energies. These observations were carried out as a part of the Yale/SMARTS program during 2008-2010 that has followed the variations in emission of the bright Fermi-LAT-monitored blazars in the southern sky with closely-spaced observations at BVRJK bands. We find the optical/IR time variability properties of these blazars to be remarkably similar to those at the gamma-ray energies. The power spectral density (PSD) functions of the R-band variability of all six blazars are fit well by simple power-law functions with negative slope such that there is higher amplitude variability on longer timescales. No clear break is identified in the PSD of any of the sources. The average slope of the PSD of R-band variability of these blazars is similar to what was found by the Fermi team for the gamma-ray variability of a larger sample of bright blazars. This is consistent with leptonic models where the optical-IR and gamma-ray emission is generated by the same population of electrons through synchrotron and inverse-Compton processes, respectively. The prominent flares present in the optical-IR as well as the gamma-ray light curves of these blazars are predominantly symmetric, i.e., have similar rise and decay timescales, indicating that the long-term variability is dominated by the crossing time of radiation or a disturbance through the emission region rather than by the acceleration or energy-loss timescales of the radiating electrons. In the blazar 3C 454.3, which has the highest-quality light curves, the location of a large gamma-ray outburst during 2009 December is consistent with being in the jet at ~18 pc from the central engine. This poses strong constraints on the models of high energy emission in the jets of blazars.

We present multiwavelength data for twelve blazars observed from 2008-2010 as part of an ongoing optical-infrared photometric monitoring project. Sources were selected to be bright, southern (dec < 20 deg) blazars observed by the Fermi Gamma-Ray Space Telescope, with daily and weekly gamma-ray fluxes made available from the start of the Fermi mission. Light curves are presented for the twelve blazars in BVRJK at near-daily cadence. We find that optical and infrared fluxes are well correlated in all sources. Gamma-ray bright flat spectrum radio quasars (FSRQs) in our sample have optical/infrared emission correlated with gamma-rays consistent with inverse Compton-scattering models for GeV emission. In FSRQs, the variability amplitude decreases towards optical/IR wavelengths, consistent with the presence of a thermal emission component from the accretion disk varying on significantly longer timescales than the jet synchrotron emission. In BL Lac objects, variability is mainly constant across wavelengths, consistent with a weak or radiatively inefficient disk. FSRQs have redder optical-infrared colors when they are brighter, while BL Lac objects show no such trend. Several objects show complicated color-magnitude behavior: AO 0235+164 appears in two different states depending on whether it is gamma-ray bright or not. OJ 287 and 3C 279 show some hysteresis tracks in their color-magnitude diagrams. Individual flares may be achromatic or otherwise depart from the trend, suggesting different jet components becoming important at different times. We present a time-dependent spectral energy distribution of the bright FSRQ 3C 454.3 during its December 2009 flare, which is well fit by an external Compton model in the bright state, although day to day changes pose challenges to a simple one-zone model. All data from the SMARTS monitoring program are publicly available on our website.

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

The blazar 3C 454.3 was revealed by the Fermi Gamma-ray Space Telescope to be in an exceptionally high flux state in July 2008. Accordingly, we performed a multi-wavelength monitoring campaign on this blazar using IR and optical observations from the SMARTS telescopes, optical, UV and X-ray data from the Swift satellite, and public-release gamma-ray data from Fermi. We find an excellent correlation between the IR, optical, UV and gamma-ray light curves, with a time lag of less than one day. The amplitude of the infrared variability is comparable to that in gamma-rays, and larger than at optical or UV wavelengths. The X-ray flux is not strongly correlated with either the gamma-rays or longer wavelength data. These variability characteristics find a natural explanation in the external Compton model, in which electrons with Lorentz factor gamma~10^(3-4) radiate synchrotron emission in the infrared-optical and also scatter accretion disk or emission line photons to gamma-ray energies, while much cooler electrons (gamma~10^(1-2)) produce X-rays by scattering synchrotron or other ambient photons.

We present the discovery of a kinematically-cold stellar population along the SE minor axis of the Andromeda galaxy (M31) that is likely the forward continuation of M31's giant southern stream. This discovery was made in the course of an on-going spectroscopic survey of red giant branch (RGB) stars in M31 using the DEIMOS instrument on the Keck II 10-m telescope. Stellar kinematics are investigated in eight fields located 9-30 kpc from M31's center (in projection). A likelihood method based on photometric and spectroscopic diagnostics is used to isolate confirmed M31 RGB stars from foreground Milky Way dwarf stars: for the first time, this is done without using radial velocity as a selection criterion, allowing an unbiased study of M31's stellar kinematics. The radial velocity distribution of the 1013 M31 RGB stars shows evidence for the presence of two components. The broad (hot) component has a velocity dispersion of 129 km/s and presumably represents M31's virialized spheroid. A significant fraction (19%) of the population is in a narrow (cold) component centered near M31's systemic velocity with a velocity dispersion that decreases with increasing radial distance, from 55.5 km/s at R_proj=12 kpc to 10.6 km/s at R_proj=18 kpc. The spatial and velocity distribution of the cold component matches that of the "Southeast shelf" predicted by the Fardal et al. (2007) orbital model of the progenitor of the giant southern stream. The metallicity distribution of the cold component matches that of the giant southern stream, but is about 0.2 dex more metal rich on average than that of the hot spheroidal component. We discuss the implications of our discovery on the interpretation of the intermediate-age spheroid population found in this region in recent ultra-deep HST imaging studies.