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.

The ESSENCE survey discovered 213 Type Ia supernovae at redshifts 0.1 < z < 0.81 between 2002 and 2008. We present their R and I-band photometry, measured from images obtained using the MOSAIC II camera at the CTIO 4 m Blanco telescope, along with rapid-response spectroscopy for each object. We use our spectroscopic follow-up observations to determine an accurate, quantitative classification and a precise redshift. Through an extensive calibration program we have improved the precision of the CTIO Blanco natural photometric system. We use several empirical metrics to measure our internal photometric consistency and our absolute calibration of the survey. We assess the effect of various potential sources of systematic bias on our measured fluxes, and we estimate that the dominant term in the systematic error budget from the photometric calibration on our absolute fluxes is ~1%.

We report the discovery and orbit determination of 14 trans-Neptunian objects (TNOs) from the ESSENCE Supernova Survey difference imaging dataset. Two additional objects discovered in a similar search of the SDSS-II Supernova Survey database were recovered in this effort. ESSENCE repeatedly observed fields far from the Solar System ecliptic (-21 deg < beta < -5 deg), reaching limiting magnitudes per observation of I~23.1 and R~23.7. We examine several of the newly detected objects in detail, including 2003 UC_414 which orbits entirely between Uranus and Neptune and lies very close to a dynamical region that would make it stable for the lifetime of the Solar System. 2003 SS_422 and 2007 TA_418 have high eccentricities and large perihelia, making them candidate members of an outer class of trans-Neptunian objects. We also report a new member of the ''extended'' or ''detached'' scattered disk, 2004 VN_112, and verify the stability of its orbit using numerical simulations. This object would have been visible to ESSENCE for only ~2% of its orbit, suggesting a vast number of similar objects across the sky. We emphasize that off-ecliptic surveys are optimal for uncovering the diversity of such objects, which in turn will constrain the history of gravitational influences that shaped our early Solar System.

We present the first large-scale effort of creating composite spectra of high-redshift type Ia supernovae (SNe Ia) and comparing them to low-redshift counterparts. Through the ESSENCE project, we have obtained 107 spectra of 88 high-redshift SNe Ia with excellent light-curve information. In addition, we have obtained 397 spectra of low-redshift SNe through a multiple-decade effort at Lick and Keck Observatories, and we have used 45 UV spectra obtained by HST/IUE. The low-redshift spectra act as a control sample when comparing to the ESSENCE spectra. In all instances, the ESSENCE and Lick composite spectra appear very similar. The addition of galaxy light to the Lick composite spectra allows a nearly perfect match of the overall spectral-energy distribution with the ESSENCE composite spectra, indicating that the high-redshift SNe are more contaminated with host-galaxy light than their low-redshift counterparts. This is caused by observing objects at all redshifts with the same slit width, which corresponds to different projected distances. After correcting for the galaxy-light contamination, subtle differences in the spectra remain. We have estimated the systematic errors when using current spectral templates for K-corrections to be ~0.02 mag. The variance in the composite spectra give an estimate of the intrinsic variance in low-redshift maximum-light SN spectra of ~3% in the optical and growing toward the UV. The difference between the maximum light low and high-redshift spectra constrain SN evolution between our samples to be < 10% in the rest-frame optical.

We present constraints on the dark energy equation-of-state parameter, w=P/(rho c^2), using 60 Type Ia supernovae (SNe Ia) from the ESSENCE supernova survey. We derive a set of constraints on the nature of the dark energy assuming a flat Universe. By including constraints on (Omega_M, w) from baryon acoustic oscillations, we obtain a value for a static equation-of-state parameter w=-1.05^{+0.13}_{-0.12} (stat; 1 sigma) +- 0.11 (sys) and Omega_M=0.274^{+0.033}_{-0.020} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.96. These results are consistent with those reported by the SuperNova Legacy Survey in a similar program measuring supernova distances and redshifts. We evaluate sources of systematic error that afflict supernova observations and present Monte Carlo simulations that explore these effects. Currently, the largest systematic currently with the potential to affect our measurements is the treatment of extinction due to dust in the supernova host galaxies. Combining our set of ESSENCE SNe Ia with the SuperNova Legacy Survey SNe Ia, we obtain a joint constraint of w=-1.07^{+0.09}_{-0.09} (stat; 1 sigma) +- 0.12 (sys), Omega_M=0.267^{+0.028}_{-0.018} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.91. The current SN Ia data are fully consistent with a cosmological constant.

