The Legacy Surveys, a combination of three ground-based imaging surveys, have mapped 16,000 deg$^2$ in three optical bands ($g$, $r$, and $z$) to a depth 1--$2$~mag deeper than the Sloan Digital Sky Survey (SDSS). Our work addresses one of the major challenges of wide-field imaging surveys conducted at ground-based observatories: the varying depth that results from varying observing conditions at Earth-bound sites. To mitigate these effects, two of the Legacy Surveys (the Dark Energy Camera Legacy Survey, or DECaLS; and the Mayall $z$-band Legacy Survey, or MzLS) employed a unique strategy to dynamically adjust the exposure times as rapidly as possible in response to the changing observing conditions. We present the tiling and observing strategies used by these surveys. We demonstrate that the tiling and dynamic observing strategies jointly result in a more uniform-depth survey that has higher efficiency for a given total observing time compared with the traditional approach of using fixed exposure times.

We present the discoveries of WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5, two low-temperature (1200$-$1400 K), high proper motion T-type subdwarfs. Both objects were discovered via their high proper motion ($>$0.5 arcsec yr$^{-1}$); WISEA J181006.18-101000.5 as part of the NEOWISE proper motion survey and WISEA J041451.67-585456.7 as part of the citizen science project Backyard Worlds; Planet 9. We have confirmed both as brown dwarfs with follow-up near-infrared spectroscopy. Their spectra and near-infrared colors are unique amongst known brown dwarfs, with some colors consistent with L-type brown dwarfs and other colors resembling those of the latest-type T dwarfs. While no forward model consistently reproduces the features seen in their near-infrared spectra, the closest matches suggest very low metallicities ([Fe/H] $\leq$ -1), making these objects likely the first examples of extreme subdwarfs of the T spectral class (esdT). WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5 are found to be part of a small population of objects that occupy the "substellar transition zone," and have the lowest masses and effective temperatures of all objects in this group.

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

We search in the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys for new strong lensing systems by using deep residual neural networks, building on previous work presented in Huang et al. These surveys together cover approximately one third of the sky visible from the northern hemisphere, reaching a z-band magnitude of ~22.5. We compile a training sample that consists of known lensing systems as well as non-lenses in the Legacy Surveys and the Dark Energy Survey. After applying our trained neural networks to the survey data, we visually inspect and rank images with probabilities above a threshold. Here we present 1014 new strong lens candidates.

We present the discovery of WISEA J083011.95+283716.0, the first Y dwarf candidate identified through the Backyard Worlds: Planet 9 citizen science project. We identified this object as a red, fast-moving source with a faint $W2$ detection in multi-epoch \textit{AllWISE} and unWISE images. We have characterized this object with Spitzer Space Telescope and \textit{Hubble Space Telescope} follow-up imaging. With mid-infrared detections in \textit{Spitzer}'s \emph{ch1} and \emph{ch2} bands and flux upper limits in Hubble Space Telescope $F105W$ and $F125W$ filters, we find that this object is both very faint and has extremely red colors ($ch1-ch2 = 3.25\pm0.23$ mag, $F125W-ch2 \geq 9.36$ mag), consistent with a T$_{eff}\sim300$ K source, as estimated from the known Y dwarf population. A preliminary parallax provides a distance of $11.1^{+2.0}_{-1.5}$ pc, leading to a slightly warmer temperature of $\sim350$ K. The extreme faintness and red Hubble Space Telescope and Spitzer Space Telescope colors of this object suggest it may be a link between the broader Y dwarf population and the coldest known brown dwarf WISE J0855$-$0714, and highlight our limited knowledge of the true spread of Y dwarf colors. We also present four additional Backyard Worlds: Planet 9 late-T brown dwarf discoveries within 30 pc.

