Accurate photometric redshift (photo-$z$) estimates are essential to the cosmological science goals of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). In this work we use simulated photometry for mock galaxy catalogs to explore how LSST photo-$z$ estimates can be improved by the addition of near-infrared (NIR) and/or ultraviolet (UV) photometry from the Euclid, WFIRST, and/or CASTOR space telescopes. Generally, we find that deeper optical photometry can reduce the standard deviation of the photo-$z$ estimates more than adding NIR or UV filters, but that additional filters are the only way to significantly lower the fraction of galaxies with catastrophically under- or over-estimated photo-$z$. For Euclid, we find that the addition of ${JH}$ $5{\sigma}$ photometric detections can reduce the standard deviation for galaxies with $z>1$ ($z>0.3$) by ${\sim}20\%$ (${\sim}10\%$), and the fraction of outliers by ${\sim}40\%$ (${\sim}25\%$). For WFIRST, we show how the addition of deep ${YJHK}$ photometry could reduce the standard deviation by ${\gtrsim}50\%$ at $z>1.5$ and drastically reduce the fraction of outliers to just ${\sim}2\%$ overall. For CASTOR, we find that the addition of its ${UV}$ and $u$-band photometry could reduce the standard deviation by ${\sim}30\%$ and the fraction of outliers by ${\sim}50\%$ for galaxies with $z<0.5$. We also evaluate the photo-$z$ results within sky areas that overlap with both the NIR and UV surveys, and when spectroscopic training sets built from the surveys' small-area deep fields are used.

Earth's temporarily-captured orbiters (TCOs) are a sub-population of near-Earth objects (NEOs). TCOs can provide constraints for NEO population models in the 1--10-metre-diameter range, and they are outstanding targets for in situ exploration of asteroids due to a low requirement on $\Delta v$. So far there has only been a single serendipitous discovery of a TCO. Here we assess in detail the possibility of their discovery with the upcoming Large Synoptic Survey Telescope (LSST), previously identified as the primary facility for such discoveries. We simulated observations of TCOs by combining a synthetic TCO population with an LSST survey simulation. We then assessed the detection rates, detection linking and orbit computation, and sources for confusion. Typical velocities of detectable TCOs will range from 1$^{\circ}$/day to 50$^{\circ}$/day, and typical apparent $V$ magnitudes from 21 to 23. Potentially-hazardous asteroids have observational characteristics similar to TCOs, but the two populations can be distinguished based on their orbits with LSST data alone. We predict that a TCO can be discovered once every year with the baseline moving-object processing system (MOPS). The rate can be increased to one TCO discovery every two months if tools complementary to the baseline MOPS are developed for the specific purpose of discovering these objects.

We present the discovery of an optical transient (OT) in Messier 51, designated M51 OT2019-1 (also ZTF19aadyppr, AT 2019abn, ATLAS19bzl), by the Zwicky Transient Facility (ZTF). The OT rose over 15 days to an observed luminosity of $M_r=-13$ (${\nu}L_{\nu}=9\times10^6~L_{\odot}$), in the luminosity gap between novae and typical supernovae (SNe). Spectra during the outburst show a red continuum, Balmer emission with a velocity width of $\approx400$ km s$^{-1}$, Ca II and [Ca II] emission, and absorption features characteristic of an F-type supergiant. The spectra and multiband light curves are similar to the so-called "SN impostors" and intermediate-luminosity red transients (ILRTs). We directly identify the likely progenitor in archival Spitzer Space Telescope imaging with a $4.5~\mu$m luminosity of $M_{[4.5]}\approx-12.2$ and a $[3.6]-[4.5]$ color redder than 0.74 mag, similar to those of the prototype ILRTs SN 2008S and NGC 300 OT2008-1. Intensive monitoring of M51 with Spitzer further reveals evidence for variability of the progenitor candidate at [4.5] in the years before the OT. The progenitor is not detected in pre-outburst Hubble Space Telescope optical and near-IR images. The optical colors during outburst combined with spectroscopic temperature constraints imply a higher reddening of $E(B-V)\approx0.7$ mag and higher intrinsic luminosity of $M_r\approx-14.9$ (${\nu}L_{\nu}=5.3\times10^7~L_{\odot}$) near peak than seen in previous ILRT candidates. Moreover, the extinction estimate is higher on the rise than on the plateau, suggestive of an extended phase of circumstellar dust destruction. These results, enabled by the early discovery of M51 OT2019-1 and extensive pre-outburst archival coverage, offer new clues about the debated origins of ILRTs and may challenge the hypothesis that they arise from the electron-capture induced collapse of extreme asymptotic giant branch stars.

