The International Virtual Observatory Alliance (IVOA) held its bi-annual Interoperability Meetings in May 2019, and in October 2019 following the ADASS XXIX conference. We provide a brief report on the status of the IVOA and the activities of the Interoperability Meetings.

Many astronomy data centres still work on filesystems. Industry has moved on; current practice in computing infrastructure is to achieve Big Data scalability using object stores rather than POSIX file systems. This presents us with opportunities for portability and reuse of software underlying processing and archive systems but it also causes problems for legacy implementations in current data centers.

The International Virtual Observatory Alliance (IVOA) held its bi-annual Interoperability Meeting over two and half days prior to the ADASS 2018 conference. We provide a brief report on the status of the IVOA and the activities of the Interoperability Meeting held in College Park.

Writing generic data processing pipelines requires that the algorithmic code does not ever have to know about data formats of files, or the locations of those files. At LSST we have a software system known as "the Data Butler," that abstracts these details from the software developer. Scientists can specify the dataset they want in terms they understand, such as filter, observation identifier, date of observation, and instrument name, and the Butler translates that to one or more files which are read and returned to them as a single Python object. Conversely, once they have created a new dataset they can give it back to the Butler, with a label describing its new status, and the Butler can write it in whatever format it has been configured to use. All configuration is in YAML and supports standard defaults whilst allowing overrides.

The Large Synoptic Survey Telescope (LSST) is an ambitious astronomical survey with a similarly ambitious Data Management component. Data Management for LSST includes processing on both nightly and yearly cadences to generate transient alerts, deep catalogs of the static sky, and forced photometry light-curves for billions of objects at hundreds of epochs, spanning at least a decade. The algorithms running in these pipelines are individually sophisticated and interact in subtle ways. This paper provides an overview of those pipelines, focusing more on those interactions than the details of any individual algorithm.

The JCMT Gould Belt Survey was one of the first Legacy Surveys with the James Clerk Maxwell Telescope in Hawaii, mapping 47 square degrees of nearby (< 500 pc) molecular clouds in both dust continuum emission at 850 $\mu$m and 450 $\mu$m, as well as a more-limited area in lines of various CO isotopologues. While molecular clouds and the material that forms stars have structures on many size scales, their larger-scale structures are difficult to observe reliably in the submillimetre regime using ground-based facilities. In this paper, we quantify the extent to which three subsequent data-reduction methods employed by the JCMT GBS accurately recover emission structures of various size scales, in particular, dense cores which are the focus of many GBS science goals. With our current best data-reduction procedure, we expect to recover 100% of structures with Gaussian sigma sizes of $\le$30" and intensity peaks of at least five times the local noise for isolated peaks of emission. The measured sizes and peak fluxes of these compact structures are reliable (within 15% of the input values), but source recovery and reliability both decrease significantly for larger emission structures and for fainter peaks. Additional factors such as source crowding have not been tested in our analysis. The most recent JCMT GBS data release includes pointing corrections, and we demonstrate that these tend to decrease the sizes and increase the peak intensities of compact sources in our dataset, mostly at a low level (several percent), but occasionally with notable improvement.

(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.

Since its inception in the early 2000, the Virtual Observatory (VO), developed as a collaboration of many national and international projects, has become a major factor in the discovery and dissemination of astronomical information worldwide. The IVOA has been coordinating all these efforts worldwide to ensure a common VO framework that enables transparent access to and interoperability of astronomy resources (data and software) around the world. The VO is not a magic solution to all astronomy data management challenges but it does bring useful solutions in many areas borne out by the fact that VO interfaces are broadly found in astronomy major data centres and projects worldwide. Astronomy data centres have been building VO services on top of their existing data services to increase interoperability with other VO-compliant data resources to take advantage of the continuous and increasing development of VO applications. VO applications have made multi-instrument and multi-wavelength science, a difficult and fruitful part of astronomy, somewhat easier. More recently, several major new astronomy projects have been directly adopting VO standards to build their data management infrastructure, giving birth to VO built-in archives. Embracing the VO framework from the beginning brings the double gain of not needing to reinvent the wheel and ensuring from the start interoperability with other astronomy VO resources. Some of the IVOA standards are also starting to be used by neighbour disciplines like planetary sciences. There is still quite a lot to be done on the VO, in particular tackling the upcoming big data challenge and how to find interoperable solutions to the new data analysis paradigm of bringing and running the software close to the data.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

The unprecedented volume and rate of transient events that will be discovered by the Large Synoptic Survey Telescope (LSST) demands that the astronomical community update its followup paradigm. Alert-brokers -- automated software system to sift through, characterize, annotate and prioritize events for followup -- will be critical tools for managing alert streams in the LSST era. The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is one such broker. In this work, we develop a machine learning pipeline to characterize and classify variable and transient sources only using the available multiband optical photometry. We describe three illustrative stages of the pipeline, serving the three goals of early, intermediate and retrospective classification of alerts. The first takes the form of variable vs transient categorization, the second, a multi-class typing of the combined variable and transient dataset, and the third, a purity-driven subtyping of a transient class. While several similar algorithms have proven themselves in simulations, we validate their performance on real observations for the first time. We quantitatively evaluate our pipeline on sparse, unevenly sampled, heteroskedastic data from various existing observational campaigns, and demonstrate very competitive classification performance. We describe our progress towards adapting the pipeline developed in this work into a real-time broker working on live alert streams from time-domain surveys.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

The LSST software systems make extensive use of Python, with almost all of it initially being developed solely in Python 2. Since LSST will be commissioned when Python 2 is end-of-lifed it is critical that we have all our code support Python 3 before commissioning begins. Over the past year we have made significant progress in migrating the bulk of the code from the Data Management system onto Python 3. This paper presents our migration methodology, and the current status of the port, with our eventual aim to be running completely on Python 3 by early 2018. We also discuss recent modernizations to our Python codebase.

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

We present the first data release of the James Clerk Maxwell Telescope (JCMT) Plane Survey (JPS), the JPS Public Release 1 (JPSPR1). JPS is an 850-um continuum survey of six fields in the northern inner Galactic Plane in a longitude range of l=7-63, made with the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2). This first data release consists of emission maps of the six JPS regions with an average pixel-to-pixel noise of 7.19 mJy beam^-1, when smoothed over the beam, and a compact-source catalogue containing 7,813 sources. The 95 per cent completeness limits of the catalogue are estimated at 0.04 Jy beam^-1 and 0.3 Jy for the peak and integrated flux densities, respectively. The emission contained in the compact-source catalogue is 42 +- 5 per cent of the total and, apart from the large-scale (greater than 8') emission, there is excellent correspondence with features in the 500-um Herschel maps. We find that, with two-dimensional matching, 98 +- 2 per cent of sources within the fields centred at l=20, 30, 40 and 50 are associated with molecular clouds, with 91 +- 3 per cent of the l=30 and 40 sources associated with dense molecular clumps. Matching the JPS catalogue to Herschel 70-um sources, we find that 38 +- 1 per cent of sources show evidence of ongoing star formation. The images and catalogue will be a valuable resource for studies of star formation in the Galaxy and the role of environment and spiral arms in the star formation process.

The SCUBA-2 Ambitious Sky Survey (SASSy) is composed of shallow 850-$\umu$m imaging using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) on the James Clerk Maxwell Telescope. Here we describe the extraction of a catalogue of beam-sized sources from a roughly $120\,{\rm deg}^2$ region of the Galactic plane mapped uniformly (to an rms level of about 40\,mJy), covering longitude 120\degr\,$<$\,\textit{l}\,$<$\,140\degr\ and latitude $\abs{\textit{b}}$\,$<$\,2.9\degr. We used a matched-filtering approach to increase the signal-to-noise (S/N) ratio in these noisy maps and tested the efficiency of our extraction procedure through estimates of the false discovery rate, as well as by adding artificial sources to the real images. The primary catalogue contains a total of 189 sources at 850\,$\umu$m, down to a S/N threshold of approximately 4.6. Additionally, we list 136 sources detected down to ${\rm S/N}=4.3$, but recognise that as we go lower in S/N, the reliability of the catalogue rapidly diminishes. We perform follow-up observations of some of our lower significance sources through small targeted SCUBA-2 images, and list 265 sources detected in these maps down to ${\rm S/N}=5$. This illustrates the real power of SASSy: inspecting the shallow maps for regions of 850-$\umu$m emission and then using deeper targeted images to efficiently find fainter sources. We also perform a comparison of the SASSy sources with the Planck Catalogue of Compact Sources and the \textit{IRAS} Point Source Catalogue, to determine which sources discovered in this field might be new, and hence potentially cold regions at an early stage of star formation.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint effort of NOAO and the Department of Computer Science at the University of Arizona to build prototype software to process alerts from time-domain surveys, especially LSST, to identify those alerts that must be followed up immediately. Value is added by annotating incoming alerts with existing information from previous surveys and compilations across the electromagnetic spectrum and from the history of past alerts. Comparison against a knowledge repository of properties and features of known or predicted kinds of variable phenomena is used for categorization. The architecture and algorithms being employed are described.

