Montage: a grid portal and software toolkit for science-grade astronomical image mosaicking

Montage is a portable software toolkit for constructing custom, science-grade mosaics by composing multiple astronomical images. The mosaics constructed by Montage preserve the astrometry (position) and photometry (intensity) of the sources in the input images. The mosaic to be constructed is specified by the user in terms of a set of parameters, including dataset and wavelength to be used, location and size on the sky, coordinate system and projection, and spatial sampling rate. Many astronomical datasets are massive, and are stored in distributed archives that are, in most cases, remote with respect to the available computational resources. Montage can be run on both single- and multi-processor computers, including clusters and grids. Standard grid tools are used to run Montage in the case where the data or computers used to construct a mosaic are located remotely on the Internet. This paper describes the architecture, algorithms, and usage of Montage as both a software toolkit and as a grid portal. Timing results are provided to show how Montage performance scales with number of processors on a cluster computer. In addition, we compare the performance of two methods of running Montage in parallel on a grid.Comment: 16 pages, 11 figure

Grids and the Virtual Observatory

We consider several projects from astronomy that benefit from the Grid paradigm and associated technology, many of which involve either massive datasets or the federation of multiple datasets. We cover image computation (mosaicking, multi-wavelength images, and synoptic surveys); database computation (representation through XML, data mining, and visualization); and semantic interoperability (publishing, ontologies, directories, and service descriptions)

The ATLAS-SPT Radio Survey of Cluster Galaxies

Using a high-performance computing cluster to mosaic 4,787 pointings, we have imaged the 100 sq. deg. South Pole Telescope (SPT) deep-field at 2.1 GHz using the Australian Telescope Compact Array to an rms of 80 $\mu$Jy and a resolution of 8". Our goal is to generate an independent sample of radio-selected galaxy clusters to study how the radio properties compare with cluster properties at other wavelengths, over a wide range of redshifts in order to construct a timeline of their evolution out to $z \sim 1.3$. A preliminary analysis of the source catalogue suggests there is no spatial correlation between the clusters identified in the SPT-SZ catalogue and our wide-angle tail galaxies.Comment: 9 pages, 4 figures. Submitted to Proceedings of Science for "The many facets of extragalactic radio surveys: towards new scientific challenges", Bologna, Italy 20-23 October 2015 (EXTRA-RADSUR2015

Surveying the Agents of Galaxy Evolution in the Tidally Stripped, Low Metallicity Small Magellanic Cloud (SAGE-SMC). I. Overview

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 deg^2) including the body, wing, and tail in seven bands from 3.6 to 160 μm using IRAC and MIPS on the Spitzer Space Telescope. The data were reduced and mosaicked, and the point sources were measured using customized routines specific for large surveys. We have made the resulting mosaics and point-source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 μm/24 μm ratio is somewhat lower than the average and the 70 μm/160 μm ratio is somewhat higher than the average. The global infrared spectral energy distribution (SED) shows that the SMC has approximately 1/3 the aromatic emission/polycyclic aromatic hydrocarbon abundance of most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average SED of H II/star formation regions is compared to the equivalent Large Magellanic Cloud average H II/star formation region SED. These preliminary results will be expanded in detail in subsequent papers

PINGS: the PPAK IFS Nearby Galaxies Survey

We present the PPAK Integral Field Spectroscopy (IFS) Nearby Galaxies Survey: PINGS, a 2-dimensional spectroscopic mosaicking of 17 nearby disk galaxies in the optical wavelength range. This project represents the first attempt to obtain continuous coverage spectra of the whole surface of a galaxy in the nearby universe. The final data set comprises more than 50000 individual spectra, covering in total an observed area of nearly 80 arcmin^2. In this paper we describe the main astrophysical issues to be addressed by the PINGS project, we present the galaxy sample and explain the observing strategy, the data reduction process and all uncertainties involved. Additionally, we give some scientific highlights extracted from the first analysis of the PINGS sample.Comment: Accepted for publication in MNRAS, 26 pages, 14 figures (some in low resolution), 3 table

Methodology for the metric restoration of the historical cartography applied to Francisco Coello's cartografic series of the Royal Site of Aranjuez