We introduce AXS (Astronomy eXtensions for Spark), a scalable open-source astronomical data analysis framework built on Apache Spark, a widely used industry-standard engine for big data processing. Building on capabilities present in Spark, AXS aims to enable querying and analyzing almost arbitrarily large astronomical catalogs using familiar Python/AstroPy concepts, DataFrame APIs, and SQL statements. We achieve this by i) adding support to Spark for efficient on-line positional cross-matching and ii) supplying a Python library supporting commonly-used operations for astronomical data analysis. To support scalable cross-matching, we developed a variant of the ZONES algorithm (Gray et al. 2004) capable of operating in distributed, shared-nothing architecture. We couple this to a data partitioning scheme that enables fast catalog cross-matching and handles the data skew often present in deep all-sky data sets. The cross-match and other often-used functionalities are exposed to the end users through an easy-to-use Python API. We demonstrate AXS' technical and scientific performance on SDSS, ZTF, Gaia DR2, and AllWise catalogs. Using AXS we were able to perform on-the-fly cross-match of Gaia DR2 (1.8 billion rows) and AllWise (900 million rows) data sets in ~ 30 seconds. We discuss how cloud-ready distributed systems like AXS provide a natural way to enable comprehensive end-user analyses of large datasets such as LSST.

We present new mass estimates and cumulative mass profiles (CMPs) with Bayesian credible regions for the Milky Way (MW) Galaxy, given the kinematic data of globular clusters as provided by (1) the $\textit{Gaia}$ DR2 collaboration and the HSTPROMO team, and (2) the new catalog in Vasiliev (2019). We use globular clusters beyond 15kpc to estimate the CMP of the MW, assuming a total gravitational potential model $\Phi(r) = \Phi_{\circ}r^{-\gamma}$, which approximates an NFW-type potential at large distances when $\gamma=0.5$. We compare the resulting CMPs given data sets (1) and (2), and find the results to be nearly identical. The median estimate for the total mass is $M_{200}= 0.70 \times 10^{12} M_{\odot}$ and the $50\%$ Bayesian credible interval is $(0.62, 0.81)\times10^{12}M_{\odot}$. However, because the Vasiliev catalog contains more complete data at large $r$, the MW total mass is slightly more constrained by these data. In this work, we also supply instructions for how to create a CMP for the MW with Bayesian credible regions, given a model for $M(<r)$ and samples drawn from a posterior distribution. With the CMP, we can report median estimates and $50\%$ Bayesian credible regions for the MW mass within any distance (e.g., $M(r=25\text{kpc})= 0.26~(0.20, 0.36)\times10^{12}M_{\odot}$, $M(r=50\text{kpc})= 0.37~(0.29, 0.51) \times10^{12}M_{\odot}$, $M(r=100\text{kpc}) = 0.53~(0.41, 0.74) \times10^{12}M_{\odot}$, etc.), making it easy to compare our results directly to other studies.

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The realtime pipeline employs a novel image-differencing algorithm, optimized for the detection of point source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. The reconstructed astrometric accuracy per science image with respect to Gaia is typically 45 to 85 milliarcsec. This is the RMS per axis on the sky for sources extracted with photometric S/N >= 10. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae, rare flux transients, variable stars, eclipsing binaries, variability from active galactic nuclei, counterparts to gravitational wave sources, a more complete census of Type Ia supernovae, and Solar System objects.

A foundational goal of the Large Synoptic Survey Telescope (LSST) is to map the Solar System small body populations that provide key windows into understanding of its formation and evolution. This is especially true of the populations of the Outer Solar System -- objects at the orbit of Neptune $r > 30$AU and beyond. In this whitepaper, we propose a minimal change to the LSST cadence that can greatly enhance LSST's ability to discover faint distant Solar System objects across the entire wide-fast-deep (WFD) survey area. Specifically, we propose that the WFD cadence be constrained so as to deliver least one sequence of $\gtrsim 10$ visits per year taken in a $\sim 10$ day period in any combination of $g, r$, and $i$ bands. Combined with advanced shift-and-stack algorithms (Whidden et al. 2019) this modification would enable a nearly complete census of the outer Solar System to $\sim 25.5$ magnitude, yielding $4-8$x more KBO discoveries than with single-epoch baseline, and enabling rapid identification and follow-up of unusual distant Solar System objects in $\gtrsim 5$x greater volume of space. These increases would enhance the science cases discussed in Schwamb et al. (2018) whitepaper, including probing Neptune's past migration history as well as discovering hypothesized planet(s) beyond the orbit of Neptune (or at least placing significant constraints on their existence).

We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over 10^10 possible asteroid trajectories in stacks of the order 10-15 4K x 4K images in under a minute using a single consumer grade GPU. We apply this algorithm to data from the 2015 campaign of the High Cadence Transient Survey (HiTS) obtained with the Dark Energy Camera (DECam). We find 39 previously unknown Kuiper Belt Objects in the 150 square degrees of the survey. Comparing these asteroids to an existing model for the inclination distribution of the Kuiper Belt we demonstrate that we recover a KBO population above our detection limit consistent with previous studies. Software used in this analysis is made available as an open source package.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

