Radio relics are patches of diffuse synchrotron radio emission that trace shock waves. Relics are thought to form when intra-cluster medium electrons are accelerated by cluster merger induced shock waves through the diffusive shock acceleration mechanism. In this paper, we present observations spanning 150 MHz to 30 GHz of the `Sausage' and `Toothbrush' relics from the Giant Metrewave and Westerbork telescopes, the Karl G. Jansky Very Large Array, the Effelsberg telescope, the Arcminute Microkelvin Imager and Combined Array for Research in Millimeter-wave Astronomy. We detect both relics at 30 GHz, where the previous highest frequency detection was at 16 GHz. The integrated radio spectra of both sources clearly steepen above 2 GHz, at the >6$\sigma$ significance level, supports the spectral steepening previously found in the `Sausage' and the Abell 2256 relic. Our results challenge the widely adopted simple formation mechanism of radio relics and suggest more complicated models have to be developed that, for example, involve re-acceleration of aged seed electrons.

We use a Bayesian software package to analyze CARMA-8 data towards 19 unconfirmed Planck SZ-cluster candidates from Rodriguez-Gonzalvez et al. (2015), that are associated with significant overdensities in WISE. We used two cluster parameterizations, one based on a (fixed shape) generalized-NFW pressure profile and another based on a beta-gas-density profile (with varying shape parameters) to obtain parameter estimates for the nine CARMA-8 SZ-detected clusters. We find our sample is comprised of massive, Y_{500}=0.0010 \pm 0.0015 arcmin^2, relatively compact, theta_{500}= 3.9 \pm 2.0 arcmin systems. Results from the beta model show that our cluster candidates exhibit a heterogeneous set of brightness-temperature profiles. Comparison of Planck and CARMA-8 measurements showed good agreement in Y_{500} and an absence of obvious biases. We estimated the total cluster mass M_{500} as a function of z for one of the systems; at the preferred photometric redshift of 0.5, the derived mass, M_{500} \approx 0.8 \pm 0.2 \times 10^{15} Msun. Spectroscopic Keck/MOSFIRE data confirmed a galaxy member of one of our cluster candidates to be at z=0.565. Applying a Planck prior in Y_{500} to the CARMA-8 results reduces uncertainties for both parameters by a factor >4, relative to the independent Planck or CARMA-8 measurements. We here demonstrate a powerful technique to find massive clusters at intermediate z \gtrsim 0.5 redshifts using a cross-correlation between Planck and WISE data, with high-resolution follow-up with CARMA-8. We also use the combined capabilities of Planck and CARMA-8 to obtain a dramatic reduction by a factor of several, in parameter uncertainties.

We present 1-2 arcmin spatial resolution CARMA-8 31-GHz observations towards 19 unconfirmed Planck cluster candidates, selected to have significant galaxy overdensities from the WISE early data release and thought to be at z>1 from the WISE colors of the putative brightest cluster galaxy (BCG). We find a Sunyaev-Zeldovich (SZ) detection in the CARMA-8 data towards 9 candidate clusters, where one detection is considered tentative. For each cluster candidate we present CARMA-8 maps, a study of their radio-source environment and we assess the reliability of the SZ detection. The CARMA SZ detections appear to be SZ-bright, with the mean, primary-beam-corrected peak flux density of the decrement being -2.9mJy/beam with a standard deviation of 0.8, and are typically offset from the Planck position by approximately 80 arcsec. Using archival imaging data in the vicinity of the CARMA SZ centroids, we present evidence that one cluster matches Abell 586-a known z~0.2 cluster; four candidate clusters are likely to have 0.3<z<0.7; and, for the remaining 4, the redshift information is inconclusive. We also argue that the sensitivity limits resulting from the cross-correlation between Planck and WISE makes it challenging to use our selection criterion to identify clusters at z > 1.

We have obtained deep SZ observations towards 15 of the apparently hottest XMM Cluster Survey (XCS) clusters that can be observed with the Arcminute Microkelvin Imager (AMI). We use a Bayesian analysis to quantify the significance of our SZ detections. We detect the SZ effect at high significance towards three of the clusters and at lower significance for a further two clusters. Towards the remaining ten clusters, no clear SZ signal was measured. We derive cluster parameters using the XCS mass estimates as a prior in our Bayesian analysis. For all AMI-detected clusters, we calculate large-scale mass and temperature estimates while for all undetected clusters we determine upper limits on these parameters. We find that the large- scale mean temperatures derived from our AMI SZ measurements (and the upper limits from null detections) are substantially lower than the XCS-based core-temperature estimates. For clusters detected in the SZ, the mean temperature is, on average, a factor of 1.4 lower than temperatures from the XCS. For clusters undetected in SZ, the average 68% upper limit on the mean temperature is a factor of 1.9 below the XCS temperature.

