An attitude determination system for balloon-borne experiments is presented. The system provides pointing information in azimuth and elevation for instruments flying on stratospheric balloons over Antarctica. In-flight attitude is given by the real-time combination of readings from star cameras, a magnetometer, sun sensors, GPS, gyroscopes, tilt sensors and an elevation encoder. Post-flight attitude reconstruction is determined from star camera solutions, interpolated by the gyroscopes using an extended Kalman Filter. The multi-sensor system was employed by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol), an experiment that measures polarized thermal emission from interstellar dust clouds. A similar system was designed for the upcoming flight of SPIDER, a Cosmic Microwave Background polarization experiment. The pointing requirements for these experiments are discussed, as well as the challenges in designing attitude reconstruction systems for high altitude balloon flights. In the 2010 and 2012 BLASTPol flights from McMurdo Station, Antarctica, the system demonstrated an accuracy of <5' rms in-flight, and <5" rms post-flight.

We present the thermal model of the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This instrument was successfully flown in two circumpolar flights from McMurdo, Antarctica in 2010 and 2012. During these two flights, BLASTPol obtained unprecedented information about the magnetic field in molecular clouds through the measurement of the polarized thermal emission of interstellar dust grains. The thermal design of the experiment addresses the stability and control of the payload necessary for this kind of measurement. We describe the thermal modeling of the payload including the sun-shielding strategy. We present the in-flight thermal performance of the instrument and compare the predictions of the model with the temperatures registered during the flight. We describe the difficulties of modeling the thermal behavior of the balloon-borne platform and establish landmarks that can be used in the design of future balloon-borne instruments.

We present the second generation BLASTbus electronics. The primary purposes of this system are detector readout, attitude control, and cryogenic housekeeping, for balloon-borne telescopes. Readout of neutron transmutation doped germanium (NTD-Ge) bolometers requires low noise and parallel acquisition of hundreds of analog signals. Controlling a telescope's attitude requires the capability to interface to a wide variety of sensors and motors, and to use them together in a fast, closed loop. To achieve these different goals, the BLASTbus system employs a flexible motherboard-daughterboard architecture. The programmable motherboard features a digital signal processor (DSP) and field-programmable gate array (FPGA), as well as slots for three daughterboards. The daughterboards provide the interface to the outside world, with versions for analog to digital conversion, and optoisolated digital input/output. With the versatility afforded by this design, the BLASTbus also finds uses in cryogenic, thermometry, and power systems. For accurate timing control to tie everything together, the system operates in a fully synchronous manner. BLASTbus electronics have been successfully deployed to the South Pole, and flown on stratospheric balloons.

Turbulence and magnetic fields are expected to be important for regulating molecular cloud formation and evolution. However, their effects on subparsec to 100 parsec scales, leading to the formation of starless cores, is not well understood. We investigate the prestellar core structure morphologies obtained from analysis of the Herschel-SPIRE 350 $\mu$m maps of the Lupus I cloud. This distribution is first compared on a statistical basis to the large scale shape of the main filament. We find the distribution of the elongation position angle of the cores to be consistent with a random distribution, which means no specific orientation of the morphology of the cores is observed with respect to a large-scale filament shape model for Lupus I, or relative to a large-scale bent filament model. This distribution is also compared to the mean orientation of the large-scale magnetic fields probed at 350 $\mu$m with the Balloon-borne Large Aperture Telescope for Polarimetry (BLASTPol) during its 2010 campaign. Here again we do not find any correlation between the core morphology distribution and the average orientation of the magnetic fields on parsec scales. Our main conclusion is that the local filament dynamics - including secondary filaments that often run orthogonally to the primary filament - and possibly small-scale variations in the local magnetic field direction, could be the dominant factors for explaining the final orientation of each core.

A cryogenic achromatic half-wave plate (HWP) for submillimetre astronomical polarimetry has been designed, manufactured, tested, and deployed in the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). The design is based on the five-slab Pancharatnam recipe and it works in the wavelength range 200-600 micron, making it the broadest-band HWP built to date at (sub)millimetre wavelengths. The frequency behaviour of the HWP has been fully characterised at room and cryogenic temperatures with incoherent radiation from a polarising Fourier transform spectrometer. We develop a novel empirical model, complementary to the physical and analytical ones available in the literature, that allows us to recover the HWP Mueller matrix and phase shift as a function of frequency and extrapolated to 4K. We show that most of the HWP non-idealities can be modelled by quantifying one wavelength-dependent parameter, the position of the HWP equivalent axes, which is then readily implemented in a map-making algorithm. We derive this parameter for a range of spectral signatures of input astronomical sources relevant to BLASTPol, and provide a benchmark example of how our method can yield improved accuracy on measurements of the polarisation angle on the sky at submillimetre wavelengths.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was created by adding polarimetric capability to the BLAST experiment that was flown in 2003, 2005, and 2006. BLASTPol inherited BLAST's 1.8 m primary and its Herschel/SPIRE heritage focal plane that allows simultaneous observation at 250, 350, and 500 {\mu}m. We flew BLASTPol in 2010 and again in 2012. Both were long duration Antarctic flights. Here we present polarimetry of the nearby filamentary dark cloud Lupus I obtained during the 2010 flight. Despite limitations imposed by the effects of a damaged optical component, we were able to clearly detect submillimeter polarization on degree scales. We compare the resulting BLASTPol magnetic field map with a similar map made via optical polarimetry (The optical data were published in 1998 by J. Rizzo and collaborators.). The two maps partially overlap and are reasonably consistent with one another. We compare these magnetic field maps to the orientations of filaments in Lupus I, and we find that the dominant filament in the cloud is approximately perpendicular to the large-scale field, while secondary filaments appear to run parallel to the magnetic fields in their vicinities. This is similar to what is observed in Serpens South via near-IR polarimetry, and consistent with what is seen in MHD simulations by F. Nakamura and Z. Li.