In this paper we present and characterize a nearest-neighbors color-matching photometric redshift estimator that features a direct relationship between the precision and accuracy of the input magnitudes and the output photometric redshifts. This aspect makes our estimator an ideal tool for evaluating the impact of changes to LSST survey parameters that affect the measurement errors of the photometry, which is the main motivation of our work (i.e., it is not intended to provide the "best" photometric redshifts for LSST data). We show how the photometric redshifts will improve with time over the 10-year LSST survey and confirm that the nominal distribution of visits per filter provides the most accurate photo-$z$ results. The LSST survey strategy naturally produces observations over a range of airmass, which offers the opportunity of using an SED- and $z$-dependent atmospheric affect on the observed photometry as a color-independent redshift indicator. We show that measuring this airmass effect and including it as a prior has the potential to improve the photometric redshifts and can ameliorate extreme outliers, but that it will only be adequately measured for the brightest galaxies, which limits its overall impact on LSST photometric redshifts. We furthermore demonstrate how this airmass effect can induce a bias in the photo-$z$ results, and caution against survey strategies that prioritize high-airmass observations for the purpose of improving this prior. Ultimately, we intend for this work to serve as a guide for the expectations and preparations of the LSST science community with regards to the minimum quality of photo-$z$ as the survey progresses.

Using the most recent prototypes, design, and as-built system information, we test and quantify the capability of the Large Synoptic Survey Telescope (LSST) to discover Potentially Hazardous Asteroids (PHAs) and Near-Earth Objects (NEOs). We empirically estimate an expected upper limit to the false detection rate in LSST image differencing, using measurements on DECam data and prototype LSST software and find it to be about $450$~deg$^{-2}$. We show that this rate is already tractable with current prototype of the LSST Moving Object Processing System (MOPS) by processing a 30-day simulation consistent with measured false detection rates. We proceed to evaluate the performance of the LSST baseline survey strategy for PHAs and NEOs using a high-fidelity simulated survey pointing history. We find that LSST alone, using its baseline survey strategy, will detect $66\%$ of the PHA and $61\%$ of the NEO population objects brighter than $H=22$, with the uncertainty in the estimate of $\pm5$ percentage points. By generating and examining variations on the baseline survey strategy, we show it is possible to further improve the discovery yields. In particular, we find that extending the LSST survey by two additional years and doubling the MOPS search window increases the completeness for PHAs to $86\%$ (including those discovered by contemporaneous surveys) without jeopardizing other LSST science goals ($77\%$ for NEOs). This equates to reducing the undiscovered population of PHAs by additional $26\%$ ($15\%$ for NEOs), relative to the baseline survey.

The Large Synoptic Survey Telescope (LSST), the next-generation optical imaging survey sited at Cerro Pachon in Chile, will provide an unprecedented database of astronomical measurements. The LSST design, with an 8.4m (6.7m effective) primary mirror, a 9.6 sq. deg. field of view, and a 3.2 Gigapixel camera, will allow about 10,000 sq. deg. of sky to be covered twice per night, every three to four nights on average, with typical 5-sigma depth for point sources of $r$=24.5 (AB). With over 800 observations in $ugrizy$ bands over a 10-year period, these data will enable a deep stack reaching $r$=27.5 (about 5 magnitudes deeper than SDSS) and faint time-domain astronomy. The measured properties of newly discovered and known astrometric and photometric transients will be publicly reported within 60 sec after observation. The vast database of about 30 trillion observations of 40 billion objects will be mined for the unexpected and used for precision experiments in astrophysics. In addition to a brief introduction to LSST, we discuss a number of astro-statistical challenges that need to be overcome to extract maximum information and science results from LSST dataset.

The Large Synoptic Survey Telescope (LSST) is a large-aperture, wide-field, ground-based survey system that will image the sky in six optical bands from 320 to 1050 nm, uniformly covering approximately $18,000$deg$^2$ of the sky over 800 times. The LSST is currently under construction on Cerro Pach\'on in Chile, and expected to enter operations in 2022. Once operational, the LSST will explore a wide range of astrophysical questions, from discovering "killer" asteroids to examining the nature of Dark Energy. The LSST will generate on average 15 TB of data per night, and will require a comprehensive Data Management system to reduce the raw data to scientifically useful catalogs and images with minimum human intervention. These reductions will result in a real-time alert stream, and eleven data releases over the 10-year duration of LSST operations. To enable this processing, the LSST project is developing a new, general-purpose, high-performance, scalable, well documented, open source data processing software stack for O/IR surveys. Prototypes of this stack are already capable of processing data from existing cameras (e.g., SDSS, DECam, MegaCam), and form the basis of the Hyper-Suprime Cam (HSC) Survey data reduction pipeline.