We describe the implementation and optimization of the ESSENCE supernova survey, which we have undertaken to measure the equation of state parameter of the dark energy. We present a method for optimizing the survey exposure times and cadence to maximize our sensitivity to the dark energy equation of state parameter w=P/rho c^2 for a given fixed amount of telescope time. For our survey on the CTIO 4m telescope, measuring the luminosity distances and redshifts for supernovae at modest redshifts (z~0.5 +- 0.2) is optimal for determining w. We describe the data analysis pipeline based on using reliable and robust image subtraction to find supernovae automatically and in near real-time. Since making cosmological inferences with supernovae relies crucially on accurate measurement of their brightnesses, we describe our efforts to establish a thorough calibration of the CTIO 4m natural photometric system. In its first four years, ESSENCE has discovered and spectroscopically confirmed 102 type Ia SNe, at redshifts from 0.10 to 0.78, identified through an impartial, effective methodology for spectroscopic classification and redshift determination. We present the resulting light curves for the all type Ia supernovae found by ESSENCE and used in our measurement of w, presented in Wood-Vasey et al, 2007.

We present a set of 11 type Ia supernova (SN Ia) lightcurves with dense, pre-maximum sampling. These supernovae (SNe), in galaxies behind the Large Magellanic Cloud (LMC), were discovered by the SuperMACHO survey. The SNe span a redshift range of z = 0.11 - 0.35. Our lightcurves contain some of the earliest pre-maximum observations of SNe Ia to date. We also give a functional model that describes the SN Ia lightcurve shape (in our VR-band). Our function uses the "expanding fireball" model of Goldhaber et al. (1998) to describe the rising lightcurve immediately after explosion but constrains it to smoothly join the remainder of the lightcurve. We fit this model to a composite observed VR-band lightcurve of three SNe between redshifts of 0.135 to 0.165. These SNe have not been K-corrected or adjusted to account for reddening. In this redshift range, the observed VR-band most closely matches the rest frame V-band. Using the best fit to our functional description of the lightcurve, we find the time between explosion and observed VR-band maximum to be 17.6+-1.3(stat)+-0.07(sys) rest-frame days for a SN Ia with a VR-band Delta m_{-10} of 0.52mag. For the redshifts sampled, the observed VR-band time-of-maximum brightness should be the same as the rest-frame V-band maximum to within 1.1 rest-frame days.

Using archival data of low-redshift (z < 0.01) Type Ia supernovae (SN Ia) and recent observations of high-redshift (0.16 < z <0.64; Matheson et al. 2005) SN Ia, we study the "uniformity'' of the spectroscopic properties of nearby and distant SN Ia. We find no difference in the measures we describe here. In this paper, we base our analysis solely on line-profile morphology, focusing on measurements of the velocity location of maximum absorption (vabs) and peak emission (vpeak). We find that the evolution of vabs and vpeak for our sample lines (Ca II 3945, Si II 6355, and S II 5454, 5640) is similar for both the low- and high-redshift samples. We find that vabs for the weak S II 5454, 5640 lines, and vpeak for S II 5454, can be used to identify fast-declining [dm15 > 1.7] SN Ia, which are also subluminous. In addition, we give the first direct evidence in two high-z SN Ia spectra of a double-absorption feature in Ca II 3945, an event also observed, though infrequently, in low-redshift SN Ia spectra (6/22 SN Ia in our local sample). We report for the first time the unambiguous and systematic intrinsic blueshift of peak emission of optical P-Cygni line profiles in Type Ia spectra, by as much as 8000 km/s. All the high-z SN Ia analyzed in this paper were discovered and followed up by the ESSENCE collaboration, and are now publicly available.