The LSST data management science pipelines software consists of more than 100,000 lines of Python 2 code. LSST operations will begin after support for Python 2 has been dropped by the Python community in 2020, and we must therefore plan to migrate the codebase to Python 3. During the transition period we must also support our community of active Python 2 users and this complicates the porting significantly. We have decided to use the Python future package as the basis for our port to enable support for Python 2 and Python 3 simultaneously, whilst developing with a mindset more suited to Python 3. In this paper we report on the current status of the port and the difficulties that have been encountered.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $\sigma_{450}=1.2$ mJy/beam and $\sigma_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

The Orion A molecular cloud is one of the most well-studied nearby star-forming regions, and includes regions of both highly clustered and more dispersed star formation across its full extent. Here, we analyze dense, star-forming cores identified in the 850 {\mu}m and 450 {\mu}m SCUBA-2 maps from the JCMT Gould Belt Legacy Survey. We identify dense cores in a uniform manner across the Orion A cloud and analyze their clustering properties. Using two independent lines of analysis, we find evidence that clusters of dense cores tend to be mass segregated, suggesting that stellar clusters may have some amount of primordial mass segregation already imprinted in them at an early stage. We also demonstrate that the dense core clusters have a tendency to be elongated, perhaps indicating a formation mechanism linked to the filamentary structure within molecular clouds.

We present a comparison of SCUBA-2 850-$\mu$m and Herschel 70--500-$\mu$m observations of the L1495 filament in the Taurus Molecular Cloud with the goal of characterising the SCUBA-2 Gould Belt Survey (GBS) data set. We identify and characterise starless cores in three data sets: SCUBA-2 850-$\mu$m, Herschel 250-$\mu$m, and Herschel 250-$\mu$m spatially filtered to mimic the SCUBA-2 data. SCUBA-2 detects only the highest-surface-brightness sources, principally detecting protostellar sources and starless cores embedded in filaments, while Herschel is sensitive to most of the cloud structure, including extended low-surface-brightness emission. Herschel detects considerably more sources than SCUBA-2 even after spatial filtering. We investigate which properties of a starless core detected by Herschel determine its detectability by SCUBA-2, and find that they are the core's temperature and column density (for given dust properties). For similar-temperature cores, such as those seen in L1495, the surface brightnesses of the cores are determined by their column densities, with the highest-column-density cores being detected by SCUBA-2. For roughly spherical geometries, column density corresponds to volume density, and so SCUBA-2 selects the densest cores from a population at a given temperature. This selection effect, which we quantify as a function of distance, makes SCUBA-2 ideal for identifying those cores in Herschel catalogues that are closest to forming stars. Our results can now be used by anyone wishing to use the SCUBA-2 GBS data set.

We investigate the properties of the galaxies selected from the deepest 850-micron survey undertaken to date with SCUBA-2 on the JCMT. This deep 850-micron imaging was taken in parallel with deep 450-micron imaging in the very best observing conditions as part of the SCUBA-2 Cosmology Legacy Survey. A total of 106 sources were uncovered at 850 microns from ~150, sq. arcmin in the centre of the COSMOS/UltraVISTA/CANDELS field, imaged to a typical rms depth of ~0.25 mJy. We utilise the wealth of available deep multi-frequency data to establish the complete redshift distribution for this sample, yielding <z> = 2.38 +- 0.09, a mean redshift comparable with that derived for all but the brightest previous sub-mm samples. We have also been able to establish the stellar masses of the majority of the galaxy identifications, enabling us to explore their location on the star-formation-rate:stellar-mass (SFR:M*) plane. Crucially, our new deep sample reaches flux densities equivalent to SFR ~ 100 Msun/yr, enabling us to confirm that sub-mm galaxies form the high-mass end of the `main sequence' (MS) of star-forming galaxies at z > 1.5 (with a mean specific SFR of sSFR = 2.25 +- 0.19 /Gyr at z ~ 2.5). Our results are consistent with no significant flattening of the MS towards high masses at these redshifts, suggesting that reports of such flattening possibly arise from under-estimates of dust-enshrouded star-formation activity in massive star-forming galaxies. However, our findings add to the growing evidence that average sSFR rises only slowly at high redshift, resulting in log(sSFR) being an apparently simple linear function of the age of the Universe.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

Ly-alpha blobs (LABs) offer insight into the complex interface between galaxies and their circumgalactic medium. Whilst some LABs have been found to contain luminous star-forming galaxies and active galactic nuclei that could potentially power the Ly-alpha emission, others appear not to be associated with obvious luminous galaxy counterparts. It has been speculated that LABs may be powered by cold gas streaming on to a central galaxy, providing an opportunity to directly observe the `cold accretion' mode of galaxy growth. Star-forming galaxies in LABs could be dust obscured and therefore detectable only at longer wavelengths. We stack deep SCUBA-2 observations of the SSA22 field to determine the average 850um flux density of 34 LABs. We measure S_850 = 0.6 +/- 0.2mJy for all LABs, but stacking the LABs by size indicates that only the largest third (area > 1794 kpc^2) have a mean detection, at 4.5 sigma, with S_850 = 1.4 +/- 0.3mJy. Only two LABs (1 and 18) have individual SCUBA-2 > 3.5 sigma detections at a depth of 1.1mJy/beam. We consider two possible mechanisms for powering the LABs and find that central star formation is likely to dominate the emission of Ly-alpha, with cold accretion playing a secondary role.

In view of increased interest in object-oriented systems for describing coordinate information, we present a description of the data model used by the Starlink AST library. AST provides a comprehensive range of facilities for attaching world co-ordinate systems to astronomical data, and for retrieving and interpreting that information in a variety of formats, including FITS-WCS. AST is a mature system that has been in use for more than 17 years, and may consequently be useful as a means of informing development of similar systems in the future.

The JCMT Gould Belt Legacy Survey obtained SCUBA-2 observations of dense cores within three sub-regions of Orion B: LDN 1622, NGC 2023/2024, and NGC 2068/2071, all of which contain clusters of cores. We present an analysis of the clustering properties of these cores, including the two-point correlation function and Cartwright's Q parameter. We identify individual clusters of dense cores across all three regions using a minimal spanning tree technique, and find that in each cluster, the most massive cores tend to be centrally located. We also apply the independent M-Sigma technique and find a strong correlation between core mass and the local surface density of cores. These two lines of evidence jointly suggest that some amount of mass segregation in clusters has happened already at the dense core stage.