We present CARMA observations of a massive galaxy cluster discovered in the AMI blind SZ survey. Without knowledge of the cluster redshift a Bayesian analysis of the AMI, CARMA and joint AMI & CARMA uv-data is used to quantify the detection significance and parameterise both the physical and observational properties of the cluster whilst accounting for the statistics of primary CMB anisotropies, receiver noise and radio sources. The joint analysis of the AMI & CARMA uv-data was performed with two parametric physical cluster models: the {\beta}-model; and the model described in Olamaie et al. 2012 with the pressure profile fixed according to Arnaud et al. 2010. The cluster mass derived from these different models is comparable but our Bayesian evidences indicate a preference for the {\beta}-profile which we, therefore, use throughout our analysis. From the CARMA data alone we obtain a Bayesian probability of detection ratio of 12.8:1 when assuming that a cluster exists within our search area; alternatively assuming that Jenkins et al. 2001 accurately predicts the number of clusters as a function of mass and redshift, the Bayesian probability of detection is 0.29:1. From the analysis of the AMI or AMI & CARMA data the probability of detection ratio exceeds 4.5x10^3:1. Performing a joint analysis of the AMI & CARMA data with a physical cluster model we derive the total mass internal to r200 as MT,200 = 4.1x10^14Msun. Using a phenomenological {\beta}-model to quantify the temperature decrement as a function of angular distance we find a central SZ temperature decrement of 170{\mu}K in the AMI & CARMA data. The SZ decrement in the CARMA data is weaker than expected and we speculate that this is a consequence of the cluster morphology. In a forthcoming study we will assess the impact of cluster morphology on the SZ decrements that are observed with interferometers such as AMI and CARMA.

We present 16-GHz Sunyaev-Zel'dovich observations using the Arcminute Microkelvin Imager (AMI) and subsequent Bayesian analysis of six galaxy clusters at redshift ($z \approx 1$) chosen from an X-ray and Infrared selected sample from Culverhouse et al. (2010). In the subsequent analysis we use two cluster models, an isothermal \beta-model and a Dark Matter GNFW (DM-GNFW) model in order to derive a formal detection probability and the cluster parameters. We detect two clusters (CLJ1415+3612 & XMJ0830+5241) and measure their total masses out to a radius of 200 $\times$ the critical density at the respective cluster's redshift. For CLJ1415+3612 and XMJ0830+5241, we find M_{\mathrm{T},200} for each model, which agree with each other for each cluster. We also present maps before and after source subtraction of the entire sample and provide 1D and 2D posterior marginalised probability distributions for each fitted cluster profile parameter of the detected clusters. Using simulations which take into account the measured source environment from the AMI Large Array (LA), source confusion noise, CMB primordials, instrument noise, we estimate from low-radius X-ray data from Culverhouse et al. (2010), the detectability of each cluster in the sample and compare it with the result from the Small Array (SA) data. Furthermore, we discuss the validity of the assumptions of isothermality and constant gas mass fraction. We comment on the bias that these small-radius estimates introduce to large-radius SZ predictions. In addition, we follow-up the two detections with deep, single-pointed LA observations. We find a 3 sigma tentative decrement toward CLJ1415+3612 at high-resolution and a 5 sigma high-resolution decrement towards XMJ0830+5241.

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

