We present a model to predict the luminosity function for radio galaxies and their linear size distribution at any redshift. The model takes a black hole mass function and Eddington ratio distribution as input and tracks the evolution of radio sources, taking into account synchrotron, adiabatic and inverse Compton energy losses. We first test the model at z = 2 where plenty of radio data is available and show that the radio luminosity function (RLF) is consistent with observations. We are able to reproduce the break in luminosity function that separates locally the FRI and FRII radio sources. Our prediction for linear size distribution at z = 2 matches the observed distribution too. We then use our model to predict a RLF and linear size distribution at z = 6, as this is the epoch when radio galaxies can be used as probes of reionisation. We demonstrate that higher inverse Compton losses lead to shorter source lifetimes and smaller sizes at high redshifts. The predicted sizes are consistent with the generally observed trend with redshift. We evolve the z = 2 RLF based on observed quasar space densities at high redshifts, and show that our RLF prediction at z = 6 is consistent. Finally, we predict the detection of 0.63, 0.092 and 0.0025 z>=6 sources per sq. degree at flux density limits of 0.1, 0.5 and 3.5 mJy. We assess the trade-off between coverage area and depth and show that LOFAR surveys with flux density limits of 0.1 and 0.5 mJy would are the most efficient at detecting a large number of z>=6 radio sources.