CatWISE is a program to catalog sources selected from combined ${\it WISE}$ and ${\it NEOWISE}$ all-sky survey data at 3.4 and 4.6 $\mu$m (W1 and W2). The CatWISE Preliminary Catalog consists of 900,849,014 sources measured in data collected from 2010 to 2016. This dataset represents four times as many exposures and spans over ten times as large a time baseline as that used for the AllWISE Catalog. CatWISE adapts AllWISE software to measure the sources in coadded images created from six-month subsets of these data, each representing one coverage of the inertial sky, or epoch. The catalog includes the measured motion of sources in 8 epochs over the 6.5 year span of the data. From comparison to ${\it Spitzer}$, the SNR=5 limits in magnitudes in the Vega system are W1=17.67 and W2=16.47, compared to W1=16.96 and W2=16.02 for AllWISE. From comparison to ${\it Gaia}$, CatWISE positions have typical accuracies of 50 mas for stars at W1=10 mag and 275 mas for stars at W1=15.5 mag. Proper motions have typical accuracies of 10 mas yr$^{-1}$ and 30 mas yr$^{-1}$ for stars with these brightnesses, an order of magnitude better than from AllWISE. The catalog is available in the WISE/NEOWISE Enhanced and Contributed Products area of the NASA/IPAC Infrared Science Archive.

We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging dataset, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The LRG sample contains 2.7 million objects over the redshift range of $0.4 < z < 0.9$ over 5655 sq. degrees. We have developed new photometric redshift (photo-$z$) estimates using the Legacy Survey DECam and WISE photometry, with $\sigma_{\mathrm{NMAD}} = 0.02$ precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a 5-parameter HOD model fit using the MultiDark halo catalog, we find that there is little evolution in HOD parameters except at the highest-redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of MCMC fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogs, and present a new technique for finding best-fit parameters in this situation. Accompanying this paper we have released the PRLS (Photometric Redshifts for the Legacy Surveys) catalog of photo-$z$'s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 dataset. This catalog provides accurate photometric redshifts for objects with $z < 21$ over more than 16,000 square degrees of sky.

We present follow-up $Spitzer$ observations at 3.6$\mu$m (ch1) and 4.5$\mu$m (ch2) of CWISEP J144606.62$-$231717.8, one of the coldest known brown dwarfs in the solar neighborhood. This object was found by mining the Wide-field Infrared Survey Explorer ($WISE$) and $NEOWISE$ data via the CatWISE Preliminary Catalog by Meisner et al. (2019b), where an initial $Spitzer$ color of ch1$-$ch2 = 3.71$\pm$0.44 mag was reported, implying it could be one of the reddest, and hence coldest, known brown dwarfs. Additional $Spitzer$ data presented here allows us to revise its color to ch1$-$ch2 = 2.986$\pm$0.048 mag, which makes CWISEP J144606.62$-$231717.8 the 5th reddest brown dwarf ever observed. A preliminary trigonometric parallax measurement, based on a combination of $WISE$ and $Spitzer$ astrometry, places this object at a distance of 10.1$^{+1.7}_{-1.3}$ pc. Based on our improved $Spitzer$ color and preliminary parallax, CWISEP J144606.62$-$231717.8 has a $T_{\rm eff}$ in the 310$-$360 K range. Assuming an age of 0.5$-$13 Gyr, this corresponds to a mass between 2 and 20 $M_{\rm Jup}$.

We present Spitzer 3.6$\mu$m and 4.5$\mu$m follow-up of 170 candidate extremely cool brown dwarfs newly discovered via the combination of WISE and NEOWISE imaging at 3$-$5$\mu$m. CatWISE, a joint analysis of archival WISE and NEOWISE data, has improved upon the motion measurements of AllWISE by leveraging a $>$10$\times$ time baseline enhancement, from 0.5 years (AllWISE) to 6.5 years (CatWISE). As a result, CatWISE motion selection has yielded a large sample of previously unrecognized brown dwarf candidates, many of which have archival detections exclusively in the WISE 4.6$\mu$m (W2) channel, suggesting that they could be both exceptionally cold and nearby. Where these objects go undetected in WISE W1 (3.4$\mu$m), Spitzer can provide critically informative detections at 3.6$\mu$m. Of our motion-confirmed discoveries, seventeen have a best-fit Spitzer [3.6]$-$[4.5] color most consistent with spectral type Y. CWISEP J144606.62$-$231717.8 ($\mu \approx 1.3''$/yr) is likely the reddest, and therefore potentially coldest, member of our sample with a very uncertain [3.6]$-$[4.5] color of 3.71 $\pm$ 0.44 magnitudes. We also highlight our highest proper motion discovery, WISEA J153429.75$-$104303.3, with $\mu \approx 2.7''$/yr. Given that the prior list of confirmed and presumed Y dwarfs consists of just 27 objects, the Spitzer follow-up presented in this work has substantially expanded the sample of identified Y dwarfs. Our new discoveries thus represent significant progress toward understanding the bottom of the substellar mass function, investigating the diversity of the Y dwarf population, and selecting optimal brown dwarf targets for JWST spectroscopy.