We present 16-GHz observations using the Arcminute Microkelvin Imager (AMI) of 11 clusters with 7 x 10^{37}W < L_X < 11 x 10^{37}W (h_{50}=1.0) selected from the Local Cluster Substructure Survey (LoCuSS) and compare them to X-ray data. We use a fast, Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model and obtain robust cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy. Our analysis fits a spherical, isothermal beta-model to our data and assumes the cluster follows the theoretical mass-temperature relation. Large-scale cluster parameters internal to r_{500} are derived under the assumption of hydrostatic equilibrium. Posterior distributions for the large-scale parameters of 8 of our clusters are given; SZ effects towards Abell 1704 and Zw0857.9+2107 were not detected and our spherical beta-profile was found to be an inadequate fit to the decrement on our map for Abell 2409.

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our target sensitivity of 100{\mu}Jy/beam. In follow-up deep pointed observations the SZ effect is detected with a maximum peak decrement greater than 8 \times the thermal noise. No corresponding emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palomar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent is large (\approx 10') and complex; our analysis favours a model containing two clusters rather than a single cluster. Our Bayesian analysis is currently limited to modelling each cluster with an ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellipsoid to the deeper candidate we find the following. (a) Assuming that the Evrard et al. (2002) approximation to Press & Schechter (1974) correctly gives the number density of clusters as a function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of the AMI detection of this cluster is 7.9 \times 10^4:1; alternatively assuming Jenkins et al. (2001) as the true prior, the formal Bayesian probability ratio of detection is 2.1 \times 10^5:1. (b) The cluster mass is MT,200 = 5.5+1.2\times 10^14h-1M\odot. (c) Abandoning a physical model with num- -1.3 70 ber density prior and instead simply modelling the SZ decrement using a phenomenological {\beta}-model of temperature decrement as a function of angular distance, we find a central SZ temperature decrement of -295+36 {\mu}K - this allows for CMB primary anisotropies, receiver -15 noise and radio sources. We are unsure if the cluster system we observe is a merging system or two separate clusters.

We present 16 GHz (1.9 cm) deep radio continuum observations made with the Arcminute Microkelvin Imager (AMI) of a sample of low-mass young stars driving jets. We combine these new data with archival information from an extensive literature search to examine spectral energy distributions (SEDs) for each source and calculate both the radio and sub-mm spectral indices in two different scenarios: (1) fixing the dust temperature (Td) according to evolutionary class; (2) allowing Td to vary. We use the results of this analysis to place constraints on the physical mechanisms responsible for the radio emission. From AMI data alone, as well as from model fitting to the full SED in both scenarios, we find that 80 per cent of the objects in this sample have spectral indices consistent with free-free emission. We find an average spectral index in both Td scenarios consistent with free-free emission. We examine correlations of the radio luminosity with bolometric luminosity, envelope mass, and outflow force and find that these data are consistent with the strong correlation with envelope mass seen in lower luminosity samples. We examine the errors associated with determining the radio luminosity and find that the dominant source of error is the uncertainty on the opacity index, beta. We examine the SEDs for variability in these young objects, and find evidence for possible radio flare events in the histories of L1551 IRS 5 and Serpens SMM 1.

Most Sunyaev--Zel'dovich (SZ) and X-ray analyses of galaxy clusters try to constrain the cluster total mass and/or gas mass using parameterised models and assumptions of spherical symmetry and hydrostatic equilibrium. By numerically exploring the probability distributions of the cluster parameters given the simulated interferometric SZ data in the context of Bayesian methods, and assuming a beta-model for the electron number density we investigate the capability of this model and analysis to return the simulated cluster input quantities via three rameterisations. In parameterisation I we assume that the T is an input parameter. We find that parameterisation I can hardly constrain the cluster parameters. We then investigate parameterisations II and III in which fg(r200) replaces temperature as a main variable. In parameterisation II we relate M_T(r200) and T assuming hydrostatic equilibrium. We find that parameterisation II can constrain the cluster physical parameters but the temperature estimate is biased low. In parameterisation III, the virial theorem replaces the hydrostatic equilibrium assumption. We find that parameterisation III results in unbiased estimates of the cluster properties. We generate a second simulated cluster using a generalised NFW (GNFW) pressure profile and analyse it with an entropy based model to take into account the temperature gradient in our analysis and improve the cluster gas density distribution. This model also constrains the cluster physical parameters and the results show a radial decline in the gas temperature as expected. The mean cluster total mass estimates are also within 1 sigma from the simulated cluster true values. However, we find that for at least interferometric SZ analysis in practice at the present time, there is no differences in the AMI visibilities between the two models. This may of course change as the instruments improve.

We present deep 1.8 cm (16 GHz) radio continuum imaging of seven young stellar objects in the Taurus molecular cloud. These objects have previously been extensively studied in the sub-mm to NIR range and their SEDs modelled to provide reliable physical and geometrical parametres.We use this new data to constrain the properties of the long-wavelength tail of the greybody spectrum, which is expected to be dominated by emission from large dust grains in the protostellar disk. We find spectra consistent with the opacity indices expected for such a population, with an average opacity index of beta = 0.26+/-0.22 indicating grain growth within the disks. We use spectra fitted jointly to radio and sub-mm data to separate the contributions from thermal dust and radio emission at 1.8 cm and derive disk masses directly from the cm-wave dust contribution. We find that disk masses derived from these flux densities under assumptions consistent with the literature are systematically higher than those calculated from sub-mm data, and meet the criteria for giant planet formation in a number of cases.

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Serpens molecular cloud by the Spitzer c2d programme at a wavelength of 1.8cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). These observations have a resolution of ~30arcsec and an average sensitivity of 19microJy/beam. The targets are predominantly Class I sources, and we find the detection rate for Class I objects in this sample to be low (18%) compared to that of Class 0 objects (67%), consistent with previous works. For detected objects we examine correlations of radio luminosity with bolometric luminosity and envelope mass and find that these data support correlations found by previous samples, but do not show any indiction of the evolutionary divide hinted at by similar data from the Perseus molecular cloud when comparing radio luminosity with envelope mass. We conclude that envelope mass provides a better indicator for radio luminosity than bolometric luminosity, based on the distribution of deviations from the two correlations. Combining these new data with archival 3.6cm flux densities we also examine the spectral indices of these objects and find an average spectral index of 0.53+/-1.14, consistent with the canonical value for a partially optically thick spherical or collimated stellar wind. However, we caution that possible inter-epoch variability limits the usefulness of this value, and such variability is supported by our identification of a possible flare in the radio history of Serpens SMM 1.