At bright radio powers ($P_{\rm 1.4 GHz} > 10^{25}$ W/Hz) the space density of the most powerful sources peaks at higher redshift than that of their weaker counterparts. This paper establishes whether this luminosity-dependent evolution persists for sources an order of magnitude fainter than those previously studied, by measuring the steep--spectrum radio luminosity function (RLF) across the range $10^{24} < P_{\rm 1.4 GHz} < 10^{28}$ W/Hz, out to high redshift. A grid-based modelling method is used, in which no assumptions are made about the RLF shape and high-redshift behaviour. The inputs to the model are the same as in Rigby et al. (2011): redshift distributions from radio source samples, together with source counts and determinations of the local luminosity function. However, to improve coverage of the radio power vs. redshift plane at the lowest radio powers, a new faint radio sample is introduced. This covers 0.8 sq. deg., in the Subaru/XMM-Newton Deep Field, to a 1.4 GHz flux density limit of $S_{\rm 1.4 GHz} \geq 100~\mu$Jy, with 99% redshift completeness. The modelling results show that the previously seen high-redshift declines in space density persist to $P_{\rm 1.4 GHz} < 10^{25}$ W/Hz. At $P_{\rm 1.4 GHz} > 10^{26}$ W/Hz the redshift of the peak space density increases with luminosity, whilst at lower radio luminosities the position of the peak remains constant within the uncertainties. This `cosmic downsizing' behaviour is found to be similar to that seen at optical wavelengths for quasars, and is interpreted as representing the transition from radiatively efficient to inefficient accretion modes in the steep-spectrum population. This conclusion is supported by constructing simple models for the space density evolution of these two different radio galaxy classes; these are able to successfully reproduce the observed variation in peak redshift.

We present results from a narrow-band survey of the field around the high redshift radio galaxy MRC 2104-242. We have selected Halpha emitters in a 7sq.arcmin field and compared the measured number density with that of a field sample at similar redshift. We find that MRC 2104-242 lies in an overdensity of galaxies that is 8.0 +/- 0.8 times the average density of a blank field, suggesting it resides in a large-scale structure that may eventually collapse to form a massive cluster. We find that there is more dust obscured star formation in the protocluster galaxies than in similarly selected control field galaxies and there is tentative evidence of a higher fraction of starbursting galaxies in the denser environment. However, on average we do not find a difference between the star formation rate (SFR)-mass relations of the protocluster and field galaxies and so conclude that the SFR of these galaxies at z~2.5 is governed predominantly by galaxy mass and not the host environment. We also find that the stellar mass distribution of the protocluster galaxies is skewed towards higher masses and there is a significant lack of galaxies at M < 10^10Msun within our small field of view. Based on the level of overdensity we expect to find ~22 star forming galaxies below 10^10Msun in the protocluster and do not detect any. This lack of low mass galaxies affects the level of overdensity which we detect. If we only consider high mass (M > 10^10.5Msun) galaxies, the density of the protocluster field increases to ~55 times the control field density.

This paper presents the first results of a far-infrared search for protocluster-associated galaxy overdensities using the SPIRE instrument on-board the {\it Herschel} Space Observatory. Large ($\sim$400 arcmin$^{2}$) fields surrounding 26 powerful high-redshift radio galaxies ($2.0 < z < 4.1$; $L_{\rm 500 MHz} > 10^{28.5}$ WHz$^{-1}$) are mapped at 250, 350 and 500\mic to give a unique wide-field sample. On average the fields have a higher than expected, compared to blank fields, surface density of 500\mic sources within 6 comoving Mpc of the radio galaxy. The analysis is then restricted to potential protocluster members only, which are identified using a far-infrared colour selection; this reveals significant overdensities of galaxies in 2 fields, neither of which are previously known protoclusters. The probability of finding 2 overdensities of this size by chance, given the number of fields observed is $5 \times 10^{-4}$. Overdensities here exist around radio galaxies with $L_{\rm 500 MHz} \gtrsim 10^{29}$ WHz$^{-1}$ and $z < 3$. The radial extent of the average far-infrared overdensity is found to be $\sim$6 comoving Mpc. Comparison with predictions from numerical simulations shows that the overdensities are consistent with having masses $> 10^{14}$Msolar. However, the large uncertainty in the redshift estimation means that it is possible that these far-infrared overdensities consist of several structures across the redshift range searched.

We present an analysis of CO molecular gas tracers in a sample of 500{\mu}m-selected Herschel-ATLAS galaxies at z<0.05 (cz<14990km/s). Using 22-500{\mu}m photometry from WISE, IRAS and Herschel, with HI data from the literature, we investigate correlations between warm and cold dust, and tracers of the gas in different phases. The correlation between global CO(3-2) line fluxes and FIR-submillimetre fluxes weakens with increasing IR wavelength ({\lambda}>60{\mu}m), as a result of colder dust being less strongly associated with dense gas. Conversely, CO(2-1) and HI line fluxes both appear to be better correlated with longer wavelengths, suggesting that cold dust is more strongly associated with diffuse atomic and molecular gas phases, consistent with it being at least partially heated by radiation from old stellar populations. The increased scatter at long wavelengths implies that submillimetre fluxes are a poorer tracer of SFR. Fluxes at 22 and 60{\mu}m are also better correlated with diffuse gas tracers than dense CO(3-2), probably due to very-small-grain emission in the diffuse interstellar medium, which is not correlated with SFR. The FIR/CO luminosity ratio and the dust mass/CO luminosity ratio both decrease with increasing luminosity, as a result of either correlations between mass and metallicity (changing CO/H2) or between CO luminosity and excitation [changing CO(3-2)/CO(1-0)].

We have constructed a sample of radio-loud objects with optical spectroscopy from the Galaxy and Mass Assembly (GAMA) project over the Herschel-ATLAS Phase 1 fields. Classifying the radio sources in terms of their optical spectra, we find that strong-emission-line sources (`high-excitation radio galaxies') have, on average, a factor ~4 higher 250-micron Herschel luminosity than weak-line (`low-excitation') radio galaxies and are also more luminous than magnitude-matched radio-quiet galaxies at the same redshift. Using all five H-ATLAS bands, we show that this difference in luminosity between the emission-line classes arises mostly from a difference in the average dust temperature; strong-emission-line sources tend to have comparable dust masses to, but higher dust temperatures than, radio galaxies with weak emission lines. We interpret this as showing that radio galaxies with strong nuclear emission lines are much more likely to be associated with star formation in their host galaxy, although there is certainly not a one-to-one relationship between star formation and strong-line AGN activity. The strong-line sources are estimated to have star-formation rates at least a factor 3-4 higher than those in the weak-line objects. Our conclusion is consistent with earlier work, generally carried out using much smaller samples, and reinforces the general picture of high-excitation radio galaxies as being located in lower-mass, less evolved host galaxies than their low-excitation counterparts.

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations.

We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the submm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into submm imaging at 250, 350 and 500{\mu}m to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250{\mu}m background (4% at 350{\mu}m; 3% at 500{\mu}m), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250{\mu}m/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Submm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years.