We present the $^{13}$CO/C$^{18}$O (J=3-2) Heterodyne Inner Milky Way Plane Survey (CHIMPS) which has been carried out using the Heterodyne Array Receiver Program on the 15 m James Clerk Maxwell Telescope (JCMT) in Hawaii. The high-resolution spectral survey currently covers |b| < 0.5 deg and 28 < l < 46 deg, with an angular resolution of 15 arcsec in 0.5 km/s velocity channels. The spectra have a median rms of $\sim$ 0.6 K at this resolution, and for optically thin gas at an excitation temperature of 10 K, this sensitivity corresponds to column densities of $N_{\mathrm{H}_{2}} \sim 3 \times 10^{20}\,$cm$^{-2}$ and $N_{\mathrm{H}_{2}} \sim 4 \times 10^{21}\,$cm$^{-2}$ for $^{13}$CO and C$^{18}$O, respectively. The molecular gas that CHIMPS traces is at higher column densities and is also more optically thin than in other publicly available CO surveys due to its rarer isotopologues, and thus more representative of the three-dimensional structure of the clouds. The critical density of the J=3-2 transition of CO is $\gtrsim 10^{4}$ cm$^{-3}$ at temperatures of $\leq 20$ K, and so the higher density gas associated with star formation is well traced. These data complement other existing Galactic plane surveys, especially the JCMT Galactic Plane Survey which has similar spatial resolution and column density sensitivity, and the Herschel infrared Galactic Plane Survey. In this paper, we discuss the observations, data reduction and characteristics of the survey, presenting integrated emission maps for the region covered. Position-velocity diagrams allow comparison with Galactic structure models of the Milky Way, and while we find good agreement with a particular four arm model, there are some significant deviations.

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 first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

The Large Synoptic Survey Telescope (LSST) is an 8-m optical ground-based telescope being constructed on Cerro Pachon in Chile. LSST will survey half the sky every few nights in six optical bands. The data will be transferred to NCSA and within 60 seconds they will be reduced using difference imaging techniques and detected transients will be announced to the community in the VOEvent format. Annual data releases will be made from all the data during the 10-year mission, with unprecedented depth of coadds and time resolution of catalogs for such a large region of sky. In this paper we present the current status of the data processing software, and describe how to obtain it.

Performing ground-based submillimetre observations is a difficult task as the measurements are subject to absorption and emission from water vapour in the Earth's atmosphere and time variation in weather and instrument stability. Removing these features and other artifacts from the data is a vital process which affects the characteristics of the recovered astronomical structure we seek to study. In this paper, we explore two data reduction methods for data taken with the Submillimetre Common-User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). The JCMT Legacy Reduction 1 (JCMT LR1) and The Gould Belt Legacy Survey Legacy Release 1 (GBS LR1) reduction both use the same software, Starlink, but differ in their choice of data reduction parameters. We find that the JCMT LR1 reduction is suitable for determining whether or not compact emission is present in a given region and the GBS LR1 reduction is tuned in a robust way to uncover more extended emission, which better serves more in-depth physical analyses of star-forming regions. Using the GBS LR1 method, we find that compact sources are recovered well, even at a peak brightness of only 3 times the noise, whereas the reconstruction of larger objects requires much care when drawing boundaries around the expected astronomical signal in the data reduction process. Incorrect boundaries can lead to false structure identification or it can cause structure to be missed. In the JCMT LR1 reduction, the extent of the true structure of objects larger than a point source is never fully recovered.

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

With the advent of modern multi-detector heterodyne instruments that can result in observations generating thousands of spectra per minute it is no longer feasible to reduce these data as individual spectra. We describe the automated data reduction procedure used to generate baselined data cubes from heterodyne data obtained at the James Clerk Maxwell Telescope. The system can automatically detect baseline regions in spectra and automatically determine regridding parameters, all without input from a user. Additionally it can detect and remove spectra suffering from transient interference effects or anomalous baselines. The pipeline is written as a set of recipes using the ORAC-DR pipeline environment with the algorithmic code using Starlink software packages and infrastructure. The algorithms presented here can be applied to other heterodyne array instruments and have been applied to data from historical JCMT heterodyne instrumentation.

The General Single-Dish Data format (GSDD) was developed in the mid-1980s as a data model to support centimeter, millimeter and submillimeter instrumentation at NRAO, JCMT, the University of Arizona and IRAM. We provide an overview of the GSDD requirements and associated data model, discuss the implementation of the resultant file formats, describe its usage in the observatories and provide a retrospective on the format.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard is, however, showing its age. Developed in the late 1970s the FITS authors made a number of implementation choices for the format that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not appreciate the challenges which we would be facing today in astronomical computing. Difficulties we now face include, but are not limited to, having to address the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets and the need to capture significantly more complex metadata and data relationships. There are members of the community today who find some (or all) of these limitations unworkable, and have decided to move ahead with storing data in other formats. This reaction should be taken as a wakeup call to the FITS community to make changes in the FITS standard, or to see its usage fall. In this paper we detail some selected important problems which exist within the FITS standard today. It is not our intention to prescribe specific remedies to these issues; rather, we hope to call attention of the FITS and greater astronomical computing communities to these issues in the hopes that it will spur action to address them.

The Starlink Hierarchical Data System has been a very successful niche astronomy file format and library for over 30 years. Development of the library was frozen ten years ago when funding for Starlink was stopped and almost no-one remains who understands the implementation details. To ensure the long-term sustainability of the Starlink application software and to make the extensible N-Dimensional Data Format accessible to a broader range of users, we propose to re-implement the HDS library application interface as a layer on top of the Hierarchical Data Format version 5. We present an overview of the new implementation of version 5 of the HDS file format and describe differences between the expectations of the HDS and HDF5 library interfaces. We finish by comparing the old and new HDS implementations by looking at a comparison of file sizes and by comparing performance benchmarks.

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard, however, is showing its age. Developed in the late 1970s, the FITS authors made a number of implementation choices that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not anticipate the challenges which we are facing today in astronomical computing. Difficulties we now face include, but are not limited to, addressing the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets, and capturing significantly more complex metadata and data relationships. There are members of the community today who find some or all of these limitations unworkable, and have decided to move ahead with storing data in other formats. If this fragmentation continues, we risk abandoning the advantages of broad interoperability, and ready archivability, that the FITS format provides for astronomy. In this paper we detail some selected important problems which exist within the FITS standard today. These problems may provide insight into deeper underlying issues which reside in the format and we provide a discussion of some lessons learned. It is not our intention here to prescribe specific remedies to these issues; rather, it is to call attention of the FITS and greater astronomical computing communities to these problems in the hope that it will spur action to address them.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

We present, as a case study, a description of the partnership between an observatory (JCMT) and a data centre (CADC) that led to the development of the JCMT Science Archive (JSA). The JSA is a successful example of a service designed to use Virtual Observatory (VO) technologies from the start. We describe the motivation, process and lessons learned from this approach.

The future of astronomy is inextricably entwined with the care and feeding of astronomical data products. Community standards such as FITS and NDF have been instrumental in the success of numerous astronomy projects. Their very success challenges us to entertain pragmatic strategies to adapt and evolve the standards to meet the aggressive data-handling requirements of facilities now being designed and built. We discuss characteristics that have made standards successful in the past, as well as desirable features for the future, and an open discussion follows.

CCAT will be a large sub-millimeter telescope to be built near the ALMA site in northern Chile. The telescope must support a varied set of instrumentation including large format KID cameras, a large heterodyne array and a KID-based direct detection multi-object spectrometer. We are developing a new raw data model based on HDF5 that can cope with the expected data rates of order Gbit/s and is flexible enough to hold raw data from all planned instruments.

CCAT is a large submillimetre telescope to be built near the ALMA site in northern Chile. A large-format KID camera, with up to 48,000 detectors at a single waveband sampled at about 1 kHz, will have a data rate about 50 times larger than SCUBA-2, the largest existing submillimetre camera. Creating a map from this volume of data will be a challenge, both in terms of memory and processing time required. We investigate how to extend SMURF, the iterative map-maker used for reducing SCUBA-2 observations, to a distributed-node parallel system, and estimate how the processing time scales with the number of nodes in the system.

The extensible N-Dimensional Data Format (NDF) was designed and developed in the late 1980s to provide a data model suitable for use in a variety of astronomy data processing applications supported by the UK Starlink Project. Starlink applications were used extensively, primarily in the UK astronomical community, and form the basis of a number of advanced data reduction pipelines today. This paper provides an overview of the historical drivers for the development of NDF and the lessons learned from using a defined hierarchical data model for many years in data reduction software, data pipelines and in data acquisition systems.