The submillimeter opacity of dust in the diffuse Galactic interstellar medium (ISM) has been quantified using a pixel-by-pixel correlation of images of continuum emission with a proxy for column density. We used three BLAST bands at 250, 350, and 500 \mu m and one IRAS at 100 \mu m. The proxy is the near-infrared color excess, E(J-Ks), obtained from 2MASS. Based on observations of stars, we show how well this color excess is correlated with the total hydrogen column density for regions of moderate extinction. The ratio of emission to column density, the emissivity, is then known from the correlations, as a function of frequency. The spectral distribution of this emissivity can be fit by a modified blackbody, whence the characteristic dust temperature T and the desired opacity \sigma_e(1200) at 1200 GHz can be obtained. We have analyzed 14 regions near the Galactic plane toward the Vela molecular cloud, mostly selected to avoid regions of high column density (N_H > 10^{22} cm^-2) and small enough to ensure a uniform T. We find \sigma_e(1200) is typically 2 to 4 x 10^{-25} cm^2/H and thus about 2 to 4 times larger than the average value in the local high Galactic latitude diffuse atomic ISM. This is strong evidence for grain evolution. There is a range in total power per H nucleon absorbed (re-radiated) by the dust, reflecting changes in the interstellar radiation field and/or the dust absorption opacity. These changes affect the equilibrium T, which is typically 15 K, colder than at high latitudes. Our analysis extends, to higher opacity and lower T, the trend of increasing opacity with decreasing T that was found at high latitudes. The recognition of changes in the emission opacity raises a cautionary flag because all column densities deduced from dust emission maps, and the masses of compact structures within them, depend inversely on the value adopted.