We report the discovery of WISE2150-7520AB (W2150AB): a widely separated (~ 341 AU) very low mass L1 + T8 co-moving system. The system consists of the previously known L1 primary 2MASS J21501592-7520367 and a newly discovered T8 secondary found at position 21:50:18.99 -75:20:54.6 (MJD=57947) using Wide-field Infrared Survey Explorer (WISE) data via the Backyard Worlds: Planet 9 citizen science project. We present Spitzer ch1 and ch2 photometry (ch1-ch2= 1.41 +/-0.04 mag) of the secondary and FIRE prism spectra of both components. The sources show no peculiar spectral or photometric signatures indicating that each component is likely field age. Using all observed data and the Gaia DR2 parallax of 41.3593 +/- 0.2799 mas for W2150A we deduce fundamental parameters of log(Lbol/Lsun)=-3.69 +/- 0.01, Teff=2118 +/- 62 K, and an estimated mass=72 +/- 12 MJup for the L1 and log(Lbol/Lsun)=-5.64 +/- 0.02, Teff=719 +/- 61 K, and an estimated mass=34 +/- 22 MJup for the T8. At a physical separation of ~341 AU this system has Ebin = 10^41 erg making it the lowest binding energy system of any pair with Mtot < 0.1 Msun not associated with a young cluster. It is equivalent in estimated mass ratio, Ebin, and physical separation to the ~ 2 Myr M7.25 + M8.25 binary brown dwarf 2MASS J11011926-7732383AB (2M1101AB) found in the Chameleon star forming region. W2150AB is the widest companion system yet observed in the field where the primary is an L dwarf or later.

We present full-sky coadded maps created by uniformly combining the first five years of Wide-field Infrared Survey Explorer (WISE) and NEOWISE imaging at 3.4 microns (W1) and 4.6 microns (W2). By incorporating both pre-hibernation WISE exposures from 2010-2011 and the first four years (2013-2017) of post-hibernation exposures from the NEOWISE-Reactivation mission, we are able to provide W1/W2 coadds that span a 15 times longer time baseline and are substantially deeper than the standard AllWISE data products. Our new five-year "full-depth" coadds are now the deepest ever all-sky maps at 3-5 microns, permitting detection of sources ~2 times (~0.7 mag) fainter than AllWISE at 5 sigma significance. We additionally present an updated set of "time-resolved" W1/W2 coadds, which separately stack each of ~10 sky passes at each inertial sky location, enabling motion and variability measurements for faint infrared sources over a long ~7.5 year time baseline. We highlight new processing improvements relative to our previous "unWISE" coadd releases, focusing on astrometric calibration and artifact flagging. The deep WISE stacks presented here are already being used to perform target selection for the Dark Energy Spectroscopic Instrument (DESI), and our full-sky coadded WISE/NEOWISE products will be key precursor data sets for upcoming wide-field infrared missions including SPHEREx and NEOCam.