ORAC-DR is a general purpose data reduction pipeline system designed to be instrument and observatory agnostic. The pipeline works with instruments as varied as infrared integral field units, imaging arrays and spectrographs, and sub-millimeter heterodyne arrays & continuum cameras. This paper describes the architecture of the pipeline system and the implementation of the core infrastructure. We finish by discussing the lessons learned since the initial deployment of the pipeline system in the late 1990s.

We describe a search for submillimeter emission in the vicinity of one of the most distant, luminous galaxies known, HerMES FLS3 at z=6.34, exploiting it as a signpost to a potentially biased region of the early Universe, as might be expected in hierarchical structure formation models. Imaging to the confusion limit with the innovative, wide-field submillimeter bolometer camera, SCUBA-2, we are sensitive to colder and/or less luminous galaxies in the surroundings of HFLS3. We use the Millennium Simulation to illustrate that HFLS3 may be expected to have companions if it is as massive as claimed, but find no significant evidence from the surface density of SCUBA-2 galaxies in its vicinity, or their colors, that HFLS3 marks an over-density of dusty, star-forming galaxies. We cannot rule out the presence of dusty neighbours with confidence, but deeper 450-um imaging has the potential to more tightly constrain the redshifts of nearby galaxies, at least one of which likely lies at z>~5. If associations with HFLS3 can be ruled out, this could be taken as evidence that HFLS3 is less biased than a simple extrapolation of the Millennium Simulation may imply. This could suggest either that it represents a rare short-lived, but highly luminous, phase in the evolution of an otherwise typical galaxy, or that this system has suffered amplification due to a foreground gravitational lens and so is not as intrinsically luminous as claimed.

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint project of the National Optical Astronomy Observatory and the Department of Computer Science at the University of Arizona. The goal is to build the software infrastructure necessary to process and filter alerts produced by time-domain surveys, with the ultimate source of such alerts being the Large Synoptic Survey Telescope (LSST). The ANTARES broker will add value to alerts by annotating them with information from external sources such as previous surveys from across the electromagnetic spectrum. In addition, the temporal history of annotated alerts will provide further annotation for analysis. These alerts will go through a cascade of filters to select interesting candidates. For the prototype, `interesting' is defined as the rarest or most unusual alert, but future systems will accommodate multiple filtering goals. The system is designed to be flexible, allowing users to access the stream at multiple points throughout the process, and to insert custom filters where necessary. We describe the basic architecture of ANTARES and the principles that will guide development and implementation.

We present statistically significant detections at 850um of the Lyman Break Galaxy (LBG) population at z=3, 4, and 5 using data from the Submillimetre Common User Bolometer Array 2 (SCUBA-2) Cosmology Legacy Survey (S2CLS) in the United Kingdom Infrared Deep Sky Survey Ultra Deep Survey (UKIDSS-UDS) field. We employ a stacking technique to probe beneath the survey limit to measure the average 850um flux density of LBGs at z=3, 4, and 5 with typical ultraviolet luminosities of L(1700A)~10^29 erg/s/Hz. We measure 850um flux densities of (0.25 +/- 0.03, (0.41 +/- 0.06), and (0.88 +/- 0.23) mJy respectively, and find that they contribute at most 20 per cent to the cosmic far-infrared background at 850um. Fitting an appropriate range of spectral energy distributions to the z=3, 4, and 5 LBG stacked 24-850um fluxes, we derive infrared (IR) luminosities of L(8-1000um)~3.2, 5.5, and 11.0x10^11 Lsun (corresponding to star formation rates of ~50-200 Msun/yr) respectively. We find that the evolution in the IR luminosity density of LBGs is broadly consistent with model predictions for the expected contribution of luminous IR galaxy (LIRG) to ultraluminous IR galaxy (ULIRG) type systems at these epochs. We also see a strong positive correlation between stellar mass and IR luminosity. Our data are consistent with the main sequence of star formation showing little or no evolution from z=3 to 5. We have also confirmed that, for a fixed mass, the reddest LBGs (UV slope Beta -> 0) are indeed redder due to dust extinction, with SFR(IR)/SFR(UV) increasing by approximately an order of magnitude over -2<Beta<0 such that SFR(IR)/SFR(UV)~20 for the reddest LBGs. Furthermore, the most massive LBGs also tend to have higher obscured-to-unobscured ratio, hinting at a variation in the obscuration properties across the mass range.

We report the detection at 850um of the central source in SSA22-LAB1, the archetypal Lyman-alpha Blob (LAB), a 100kpc-scale radio-quiet emission-line nebula at z=3.1. The flux density of the source, $S_{850}=4.6\pm1.1$mJy implies the presence of a galaxy, or group of galaxies, with a total luminosity of $L_{\rm IR}\approx10^{12}L_\odot$. The position of an active source at the center of a ~50kpc-radius ring of linearly polarized Ly-alpha emission detected by Hayes et al. (2011) suggests that the central source is leaking Ly-alpha photons preferentially in the plane of the sky, which undergo scattering in HI clouds at large galactocentric radius. The Ly-alpha morphology around the submillimeter detection is reminiscent of biconical outflow, and the average Ly-alpha line profiles of the two `lobes' are dominated by a red peak, expected for a resonant line emerging from a medium with a bulk velocity gradient that is outflowing relative to the line center. Taken together, these observations provide compelling evidence that the central active galaxy (or galaxies) is responsible for a large fraction of the extended Ly-alpha emission and morphology. Less clear is the history of the cold gas in the circumgalactic medium being traced by Ly-alpha: is it mainly pristine material accreting into the halo that has not yet been processed through an interstellar medium (ISM), now being blown back as it encounters an outflow, or does it mainly comprise gas that has been swept-up within the ISM and expelled from the galaxy?

SPLAT-VO is a powerful graphical tool for displaying, comparing, modifying and analyzing astronomical spectra, as well as searching and retrieving spectra from services around the world using Virtual Observatory (VO) protocols and services. The development of SPLAT-VO started in 1999, as part of the Starlink StarJava initiative, sometime before that of the VO, so initial support for the VO was necessarily added once VO standards and services became available. Further developments were supported by the Joint Astronomy Centre, Hawaii until 2009. Since end of 2011 development of SPLAT-VO has been continued by the German Astrophysical Virtual Observatory, and the Astronomical Institute of the Academy of Sciences of the Czech Republic. From this time several new features have been added, including support for the latest VO protocols, along with new visualization and spectra storing capabilities. This paper presents the history of SPLAT-VO, it's capabilities, recent additions and future plans, as well as a discussion on the motivations and lessons learned up to now.

CCAT will be a 25m diameter sub-millimeter telescope capable of operating in the 0.2 to 2.1mm wavelength range. It will be located at an altitude of 5600m on Cerro Chajnantor in northern Chile near the ALMA site. The anticipated first generation instruments include large format (60,000 pixel) kinetic inductance detector (KID) cameras, a large format heterodyne array and a direct detection multi-object spectrometer. The paper describes the architecture of the CCAT software and the development strategy.

Current ground-based submillimeter instruments (e.g. SCUBA-2, SHARC-2 and LABOCA) have hundreds to thousands of detectors, sampled at tens to hundreds of hertz, generating up to hundreds of gigabytes per night. Since noise is correlated between detectors and in time, due to atmospheric signals and temperature oscillations, naive map-making is not applicable. In addition, the size of the data sets makes direct likelihood based inversion techniques intractable. As a result, the data reduction approach for most current submm cameras is to adopt iterative methods in order to separate noise from sky signal, and hence effectively produce astronomical images. We investigate how today's map-makers scale to the next generation of instruments, which will have tens of thousands of detectors sampled at thousands of hertz, leading to data sets of challenging size. We propose strategies for reducing such large data sets.

CCAT is a sub-millimeter telescope to be built on Cerro Chajnantor in Chile near the ALMA site. The remote location means that all observing will be done by remote observers with the future goal of fully autonomous observing using a dynamic scheduler. The fully autonomous observing mode provides a natural means for accepting transient alert notifications for immediate follow up.