We present follow-up observations of two galaxy clusters detected blindly via the Sunyaev-Zel'dovich (SZ) effect and released in the Planck Early Release Compact Source Catalogue. We use the Arcminute Microkelvin Imager, a dual-array 14-18 GHz radio interferometer. After radio source subtraction, we find a SZ decrement of integrated flux density -1.08+/-0.10 mJy toward PLCKESZ G121.11+57.01, and improve the position measurement of the cluster, finding the centre to be RA 12 59 36.4, Dec +60 04 46.8, to an accuracy of 20 arcseconds. The region of PLCKESZ G115.71+17.52 contains strong extended emission, so we are unable to confirm the presence of this cluster via the SZ effect.

We have observed an area of approximatley 27 deg^2 to an rms noise level of less than 0.2 mJy at 15.7 GHz, using the Arcminute Microkelvin Imager Large Array. These observations constitute the most sensitive radio-source survey of any extent (greater than approximately 0.2 deg^2) above 1.4 GHz. This paper presents the techniques employed for observing, mapping and source extraction. We have used a systematic procedure for extracting information and producing source catalogues, from maps with varying noise and uv-coverage. We have performed simulations to test our mapping and source-extraction procedures, and developed methods for identifying extended, overlapping and spurious sources in noisy images. In an accompanying paper, AMI Consortium: Davies et al. 2010, the first results from the 10C survey, including the deep 15.7-GHz source count, are presented.

The first results from the Tenth Cambridge (10C) Survey of Radio Sources, carried out using the AMI Large Array (LA) at an observing frequency of 15.7 GHz, are presented. The survey fields cover an area of approximately 27 sq. degrees to a flux-density completeness of 1 mJy. Results for some deeper areas, covering approximately 12 sq. degrees, wholly contained within the total areas and complete to 0.5 mJy, are also presented. The completeness for both areas is estimated to be at least 93 per cent. The source catalogue contains 1897 entries and is available at www.mrao.cam.ac.uk/surveys/10C. It has been combined with that of the 9C Survey to calculate the 15.7-GHz source counts. A broken power law is found to provide a good parameterisation of the differential count between 0.5 mJy and 1 Jy. The measured count has been compared to that predicted by de Zotti et al. (2005). The model displays good agreement with the data at the highest flux densities but under-predicts the integrated count between 0.5 mJy and 1 Jy by about 30 per cent. Entries from the source catalogue have been matched to those contained in the catalogues of NVSS and FIRST (both of which have observing frequencies of 1.4 GHz). This matching provides evidence for a shift in the typical 1.4-to-15.7-GHz spectral index of the 15.7-GHz-selected source population with decreasing flux density towards sub-mJy levels - the spectra tend to become less steep. Automated methods for detecting extended sources have been applied to the data; approximately 5 per cent of the sources are found to be extended relative to the LA synthesised beam of approximately 30 arcsec. Investigations using higher-resolution data showed that most of the genuinely extended sources at 16 GHz are classical doubles, although some nearby galaxies and twin-jet sources were also identified.