We present the discovery of an extremely cold, nearby brown dwarf in the solar neighborhood, found in the CatWISE catalog (Eisenhardt et al., in prep.). Photometric follow-up with Spitzer reveals that the object, CWISEP J193518.59-154620.3, has ch1$-$ch2 = 3.24$\,\pm\,$0.31 mag, making it one of the reddest brown dwarfs known. Using the Spitzer photometry and the polynomial relations from Kirkpatrick et al. (2019) we estimate an effective temperature in the $\sim$270--360 K range, and a distance estimate in the 5.6$-$10.9 pc range. We combined the WISE, NEOWISE, and Spitzer data to measure a proper motion of $\mu_\alpha \cos \delta = 337\pm69$ mas yr$^{-1}$, $\mu_\delta = -50\pm97$ mas yr$^{-1}$, which implies a relatively low tangential velocity in the range 7$-$22 km s$^{-1}$.

We perform a semi-automated search for strong gravitational lensing systems in the 9,000 deg$^2$ Dark Energy Camera Legacy Survey (DECaLS), part of the DESI Legacy Imaging Surveys (Dey et al.). The combination of the depth and breadth of these surveys are unparalleled at this time, making them particularly suitable for discovering new strong gravitational lensing systems. We adopt the deep residual neural network architecture (He et al.) developed by Lanusse et al. for the purpose of finding strong lenses in photometric surveys. We compile a training set that consists of known lensing systems in the Legacy Surveys and DES as well as non-lenses in the footprint of DECaLS. In this paper we show the results of applying our trained neural network to the cutout images centered on galaxies typed as ellipticals (Lang et al.) in DECaLS. The images that receive the highest scores (probabilities) are visually inspected and ranked. Here we present 335 candidate strong lensing systems, identified for the first time.

Infrared excesses due to dusty disks have been observed orbiting white dwarfs with effective temperatures between 7200 K and 25000 K, suggesting that the rate of tidal disruption of minor bodies massive enough to create a coherent disk declines sharply beyond 1~Gyr after white dwarf formation. We report the discovery that the candidate white dwarf LSPM J0207+3331, via the Backyard Worlds: Planet 9 citizen science project and Keck Observatory follow-up spectroscopy, is hydrogen-dominated with a luminous compact disk (L$_{\rm IR}$/L$_{\star}$=14%) and an effective temperature nearly 1000K cooler than any known white dwarf with an infrared excess. The discovery of this object places the latest time for large scale tidal disruption events to occur at $\sim$3 Gyr past the formation of the host white dwarf, making new demands of dynamical models for planetesimal perturbation and disruption around post main sequence planetary systems. Curiously, the mid-IR photometry of the disk cannot be fully explained by a geometrically thin, optically thick dust disk as seen for other dusty white dwarfs, but requires a second ring of dust near the white dwarf's Roche radius. In the process of confirming this discovery, we found that careful measurements of WISE source positions can reveal when infrared excesses for white dwarfs are co-moving with their hosts, helping distinguish them from confusion noise.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We present the unWISE Catalog, containing the positions and fluxes of roughly two billion objects observed by the Wide-field Infrared Survey Explorer (WISE) over the full sky. The unWISE Catalog has two advantages over the existing WISE catalog (AllWISE): first, it is based on significantly deeper imaging, and second, it features improved modeling of crowded regions. The deeper imaging used in the unWISE Catalog comes from the coaddition of all publicly available 3$-$5 micron WISE imaging, including that from the ongoing NEOWISE-Reactivation mission, thereby increasing the total exposure time by a factor of 5 relative to AllWISE. At these depths, even at high Galactic latitudes many sources are blended with their neighbors; accordingly, the unWISE analysis simultaneously fits thousands of sources to obtain accurate photometry. Our new catalog detects sources at 5-sigma roughly 0.7 magnitudes fainter than the AllWISE catalog and more accurately models millions of faint sources in the Galactic plane, enabling a wealth of Galactic and extragalactic science. In particular, relative to AllWISE, unWISE doubles the number of galaxies detected between redshifts 0 and 1 and triples the number between redshifts 1 and 2, cataloging more than half a billion galaxies over the whole sky.

Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear more likely than changes in obscuration. Here we report on four epochs of spectroscopy of SDSS J110057.70-005304.5, a quasar at a redshift of $z=0.378$ whose UV continuum and broad hydrogen emission lines have faded, and then returned over the past $\approx$20 years. The change in this quasar was initially identified in the infrared, and an archival spectrum from 2010 shows an intermediate phase of the transition during which the flux below rest-frame $\approx$3400\AA\ has decreased by close to an order of magnitude. This combination is unique compared to previously published examples of changing-look quasars, and is best explained by dramatic changes in the innermost regions of the accretion disk. The optical continuum has been rising since mid-2016, leading to a prediction of a rise in hydrogen emission line flux in the next year. Increases in the infrared flux are beginning to follow, delayed by a $\sim$3 year observed timescale. If our model is confirmed, the physics of changing-look quasars are governed by processes at the innermost stable circular orbit (ISCO) around the black hole, and the structure of the innermost disk. The easily identifiable and monitored changing-look quasars would then provide a new probe and laboratory of the nuclear central engine.