CCAT will be a 25-meter telescope for sub millimeter astronomy located at 5600m altitude on Cerro Chajnantor in northern Chile. CCAT will combine high sensitivity, a wide field of view, and a broad wavelength range (0.35 to 2.1mm) to provide an unprecedented capability for deep, large-area multicolor submillimeter surveys. It is planned to have a suite of instruments including large format KID cameras, a large heterodyne array and a KID-based direct detection multi-object spectrometer. The remote location drives a desire for fully autonomous observing coupled with data reduction pipelines and fast feedback to principal investigators.

We mapped five massive star forming regions with the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT). Temperature and column density maps are obtained from the SCUBA-2 450 and 850 $\mu$m images. Most of the dense clumps we find have central temperatures below 20 K with some as cold as 8 K, suggesting that they have no internal heating due to the presence of embedded protostars. This is surprising, because at the high densities inferred from these images and at these low temperatures such clumps should be unstable, collapsing to form stars and generating internal heating. The column densities at the clump centres exceed 10$^{23}$ cm$^{-2}$, and the derived peak visual extinction values are from 25-500 mag for $\beta$ = 1.5-2.5, indicating highly opaque centres. The observed cloud gas masses range from $\sim$ 10 to 10$^{3}$ M$_{\odot}$. The outer regions of the clumps follow an $r^{-2.36\pm0.35}$ density distribution and this power-law structure is observed outside of typically 10$^{4}$ AU. All these findings suggest that these clumps are high-mass starless clumps and most likely contain high-mass starless cores.

MS$\,$0451.6$-$0305 is a rich galaxy cluster whose strong lensing is particularly prominent at submm wavelengths. We combine new SCUBA-2 data with imaging from Herschel SPIRE and PACS and HST in order to try to understand the nature of the sources being lensed. In the region of the "giant submm arc," we uncover seven multiply imaged galaxies (up from the previously known three), of which six are found to be at a redshift of $z\sim2.9$, and possibly constitute an interacting system. Using a novel forward-modelling approach, we are able to simultaneously deblend and fit SEDs to the individual galaxies that contribute to the giant submm arc, constraining their dust temperatures, far infrared luminosities and star formation rates. The submm arc first identified by SCUBA can now be seen to be composed of at least five distinct sources, four of these within the galaxy group at $z\sim2.9$. The total unlensed luminosity for this galaxy group is $(3.1\pm0.3) \times 10^{12}\,\mathrm{L}_\odot$, which gives an unlensed star formation rate of $(450\pm50)$ M$_\odot$ yr$^{-1}$. From the properties of this system, we see no evidence of evolution towards lower temperatures in the dust temperature versus far-infrared luminosity relation for high redshift galaxies.

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

The Submillimetre Common User Bolometer Array 2 (SCUBA-2) is an instrument operating on the 15-m James Clerk Maxwell Telescope, nominally consisting of 5120 bolometers in each of two simultaneous imaging bands centred over 450 and 850 um. The camera is operated by scanning across the sky and recording data at a rate of 200 Hz. As the largest of a new generation of multiplexed kilopixel bolometer cameras operating in the (sub)millimetre, SCUBA-2 data analysis represents a significant challenge. We describe the production of maps using the Sub-Millimetre User Reduction Facility (SMURF) in which we have adopted a fast, iterative approach to map-making that enables data reduction on single, modern, high-end desktop computers, with execution times that are typically shorter than the observing times. SMURF is used in an automated setting, both at the telescope for real-time feedback to observers, as well as for the production of science products for the JCMT Science Archive at the Canadian Astronomy Data Centre. Three detailed case studies are used to: (i) explore convergence properties of the map-maker using simple prior constraints (Uranus -- a point source); (ii) achieve the white-noise limit for faint point-source studies (extragalactic blank-field survey of the Lockman Hole); and (iii) demonstrate that our strategy is capable of recovering angular scales comparable to the size of the array footprint (approximately 5 arcmin) for bright extended sources (star-forming region M17).

We present Herschel observations from the Herschel Gould Belt Survey and SCUBA-2 science verification observations from the JCMT Gould Belt Survey of the B1 clump in the Perseus molecular cloud. We determined the dust emissivity index using four different techniques to combine the Herschel PACS+SPIRE data at 160 - 500 microns with the SCUBA-2 data at 450 microns and 850 microns. Of our four techniques, we found the most robust method was to filter-out the large-scale emission in the Herschel bands to match the spatial scales recovered by the SCUBA-2 reduction pipeline. Using this method, we find beta ~ 2 towards the filament region and moderately dense material and lower beta values (beta > 1.6) towards the dense protostellar cores, possibly due to dust grain growth. We find that beta and temperature are more robust with the inclusion of the SCUBA-2 data, improving estimates from Herschel data alone by factors of ~ 2 for beta and by ~ 40% for temperature. Furthermore, we find core mass differences of < 30% compared to Herschel-only estimates with an adopted beta = 2, highlighting the necessity of long wavelength submillimeter data for deriving accurate masses of prestellar and protostellar cores.

SCUBA-2 is an innovative 10000 pixel bolometer camera operating at submillimetre wavelengths on the James Clerk Maxwell Telescope (JCMT). The camera has the capability to carry out wide-field surveys to unprecedented depths, addressing key questions relating to the origins of galaxies, stars and planets. With two imaging arrays working simultaneously in the atmospheric windows at 450 and 850 microns, the vast increase in pixel count means that SCUBA-2 maps the sky 100-150 times faster than the previous SCUBA instrument. In this paper we present an overview of the instrument, discuss the physical characteristics of the superconducting detector arrays, outline the observing modes and data acquisition, and present the early performance figures on the telescope. We also showcase the capabilities of the instrument via some early examples of the science SCUBA-2 has already undertaken. In February 2012, SCUBA-2 began a series of unique legacy surveys for the JCMT community. These surveys will take 2.5 years and the results are already providing complementary data to the shorter wavelength, shallower, larger-area surveys from Herschel. The SCUBA-2 surveys will also provide a wealth of information for further study with new facilities such as ALMA, and future telescopes such as CCAT and SPICA.

SCUBA-2 is a 10000-bolometer submillimetre camera on the James Clerk Maxwell Telescope (JCMT). The instrument commissioning was completed in September 2011, and full science operations began in October 2011. To harness the full potential of this powerful new astronomical tool, the instrument calibration must be accurate and well understood. To this end, the algorithms for calculating the line-of-sight opacity have been improved, and the derived atmospheric extinction relationships at both wavebands of the SCUBA-2 instrument are presented. The results from over 500 primary and secondary calibrator observations have allowed accurate determination of the flux conversion factors (FCF) for the 850 and 450 micron arrays. Descriptions of the instrument beam-shape and photometry methods are presented. The calibration factors are well determined, with relative calibration accuracy better than 5 per cent at 850 microns and 10 per cent at 450 microns, reflecting the success of the derived opacity relations as well as the stability of the performance of the instrument over several months. The sample-size of the calibration observations and accurate FCFs have allowed the determination of the 850 and 450 micron fluxes of several well-known submillimetre sources, and these results are compared with previous measurements from SCUBA.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

Recent developments in the AST library are described, including a Python interface, support for the FITS-WCS "-TAB" system for storing tabular co-ordinate information, and extended support for representing distortions in spatial projections, using several schemes in common use (IRAF TNX/ZPX, Spitzer SIP, NOAO TPV and SCAMP).

SCUBA-2 is the largest submillimetre wide-field bolometric camera ever built. This 43 square arc-minute field-of-view instrument operates at two wavelengths (850 and 450 microns) and has been installed on the James Clerk Maxwell Telescope on Mauna Kea, Hawaii. SCUBA-2 has been successfully commissioned and operational for general science since October 2011. This paper presents an overview of the on-sky performance of the instrument during and since commissioning in mid-2011. The on-sky noise characteristics and NEPs of the 450 and 850 micron arrays, with average yields of approximately 3400 bolometers at each wavelength, will be shown. The observing modes of the instrument and the on-sky calibration techniques are described. The culmination of these efforts has resulted in a scientifically powerful mapping camera with sensitivities that allow a square degree of sky to be mapped to 10 mJy/beam rms at 850 micron in 2 hours and 60 mJy/beam rms at 450 micron in 5 hours in the best weather.