Hojas Kilométricas (Kilometric Sheets). Specifically, the study focuses on those sheets referring to the city centre and surrounding area of the Royal Site of Aranjuez, a town in the south of the Province of Madrid. The aim of this study is to restore the actual size and measurements of scanned images of the Hojas Kilométricas. This would allow us, among other things, to reestablish both the format and scale of the original plans. To achieve this goal it is necessary to rectify and then georeference these images, i.e. assign them a geographic reference system. This procedure is essential in the overlaying and comparison of the Hojas Kilométricas of the Royal Site with other historical cartography as well as other sources related to the same area from different time periods. Subsequent research would allow us, for example, to reconstruct the time-evolution of the urban area, to spot new construction and to pinpoint the locations of any altered or missing buildings or architectural features. In addition, this would allow us to develop and integrate databases for GIS models applicable to the management of our cultural heritage

Andromeda Optical and Infrared Disk Survey I. New Insights in Wide-Field Near-IR Surface Photometry

We present wide-field near-infrared J and Ks images of the Andromeda Galaxy taken with WIRCam on the Canada-France-Hawaii Telescope (CFHT) as part of the Andromeda Optical and Infrared Disk Survey (ANDROIDS). This data set allows simultaneous observations of resolved stars and NIR surface brightness across M31's entire bulge and disk (within R=22 kpc). The primary concern of this work is the development of NIR observation and reduction methods to recover a uniform surface brightness map across the 3x1 degree disk of M31. This necessitates sky-target nodding across 27 WIRCam fields. Two sky-target nodding strategies were tested, and we find that strictly minimizing sky sampling latency does not maximize sky subtraction accuracy, which is at best 2% of the sky level. The mean surface brightness difference between blocks in our mosaic can be reduced from 1% to 0.1% of the sky brightness by introducing scalar sky offsets to each image. The true surface brightness of M31 can be known to within a statistical zeropoint of 0.15% of the sky level (0.2 mag arcsec sq. uncertainty at R=15 kpc). Surface brightness stability across individual WIRCam frames is limited by both WIRCam flat field evolution and residual sky background shapes. To overcome flat field variability of order 1% over 30 minutes, we find that WIRCam data should be calibrated with real-time sky flats. Due either to atmospheric or instrumental variations, the individual WIRCam frames have typical residual shapes with amplitudes of 0.2% of the sky after real-time flat fielding and median sky subtraction. We present our WIRCam reduction pipeline and performance analysis here as a template for future near-infrared observers needing wide-area surface brightness maps with sky-target nodding, and give specific recommendations for improving photometry of all CFHT/WIRCam programs. (Abridged)Comment: 21 pages, 28 figures. Submitted to A

A Spitzer/IRAC Search for Substellar Companions of the Debris Disk Star epsilon Eridani

We have used the InfraRed Array Camera (IRAC) onboard the Spitzer Space telescope to search for low mass companions of the nearby debris disk star epsilon Eridani. The star was observed in two epochs 39 days apart, with different focal plane rotation to allow the subtraction of the instrumental Point Spread Function, achieving a maximum sensitivity of 0.01 MJy/sr at 3.6 and 4.5 um, and 0.05 MJy/sr at 5.8 and 8.0 um. This sensitivity is not sufficient to directly detect scattered or thermal radiation from the epsilon Eridani debris disk. It is however sufficient to allow the detection of Jovian planets with mass as low as 1 MJ in the IRAC 4.5 um band. In this band, we detected over 460 sources within the 5.70 arcmin field of view of our images. To test if any of these sources could be a low mass companion to epsilon Eridani, we have compared their colors and magnitudes with models and photometry of low mass objects. Of the sources detected in at least two IRAC bands, none fall into the range of mid-IR color and luminosity expected for cool, 1 Gyr substellar and planetary mass companions of epsilon Eridani, as determined by both models and observations of field M, L and T dwarf. We identify three new sources which have detections at 4.5 um only, the lower limit placed on their [3.6]-[4.5] color consistent with models of planetary mass objects. Their nature cannot be established with the currently available data and a new observation at a later epoch will be needed to measure their proper motion, in order to determine if they are physically associated to epsilon Eridani.Comment: 36 pages, to be published on The Astrophysical Journal, vol. 647, August 200