We have used the first ~3 years of 3.4 micron (W1) and 4.6 micron (W2) observations from the WISE and NEOWISE missions to create a full-sky set of time-resolved coadds. As a result of the WISE survey strategy, a typical sky location is visited every six months and is observed during 12 or more exposures per visit, with these exposures spanning a ~1 day time interval. We have stacked the exposures within such ~1 day intervals to produce one coadd per band per visit -- that is, one coadd every six months at a given position on the sky in each of W1 and W2. For most parts of the sky we have generated six epochal coadds per band, with one visit during the fully cryogenic WISE mission, one visit during NEOWISE, and then, after a 33 month gap, four more visits during the NEOWISE-Reactivation mission phase. These coadds are suitable for studying long-timescale mid-infrared variability and measuring motions to ~1.3 magnitudes fainter than the single-exposure detection limit. In most sky regions, our coadds span a 5.5 year time period and therefore provide a >10x enhancement in time baseline relative to that available for the AllWISE catalog's apparent motion measurements. As such, the signature application of these new coadds is expected to be motion-based identification of relatively faint brown dwarfs, especially those cold enough to remain undetected by Gaia.

We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 \mu$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$\beta$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios.

We have uniformly reprocessed ~140 terabytes of WISE and NEOWISE exposures to create the deepest ever full-sky maps at 3.4 microns (W1) and 4.6 microns (W2). Our coadds include ~4 years of observations and therefore feature ~4 times greater integer frame coverage than the AllWISE Atlas stacks. Our new, publicly available maps should find a wide range of applications, and in particular will enable the selection of luminous red galaxy and quasar targets for the Dark Energy Spectroscopic Instrument (DESI).

The recent 'Planet Nine' hypothesis has led to many observational and archival searches for this giant planet proposed to orbit the Sun at hundreds of astronomical units. While trans-Neptunian object searches are typically conducted in the optical, models suggest Planet Nine could be self-luminous and potentially bright enough at ~3-5 microns to be detected by the Wide-field Infrared Survey Explorer (WISE). We have previously demonstrated a Planet Nine search methodology based on time-resolved WISE coadds, allowing us to detect moving objects much fainter than would be possible using single-frame extractions. In the present work, we extend our 3.4 micron (W1) search to cover more than three quarters of the sky and incorporate four years of WISE observations spanning a seven year time period. This represents the deepest and widest-area WISE search for Planet Nine to date. We characterize the spatial variation of our survey's sensitivity and rule out the presence of Planet Nine in the parameter space searched at W1 < 16.7 in high Galactic latitude regions (90% completeness).

The DECam Plane Survey is a five-band optical and near-infrared survey of the southern Galactic plane with the Dark Energy Camera at Cerro Tololo. The survey is designed to reach past the main-sequence turn-off at the distance of the Galactic center through a reddening E(B-V) of 1.5 mag. Typical single-exposure depths are 23.7, 22.8, 22.3, 21.9, and 21.0 mag in the grizY bands, with seeing around 1 arcsecond. The footprint covers the Galactic plane with |b| < 4 degrees, 5 degrees > l > -120 degrees. The survey pipeline simultaneously solves for the positions and fluxes of tens of thousands of sources in each image, delivering positions and fluxes of roughly two billion stars with better than 10 mmag precision. Most of these objects are highly reddened and deep in the Galactic disk, probing the structure and properties of the Milky Way and its interstellar medium. The full survey is publicly available.