For more than a decade the Joint Astronomy Centre has been developing software tools to simplify observing and make it possible to use the telescopes in many different operational modes. In order to support remote operations the data handling systems need to be in place to allow observation preparation, flexible queue scheduling, data quality pipelines and science archives all to be connected in a data-driven environment. We discuss the history of these developments at UKIRT and JCMT and how the decision to combine software development at both telescopes led each to get features that they could not have justified if they were treated independently.

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

SCUBA-2 is the largest submillimetre array camera in the world and was commissioned on the James Clerk Maxwell Telescope (JCMT) with two arrays towards the end of 2009. A period of shared-risks observing was then completed and the full planned complement of 8 arrays, 4 at 850 microns and 4 at 450 microns, are now installed and ready to be commissioned. SCUBA-2 has 10,240 bolometers, corresponding to a data rate of 8 MB/s when sampled at the nominal rate of 200 Hz. The pipeline produces useful maps in near real time at the telescope and often publication quality maps in the JCMT Science Archive (JSA) hosted at the Canadian Astronomy Data Centre (CADC).

Commissioning of SCUBA-2 included a program of skydips and observations of calibration sources intended to be folded into regular observing as standard methods of source flux calibration and to monitor the atmospheric opacity and stability. During commissioning, it was found that these methods could also be utilised to characterise the fundamental instrument response to sky noise and astronomical signals. Novel techniques for analysing on-sky performance and atmospheric conditions are presented, along with results from the calibration observations and skydips.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We present a catalogue of Molecular Hydrogen emission-line Objects (MHOs) in outflows from young stars, most of which are embedded. All objects are identified in the near-infrared lines of molecular hydrogen, all reside in the Milky Way, and all are associated with jets or molecular outflows from young stars. Objects in both low and high-mass star forming regions are included. This catalogue complements the existing database of Herbig-Haro objects; indeed, for completeness, HH objects that are detected in H2 emission are included in the MHO catalogue.

Calibrated data for 143 flat-spectrum extragalactic radio sources are presented at a wavelength of 850 microns covering a five-year period from April 2000. The data, obtained at the James Clerk Maxwell Telescope using the SCUBA camera in pointing mode, were analysed using an automated pipeline process based on the Observatory Reduction and Acquisition Control - Data Reduction (ORAC-DR) system. This paper describes the techniques used to analyse and calibrate the data, and presents the database of results along with a representative sample of the better-sampled lightcurves. A re-analysis of previously published data from 1997 to 2000 is also presented. The combined catalogue, comprising 10493 flux density measurements, provides a unique and valuable resource for studies of extragalactic radio sources.

This paper describes a new Heterodyne Array Receiver Programme (HARP) and Auto-Correlation Spectral Imaging System (ACSIS) that have recently been installed and commissioned on the James Clerk Maxwell Telescope (JCMT). The 16-element focal-plane array receiver, operating in the submillimetre from 325 to 375 GHz, offers high (three-dimensional) mapping speeds, along with significant improvements over single-detector counterparts in calibration and image quality. Receiver temperatures are $\sim$120 K across the whole band and system temperatures of $\sim$300K are reached routinely under good weather conditions. The system includes a single-sideband filter so these are SSB figures. Used in conjunction with ACSIS, the system can produce large-scale maps rapidly, in one or more frequency settings, at high spatial and spectral resolution. Fully-sampled maps of size 1 square degree can be observed in under 1 hour. The scientific need for array receivers arises from the requirement for programmes to study samples of objects of statistically significant size, in large-scale unbiased surveys of galactic and extra-galactic regions. Along with morphological information, the new spectral imaging system can be used to study the physical and chemical properties of regions of interest. Its three-dimensional imaging capabilities are critical for research into turbulence and dynamics. In addition, HARP/ACSIS will provide highly complementary science programmes to wide-field continuum studies, and produce the essential preparatory work for submillimetre interferometers such as the SMA and ALMA.