We present the dust properties and star-formation histories of local submillimetre-selected galaxies in Herschel-ATLAS, classified by optical morphology. The early-type galaxies (ETGs) that are detected contain as much dust as typical spirals, and form a unique sample that has been blindly selected at submillimetre wavelengths. Comparing H-ATLAS galaxies to a control sample of optically selected galaxies, we find 5.5% of luminous ETGs are detected in H-ATLAS. The H-ATLAS ETGs contain a significant mass of cold dust: the mean dust mass is 5.5x10^7 Msun, with individual galaxies ranging from 9x10^5-4x10^8 Msun. This is comparable to that of spirals in our sample, and is an order of magnitude more dust than that found for the control ETGs, which have a median dust mass inferred from stacking of (0.8-4.0)x10^6 Msun. The ETGs detected in H-ATLAS have bluer NUV-r colours, higher specific star-formation rates and younger stellar populations than ETGs which are optically selected, and may be transitioning from the blue cloud to the red sequence. We also find that H-ATLAS and control ETGs inhabit similar low-density environments. We conclude that the dust in H-ATLAS and control ETGs cannot be solely from stellar sources, and a large contribution from dust formed in the ISM or external sources is required. Alternatively, dust destruction may not be as efficient as predicted. We also explore the properties of the most passive spiral galaxies in our sample with SSFR<10^-11/yr. We find these passive spirals have lower dust-to-stellar mass ratios, higher stellar masses and older stellar population ages than normal spirals. The passive spirals inhabit low density environments similar to those of the normal spiral galaxies in our sample. This shows that the processes which turn spirals passive do not occur solely in the intermediate density environments of group and cluster outskirts. (Abridged)

We present a technique to identify optical counterparts of 250 um-selected sources from the Herschel-ATLAS survey. Of the 6621 250 um > 32 mJy sources in our science demonstration catalogue we find that ~60 percent have counterparts brighter than r=22.4 mag in the Sloan Digital Sky Survey. Applying a likelihood ratio technique we are able to identify 2423 of the counterparts with a reliability R > 0.8. This is approximately 37 percent of the full 250 micron catalogue. We have estimated photometric redshifts for each of these 2423 reliable counterparts, while 1099 also have spectroscopic redshifts collated from several different sources, including the GAMA survey. We estimate the completeness of identifying counterparts as a function of redshift, and present evidence that 250 um-selected Herschel-ATLAS galaxies have a bimodal redshift distribution. Those with reliable optical identifications have a redshift distribution peaking at z ~ 0.25 +/- 0.05, while sub-mm colours suggest that a significant fraction with no counterpart above the r-band limit have z > 1. We also suggest a method for selecting populations of strongly-lensed high redshift galaxies. Our identifications are matched to UV--NIR photometry from the GAMA survey, and these data are available as part of the Herschel-ATLAS public data release.