We have reprocessed over 100 terabytes of single-exposure WISE/NEOWISE images to create the deepest ever full-sky maps at 3-5 microns. We incorporate all publicly available W1 and W2 imaging - a total of ~8 million exposures in each band - from ~37 months of observations spanning 2010 January to 2015 December. Our coadds preserve the native WISE resolution and feature depth of coverage ~3 times greater than that of the AllWISE Atlas stacks. Our coadds are designed to enable deep forced photometry, in particular for the Dark Energy Camera Legacy Survey (DECaLS) and Mayall z-Band Legacy Survey (MzLS), both of which are being used to select targets for the Dark Energy Spectroscopic Instrument (DESI). We describe newly introduced processing steps aimed at leveraging added redundancy to remove artifacts, with the intent of facilitating uniform target selection and searches for rare/exotic objects (e.g. high-redshift quasars and distant galaxy clusters). Forced photometry depths achieved with these coadds extend 0.56 (0.46) magnitudes deeper in W1 (W2) than is possible with only pre-hibernation WISE imaging.

The Wide-field Infrared Survey Explorer (WISE) is a powerful tool for finding nearby brown dwarfs and searching for new planets in the outer solar system, especially with the incorporation of NEOWISE and NEOWISE-Reactivation data. So far, searches for brown dwarfs in WISE data have yet to take advantage of the full depth of the WISE images. To efficiently search this unexplored space via visual inspection, we have launched a new citizen science project, called "Backyard Worlds: Planet 9," which asks volunteers to examine short animations composed of difference images constructed from time-resolved WISE coadds. We report the discovery of the first new substellar object found by this project, WISEA J110125.95+540052.8, a T5.5 brown dwarf located approximately 34 pc from the Sun with a total proper motion of $\sim$0.7 as yr$^{-1}$. WISEA J110125.95+540052.8 has a WISE $W2$ magnitude of $W2=15.37 \pm 0.09$, this discovery demonstrates the ability of citizen scientists to identify moving objects via visual inspection that are 0.9 magnitudes fainter than the $W2$ single-exposure sensitivity, a threshold that has limited prior motion-based brown dwarf searches with WISE.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

A distant, as yet unseen ninth planet has been invoked to explain various observations of the outer solar system. While such a 'Planet Nine', if it exists, is most likely to be discovered via reflected light in the optical, it may emit much more strongly at 3$-$5$\mu$m than simple blackbody predictions would suggest, depending on its atmospheric properties (Fortney et al. 2016). As a result, Planet Nine may be detectable at 3.4$\mu$m with WISE, but single exposures are too shallow except at relatively small distances ($d_9 \lesssim 430$ AU). We develop a method to search for Planet Nine far beyond the W1 single-exposure sensitivity, to distances as large as 800 AU, using inertial coadds of W1 exposures binned into $\sim$1 day intervals. We apply our methodology to $\sim$2000 square degrees of sky identified by Holman & Payne (2016) as a potentially likely Planet Nine location, based on the Fienga et al. (2016) Cassini ranging analysis. We do not detect a plausible Planet Nine candidate, but are able to derive a detailed completeness curve, ruling out its presence within the parameter space searched at $W1 < 16.66$ (90% completeness). Our method uses all publicly available W1 imaging, spanning 2010 January to 2015 December, and will become more sensitive with future NEOWISE-Reactivation releases of additional W1 exposures. We anticipate that our method will be applicable to the entire high Galactic latitude sky, and we will extend our search to that full footprint in the near future.

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