We present 13.9-18.2 GHz observations of the Sunyaev-Zel'dovich (SZ) effect towards Abell 2146 using the Arcminute Microkelvin Imager (AMI). The cluster is detected with a peak SNR ratio of 13 sigma in the radio source subtracted map. Comparison of the SZ and X-ray images suggests that they both have extended regions which lie approximately perpendicular to one another, with their emission peaks significantly displaced. These features indicate non-uniformities in the distributions of the gas temperature and pressure, indicative of a cluster merger. We use a Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model to obtain cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy; the probability of SZ + CMB primordial structure + radio sources + receiver noise to CMB + radio sources + receiver noise is 3 x 10^{6}:1. We compare the results from three different cluster models. Our preferred model exploits the observation that the gas fractions do not appear to vary greatly between clusters. Given the relative masses of the two merging systems in Abell 2146, the mean gas temperature can be deduced from the virial theorem (assuming all of the kinetic energy is in the form of internal gas energy) without being affected significantly by the merger event, provided the primary cluster was virialized before the merger. In this model we fit a simple spherical isothermal beta-model, despite the inadequacy of this model for a merging system like Abell 2146, and assume the cluster follows the mass-temperature relation of a virialized, singular, isothermal sphere. We note that this model avoids inferring large-scale cluster parameters internal to r_200 under the widely used assumption of hydrostatic equilibrium. We find that at r_200 M_T= 4.1 \pm 0.5 x 10^{14} h^{-1}M_sun and T=4.5 \pm 0.5 keV.

We present observations from the Small Array of the Arcminute Microkelvin Imager (AMI) of eight high X-ray luminosity galaxy cluster systems selected from the Local Cluster Substructure Survey (LoCuSS) sample.We detect the Sunyaev-Zel'dovich (SZ) effect in seven of these clusters. With the assumptions that galaxy clusters are isothermal, have a density profile described by a spherical b -model and obey the theoretical M-T relation, we are able to derive cluster parameters at r200 from our SZ data. With the additional assumption of hydrostatic equilibrium we are able to derive parameters at r500. We present posterior probability distributions for cluster parameters such as mass, radius and temperature (TSZ, MT). Combining our sample with that of AMI Consortium: Rodr'\iguez-Gonz\'alvez et al. (2011) and using large-radius X-ray temperature estimates (TX) from Chandra and Suzaku observations, we find that there is reasonable correspondence between TX and TSZ,MT values at low TX, but that for clusters with TX above around 6keV the correspondence breaks down with TX exceeding TSZ, MT; we stress that this finding is based on just ten clusters.

We present observations using the Small Array of the Arcminute Microkelvin Imager (AMI; 14-18 GHz) of four Abell and three MACS clusters spanning 0.171-0.686 in redshift. We detect Sunyaev-Zel'dovich (SZ) signals in five of these without any attempt at source subtraction, although strong source contamination is present. With radio-source measurements from high-resolution observations, and under the assumptions of spherical $\beta$-model, isothermality and hydrostatic equilibrium, a Bayesian analysis of the data in the visibility plane detects extended SZ decrements in all seven clusters over and above receiver noise, radio sources and primary CMB imprints. Bayesian evidence ratios range from 10^{11}:1 to 10^{43}:1 for six of the clusters and 3000:1 for one with substantially less data than the others. We present posterior probability distributions for, e.g., total mass and gas fraction averaged over radii internal to which the mean overdensity is 1000, 500 and 200, r_200 being the virial radius. Reaching r_200 involves some extrapolation for the nearer clusters but not for the more-distant ones. We find that our estimates of gas fraction are low (compared with most in the literature) and decrease with increasing radius. These results appear to be consistent with the notion that gas temperature in fact falls with distance (away from near the cluster centre) out to the virial radius.

We report new cm-wave measurements at five frequencies between 15 and 18GHz of the continuum emission from the reportedly anomalous "region 4" of the nearby galaxy NGC6946. We find that the emission in this frequency range is significantly in excess of that measured at 8.5GHz, but has a spectrum from 15-18GHz consistent with optically thin free-free emission from a compact HII region. In combination with previously published data we fit four emission models containing different continuum components using the Bayesian spectrum analysis package radiospec. These fits show that, in combination with data at other frequencies, a model with a spinning dust component is slightly preferred to those that possess better-established emission mechanisms.