We present a method to infer reddenings and distances to stars, based only on their broad-band photometry, and show how this method can be used to produce a three-dimensional dust map of the Galaxy. Our method samples from the full probability density function of distance, reddening and stellar type for individual stars, as well as the full uncertainty in reddening as a function of distance in the 3D dust map. We incorporate prior knowledge of the distribution of stars in the Galaxy and the detection limits of the survey. For stars in the Pan-STARRS 1 (PS1) 3 pi survey, we demonstrate that our reddening estimates are unbiased, and accurate to ~0.13 mag in E(B-V) for the typical star. Based on comparisons with mock catalogs, we expect distances for main-sequence stars to be constrained to within ~20% - 60%, although this range can vary, depending on the reddening of the star, the precise stellar type and its position on the sky. A further paper will present a 3D map of dust over the three quarters of the sky surveyed by PS1. Both the individual stellar inferences and the 3D dust map will enable a wealth of Galactic science in the plane. The method we present is not limited to the passbands of the PS1 survey, but may be extended to incorporate photometry from other surveys, such as 2MASS, SDSS (where available), and in the future, LSST and Gaia.

We present a sample of ~5,000 RR Lyrae stars selected from the recalibrated LINEAR dataset and detected at heliocentric distances between 5 kpc and 30 kpc over ~8,000 deg^2 of sky. The coordinates and light curve properties, such as period and Oosterhoff type, are made publicly available. We find evidence for the Oosterhoff dichotomy among field RR Lyrae stars, with the ratio of the type II and I subsamples of about 1:4. The number density distribution of halo RRab stars as a function of galactocentric distance can be described as an oblate ellipsoid with the axis ratio q=0.63 and with either a single or a double power law with a power-law index in the range -2 to -3. Using a group-finding algorithm EnLink, we detected seven candidate halo groups, only one of which is statistically spurious. Three of these groups are near globular clusters (M53/NGC 5053, M3, M13), and one is near a known halo substructure (Virgo Stellar Stream); the remaining three groups do not seem to be near any known halo substructures or globular clusters, and seem to have a higher ratio of Oosterhoff type II to Oosterhoff type I RRab stars than what is found in the halo. The extended morphology and the position (outside the tidal radius) of some of the groups near globular clusters is suggestive of tidal streams possibly originating from globular clusters. Spectroscopic followup of detected halo groups is encouraged.

We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, $R_V$, can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes. Our results are in good agreement with the extinction normalization given by the Schlegel et al. (1998, SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with Schlafly et al. (2010). The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses support the models by Fitzpatrick (1999) and Cardelli et al. (1989). For the latter, we find an $R_V=3.0\pm0.1$(random) $\pm0.1$(systematic) over most of the high-latitude sky. At low Galactic latitudes (|b|<5), we demonstrate that the SFD map cannot be reliably used to correct for extinction as most stars are embedded in dust, rather than behind it. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed "dusty parallax relation", which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format.

We present measurements of reddening due to dust using the colors of stars in the Sloan Digital Sky Survey (SDSS). We measure the color of main sequence turn-off stars by finding the "blue tip" of the stellar locus: the prominent blue edge in the distribution of stellar colors. The method is sensitive to color changes of order 18, 12, 7, and 8 mmag of reddening in the colors u-g, g-r, r-i, and i-z, respectively, in regions measuring 90' by 14'. We present maps of the blue tip colors in each of these bands over the entire SDSS footprint, including the new dusty southern Galactic cap data provided by the SDSS-III. The results disfavor the best fit O'Donnell (1994) and Cardelli et al. (1989) reddening laws, but are well described by a Fitzpatrick (1999) reddening law with R_V = 3.1. The SFD dust map is found to trace the dust well, but overestimates reddening by factors of 1.4, 1.0, 1.2, and 1.4 in u-g, g-r, r-i, and i-z, largely due to the adopted reddening law. In select dusty regions of the sky, we find evidence for problems in the SFD temperature correction. A dust map normalization difference of 15% between the Galactic north and south sky may be due to these dust temperature errors.

We have produced a new conformal map of the universe illustrating recent discoveries, ranging from Kuiper belt objects in the Solar system, to the galaxies and quasars from the Sloan Digital Sky Survey. This map projection, based on the logarithm map of the complex plane, preserves shapes locally, and yet is able to display the entire range of astronomical scales from the Earth's neighborhood to the cosmic microwave background. The conformal nature of the projection, preserving shapes locally, may be of particular use for analyzing large scale structure. Prominent in the map is a Sloan Great Wall of galaxies 1.37 billion light years long, 80% longer than the Great Wall discovered by Geller and Huchra and therefore the largest observed structure in the universe.