The Near Earth Object Wide-field Infrared Survey Explorer (NEOWISE) Reactivation mission released data from its first full year of observations in 2015. This data set includes ~2.5 million exposures in each of W1 and W2, effectively doubling the amount of WISE imaging available at 3.4 and 4.6 microns relative to the AllWISE release. We have created the first ever full-sky set of coadds combining all publicly available W1 and W2 exposures from both the AllWISE and NEOWISE-Reactivation (NEOWISER) mission phases. We employ an adaptation of the unWISE image coaddition framework (Lang 2014), which preserves the native WISE angular resolution and is optimized for forced photometry. By incorporating two additional scans of the entire sky, we not only improve the W1/W2 depths, but also largely eliminate time-dependent artifacts such as off-axis scattered moonlight. We anticipate that our new coadds will have a broad range of applications, including target selection for upcoming spectroscopic cosmology surveys, identification of distant/massive galaxy clusters, and discovery of high-redshift quasars. In particular, our full-depth AllWISE+NEOWISER coadds will be an important input for the Dark Energy Spectroscopic Instrument (DESI) selection of luminous red galaxy and quasar targets. Our full-depth W1/W2 coadds are already in use within the DECam Legacy Survey (DECaLS) and Mayall z-band Legacy Survey (MzLS) reduction pipelines. Much more work still remains in order to fully leverage NEOWISER imaging for astrophysical applications beyond the solar system.

The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity.

We employ an all-sky map of the anomalous microwave emission (AME) produced by component separation of the microwave sky to study correlations between the AME and Galactic dust properties. We find that while the AME is highly correlated with all tracers of dust emission, the best predictor of the AME strength is the dust radiance. Fluctuations in the AME intensity per dust radiance are uncorrelated with fluctuations in the emission from polycyclic aromatic hydrocarbons (PAHs), casting doubt on the association between AME and PAHs. The PAH abundance is strongly correlated with the dust optical depth and dust radiance, consistent with PAH destruction in low density regions. We find that the AME intensity increases with increasing radiation field strength, at variance with predictions from the spinning dust hypothesis. Finally, the temperature-dependence of the AME per dust radiance disfavors the interpretation of the AME as thermal emission. A reconsideration of other AME carriers, such as ultrasmall silicates, and other emission mechanisms, such as magnetic dipole emission, is warranted.

During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.

We apply the Finkbeiner et al. (1999) two-component thermal dust emission model to the Planck HFI maps. This parametrization of the far-infrared dust spectrum as the sum of two modified blackbodies serves as an important alternative to the commonly adopted single modified blackbody (MBB) dust emission model. Analyzing the joint Planck/DIRBE dust spectrum, we show that two-component models provide a better fit to the 100-3000 GHz emission than do single-MBB models, though by a lesser margin than found by Finkbeiner et al. (1999) based on FIRAS and DIRBE. We also derive full-sky 6.1' resolution maps of dust optical depth and temperature by fitting the two-component model to Planck 217-857 GHz along with DIRBE/IRAS 100 micron data. Because our two-component model matches the dust spectrum near its peak, accounts for the spectrum's flattening at millimeter wavelengths, and specifies dust temperature at 6.1' FWHM, our model provides reliable, high-resolution thermal dust emission foreground predictions from 100 to 3000 GHz. We find that, in diffuse sky regions, our two-component 100-217 GHz predictions are on average accurate to within 2.2%, while extrapolating the Planck Collaboration (2013a) single-MBB model systematically underpredicts emission by 18.8% at 100 GHz, 12.6% at 143 GHz and 7.9% at 217 GHz. We calibrate our two-component optical depth to reddening, and compare with reddening estimates based on stellar spectra. We find the dominant systematic problems in our temperature/reddening maps to be zodiacal light on large angular scales and the cosmic infrared background anisotropy on small angular scales.

We describe our custom processing of the entire Wide-field Infrared Survey Explorer (WISE) 12 micron imaging data set, and present a high-resolution, full-sky map of diffuse Galactic dust emission that is free of compact sources and other contaminating artifacts. The principal distinctions between our resulting co-added images and the WISE Atlas stacks are our removal of compact sources, including their associated electronic and optical artifacts, and our preservation of spatial modes larger than 1.5 degrees. We provide access to the resulting full-sky map via a set of 430 12.5 degree by 12.5 degree mosaics. These stacks have been smoothed to 15" resolution and are accompanied by corresponding coverage maps, artifact images, and bit-masks for point sources, resolved compact sources, and other defects. When combined appropriately with other mid-infrared and far-infrared data sets, we expect our WISE 12 micron co-adds to form the basis for a full-sky dust extinction map with angular resolution several times better than Schlegel et al. (1998).