We present measurements of the mean mid-infrared-to-submillimeter flux densities of massive (M\ast \approx 2 \times 10^11 Msun) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS) catalog on maps: at 24 {\mu}m (Spitzer/MIPS); 70, 100, and 160{\mu}m (Herschel/PACS); 250, 350, 500{\mu}m (BLAST); and 870{\mu}m (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star-formation rate of SFR = 63 [48, 81] Msun yr^-1 . We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disk-like and spheroid-like, according to their Sersic indices, n. We find evidence that most of the star formation is occurring in n \leq 2 (disk-like) galaxies, with median [interquartile] SFR = 122 [100,150] Msun yr^-1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14 [9,20] Msun yr^-1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.

We constrain the evolution of the rest-frame far-infrared (FIR) luminosity function out to high redshift, by combining several pieces of complementary information provided by the deep Balloon-borne Large-Aperture Submillimeter Telescope surveys at 250, 350 and 500 micron, as well as other FIR and millimetre data. Unlike most other phenomenological models, we characterise the uncertainties in our fitted parameters using Monte Carlo Markov Chains. We use a bivariate local luminosity function that depends only on FIR luminosity and 60-to-100 micron colour, along with a single library of galaxy spectral energy distributions indexed by colour, and apply simple luminosity and density evolution. We use the surface density of sources, Cosmic Infrared Background (CIB) measurements and redshift distributions of bright sources, for which identifications have been made, to constrain this model. The precise evolution of the FIR luminosity function across this crucial range has eluded studies at longer wavelengths (e.g., using SCUBA and MAMBO) and at shorter wavelengths (e.g., Spitzer), and should provide a key piece of information required for the study of galaxy evolution. Our adoption of Monte Carlo methods enables us not only to find the best-fit evolution model, but also to explore correlations between the fitted parameters. Our model-fitting approach allows us to focus on sources of tension coming from the combination of data-sets. We specifically find that our choice of parameterisation has difficulty fitting the combination of CIB measurements and redshift distribution of sources near 1 mm. Existing and future data sets will be able to dramatically improve the fits, as well as break strong degeneracies among the models. [abridged]

We identify and investigate the nature of the 20 brightest 250-micron sources detected by the BLAST within the central 150 sq. arcmin of the GOODS-South field. Aided by the available deep VLA radio imaging, reaching S_1.4 = 30 micro-Jy, we have identified radio counterparts for 17/20 of the 250-micron sources. The resulting enhanced positional accuracy of ~1 arcsec has then allowed us to exploit the deep multi-frequency imaging of GOODS-South to establish secure galaxy counterparts for the 17 radio-identified sources, and plausible galaxy candidates for the 3 radio-unidentified sources. Confusion is a serious issue for this deep BLAST 250-micron survey, due to the large size of the beam. Nevertheless, we argue that our chosen counterparts are significant, and often dominant contributors to the measured BLAST flux densities. For all of these 20 galaxies we have been able to determine spectroscopic (8) or photometric (12) redshifts. The result is the first `complete' redshift distribution for a deep 250-micron selected galaxy sample. This reveals that 250-micron surveys reaching detection limits of ~30 mJy contain not only low-redshift spirals/LIRGs, but also the extreme z~2 dust-enshrouded starburst galaxies previously discovered at sub-millimetre wavelengths. Based on their IRAC colours, we find that virtually all of the BLAST galaxy identifications appear better described as analogues of the M82 starburst galaxy, or Sc star-forming discs rather than highly obscured ULIRGs. Inspection of the LABOCA 870-micron imaging of the GOODS-South field yields detections of 7/11 of the z>1 BLAST sources, and reveals 250/870 flux-density ratios consistent with a standard 40K modified black-body fit with a dust emissivity index beta=1.5.

We have carried out the first general submillimeter analysis of the field towards GRSMC 45.46+0.05, a massive star forming region in Aquila. The deconvolved 6 deg^2 (3\degree X 2\degree) maps provided by BLAST in 2005 at 250, 350, and 500 micron were used to perform a preliminary characterization of the clump population previously investigated in the infrared, radio, and molecular maps. Interferometric CORNISH data at 4.8 GHz have also been used to characterize the Ultracompact HII regions (UCHIIRs) within the main clumps. By means of the BLAST maps we have produced an initial census of the submillimeter structures that will be observed by Herschel, several of which are known Infrared Dark Clouds (IRDCs). Our spectral energy distributions of the main clumps in the field, located at ~7 kpc, reveal an active population with temperatures of T~35-40 K and masses of ~10^3 Msun for a dust emissivity index beta=1.5. The clump evolutionary stages range from evolved sources, with extended HII regions and prominent IR stellar population, to massive young stellar objects, prior to the formation of an UCHIIR.The CORNISH data have revealed the details of the stellar content and structure of the UCHIIRs. In most cases, the ionizing stars corresponding to the brightest radio detections are capable of accounting for the clump bolometric luminosity, in most cases powered by embedded OB stellar clusters.

We present Cygnus X in a new multi-wavelength perspective based on an unbiased BLAST survey at 250, 350, and 500 micron, combined with rich datasets for this well-studied region. Our primary goal is to investigate the early stages of high mass star formation. We have detected 184 compact sources in various stages of evolution across all three BLAST bands. From their well-constrained spectral energy distributions, we obtain the physical properties mass, surface density, bolometric luminosity, and dust temperature. Some of the bright sources reaching 40 K contain well-known compact H II regions. We relate these to other sources at earlier stages of evolution via the energetics as deduced from their position in the luminosity-mass (L-M) diagram. The BLAST spectral coverage, near the peak of the spectral energy distribution of the dust, reveals fainter sources too cool (~ 10 K) to be seen by earlier shorter-wavelength surveys like IRAS. We detect thermal emission from infrared dark clouds and investigate the phenomenon of cold ``starless cores" more generally. Spitzer images of these cold sources often show stellar nurseries, but these potential sites for massive star formation are ``starless" in the sense that to date there is no massive protostar in a vigorous accretion phase. We discuss evolution in the context of the L-M diagram. Theory raises some interesting possibilities: some cold massive compact sources might never form a cluster containing massive stars; and clusters with massive stars might not have an identifiable compact cold massive precursor.

We carry out a multi-wavelength study of individual galaxies detected by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) and identified at other wavelengths, using data spanning the radio to the ultraviolet (UV). We develop a Monte Carlo method to account for flux boosting, source blending, and correlations among bands, which we use to derive deboosted far-infrared (FIR) luminosities for our sample. We estimate total star-formation rates for BLAST counterparts with z < 0.9 by combining their FIR and UV luminosities. Star formation is heavily obscured at L_FIR > 10^11 L_sun, z > 0.5, but the contribution from unobscured starlight cannot be neglected at L_FIR < 10^11 L_sun, z < 0.25. We assess that about 20% of the galaxies in our sample show indication of a type-1 active galactic nucleus (AGN), but their submillimeter emission is mainly due to star formation in the host galaxy. We compute stellar masses for a subset of 92 BLAST counterparts; these are relatively massive objects, with a median mass of ~10^11 M_sun, which seem to link the 24um and SCUBA populations, in terms of both stellar mass and star-formation activity. The bulk of the BLAST counterparts at z<1 appear to be run-of-the-mill star-forming galaxies, typically spiral in shape, with intermediate stellar masses and practically constant specific star-formation rates. On the other hand, the high-z tail of the BLAST counterparts significantly overlaps with the SCUBA population, in terms of both star-formation rates and stellar masses, with observed trends of specific star-formation rate that support strong evolution and downsizing.

We present an implementation of the iterative flux-conserving Lucy-Richardson (L-R) deconvolution method of image restoration for maps produced by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). We have analyzed its performance and convergence extensively through simulations and cross-correlations of the deconvolved images with available highresolution maps. We present new science results from two BLAST surveys, in the Galactic regions K3-50 and IC 5146, further demonstrating the benefits of performing this deconvolution. We have resolved three clumps within a radius of 4.'5 inside the star-forming molecular cloud containing K3-50. Combining the well-resolved dust emission map with available multi-wavelength data, we have constrained the Spectral Energy Distributions (SEDs) of five clumps to obtain masses (M), bolometric luminosities (L), and dust temperatures (T). The L-M diagram has been used as a diagnostic tool to estimate the evolutionary stages of the clumps. There are close relationships between dust continuum emission and both 21-cm radio continuum and 12CO molecular line emission. The restored extended large scale structures in the Northern Streamer of IC 5146 have a strong spatial correlation with both SCUBA and high resolution extinction images. A dust temperature of 12 K has been obtained for the central filament. We report physical properties of ten compact sources, including six associated protostars, by fitting SEDs to multi-wavelength data. All of these compact sources are still quite cold (typical temperature below ~ 16 K) and are above the critical Bonner-Ebert mass. They have associated low-power Young Stellar Objects (YSOs). Further evidence for starless clumps has also been found in the IC 5146 region.

We present results from a survey carried out by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) on a 9 deg^2 field near the South Ecliptic Pole at 250, 350 and 500 {\mu}m. The median 1{\sigma} depths of the maps are 36.0, 26.4 and 18.4 mJy, respectively. We apply a statistical method to estimate submillimeter galaxy number counts and find that they are in agreement with other measurements made with the same instrument and with the more recent results from Herschel/SPIRE. Thanks to the large field observed, the new measurements give additional constraints on the bright end of the counts. We identify 132, 89 and 61 sources with S/N>4 at 250, 350, 500 {\mu}m, respectively and provide a multi-wavelength combined catalog of 232 sources with a significance >4{\sigma} in at least one BLAST band. The new BLAST maps and catalogs are available publicly at http://blastexperiment.info.

We have imaged an 11.5 sq. deg. region of sky towards the South Ecliptic Pole (RA = 04h43m, Dec = -53d40m, J2000) at 24 and 70 microns with MIPS, the Multiband Imaging Photometer for Spitzer. This region is coincident with a field mapped at longer wavelengths by AKARI and the Balloon-borne Large Aperture Submillimeter Telescope. We discuss our data reduction and source extraction procedures. The median depths of the maps are 47 microJy/beam at 24 micron and 4.3 mJy/beam at 70 micron. At 24 micron, we identify 93098 point sources with signal-to-noise ratio (SNR) >5, and an additional 63 resolved galaxies; at 70 micron, we identify 891 point sources with SNR >6. From simulations, we determine a false detection rate of 1.8% (1.1%) for the 24 micron (70 micron) catalog. The 24 and 70 micron point-source catalogs are 80% complete at 230 microJy and 11 mJy, respectively. These mosaic images and source catalogs will be available to the public through the NASA/IPAC Infrared Science Archive.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLAST-Pol) is a suborbital mapping experiment designed to study the role played by magnetic fields in the star formation process. BLAST-Pol is the reconstructed BLAST telescope, with the addition of linear polarization capability. Using a 1.8 m Cassegrain telescope, BLAST-Pol images the sky onto a focal plane that consists of 280 bolometric detectors in three arrays, observing simultaneously at 250, 350, and 500 um. The diffraction-limited optical system provides a resolution of 30'' at 250 um. The polarimeter consists of photolithographic polarizing grids mounted in front of each bolometer/detector array. A rotating 4 K achromatic half-wave plate provides additional polarization modulation. With its unprecedented mapping speed and resolution, BLAST-Pol will produce three-color polarization maps for a large number of molecular clouds. The instrument provides a much needed bridge in spatial coverage between larger-scale, coarse resolution surveys and narrow field of view, and high resolution observations of substructure within molecular cloud cores. The first science flight will be from McMurdo Station, Antarctica in December 2010.

We use new large area far infrared maps ranging from 65 - 500 microns obtained with the AKARI and the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) missions to characterize the dust emission toward the Cassiopeia A supernova remnant (SNR). Using the AKARI high resolution data we find a new "tepid" dust grain population at a temperature of ~35K and with an estimated mass of 0.06 solar masses. This component is confined to the central area of the SNR and may represent newly-formed dust in the unshocked supernova ejecta. While the mass of tepid dust that we measure is insufficient by itself to account for the dust observed at high redshift, it does constitute an additional dust population to contribute to those previously reported. We fit our maps at 65, 90, 140, 250, 350, and 500 microns to obtain maps of the column density and temperature of "cold" dust (near 16 K) distributed throughout the region. The large column density of cold dust associated with clouds seen in molecular emission extends continuously from the surrounding interstellar medium to project on the SNR, where the foreground component of the clouds is also detectable through optical, X-ray, and molecular extinction. At the resolution available here, there is no morphological signature to isolate any cold dust associated only with the SNR from this confusing interstellar emission. Our fit also recovers the previously detected "hot" dust in the remnant, with characteristic temperature 100 K.

We present observations at 250, 350, and 500 um of the nearby galaxy cluster Abell 3112 (z=0.075) carried out with BLAST, the Balloon-borne Large Aperture Submillimeter Telescope. Five cluster members are individually detected as bright submillimetre sources. Their far-infrared SEDs and optical colours identify them as normal star-forming galaxies of high mass, with globally evolved stellar populations. They all have B-R colours of 1.38+/-0.08, transitional between the blue, active population and the red, evolved galaxies that dominate the cluster core. We stack to determine the mean submillimetre emission from all cluster members, which is determined to be 16.6+/-2.5, 6.1+/-1.9, and 1.5+/-1.3 mJy at 250, 350, and 500 um, respectively. Stacking analyses of the submillimetre emission of cluster members reveal trends in the mean far-infrared luminosity with respect to cluster-centric radius and Ks-band magnitude. We find that a large fraction of submillimetre emission comes from the boundary of the inner, virialized region of the cluster, at cluster-centric distances around R_500. Stacking also shows that the bulk of the submillimetre emission arises in intermediate-mass galaxies (L<L*), with Ks magnitude ~1 mag fainter than the giant ellipticals. The results and constraints obtained in this work will provide a useful reference for the forthcoming surveys to be conducted on galaxy clusters by Herschel.

The Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) carried out a 250, 350 and 500 micron survey of the galactic plane encompassing the Vela Molecular Ridge, with the primary goal of identifying the coldest dense cores possibly associated with the earliest stages of star formation. Here we present the results from observations of the Vela-D region, covering about 4 square degrees, in which we find 141 BLAST cores. We exploit existing data taken with the Spitzer MIPS, IRAC and SEST-SIMBA instruments to constrain their (single-temperature) spectral energy distributions, assuming a dust emissivity index beta = 2.0. This combination of data allows us to determine the temperature, luminosity and mass of each BLAST core, and also enables us to separate starless from proto-stellar sources. We also analyze the effects that the uncertainties on the derived physical parameters of the individual sources have on the overall physical properties of starless and proto-stellar cores, and we find that there appear to be a smooth transition from the pre- to the proto-stellar phase. In particular, for proto-stellar cores we find a correlation between the MIPS24 flux, associated with the central protostar, and the temperature of the dust envelope. We also find that the core mass function of the Vela-D cores has a slope consistent with other similar (sub)millimeter surveys.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has recently surveyed ~8.7 deg^2 centered on GOODS-South at 250, 350, and 500 microns. In Dye et al. (2009) we presented the catalogue of sources detected at 5-sigma in at least one band in this field and the probable counterparts to these sources in other wavebands. In this paper, we present the results of a redshift survey in which we succeeded in measuring redshifts for 82 of these counterparts. The spectra show that the BLAST counterparts are mostly star-forming galaxies but not extreme ones when compared to those found in the Sloan Digital Sky Survey. Roughly one quarter of the BLAST counterparts contain an active nucleus. We have used the spectroscopic redshifts to carry out a test of the ability of photometric redshift methods to estimate the redshifts of dusty galaxies, showing that the standard methods work well even when a galaxy contains a large amount of dust. We have also investigated the cases where there are two possible counterparts to the BLAST source, finding that in at least half of these there is evidence that the two galaxies are physically associated, either because they are interacting or because they are in the same large-scale structure. Finally, we have made the first direct measurements of the luminosity function in the three BLAST bands. We find strong evolution out to z=1, in the sense that there is a large increase in the space-density of the most luminous galaxies. We have also investigated the evolution of the dust-mass function, finding similar strong evolution in the space-density of the galaxies with the largest dust masses, showing that the luminosity evolution seen in many wavebands is associated with an increase in the reservoir of interstellar matter in galaxies.

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

We describe the application of a statistical method to estimate submillimeter galaxy number counts from confusion limited observations by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). Our method is based on a maximum likelihood fit to the pixel histogram, sometimes called 'P(D)', an approach which has been used before to probe faint counts, the difference being that here we advocate its use even for sources with relatively high signal-to-noise ratios. This method has an advantage over standard techniques of source extraction in providing an unbiased estimate of the counts from the bright end down to flux densities well below the confusion limit. We specifically analyse BLAST observations of a roughly 10 sq. deg. map centered on the Great Observatories Origins Deep Survey South (GOODS-S) field. We provide estimates of number counts at the three BLAST wavelengths, 250, 350, and 500 microns; instead of counting sources in flux bins we estimate the counts at several flux density nodes connected with power-laws. We observe a generally very steep slope for the counts of about -3.7 at 250 microns and -4.5 at 350 and 500 microns, over the range ~0.02-0.5 Jy, breaking to a shallower slope below about 0.015 Jy at all three wavelengths. We also describe how to estimate the uncertainties and correlations in this method so that the results can be used for model-fitting. This method should be well-suited for analysis of data from the Herschel satellite.

We detect correlations in the cosmic far-infrared background due to the clustering of star-forming galaxies in observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST, at 250, 350, and 500 microns. We perform jackknife and other tests to confirm the reality of the signal. The measured correlations are well fit by a power law over scales of 5-25 arcminutes, with Delta I/I = 15.1 +/- 1.7%. We adopt a specific model for submillimeter sources in which the contribution to clustering comes from sources in the redshift ranges 1.3 <= z <= 2.2, 1.5 <= z <= 2.7, and 1.7 <= z <= 3.2, at 250, 350, and 500 microns, respectively. With these distributions, our measurement of the power spectrum, P(k_theta), corresponds to linear bias parameters, b = 3.8 +/- 0.6, 3.9 +/- 0.6 and 4.4 +/- 0.7, respectively. We further interpret the results in terms of the halo model, and find that at the smaller scales, the simplest halo model fails to fit our results. One way to improve the fit is to increase the radius at which dark matter halos are artificially truncated in the model, which is equivalent to having some star-forming galaxies at z >= 1 located in the outskirts of groups and clusters. In the context of this model we find a minimum halo mass required to host a galaxy is log (M_min / M_sun) = 11.5 (+0.4/-0.1), and we derive effective biases $b_eff = 2.2 +/- 0.2, 2.4 +/- 0.2, and 2.6 +/- 0.2, and effective masses log (M_eff / M_sun) = 12.9 +/- 0.3, 12.8 +/- 0.2, and 12.7 +/- 0.2, at 250, 350, and 500 microns, corresponding to spatial correlation lengths of r_0 = 4.9, 5.0, and 5.2 +/- 0.7 h^-1 Mpc, respectively. Finally, we discuss implications for clustering measurement strategies with Herschel and Planck.

Over the course of two flights, the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) made resolved maps of seven nearby (<25 Mpc) galaxies at 250, 350, and 500 microns. During its June 2005 flight from Sweden, BLAST observed a single nearby galaxy, NGC 4565. During the December 2006 flight from Antarctica, BLAST observed the nearby galaxies NGC 1097, NGC 1291, NGC 1365, NGC 1512, NGC 1566, and NGC 1808. We fit physical dust models to a combination of BLAST observations and other available data for the galaxies observed by Spitzer. We fit a modified blackbody to the remaining galaxies to obtain total dust mass and mean dust temperature. For the four galaxies with Spitzer data, we also produce maps and radial profiles of dust column density and temperature. We measure the fraction of BLAST detected flux originating from the central cores of these galaxies and use this to calculate a "core fraction," an upper limit on the "AGN fraction" of these galaxies. We also find our resolved observations of these galaxies give a dust mass estimate 5-19 times larger than an unresolved observations would predict. Finally, we are able to use these data to derive a value for the dust mass absorption co-efficient of kappa = 0.29 +/-0.03 m^2 kg^-1 at 250 microns. This study is an introduction to future higher-resolution and higher-sensitivity studies to be conducted by Herschel and SCUBA-II.

We report multi-wavelength power spectra of diffuse Galactic dust emission from BLAST observations at 250, 350, and 500 microns in Galactic Plane fields in Cygnus X and Aquila. These submillimeter power spectra statistically quantify the self-similar structure observable over a broad range of scales and can be used to assess the cirrus noise which limits the detection of faint point sources. The advent of submillimeter surveys with the Herschel Space Observatory makes the wavelength dependence a matter of interest. We show that the observed relative amplitudes of the power spectra can be related through a spectral energy distribution (SED). Fitting a simple modified black body to this SED, we find the dust temperature in Cygnus X to be 19.9 +/- 1.3 K and in the Aquila region 16.9 +/- 0.7 K. Our empirical estimates provide important new insight into the substantial cirrus noise that will be encountered in forthcoming observations.

We present the 250, 350, and 500 micron detection of bright submillimeter emission in the direction of the Bullet Cluster measured by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). The 500 micron centroid is coincident with an AzTEC 1.1 mm point-source detection at a position close to the peak lensing magnification produced by the cluster. However, the 250 micron and 350 micron centroids are elongated and shifted toward the south with a differential shift between bands that cannot be explained by pointing uncertainties. We therefore conclude that the BLAST detection is likely contaminated by emission from foreground galaxies associated with the Bullet Cluster. The submillimeter redshift estimate based on 250-1100 micron photometry at the position of the AzTEC source is z_phot = 2.9 (+0.6 -0.3), consistent with the infrared color redshift estimation of the most likely IRAC counterpart. These flux densities indicate an apparent far-infrared luminosity of L_FIR = 2E13 Lsun. When the amplification due to the gravitational lensing of the cluster is removed, the intrinsic far-infrared luminosity of the source is found to be L_FIR <= 10^12 Lsun, consistent with typical luminous infrared galaxies.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has made one square degree, deep, confusion limited maps at three different bands, centered on the Great Observatories Origins Deep Survey South field. By calculating the covariance of these maps with catalogs of 24 micron sources from the Far-Infrared Deep Extragalactic Legacy Survey (FIDEL), we have determined that the total submillimeter intensities are 8.60 +/- 0.59, 4.93 +/- 0.34, and 2.27 +/- 0.20 nW m^2 sr^(-1) at 250, 350, and 500 micron, respectively. These numbers are more precise than previous estimates of the cosmic infrared background (CIB) and are consistent with 24 micron-selected galaxies generating the full intensity of the CIB. We find that the fraction of the CIB that originates from sources at z >= 1.2 increases with wavelength, with 60% from high redshift sources at 500 micron. At all BLAST wavelengths, the relative intensity of high-z sources is higher for 24 micron-faint sources than it is for 24 micron-bright sources. Galaxies identified as active galactic nuclei (AGN) by their Infrared Array Camera (IRAC) colors are 1.6-2.6 times brighter than the average population at 250-500 micron, consistent with what is found for X-ray-selected AGN. BzK-selected galaxies are found to be moderately brighter than typical 24 micron-selected galaxies in the BLAST bands. These data provide high precision constraints for models of the evolution of the number density and intensity of star forming galaxies at high redshift.

We directly measure redshift evolution in the mean physical properties (far-infrared luminosity, temperature, and mass) of the galaxies that produce the cosmic infrared background (CIB), using measurements from the Balloon-borne Large Aperture Sub-millimeter Telescope (BLAST), and Spitzer which constrain the CIB emission peak. This sample is known to produce a surface brightness in the BLAST bands consistent with the full CIB, and photometric redshifts are identified for all of the objects. We find that most of the 70 micron background is generated at z <~ 1 and the 500 micron background generated at z >~ 1. A significant growth is observed in the mean luminosity from ~ 10^9 - 10^12 L_sun, and in the mean temperature by 10 K, from redshifts 0< z < 3. However, there is only weak positive evolution in the comoving dust mass in these galaxies across the same redshift range. We also measure the evolution of the far-infrared luminosity density, and the star-formation rate history for these objects, finding good agreement with other infrared studies up to z ~1, exceeding the contribution attributed to optically-selected galaxies.

We present first results from an unbiased 50 deg^2 submillimeter Galactic survey at 250, 350, and 500 micron from the 2006 flight of the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). The map has resolution ranging from 36 arcsec to 60 arcsec in the three submillimeter bands spanning the thermal emission peak of cold starless cores. We determine the temperature, luminosity, and mass of more than one thousand compact sources in a range of evolutionary stages and an unbiased statistical characterization of the population. From comparison with C^(18)O data, we find the dust opacity per gas mass, kappa r = 0.16 cm^2 g^(-1) at 250 micron, for cold clumps. We find that 2% of the mass of the molecular gas over this diverse region is in cores colder than 14 K, and that the mass function for these cold cores is consistent with a power law with index alpha = -3.22 +/- 0.14 over the mass range 14 M_sun < M < 80 M_sun. Additionally, we infer a mass-dependent cold core lifetime of t_c(M) = 4E6 (M/20 M_sun)^(-0.9) years - longer than what has been found in previous surveys of either low or high mass cores, and significantly longer than free fall or likely turbulent decay times. This implies some form of non-thermal support for cold cores during this early stage of star formation.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 250-hour flight over Antarctica in December 2006 (BLAST06). As part of the calibration and pointing procedures, the red hypergiant star VY CMa was observed and used as the primary calibrator. Details of the overall BLAST06 calibration procedure are discussed. The 1-sigma absolute calibration is accurate to 10, 12, and 13% at the 250, 350, and 500 micron bands, respectively. The errors are highly correlated between bands resulting in much lower error for the derived shape of the 250-500 micron continuum. The overall pointing error is <5" rms for the 36, 42, and 60" beams. The performance of the optics and pointing systems is discussed.

We have identified radio and/or mid-infrared counterparts to 198 out of 350 sources detected at >=5 sigma over ~ 9 square degrees centered on the Chandra Deep Field South (CDFS) by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) at 250, 350 and 500 um. We have matched 114 of these counterparts to optical sources with previously derived photometric redshifts and fitted SEDs to the BLAST fluxes and fluxes at 70 and 160 um acquired with the Spitzer Space Telescope. In this way, we have constrained dust temperatures, total far-infrared/sub-millimeter luminosities and star formation rates for each source. Our findings show that on average, the BLAST sources lie at significantly lower redshifts and have significantly lower rest-frame dust temperatures compared to submm sources detected in surveys conducted at 850 um. We demonstrate that an apparent increase in dust temperature with redshift in our sample arises as a result of selection effects. Finally, we provide the full multi-wavelength catalog of >= 5 sigma BLAST sources contained within the complete ~ 9 square degree survey area.

Submillimetre surveys during the past decade have discovered a population of luminous, high-redshift, dusty starburst galaxies. In the redshift range 1 <= z <= 4, these massive submillimetre galaxies go through a phase characterized by optically obscured star formation at rates several hundred times that in the local Universe. Half of the starlight from this highly energetic process is absorbed and thermally re-radiated by clouds of dust at temperatures near 30 K with spectral energy distributions peaking at 100 microns in the rest frame. At 1 <= z <= 4, the peak is redshifted to wavelengths between 200 and 500 microns. The cumulative effect of these galaxies is to yield extragalactic optical and far-infrared backgrounds with approximately equal energy densities. Since the initial detection of the far-infrared background (FIRB), higher-resolution experiments have sought to decompose this integrated radiation into the contributions from individual galaxies. Here we report the results of an extragalactic survey at 250, 350 and 500 microns. Combining our results at 500 microns with those at 24 microns, we determine that all of the FIRB comes from individual galaxies, with galaxies at z >= 1.2 accounting for 70 per cent of it. As expected, at the longest wavelengths the signal is dominated by ultraluminous galaxies at z > 1.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital experiment designed to study the process of star formation in local galaxies (including the Milky Way) and in galaxies at cosmological distances. Using a 2-m Cassegrain telescope, BLAST images the sky onto a focal plane, which consists of 270 bolometric detectors split between three arrays, observing simultaneously in 30% wide bands, centered at 250, 350, and 500 microns. The diffraction-limited optical system provides a resolution of 30" at 250 microns. The pointing system enables raster-like scans with a positional accuracy of ~30", reconstructed to better than 5" rms in post-flight analysis. BLAST had two successful flights, from the Arctic in 2005, and from Antarctica in 2006, which provided the first high-resolution and large-area (~0.8-200 deg^2) submillimeter surveys at these wavelengths. As a pathfinder for the SPIRE instrument on Herschel, BLAST shares with the ESA satellite similar focal plane technology and scientific motivation. A third flight in 2009 will see the instrument modified to be polarization-sensitive (BLAST-Pol). With its unprecedented mapping speed and resolution, BLAST-Pol will provide insights into Galactic star-forming nurseries, and give the necessary link between the larger, coarse resolution surveys and the narrow, resolved observations of star-forming structures from space and ground based instruments being commissioned in the next 5 years.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 100-hour flight from northern Sweden in June 2005 (BLAST05). As part of the calibration and pointing procedures, several compact sources were mapped, including solar system, Galactic, and extragalactic targets, specifically Pallas, CRL 2688, LDN 1014, IRAS 20126+4104, IRAS 21078+5211, IRAS 21307+5049, IRAS 22134+5834, IRAS 23011+6126, K3-50, W 75N, and Mrk 231. One additional source, Arp 220, was observed and used as our primary calibrator. Details of the overall BLAST05 calibration procedure are discussed here. The BLAST observations of each compact source are described, flux densities and spectral energy distributions are reported, and these are compared with previous measurements at other wavelengths. The 250, 350, and 500 um BLAST data can provide useful constraints to the amplitude and slope of the submillimeter continuum, which in turn may be useful for the improved calibration of other submillimeter instruments.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital surveying experiment designed to study the evolutionary history and processes of star formation in local galaxies (including the Milky Way) and galaxies at cosmological distances. The BLAST continuum camera, which consists of 270 detectors distributed between 3 arrays, observes simultaneously in broad-band (30%) spectral-windows at 250, 350, and 500 microns. The optical design is based on a 2m diameter telescope, providing a diffraction-limited resolution of 30" at 250 microns. The gondola pointing system enables raster mapping of arbitrary geometry, with a repeatable positional accuracy of ~30"; post-flight pointing reconstruction to ~5" rms is achieved. The on-board telescope control software permits autonomous execution of a pre-selected set of maps, with the option of manual override. In this paper we describe the primary characteristics and measured in-flight performance of BLAST. BLAST performed a test-flight in 2003 and has since made two scientifically productive long-duration balloon flights: a 100-hour flight from ESRANGE (Kiruna), Sweden to Victoria Island, northern Canada in June 2005; and a 250-hour, circumpolar-flight from McMurdo Station, Antarctica, in December 2006.

We describe a map-making method which we have developed for the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) experiment, but which should have general application to data from other submillimeter arrays. Our method uses a Maximum Likelihood based approach, with several approximations, which allows images to be constructed using large amounts of data with fairly modest computer memory and processing requirements. This new approach, Signal And Noise Estimation Procedure Including Correlations (SANEPIC), builds upon several previous methods, but focuses specifically on the regime where there is a large number of detectors sampling the same map of the sky, and explicitly allowing for the the possibility of strong correlations between the detector timestreams. We provide real and simulated examples of how well this method performs compared with more simplistic map-makers based on filtering. We discuss two separate implementations of SANEPIC: a brute-force approach, in which the inverse pixel-pixel covariance matrix is computed; and an iterative approach, which is much more efficient for large maps. SANEPIC has been successfully used to produce maps using data from the 2005 BLAST flight.

We present the first results from a new 250, 350, and 500 micron Galactic Plane survey taken with the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) in 2005. This survey's primary goal is to identify and characterize high-mass proto-stellar objects (HMPOs). The region studied here covers 4 sq. deg near the open cluster NGC 6823 in the constellation Vulpecula (l=59). We find 60 compact sources (<60'' diameter) detected simultaneously in all three bands. Their spectral energy distributions (SEDs) are constrained through BLAST, IRAS, Spitzer MIPS, and MSX photometry, with inferred dust temperatures spanning ~12-40K assuming a dust emissivity index beta=1.5. The luminosity-to-mass ratio, a distance-independent quantity, spans ~0.2-130 L_\odot M_\odot^{-1}. Distances are estimated from coincident 13CO (1->0) velocities combined with a variety of other velocity and morphological data in the literature. In total, 49 sources are associated with a molecular cloud complex encompassing NGC 6823 (distance ~2.3kpc), 10 objects with the Perseus Arm (~8.5kpc) and one object is probably in the outer Galaxy (~14kpc). Near NGC 6823, the inferred luminosities and masses of BLAST sources span ~40-10^4 L_\odot, and ~15-700 M_\odot, respectively. The mass spectrum is compatible with molecular gas masses in other high-mass star forming regions. Several luminous sources appear to be Ultra Compact HII regions powered by early B stars. However, many of the objects are cool, massive gravitationally-bound clumps with no obvious internal radiation from a protostar, and hence excellent HMPO candidates.