We present 25 arcsecond resolution radio images of five Lynds Dark Nebulae (L675, L944, L1103, L1111 & L1246) at 16 GHz made with the Arcminute Microkelvin Imager (AMI) Large Array. These objects were previously observed with the AMI Small Array to have an excess of emission at microwave frequencies relative to lower frequency radio data. In L675 we find a flat spectrum compact radio counterpart to the 850 micron emission seen with SCUBA and suggest that it is cm-wave emission from a previously unknown deeply embedded young protostar. In the case of L1246 the cm-wave emission is spatially correlated with 8 micron emission seen with Spitzer. Since the MIR emission is present only in Spitzer band 4 we suggest that it arises from a population of PAH molecules, which also give rise to the cm-wave emission through spinning dust emission.

Observations at 14.2 to 17.9 GHz made with the AMI Small Array towards fourteen Lynds Dark Nebulae with a resolution of 2' are reported. These sources are selected from the SCUBA observations of Visser et al. (2001) as small angular diameter clouds well matched to the synthesized beam of the AMI Small Array. Comparison of the AMI observations with radio observations at lower frequencies with matched uv-plane coverage is made, in order to search for any anomalous excess emission which can be attributed to spinning dust. Possible emission from spinning dust is identified as a source within a 2' radius of the Scuba position of the Lynds dark nebula, exhibiting an excess with respect to lower frequency radio emission. We find five sources which show a possible spinning dust component in their spectra. These sources have rising spectral indices in the frequency range 14.2--17.9 GHz. Of these five one has already been reported, L1111, we report one new definite detection, L675, and three new probable detections (L944, L1103 and L1246). The relative certainty of these detections is assessed on the basis of three criteria: the extent of the emission, the coincidence of the emission with the Scuba position and the likelihood of alternative explanations for the excess. Extended microwave emission makes the likelihood of the anomalous emission arising as a consequence of a radio counterpart to a protostar or a proto-planetary disk unlikely. We use a 2' radius in order to be consistent with the IRAS identifications of dark nebulae (Parker 1988), and our third criterion is used in the case of L1103 where a high flux density at 850 microns relative to the FIR data suggests a more complicated emission spectrum.

Using the Arcminute Microkelvin Imager (AMI) at 16 GHz and the Very Small Array (VSA) at 33 GHz to make follow-up observations of sources in the New Extragalactic WMAP Point Source Catalogue, we have investigated the flux-density variability in a complete sample of 97 sources over timescales of a few months to approximately 1.5 years. We find that 53 per cent of the 93 sources, for which we have multiple observations, are variable, at the 99 per cent confidence level, above the flux density calibration uncertainties of approximately 4 per cent at 16 GHz; the fraction of sources having varied by more than 20 per cent is 15 per cent at 16 GHz and 20 per cent at 33 GHz. Not only is this common occurrence of variability at high frequency of interest for source physics, but strategies for coping with source contamination in CMB work must take this variability into account. There is no strong evidence of a correlation between variability and flux density for the sample as a whole. Using a maximum-likelihood method, we calculate the correlation in the variability at the two frequencies in a subset of sources classified as variable from both the AMI and VSA data and find the Pearson product-moment correlation coefficient to be very high (0.955 +/- 0.034). We also find the degree of variability at 16 GHz (0.202 +/- 0.028) to be very similar to that at 33 GHz (0.224 +/- 0.039). Finally, we have investigated the relationship between variability and spectral index, alpha_13.9^33.75 (where S proportional to nu^-alpha), and find a significant difference in the spectral indices of the variable sources (-0.06 +/- 0.05) and non-variable sources (0.13 +/- 0.04).

We present follow-up observations of 97 point sources from the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data, contained within the New Extragalactic WMAP Point Source (NEWPS) catalogue between declinations of -4 and +60 degrees; the sources form a flux-density-limited sample complete to 1.1 Jy (approximately 5 sigma) at 33 GHz. Our observations were made at 16 GHz using the Arcminute Microkelvin Imager (AMI) and at 33 GHz with the Very Small Array (VSA). 94 of the sources have reliable, simultaneous -- typically a few minutes apart -- observations with both telescopes. The spectra between 13.9 and 33.75 GHz are very different from those of bright sources at low frequency: 44 per cent have rising spectra (alpha < 0.0), where flux density is proportional to frequency^-alpha, and 93 per cent have spectra with alpha < 0.5; the median spectral index is 0.04. For the brighter sources, the agreement between VSA and WMAP 33-GHz flux densities averaged over sources is very good. However, for the fainter sources, the VSA tends to measure lower values for the flux densities than WMAP. We suggest that the main cause of this effect is Eddington bias arising from variability.