We present a new large-scale survey of J=3-2 12CO emission covering 4.8 square degrees around the Rosette Nebula. Approximately 2000 compact clumps are identified, with a spatially-invariant power law mass distribution index of -1.8. Most of the inner clumps show velocity gradients of 1-3 km/s/pc, directed away from the exciting nebula. The gradients decrease with distance from the central O stars, and are consistent with a photoionised gas acceleration model, assuming clump lifetimes of a few 10^5 yrs. However, in one clear case, the observed near-constant velocity gradient is difficult to explain with simple models. Most blue-shifted but very few of the red-shifted clumps are associated with dark absorbing optical globules, confirming that the dominant molecular gas motion is expansion away from the central nebula and O stars. Many clumps also lie in a molecular ring, having an expansion velocity of 30 km/s, radius 11pc, and dynamical lifetime of ~1Myr. The J=3-2/1-0 12CO line ratios of the clumps decrease with distance from the O stars, implying a gradient in their surface temperatures; the results are consistent with a simple model of clump surface heating due to the central stars. Seven high-velocity molecular flows are found in the region, with a close correspondence between these flows and embedded young clusters. These outflows are sufficiently energetic to drive local gas turbulence within each cluster. We find 14 clear examples of association between embedded young stars seen at 24um and CO clumps; these are thought to be photoevaporating molecular envelopes. The CO clumps without evidence of embedded stars tend to have lower velocity gradients, and it is suggested that the presence of the young star may extend the lifespan of the photoevaporating envelope.

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 microns has been carried out to map the regions around several 850 micron-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared luminosities and hence star-formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350 micron flux densities for 24 850 micron sources, seven of which are detected at >2.5-sigma within a 10 arcsec search radius of the 850 micron positions. When results from the literature are included the total number of 350 micron flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-infrared (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich (Mdust~9x10^8 Msun), luminous (Lfir~2x10^12 Lsun), star-forming galaxies with intrinsic dust temperatures of ~35 K and star-formation rates of ~400 Msun/yr. We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350 micron source counts, N350(>25mJy)~200-500 deg^-2.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

SCUBA-2 is an innovative 10,000 pixel submillimeter camera due to be delivered to the James Clerk Maxwell Telescope in late 2006. The camera is expected to revolutionize submillimeter astronomy in terms of the ability to carry out wide-field surveys to unprecedented depths addressing key questions relating to the origins of galaxies, stars and planets. This paper presents an update on the project with particular emphasis on the laboratory commissioning of the instrument. The assembly and integration will be described as well as the measured thermal performance of the instrument. A summary of the performance results will be presented from the TES bolometer arrays, which come complete with in-focal plane SQUID amplifiers and multiplexed readouts, and are cooled to 100mK by a liquid cryogen-free dilution refrigerator. Considerable emphasis has also been placed on the operating modes of the instrument and the common-user aspect of the user interface and data reduction pipeline. These areas will also be described in the paper.

The L1551 molecular cloud contains two small clusters of Class 0 and I protostars, as well as a halo of more evolved Class II and III YSOs, indicating a current and at least one past burst of star formation. We present here new, sensitive maps of 850 and 450 um dust emission covering most of the L1551 cloud, new CO J=2-1 data of the molecular cloud, and a new, deep, optical image of [SII] emission. No new Class 0/I YSOs were detected. Compact sub-millimetre emitters are concentrated in two sub-clusters: IRS5 and L1551NE, and the HL~Tauri group. Both stellar groups show significant extended emission and outflow/jet activity. A jet, terminating at HH 265 and with a very weak associated molecular outflow, may originate from LkHa 358, or from a binary companion to another member of the HL Tauri group. Several Herbig Haro objects associated with IRS5/NE were clearly detected in the sub-mm, as were faint ridges of emission tracing outflow cavity walls. We confirm a large-scale molecular outflow originating from NE parallel to that from IRS5, and suggest that the "hollow shell" morphology is more likely due to two interacting outflows. We confirm the presence of a prestellar core (L1551-MC) of mass 2-3 Mo north-west of IRS5. The next generation cluster may be forming in this core. The L1551 cloud appears cometary in morphology, and appears to be illuminated and eroded from the direction of Orion, perhaps explaining the multiple episodes of star formation in this cloud. The full paper (including figures) can be downloaded at http://www.jach.hawaii.edu/~gms/l1551/l1551-apj641.pdf, or viewed at http://www.jach.hawaii.edu/~gms/l1551/.

We present 850micron imaging polarimetry of two high-mass star forming regions - W48 and S152. Within these regions we have identified two candidate high-mass protostellar objects - W48W and S152SE. The submillimetre continuum emission from the candidate HMPOs is bright in comparison to the nearby HII regions. W48W is a cold dense source, with no radio or mid-infrared emission. S152SE has an IRAS source IRAS 22566+5828 in the Southwestern part of the region, which appears in the mid-infrared 8.28micron emission, but there is no radio emission. The 850micron data shows another core within the region, in the Northeast. The polarimetry is ordered and the degree of polarisation is high over the candidate HMPOs - ~6% for W48W and ~8% for S152SE. Polarimetry results of this nature indicate a strong, ordered magnetic field threading the candidate HMPOs. The magnetic field direction in both S152SE and W48W is perpendicular to the direction of elongation of the cloud which would imply collapse along the field lines. Estimates of the magnetic field strength are derived using the Chandrasekhar & Fermi method. We calculate plane of the sky field strengths of ~0.7 mG for W48W and ~0.2 mG for S152SE. We discuss the drawbacks of using the Chandrasekhar & Fermi method with a large beam size. Mass-to-flux ratios have been calculated and both clouds are found to be roughly critical.

The S235A-B star forming region has been extensively observed in the past from the radio to the near-IR, but what was happening in the immediate surroundings of the water maser, placed in between the two nebulosities, was still unclear because of insufficient resolution especially in the spectral range from the Far-IR to the mm, even though there were sound indications that new young stellar objects (YSOs) are being formed there. We present here new high resolution maps at mm wavelengths in different molecules (HCO+, C34S, H2CS, SO2 and CH3CN), as well as in the 1.2 and 3.3 mm continuum obtained with the Plateau de Bure interferometer, and JCMT observations at 450 micron and 850 micron that unambiguously reveal the presence of new YSOs placed in between the two HII regions S235A and S235B and associated with the water maser. A molecular core and an unresolved source in the mm and in the sub-mm are centred on the maser, with indication of mass infall onto the core. Two molecular bipolar outflows and a jet originate from the same position. Weak evidence is found for a molecular rotating disk perpendicular to the direction of the main bipolar outflow. The derived parameters indicate that one of the YSOs is an intermediate luminosity object (L~1000 Lsun) in a very early evolutionary phase, embedded in a molecular core of ~100 Msun, with a temperature of 30 K. The main source of energy for the YSO could come from gravitational infall, thus making of this YSO a rare example of intermediate luminosity protostar representing a link between the earliest evolutionary phases of massive stars and low mass protostars of class 0-I.

We present an extended analysis of the SCUBA observations of the Hubble Deep Field (HDF), expanding the areal coverage of the Hughes et al. 1998 study by a factor of ~1.8 and containing at least three further sources in addition to the five in that study. We also announce the public release of the reduced data products. The map is the deepest ever made in the sub-millimetre, obtained in excellent conditions (median 850um optical depth of 0.16). Two independent reductions were made, one with SURF and the other with a wholly algorithmic IDL analysis which we present in detail here. Of the three new sources, all appear to be at z>~0.9 and one is provisionally associated with an Extremely Red Object (I-K>5). There appears to be no significant cross-correlation signal between the 850um fluctuations and sources detected by ISOCAM, the VLA or Chandra, nor with Very Red Objects (I-K>4), nor quasars and quasar candidates in the HDF (notwithstanding a small number of individual weak candidate detections). This is consistent with interpretations where the 850um-selected galaxies are at higher redshifts than those currently probed by ISOCAM/VLA, and predominantly not Compton-thin AGN. There are only one or two compelling cases for the radio source being the sub-mm source. Nevertheless, most SCUBA-HDF point sources have a nearby radio source apparently well-separated from the sub-mm centroid.

A polarimeter has been built for use with the Submillimetre Common-User Bolometer Array (SCUBA), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. SCUBA is the first of a new generation of highly sensitive submillimetre cameras, and the UK/Japan Polarimeter adds a polarimetric imaging/photometry capability in the wavelength range 350 to 2000 microns. Early science results range from measuring the synchrotron polarization of the black hole candidate Sgr A* to mapping magnetic fields inferred from polarized dust emission in Galactic star-forming clouds. We describe the instrument design, performance, observing techniques and data reduction processes, along with an assessment of the current and future scientific capability.

The calibration of ground-based submillimetre observations has always been a difficult process. We discuss how to overcome the limitations imposed by the submillimetre atmosphere. Novel ways to improve line-of-sight opacity estimates are presented, resulting in tight relations between opacities at different wavelengths. The submillimetre camera SCUBA, mounted on the JCMT, is the first large-scale submillimetre array, and as such is ideal for combatting the effects of the atmosphere. For example, we find that the off-source pixels are crucial for removing sky-noise. Benefitting from several years of SCUBA operation, a database of deep SCUBA observations has been constructed to better understand the nature of sky-noise and the effects of the atmosphere on instrument sensitivity. This has revealed several results. Firstly, there is evidence for positive correlations between sky-noise and seeing and sky-noise and sky opacity. Furthermore, 850-micron and 450-micron sky-noise are clearly correlated, suggesting that 450-micron data may be used to correct 850-micron observations for sky-noise. Perhaps most important of all: if off-source bolometers are used for sky-noise removal, there is no correlation between instrument sensitivity and chop throw, for chop throws out to 180 arcsec. Understanding the effects of submillimetre seeing is also important, and we find that the JCMT beam is not significantly broadened by seeing, nor is there an obvious correlation between seeing and pointing excursions.

The Submillimetre Common User Bolometer Array (SCUBA) instrument has been operating on the James Clerk Maxwell Telescope (JCMT) since 1997. The data archive is now sufficiently large that it can be used to investigate instrumental properties and the variability of astronomical sources. This paper describes the automated calibration and reduction scheme used to process the archive data with particular emphasis on `jiggle-map' observations of compact sources. We demonstrate the validity of our automated approach at both 850- and 450-microns and apply it to several of the JCMT secondary flux calibrators. We determine light curves for the variable sources IRC+10216 and OH231.8. This automation is made possible by using the ORAC-DR data reduction pipeline, a flexible and extensible data reduction pipeline that is used on UKIRT and the JCMT.

We present 850 micron imaging polarimetry of the W51A massive star forming region performed with SCUBA on the JCMT. From the polarimetry we infer the column-averaged magnetic field direction, projected onto the plane of the sky. We find that the magnetic field geometry in the region is complicated. We compare the field geometry with 6 cm and CS J=7-6 emission and determine that the magnetic field must be relatively weak and plays a passive role, allowing itself to be shaped by pressure forces and dynamics in the ionised and neutral gases. Comparisons are drawn between our data and 1.3 mm BIMA interferometric polarimetry data, from which we conclude that the magnetic field must increase in importance as we move to smaller scales and closer to sites of active star formation.

Calibrated data for 65 flat-spectrum extragalactic radio sources are presented at a wavelength of 850 microns covering a three-year period from April 1997. The data, obtained from the James Clerk Maxwell Telescope using the SCUBA camera in pointing mode, were analysed using an automated pipeline process based on the Observatory Reduction and Acquisition Control - Data Reduction (ORAC-DR) system. This paper describes the techniques used to analyse and calibrate the data and presents the database of results along with a representative sample of the better-sampled lightcurves.

We present first results from a submillimetre continuum survey of the Galactic Centre `Central Molecular Zone' (CMZ), made with SCUBA on the James Clerk Maxwell Telescope. SCUBA's scan-map mode has allowed us to make extremely wide-field maps of thermal dust emission with unprecedented speed and sensitivity. We also discuss some issues related to the elimination of artefacts in scan-map data. Our simultaneous 850/450 micron maps have a total size of approximately 2.8 x 0.5 degrees (400 x 75 pc) elongated along the galactic plane. They cover the Sgr A region-including Sgr A*, the circumnuclear disc, and the +20 km/s and +50 km/s clouds; the area around the Pistol; Sgr B2-the brightest feature on the map; and at their Galactic Western and Eastern edges the Sgr C and Sgr D regions. There are many striking features such as filaments and shell-like structures, as well as point sources such as Sgr A* itself. The total mass in the Central Molecular Zone is greater than that revealed in previous optically-thin molecular line maps by a factor of ~3, and new details are revealed on scales down to 0.33 pc across this 400 pc wide region.

We use an 850-micron SCUBA map of the HDF to study the dust properties of optically-selected starburst galaxies at high redshift. Although the 850-micron map does not allow detection of the majority of individual sources, we show that the galaxies with the highest UV star-formation rates are detected statistically, with a flux density of about S_{850}=0.2 mJy for an apparent UV star-formation rate of 1 h^{-2} msun yr^{-1}. This level of submillimetre output indicates that the total star-forming activity is on average a factor of approximately 6 times larger than the rate inferred from the UV output of these galaxies. The general population of optical starbursts is then predicted to contribute at least 25 per cent of the 850-micron background. We carry out a power-spectrum analysis of the map, which yields some evidence for angular clustering of the background source population, but at a level lower than that seen in Lyman-break galaxies. Together with other lines of argument, particularly from the NICMOS HDF data, this suggests that the 850-micron background originates over an extremely wide range of redshifts -- perhaps 1 < z < 6.

We present SCUBA observations of the protomultiple system NGC1333/IRAS4 at 450um and 850um. The 850um map shows significant extended emission which is most probably a remnant of the initial cloud core. At 450um, the component 4A is seen to have an elongated shape suggestive of a disk. Also we confirm that in addition to the 4A and 4B system, there exists another component 4C, which appears to lie out of the plane of the system and of the extended emission. Deconvolution of the beam reveals a binary companion to IRAS4B. Simple considerations of binary dynamics suggest that this triple 4A-4BI-4BII system is unstable and will probably not survive in its current form. Thus IRAS4 provides evidence that systems can evolve from higher to lower multiplicity as they move towards the main sequence. We construct a map of spectral index from the two wavelengths, and comment on the implications of this for dust evolution and temperature differences across the map. There is evidence that in the region of component 4A the dust has evolved, probably by coagulating into larger or more complex grains. Furthermore, there is evidence from the spectral index maps that dust from this object is being entrained in its associated outflow.

SCUBA, the Submillimetre Common-User Bolometer Array, built by the Royal Observatory Edinburgh for the James Clerk Maxwell Telescope, is the most versatile and powerful of a new generation of submillimetre cameras. It combines a sensitive dual-waveband imaging array with a three-band photometer, and is sky-background limited by the emission from the Mauna Kea atmosphere at all observing wavelengths from 350 microns to 2 mm. The increased sensitivity and array size mean that SCUBA maps close to 10,000 times faster than its single-pixel predecessor (UKT14). SCUBA is a facility instrument, open to the world community of users, and is provided with a high level of user support. We give an overview of the instrument, describe the observing modes and user interface, performance figures on the telescope, and present a sample of the exciting new results that have revolutionised submillimetre astronomy.

The Submillimetre Common-User Bolometer Array (SCUBA) is a new continuum camera operating on the James Clerk Maxwell Telescope (JCMT) on Mauna Kea, Hawaii. It consists of two arrays of bolometric detectors; a 91 pixel 350/450 micron array and a 37 pixel 750/850 micron array. Both arrays can be used simultaneously and have a field-of-view of approximately 2.4 arcminutes in diameter on the sky. Ideally, performance should be limited solely by the photon noise from the sky background at all wavelengths of operation. However, observations at submillimetre wavelengths are hampered by ``sky-noise'' which is caused by spatial and temporal fluctuations in the emissivity of the atmosphere above the telescope. These variations occur in atmospheric cells that are larger than the array diameter, and so it is expected that the resultant noise will be correlated across the array and, possibly, at different wavelengths. In this paper we describe our initial investigations into the presence of sky-noise for all the SCUBA observing modes, and explain our current technique for removing it from the data.

The Submillimetre Common-User Bolometer Array (SCUBA) is one of a new generation of cameras designed to operate in the submillimetre waveband. The instrument has a wide wavelength range covering all the atmospheric transmission windows between 300 and 2000 microns. In the heart of the instrument are two arrays of bolometers optimised for the short (350/450 microns) and long (750/850 microns) wavelength ends of the submillimetre spectrum. The two arrays can be used simultaneously, giving a unique dual-wavelength capability, and have a 2.3 arc-minute field of view on the sky. Background-limited performance is achieved by cooling the arrays to below 100 mK. SCUBA has now been in active service for over a year, and has already made substantial breakthroughs in many areas of astronomy. In this paper we present an overview of the performance of SCUBA during the commissioning phase on the James Clerk Maxwell Telescope (JCMT).

Dust emission around the nearby star epsilon Eridani has been imaged using a new submillimetre camera (SCUBA at the JCMT). At 850 microns wavelength a ring of dust is seen, peaking at 60 AU from the star and with much lower emission inside 30 AU. The mass of the ring is at least 0.01 Earth masses in dust, while an upper limit of 0.4 Earth masses in molecular gas is imposed by CO observations. The total mass is comparable to the estimated amount of material, 0.04-0.3 Earth masses, in comets orbiting the Solar System. The most probable origin of the the ring structure is that it is a young analogue to the Kuiper Belt in our Solar System, and that the central region has been partially cleared by the formation of grains into planetesimals. Dust clearing around epsilon Eri is seen within the radius of Neptune's orbit, and the peak emission at 35-75 AU lies within the estimated Kuiper Belt zone of 30-100 AU radius. epsilon Eri is a main-sequence star of type K2V (0.8 Solar masses) with an estimated age of 0.5-1.0 Gyr, so this interpretation is consistent with the early history of the Solar System where heavy bombardment occurred up to approximately 0.6 Gyr. An unexpected discovery is substructure within the ring, and these asymmetries could be due to perturbations by planets.

The advent of sensitive sub-mm array cameras now allows a proper census of dust-enshrouded massive star-formation in very distant galaxies, previously hidden activity to which even the faintest optical images are insensitive. We present the deepest sub-mm survey of the sky to date, taken with the SCUBA camera on the James Clerk Maxwell Telescope and centred on the Hubble Deep Field. The high source density found in this image implies that the survey is confusion-limited below a flux density of 2 mJy. However, within the central 80 arcsec radius independent analyses yield 5 reproducible sources with S(850um) > 2 mJy which simulations indicate can be ascribed to individual galaxies. We give positions and flux densities for these, and furthermore show using multi-frequency photometric data that the brightest sources in our map lie at redshifts z~3. These results lead to integral source counts which are completely inconsistent with a no-evolution model, and imply that massive star-formation activity continues at redshifts > 2. The combined brightness of the 5 most secure sources in our map is sufficient to account for 30 - 50% of the previously unresolved sub-mm background, and we estimate statistically that the entire background is resolved at about the 0.3 mJy level. Finally we discuss possible optical identifications and redshift estimates for the brightest sources. One source appears to be associated with an extreme starburst galaxy at z~1, whilst the remaining four appear to lie in the redshift range 2 < z < 4. This implies a star-formation density over this redshift range that is at least five times higher than that inferred from the ultraviolet output of HDF galaxies.

We present results of an extensive search for molecular gas in 14 high redshift ($z>2$) radio galaxies. Radio galaxies with redshifts between 2.0 and 4.3 were observed during several sessions on two different single-dish telescopes (IRAM and JCMT) and two interferometric arrays (VLA and OVRO). No significant CO emission to a 2$\sigma$ limit of a few times 10^{10} K km s^{-1} pc^2 was detected in any of these objects. The limits on the CO emission achieved indicate that, assuming a CO-H$_2$ conversion factor similar to the Galactic value, the mass of enriched molecular gas in these galaxies is less than $10^{11} $M$_{\odot}$. This suggests that distant radio galaxies may not be forming stars at the extremely high rates that have earlier been proposed, although they may still be forming stars at rates comparable to those of local starburst galaxies.