This paper presents a new grid-based method for investigating the evolution of the steep-spectrum radio luminosity function, with the aim of quantifying the high-redshift cut-off suggested by previous work. To achieve this, the Combined EIS-NVSS Survey of Radio Sources (CENSORS) has been developed; this is a 1.4 GHz radio survey, containing 135 sources complete to a flux density of 7.2 mJy, selected from the NRAO VLA Sky Survey (NVSS) over 6 deg^2 of the ESO Imaging Survey (EIS) Patch D. The sample is currently 73% spectroscopically complete, with the remaining redshifts estimated via the K-z or I-z magnitude-redshift relation. CENSORS is combined with additional radio data from the Parkes All-Sky, Parkes Selected Regions, Hercules and VLA COSMOS samples to provide comprehensive coverage of the radio power vs. redshift plane. The redshift distributions of these samples, together with radio source count determinations, and measurements of the local luminosity function, provide the input to the fitting process. The modelling reveals clear declines, at > 3sigma significance, in comoving density at z > 0.7 for lower luminosity sources (log P = 25-26); these turnovers are still present at log P > 27, but move to z > 3, suggesting a luminosity-dependent evolution of the redshift turnover, similar to the `cosmic downsizing' seen for other AGN populations. These results are shown to be robust to the estimated redshift errors and to increases in the spectral index for the highest redshift sources. Analytic fits to the best-fitting steep spectrum grid are provided so that the results presented here can be easily accessed by the reader, as well as allowing plausible extrapolations outside of the regions covered by the input datasets

The Herschel-ATLAS is a survey of 550 square degrees with the Herschel Space Observatory in five far--infrared and submillimetre bands. The first data for the survey, observations of a field 4x4 sq. degrees in size, were taken during the Science Demonstration Phase, and reach a 5 sigma noise level of 33 mJy/beam at 250 microns. This paper describes the source extraction methods used to create the corresponding Science Demonstration Phase catalogue, which contains 6876 sources, selected at 250 microns, within ~14 sq. degrees. SPIRE sources are extracted using a new method specifically developed for Herschel data; PACS counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 microns have undergone flux density enhancements of up to a factor of ~2, due mainly to source confusion. Correction factors are determined for these effects. The SDP dataset and corresponding catalogue will be available from http://www.h-atlas.org/.

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^{\theta} (1 + z)^{\zeta}, we find {\theta} = 0.22 +/- 0.08 and {\zeta} = 1.6 +/- 0.4, although there is substantial additional uncertainty in {\zeta} of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &\mu&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

By combining Herschel-SPIRE data with archival Spitzer, HI, and CO maps, we investigate the spatial distribution of gas and dust in the two famous grand-design spirals M99 and M100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SPIRE photometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices f350/f500 and f250/f350 decrease radially consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of both galaxies.

We present Herschel-SPIRE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, supporting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investigations. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of environmental effects on galaxy evolution.

With appropriate spatial resolution, images of spiral galaxies in thermal infrared (~10 micron and beyond) often reveal a bright central component, distinct from the stellar bulge, superimposed on a disk with prominent spiral arms. ISO and Spitzer studies have shown that much of the scatter in the mid-infrared colors of spiral galaxies is related to changes in the relative importance of these two components, rather than to other modifications, such as the morphological type or star formation rate, that affect the properties of the galaxy as a whole. With the Herschel imaging capability from 70 to 500 micron, we revisit this two-component approach at longer wavelengths, to see if it still provides a working description of the brightness distribution of galaxies, and to determine its implications on the interpretation of global far-infrared properties of galaxies.

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 {\mu}m) shows a prominent CO rotational ladder from J = 4-3 to 13-12 emitted by the central region of M82. The fundamental properties of the gas are well constrained by the high J lines observed for the first time. Radiative transfer modeling of these high-S/N 12CO and 13CO lines strongly indicates a very warm molecular gas component at ~500 K and pressure of ~3x10^6 K cm^-3, in good agreement with the H_2 rotational lines measurements from Spitzer and ISO. We suggest that this warm gas is heated by dissipation of turbulence in the interstellar medium (ISM) rather than X-rays or UV flux from the straburst. This paper illustrates the promise of the SPIRE FTS for the study of the ISM of nearby galaxies.

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh-Jeans slope of dust emission observed at 70-500 microns in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 microns are primarily dependent on radius but that the ratio of 70 to 160 micron emission shows no clear dependence on surface brightness or radius. These results along with analyses of the spectral energy distributions imply that the 160-500 micron emission traces 15-30 K dust heated by evolved stars in the bulge and disc whereas the 70 micron emission includes dust heated by the active galactic nucleus and young stars in star forming regions.

M82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind.

In this paper the cosmic evolution of the space density of Fanaroff & Riley Class I (FRI) radio sources is investigated out to z ~ 1, in order to understand the origin of the differences between these and the more powerful FRIIs. High resolution radio images are presented of the best high redshift FRI candidate galaxies, drawn from two fields of the Leiden Berkeley Deep Survey, and previously defined in Rigby, Snellen & Best (2007, Paper I). Together with lower resolution radio observations (both previously published in Paper I and, for a subset of sources, also presented here) these are used to morphologically classify the sample. Sources which are clearly resolved are classified by morphology alone, whereas barely or unresolved sources were classified using a combination of morphology and flux density loss in the higher resolution data, indicative of resolved out extended emission. The space densities of the FRIs are then calculated as a function of redshift, and compared to both measurements of the local value and the behaviour of the more powerful FRIIs. The space density of FRI radio sources with luminosities (at 1.4 GHz) > 10^25 W/Hz is enhanced by a factor of 5-9 by z ~ 1, implying moderately strong evolution of this population; this enhancement is in good agreement with models of FRII evolution at the same luminosity. There are also indications that the evolution is luminosity dependent, with the lower powered sources evolving less strongly.