The proper usage of Type Ia supernovae (SNe Ia) as distance indicators has revolutionized cosmology, and added a new dominant component to the energy density of the Universe, dark energy. Following the discovery and confirmation era, the currently ongoing SNe Ia surveys aim to determine the properties of the dark energy. ESSENCE is a five year ground-based supernova survey aimed at finding and characterizing 200 SNe Ia in the redshift domain z=[0.2-0.8]. The goal of the project is to put constraints on the equation of state parameter, w, of the dark energy with an accuracy of <10%. This paper presents these ongoing efforts in the context of the current developments in observational cosmology.

We present broad-band light curves of nine supernovae ranging in redshift from 0.5 to 0.8. The supernovae were discovered as part of the ESSENCE project, and the light curves are a combination of Cerro Tololo 4-m and Hubble Space Telescope (HST) photometry. On the basis of spectra and/or light-curve fitting, eight of these objects are definitely Type Ia supernovae, while the classification of one is problematic. The ESSENCE project is a five-year endeavor to discover about 200 high-redshift Type Ia supernovae, with the goal of tightly constraining the time average of the equation-of-state parameter [w = p/(rho c^2)] of the "dark energy." To help minimize our systematic errors, all of our ground-based photometry is obtained with the same telescope and instrument. In 2003 the highest-redshift subset of ESSENCE supernovae was selected for detailed study with HST. Here we present the first photometric results of the survey. We find that all but one of the ESSENCE SNe have slowly declining light curves, and the sample is not representative of the low-redshift set of ESSENCE Type Ia supernovae. This is unlikely to be a sign of evolution in the population. We attribute the decline-rate distribution of HST events to a selection bias at the high-redshift edge of our sample and find that such a bias will infect other magnitude-limited SN Ia searches unless appropriate precautions are taken.

We present the results of spectroscopic observations of targets discovered during the first two years of the ESSENCE project. The goal of ESSENCE is to use a sample of ~200 Type Ia supernovae (SNe Ia) at moderate redshifts (0.2 < z < 0.8) to place constraints on the equation of state of the Universe. Spectroscopy not only provides the redshifts of the objects, but also confirms that some of the discoveries are indeed SNe Ia. This confirmation is critical to the project, as techniques developed to determine luminosity distances to SNe Ia depend upon the knowledge that the objects at high redshift are the same as the ones at low redshift. We describe the methods of target selection and prioritization, the telescopes and detectors, and the software used to identify objects. The redshifts deduced from spectral matching of high-redshift SNe Ia with low-redshift SNe Ia are consistent with those determined from host-galaxy spectra. We show that the high-redshift SNe Ia match well with low-redshift templates. We include all spectra obtained by the ESSENCE project, including 52 SNe Ia, 5 core-collapse SNe, 12 active galactic nuclei, 19 galaxies, 4 possibly variable stars, and 16 objects with uncertain identifications.

We present the first results from our next-generation microlensing survey, the SuperMACHO project. We are using the CTIO 4m Blanco telescope and the MOSAIC imager to carry out a search for microlensing toward the Large Magellanic Cloud (LMC). We plan to ascertain the nature of the population responsible for the excess microlensing rate seen by the MACHO project. Our observing strategy is optimized to measure the differential microlensing rate across the face of the LMC. We find this derivative to be relatively insensitive to the details of the LMC's internal structure but a strong discriminant between Galactic halo and LMC self lensing. In December 2003 we completed our third year of survey operations. 2003 also marked the first year of real-time microlensing alerts and photometric and spectroscopic followup. We have extracted several dozen microlensing candidates, and we present some preliminary light curves and related information. Similar to the MACHO project, we find SNe behind the LMC to be a significant contaminant - this background has not been completely removed from our current single-color candidate sample. Our follow-up strategy is optimized to discriminate between SNe and true microlensing.