We present low-resolution MMT Hectospec spectroscopy of 594 candidate Monoceros stream member stars. Based on strong color-magnitude diagram overdensities, we targeted three fields within the stream's footprint, with 178 deg < l < 203 deg and -25 deg < b < 25 deg. By comparing the measured iron abundances with those expected from smooth Galactic components alone, we measure, for the first time, the spectroscopic metallicity distribution function for Monoceros. We find the stream to be chemically distinct from both the thick disk and halo, with [Fe/H] = -1, and do not detect a trend in the stream's metallicity with Galactic longitude. Passing from b = +25 deg to b = -25 deg the median Monoceros metallicity trends upward by 0.1 dex, though uncertainties in modeling sample contamination by the disk and halo make this a marginal detection. In each field, we find Monoceros to have an intrinsic [Fe/H] dispersion of 0.10-0.22 dex. From the CaII K line, we measure [Ca/Fe] for a subsample of metal poor program stars with -1.1 < [Fe/H] < -0.5. In two of three fields, we find calcium deficiencies qualitatively similar to previously reported [Ti/Fe] underabundances in Monoceros and the Sagittarius tidal stream. Further, using 90 spectra of thick disk stars in the Monoceros pointings with b ~ +/-25 deg, we detect a 0.22 dex north/south metallicity asymmetry coincident with known stellar density asymmetry at R_GC ~ 12 kpc and |Z| ~ 1.7 kpc. Our median Monoceros [Fe/H] = -1.0 and its relatively low dispersion naturally fit the expectation for an appropriately luminous M_V ~ -13 dwarf galaxy progenitor.

We develop a method for separating quasars from other variable point sources using SDSS Stripe 82 light curve data for ~10,000 variable objects. To statistically describe quasar variability, we use a damped random walk model parametrized by a damping time scale, tau, and an asymptotic amplitude (structure function), SF_inf. With the aid of an SDSS spectroscopically confirmed quasar sample, we demonstrate that variability selection in typical extragalactic fields with low stellar density can deliver complete samples with reasonable purity (or efficiency, E). Compared to a selection method based solely on the slope of the structure function, the inclusion of the tau information boosts E from 60% to 75% while maintaining a highly complete sample (98%) even in the absence of color information. For a completeness of C=90%, E is boosted from 80% to 85%. Conversely, C improves from 90% to 97% while maintaining E=80% when imposing a lower limit on tau. With the aid of color selection, the purity can be further boosted to 96%, with C= 93%. Hence, selection methods based on variability will play an important role in the selection of quasars with data provided by upcoming large sky surveys, such as Pan-STARRS and the Large Synoptic Survey Telescope (LSST). For a typical (simulated) LSST cadence over 10 years and a photometric accuracy of 0.03 mag (achieved at i~22), C is expected to be 88% for a simple sample selection criterion of tau>100 days. In summary, given an adequate survey cadence, photometric variability provides an even better method than color selection for separating quasars from stars.

We have obtained high dynamic range, good natural seeing i' images of BL Lacertae objects (BL Lacs) to search for the AGN host and thus constrain the source redshift. These objects are drawn from a sample of bright flat-spectrum radio sources that are either known (via recent Fermi LAT observations) gamma-ray emitters or similar sources that might be detected in continuing gamma-ray observations. All had spectroscopic confirmation as BL Lac sources, but no redshift solution. We detected hosts for 25/49 objects. As these galaxies have been argued to be standard candles, our measured host magnitudes provide redshift estimates (ranging from 0.2--1.0). Lower bounds are established on the redshifts of non-detections. The mean of the fit redshifts (and lower limits) is higher than those of spectroscopic solutions in the radio- and gamma-ray- loud parent samples, suggesting corrections may be needed for the luminosity function and evolution of these sources.