We present the results of sixteen Swift-triggered GRB follow-up observations taken with the VERITAS telescope array from January, 2007 to June, 2009. The median energy threshold and response time of these observations was 260 GeV and 320 s, respectively. Observations had an average duration of 90 minutes. Each burst is analyzed independently in two modes: over the whole duration of the observations and again over a shorter time scale determined by the maximum VERITAS sensitivity to a burst with a t^-1.5 time profile. This temporal model is characteristic of GRB afterglows with high-energy, long-lived emission that have been detected by the Large Area Telescope (LAT) on-board the Fermi satellite. No significant VHE gamma-ray emission was detected and upper limits above the VERITAS threshold energy are calculated. The VERITAS upper limits are corrected for gamma-ray extinction by the extragalactic background light (EBL) and interpreted in the context of the keV emission detected by Swift. For some bursts the VHE emission must have less power than the keV emission, placing constraints on inverse Compton models of VHE emission.

We describe a dedicated cosmic-ray telescope that explores a new method for detecting Cerenkov radiation from high-energy primary cosmic rays and the large particle air shower they induce upon entering the atmosphere. Using a camera comprising 16 multi-anode photomultiplier tubes for a total of 256 pixels, the Track Imaging Cerenkov Experiment (TrICE) resolves substructures in particle air showers with 0.086 degree resolution. Cerenkov radiation is imaged using a novel two-part optical system in which a Fresnel lens provides a wide-field optical trigger and a mirror system collects delayed light with four times the magnification. TrICE records well-resolved cosmic-ray air showers at rates ranging between 0.01-0.1 Hz.

Multiwavelength observations of the high-frequency-peaked blazar 1ES2344+514 were performed from 2007 October to 2008 January. The campaign represents the first contemporaneous data on the object at very high energy (VHE, E >100 GeV) {\gamma}-ray, X-ray, and UV energies. Observations with VERITAS in VHE {\gamma}-rays yield a strong detection of 20 {\sigma} with 633 excess events in a total exposure of 18.1 hours live-time. A strong VHE {\gamma}-ray flare on 2007 December 7 is measured at F(>300 GeV) = (6.76 \pm 0.62) \times 10-11 ph cm-2 s-1, corresponding to 48% of the Crab Nebula flux. Excluding this flaring episode, nightly variability at lower fluxes is observed with a time-averaged mean of F(>300 GeV) = (1.06 \pm 0.09) \times 10-11 ph cm-2 s-1 (7.6% of the Crab Nebula flux). The differential photon spectrum between 390 GeV and 8.3 TeV for the time-averaged observations excluding 2007 December 7 is well described by a power law with a photon index of {\Gamma} = 2.78 \pm 0.09stat \pm 0.15syst. Over the full period of VERITAS observations contemporaneous X-ray and UV data were taken with Swift and RXTE. The measured 2-10 keV flux ranged by a factor of ~7 during the campaign. On 2007 December 8 the highest ever observed X-ray flux from 1ES 2344+514 was measured by Swift XRT at a flux of F(2-10 keV) = (6.28 \pm 0.31) \times 10-11 erg cm-2 s-1. Evidence for a correlation between the X-ray flux and VHE {\gamma}-ray flux on nightly time-scales is indicated with a Pearson correlation coefficient of r = 0.60 \pm 0.11. Contemporaneous spectral energy distributions (SEDs) of 1ES 2344+514 are presented for two distinct flux states. A one-zone synchrotron self-Compton (SSC) model describes both SEDs using parameters consistent with previous SSC modeling of 1ES 2344+514 from non-contemporaneous observations.

We report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

Cosmic-ray proton and helium spectra have been measured with the balloon-borne Cosmic Ray Energetics And Mass experiment flown for 42 days in Antarctica in the 2004-2005 austral summer season. High-energy cosmic-ray data were collected at an average altitude of ~38.5 km with an average atmospheric overburden of ~3.9 g cm$^{-2}$. Individual elements are clearly separated with a charge resolution of ~0.15 e (in charge units) and ~0.2 e for protons and helium nuclei, respectively. The measured spectra at the top of the atmosphere are represented by power laws with a spectral index of -2.66 $\pm$ 0.02 for protons from 2.5 TeV to 250 TeV and -2.58 $\pm$ 0.02 for helium nuclei from 630 GeV/nucleon to 63 TeV/nucleon. They are harder than previous measurements at a few tens of GeV/nucleon. The helium flux is higher than that expected from the extrapolation of the power law fitted to the lower-energy data. The relative abundance of protons to helium nuclei is 9.1 $\pm$ 0.5 for the range from 2.5 TeV/nucleon to 63 TeV/nucleon. This ratio is considerably smaller than the previous measurements at a few tens of GeV/nucleon.

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including {\gamma}-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

The intermediate-frequency peaked BL Lacertae (IBL) object 3C 66A is detected during 2007 - 2008 in VHE (very high energy: E > 100 GeV) gamma-rays with the VERITAS stereoscopic array of imaging atmospheric Cherenkov telescopes. An excess of 1791 events is detected, corresponding to a significance of 21.2 standard deviations (sigma), in these observations (32.8 hours live time). The observed integral flux above 200 GeV is 6% of the Crab Nebula's flux and shows evidence for variability on the time-scale of days. The measured energy spectrum is characterized by a soft power law with photon index Gamma = 4.1 +- 0.4_stat +- 0.6_sys. The radio galaxy 3C 66B is excluded as a possible source of the VHE emission.

Indirect dark matter searches with ground-based gamma-ray observatories provide an alternative for identifying the particle nature of dark matter that is complementary to that of direct search or accelerator production experiments. We present the results of observations of the dwarf spheroidal galaxies Draco, Ursa Minor, Bootes 1, and Willman 1 conducted by VERITAS. These galaxies are nearby dark matter dominated objects located at a typical distance of several tens of kiloparsecs for which there are good measurements of the dark matter density profile from stellar velocity measurements. Since the conventional astrophysical background of very high energy gamma rays from these objects appears to be negligible, they are good targets to search for the secondary gamma-ray photons produced by interacting or decaying dark matter particles. No significant gamma-ray flux above 200 GeV was detected from these four dwarf galaxies for a typical exposure of ~20 hours. The 95% confidence upper limits on the integral gamma-ray flux are in the range 0.4-2.2x10^-12 photons cm^-2s^-1. We interpret this limiting flux in the context of pair annihilation of weakly interacting massive particles and derive constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the WIMPs. The limits are obtained under conservative assumptions regarding the dark matter distribution in dwarf galaxies and are approximately three orders of magnitude above the generic theoretical prediction for WIMPs in the minimal supersymmetric standard model framework. However significant uncertainty exists in the dark matter distribution as well as the neutralino cross sections which under favorable assumptions could further lower the limits.

We report the discovery of very high energy gamma-ray emission from the direction of the SNR G54.1+0.3 using the VERITAS ground-based gamma-ray observatory. The TeV signal has an overall significance of 6.8$\sigma$ and appears point-like given the 5$^{arcminute}$ resolution of the instrument. The integral flux above 1 TeV is 2.5% of the Crab Nebula flux and significant emission is measured between 250 GeV and 4 TeV, well described by a power-law energy spectrum dN/dE $\sim$ E$^{-\Gamma}$ with a photon index $\Gamma= 2.39\pm0.23_{stat}\pm0.30_{sys}$. We find no evidence of time variability among observations spanning almost two years. Based on the location, the morphology, the measured spectrum, the lack of variability and a comparison with similar systems previously detected in the TeV band, the most likely counterpart of this new VHE gamma-ray source is the PWN in the SNR G54.1+0.3. The measured X-ray to VHE gamma-ray luminosity ratio is the lowest among all the nebulae supposedly driven by young rotation-powered pulsars, which could indicate a particle-dominated PWN.

M87 is a nearby radio galaxy that is detected at energies ranging from radio to VHE gamma-rays. Its proximity and its jet, misaligned from our line-of-sight, enable detailed morphological studies and extensive modeling at radio, optical, and X-ray energies. Flaring activity was observed at all energies, and multi-wavelength correlations would help clarify the origin of the VHE emission. In this paper, we describe a detailed temporal and spectral analysis of the VERITAS VHE gamma-ray observations of M87 in 2008 and 2009. In the 2008 observing season, VERITAS detected an excess with a statistical significance of 7.2 sigma from M87 during a joint multi-wavelength monitoring campaign conducted by three major VHE experiments along with the Chandra X-ray Observatory. In February 2008, VERITAS observed a VHE flare from M87 occurring over a 4-day timespan. The peak nightly flux above 250GeV was 7.7% of the Crab Nebula flux. M87 was marginally detected before this 4-day flare period, and was not detected afterwards. Spectral analysis of the VERITAS observations showed no significant change in the photon index between the flare and pre-flare states. Shortly after the VHE flare seen by VERITAS, the Chandra X-ray Observatory detected the flux from the core of M87 at a historical maximum, while the flux from the nearby knot HST-1 remained quiescent. Acciari et al. (2009) presented the 2008 contemporaneous VHE gamma-ray, Chandra X-ray, and VLBA radio observations which suggest the core as the most likely source of VHE emission, in contrast to the 2005 VHE flare that was simultaneous with an X-ray flare in the HST-1 knot. In 2009, VERITAS continued its monitoring of M87 and marginally detected a 4.2 sigma excess corresponding to a flux of ~1% of the Crab Nebula. No VHE flaring activity was observed in 2009.

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

The Cosmic Ray Energetics And Mass (CREAM) is a balloon-borne experiment designed to measure the composition and energy spectra of cosmic rays of charge Z = 1 to 26 up to an energy of ~ 10^15 eV. CREAM had two successful flights on long-duration balloons (LDB) launched from Mc- Murdo Station, Antarctica, in December 2004 and December 2005. CREAM-I achieves a substantial measurement redundancy by employing multiple detector systems, namely a Timing Charge Detector and a Silicon Charge Detector (SCD) for particle identification, and a Transition Radiation Detector and a sampling tungsten/scintillating-fiber ionization calorimeter (CAL) for energy measurement. In this paper, preliminary energy spectra of various elements measured with CAL/SCD during the first 42-day flight are presented.

The Cosmic Ray Energetics And Mass (CREAM) balloon experiment had two successful flights in 2004/05 and 2005/06. It was designed to perform energy measurements from a few GeV up to 1000 TeV, taking advantage of different detection techniques. The first instrument, CREAM-1, combined a transition radiation detector with a calorimeter to provide independent energy measurements of cosmicraynuclei. Each detector was calibrated with particle beams in a limited range of energies. In order to assess the absolute energy scale of the instrument and to investigate the systematic effects of each technique, a cross-calibration was performed by comparing the two independent energy estimates on selected samples of oxygen and carbon nuclei.

We report on observations of very high-energy gamma rays from the shell-type supernova remnant Cassiopeia A with the VERITAS stereoscopic array of four imaging atmospheric Cherenkov telescopes in Arizona. The total exposure time for these observations is 22 hours, accumulated between September and November of 2007. The gamma-ray source associated with the SNR Cassiopeia A was detected above 200 GeV with a statistical significance of 8.3 s.d. The estimated integral flux for this gamma-ray source is about 3% of the Crab-Nebula flux. The photon spectrum is compatible with a power law dN/dE ~ E^(-Gamma) with an index Gamma = 2.61 +/- 0.24(stat) +/- 0.2(sys). The data are consistent with a point-like source. We provide a detailed description of the analysis results, and discuss physical mechanisms that may be responsible for the observed gamma-ray emission.

We present results from an intensive VERITAS monitoring campaign of the high-frequency peaked BL Lac object 1ES 1218+304 in 2008/2009. Although 1ES 1218+304 was detected previously by MAGIC and VERITAS at a persistent level of ~6% of the Crab Nebula flux, the new VERITAS data reveal a prominent flare reaching ~20% of the Crab. While very high energy (VHE) flares are quite common in many nearby blazars, the case of 1ES 1218+304 (redshift z = 0.182) is particularly interesting since it belongs to a group of blazars that exhibit unusually hard VHE spectra considering their redshifts. When correcting the measured spectra for absorption by the extragalactic background light, 1ES 1218+304 and a number of other blazars are found to have differential photon indices less than 1.5. The difficulty in modeling these hard spectral energy distributions in blazar jets has led to a range of theoretical gamma-ray emission scenarios, one of which is strongly constrained by these new VERITAS observations. We consider the implications of the observed light curve of 1ES 1218+304, which shows day scale flux variations, for shock acceleration scenarios in relativistic jets, and in particular for the viability of kiloparsec-scale jet emission scenarios.

VERITAS is a state-of-the-art ground-based gamma-ray observatory that operates in the very high-energy (VHE) region of 100 GeV to 50 TeV. The observatory consists of an array of four 12m-diameter imaging atmospheric Cherenkov telescopes located in southern Arizona, USA. The four-telescope array has been fully operational since September 2007, and over the last two years, VERITAS has been operating with high efficiency and with excellent performance. This talk summarizes the recent results from VERITAS, including the discovery of eight new VHE gamma-ray sources.

We report the detection of very-high-energy (VHE) gamma-ray emission from supernova remnant (SNR) G106.3+2.7. Observations performed in 2008 with the VERITAS atmospheric Cherenkov gamma-ray telescope resolve extended emission overlapping the elongated radio SNR. The 7.3 sigma (pre-trials) detection has a full angular extent of roughly 0.6deg by 0.4deg. Most notably, the centroid of the VHE emission is centered near the peak of the coincident 12CO (J = 1-0) emission, 0.4deg away from the pulsar PSR J2229+6114, situated at the northern end of the SNR. Evidently the current-epoch particles from the pulsar wind nebula are not participating in the gamma-ray production. The VHE energy spectrum measured with VERITAS is well characterized by a power law dN/dE = N_0(E/3 TeV)^{-G} with a differential index of G = 2.29 +/- 0.33stat +/- 0.30sys and a flux of N_0 = (1.15 +/- 0.27stat +/- 0.35sys)x 10^{-13} cm^{-2} s^{-1} TeV^{-1}. The integral flux above 1 TeV corresponds to ~5 percent of the steady Crab Nebula emission above the same energy. We describe the observations and analysis of the object and briefly discuss the implications of the detection in a multiwavelength context.

The recent detection by the Fermi gamma-ray space telescope of high-energy gamma-rays from the radio galaxy NGC 1275 makes the observation of the very high energy (VHE: E > 100 GeV) part of its broadband spectrum particularly interesting, especially for the understanding of active galactic nuclei (AGN) with misaligned multi-structured jets. The radio galaxy NGC 1275 was recently observed by VERITAS at energies above 100 GeV for about 8 hours. No VHE gamma-ray emission was detected by VERITAS from NGC 1275. A 99% confidence level upper limit of 2.1% of the Crab Nebula flux level is obtained at the decorrelation energy of approximately 340 GeV, corresponding to 19% of the power-law extrapolation of the Fermi Large Area Telescope (LAT) result.

Although Galactic cosmic rays (protons and nuclei) are widely believed to be dominantly accelerated by the winds and supernovae of massive stars, definitive evidence of this origin remains elusive nearly a century after their discovery [1]. The active regions of starburst galaxies have exceptionally high rates of star formation, and their large size, more than 50 times the diameter of similar Galactic regions, uniquely enables reliable calorimetric measurements of their potentially high cosmic-ray density [2]. The cosmic rays produced in the formation, life, and death of their massive stars are expected to eventually produce diffuse gamma-ray emission via their interactions with interstellar gas and radiation. M 82, the prototype small starburst galaxy, is predicted to be the brightest starburst galaxy in gamma rays [3, 4]. Here we report the detection of >700 GeV gamma rays from M 82. From these data we determine a cosmic-ray density of 250 eV cm-3 in the starburst core of M 82, or about 500 times the average Galactic density. This result strongly supports that cosmic-ray acceleration is tied to star formation activity, and that supernovae and massive-star winds are the dominant accelerators.

We present evidence that the very-high-energy (VHE, E > 100 GeV) gamma-ray emission coincident with the supernova remnant IC 443 is extended. IC 443 contains one of the best-studied sites of supernova remnant/molecular cloud interaction and the pulsar wind nebula CXOU J061705.3+222127, both of which are important targets for VHE observations. VERITAS observed IC 443 for 37.9 hours during 2007 and detected emission above 300 GeV with an excess of 247 events, resulting in a significance of 8.3 standard deviations (sigma) before trials and 7.5 sigma after trials in a point-source search. The emission is centered at 06 16 51 +22 30 11 (J2000) +- 0.03_stat +- 0.08_sys degrees, with an intrinsic extension of 0.16 +- 0.03_stat +- 0.04_sys degrees. The VHE spectrum is well fit by a power law (dN/dE = N_0 * (E/TeV)^-Gamma) with a photon index of 2.99 +- 0.38_stat +- 0.3_sys and an integral flux above 300 GeV of (4.63 +- 0.90_stat +- 0.93_sys) * 10^-12 cm^-2 s^-1. These results are discussed in the context of existing models for gamma-ray production in IC 443.

We present results from a long-term monitoring campaign on the TeV binary LSI +61 303 with VERITAS at energies above 500 GeV, and in the 2-10 keV hard X-ray bands with RXTE and Swift, sampling nine 26.5 day orbital cycles between September 2006 and February 2008. The binary was observed by VERITAS to be variable, with all integrated observations resulting in a detection at the 8.8 sigma (2006/2007) and 7.3 sigma (2007/2008) significance level for emission above 500 GeV. The source was detected during active periods with flux values ranging from 5 to 20% of the Crab Nebula, varying over the course of a single orbital cycle. Additionally, the observations conducted in the 2007-2008 observing season show marginal evidence (at the 3.6 sigma significance level) for TeV emission outside of the apastron passage of the compact object around the Be star. Contemporaneous hard X-ray observations with RXTE and Swift show large variability with flux values typically varying between 0.5 and 3.0*10^-11 ergs cm^-2 s^-1 over a single orbital cycle. The contemporaneous X-ray and TeV data are examined and it is shown that the TeV sampling is not dense enough to detect a correlation between the two bands.

The VERITAS collaboration reports the detection of very-high-energy (VHE) gamma-ray emission from the high-frequency-peaked BL Lac object 1ES 1218+304 located at a redshift of z=0.182. A gamma-ray signal was detected with a statistical significance of 10.4 standard deviations (10.4 sigma) for the observations taken during the first three months of 2007, confirming the discovery of this object made by the MAGIC collaboration. The photon spectrum between ~160 GeV and ~1.8 TeV is well described by a power law with an index of Gamma = 3.08 +/- 0.34_stat +/- 0.2_sys. The integral flux is Phi(E > 200 GeV) = (12.2 +/- 2.6) X 10^-12 cm^-2 s^-1, which corresponds to ~6% of that of the Crab Nebula. The light curve does not show any evidence for VHE flux variability. Using lower limits on the density of the extragalactic background light in the near to mid-infrared we are able to limit the range of intrinsic energy spectra for 1ES 1218+304. We show that the intrinsic photon spectrum has an index that is harder than Gamma = 2.32 +/- 0.37_stat. When including constraints from the spectra of 1ES 1101-232 and 1ES 0229+200, the spectrum of 1ES 1218+304 is likely to be harder than Gamma = 1.86 +/- 0.37_stat.

Since September 2005, the Whipple 10m Gamma-ray Telescope has been operated primarily as a blazar monitor. The five Northern Hemisphere blazars that have already been detected at the Whipple Observatory, Markarian 421, H1426+428, Markarian 501, 1ES 1959+650 and 1ES 2344+514, are monitored routinely each night that they are visible. We report on the Markarian 421 observations taken from November 2005 to June 2006 in the gamma-ray, X-ray, optical and radio bands. During this time, Markarian 421 was found to be variable at all wavelengths probed. Both the variability and the correlations among different energy regimes are studied in detail here. A tentative correlation, with large spread, was measured between the X-ray and gamma-ray bands, while no clear correlation was evident among the other energy bands. In addition to this, the well-sampled spectral energy distribution of Markarian 421 (1101+384) is presented for three different activity levels. The observations of the other blazar targets will be reported separately.

The high-frequency-peaked BL-Lacertae object \objectname{1ES 0806+524}, at redshift z=0.138, was observed in the very-high-energy (VHE) gamma-ray regime by VERITAS between November 2006 and April 2008. These data encompass the two-, and three-telescope commissioning phases, as well as observations with the full four-telescope array. \objectname{1ES 0806+524} is detected with a statistical significance of 6.3 standard deviations from 245 excess events. Little or no measurable variability on monthly time scales is found. The photon spectrum for the period November 2007 to April 2008 can be characterized by a power law with photon index $3.6 \pm 1.0_{\mathrm{stat}} \pm 0.3_{\mathrm{sys}}$ between $\sim$300 GeV and $\sim$700 GeV. The integral flux above 300 GeV is $(2.2\pm0.5_{\mathrm{stat}}\pm0.4_{\mathrm{sys}})\times10^{-12}\:\mathrm{cm}^{2}\:\mathrm{s}^{-1}$ which corresponds to 1.8% of the Crab Nebula flux. Non contemporaneous multiwavelength observations are combined with the VHE data to produce a broadband spectral energy distribution that can be reasonably described using a synchrotron-self Compton model.

VERITAS, an Imaging Atmospheric Cherenkov Telescope (IACT) system for gammma-ray astronomy in the GeV-TeV range, has recently completed its first season of observations with a full array of four telescopes. A number of astrophysical gamma-ray sources have been detected, both galactic and extragalactic, including sources previously unknown at TeV energies. We describe the status of the array and some highlight results, and assess the technical performance, sensitivity and shower reconstruction capabilities.

We present new measurements of heavy cosmic-ray nuclei at high energies per- formed during the first flight of the balloon-borne cosmic-ray experiment CREAM (Cosmic-Ray Energetics And Mass). This instrument uses multiple charge detectors and a transition radiation detector to provide the first high accuracy measurements of the relative abundances of elements from boron to oxygen up to energies around 1 TeV/n. The data agree with previous measurements at lower energies and show a relatively steep decline (~E$^-0.6$ to E$^-0.5$) at high energies. They further show the source abundance of nitrogen relative to oxygen is ~10% in the TeV/n region.

We report the detection of very high-energy gamma-ray emission from the intermediate-frequency-peaked BL Lacertae object W Comae (z=0.102) by VERITAS. The source was observed between January and April 2008. A strong outburst of gamma-ray emission was measured in the middle of March, lasting for only four days. The energy spectrum measured during the two highest flare nights is fit by a power-law and is found to be very steep, with a differential photon spectral index of Gamma = 3.81 +- 0.35_stat +- 0.34_syst. The integral photon flux above 200GeV during those two nights corresponds to roughly 9% of the flux from the Crab Nebula. Quasi-simultaneous Swift observations at X-ray energies were triggered by the VERITAS observations. The spectral energy distribution of the flare data can be described by synchrotron-self-Compton (SSC) or external-Compton (EC) leptonic jet models, with the latter offering a more natural set of parameters to fit the data.

LS I +61 303 is one of only a few high-mass X-ray binaries currently detected at high significance in very high energy gamma-rays. The system was observed over several orbital cycles (between September 2006 and February 2007) with the VERITAS array of imaging air-Cherenkov telescopes. A signal of gamma-rays with energies above 300 GeV is found with a statistical significance of 8.4 standard deviations. The detected flux is measured to be strongly variable; the maximum flux is found during most orbital cycles at apastron. The energy spectrum for the period of maximum emission can be characterized by a power law with a photon index of Gamma=2.40+-0.16_stat+-0.2_sys and a flux above 300 GeV corresponding to 15-20% of the flux from the Crab Nebula.

The multiwavelength observation of the nearby radio galaxy M87 provides a unique opportunity to study in detail processes occurring in Active Galactic Nuclei from radio waves to TeV gamma-rays. Here we report the detection of gamma-ray emission above 250 GeV from M87 in spring 2007 with the VERITAS atmospheric Cherenkov telescope array and discuss its correlation with the X-ray emission. The gamma-ray emission is measured to be point-like with an intrinsic source radius less than 4.5 arcmin. The differential energy spectrum is fitted well by a power-law function: dPhi/dE=(7.4+-1.3_{stat}+-1.5_{sys})(E/TeV)^{-2.31+-0.17_{stat}+-0.2_{sys}} 10^{-9}m^{-2}s^{-1}TeV^{-1}. We show strong evidence for a year-scale correlation between the gamma-ray flux reported by TeV experiments and the X-ray emission measured by the ASM/RXTE observatory, and discuss the possible short-time-scale variability. These results imply that the gamma-ray emission from M87 is more likely associated with the core of the galaxy than with other bright X-ray features in the jet.

The Track Imaging Cerenov Experiment (TrICE) is an air Cerenkov prototype telescope designed to use multi-anode photomultiplier to acheive a high angular resolution for measuring cosmic-ray composition at TeV-PeV energies. The TrICE camera, composed of 16 Hamamatsu R8900 16-channel multi-anode photomultiplier tubes, achieves 0.086 degree angular width per pixel over 1.5 degree wide field of view. We present a description of the TrICE camera design, calibration and performance.

Gamma-ray burst (GRB) observations at very high energies (VHE, E > 100 GeV) can impose tight constraints on some GRB emission models. Many GRB afterglow models predict a VHE component similar to that seen in blazars and plerions, in which the GRB spectral energy distribution has a double-peaked shape extending into the VHE regime. VHE emission coincident with delayed X-ray flare emission has also been predicted. GRB follow-up observations have had high priority in the observing program at the Whipple 10m Gamma-ray Telescope and GRBs will continue to be high priority targets as the next generation observatory, VERITAS, comes on-line. Upper limits on the VHE emission, at late times (>~4 hours), from seven GRBs observed with the Whipple Telescope are reported here.

We report on observations of the sky region around the unidentified TeV gamma-ray source TeV J2032+4130 carried out with the Whipple Observatory 10 m atmospheric Cherenkov telescope for a total of 65.5 hrs between 2003 and 2005. The standard two-dimensional analysis developed by the Whipple collaboration for a stand-alone telescope reveals an excess in the field of view at a pre-trials significance level of 6.1 standard deviations. The measured position of this excess is alpha(2000) =20 h 32 m 27 s, delta(2000) = 41 deg 39 min 17 s. The estimated integral flux for this gamma-ray source is about 8% of the Crab-Nebula flux. The data are consistent with a point-like source. Here we present a detailed description of the standard two-dimensional analysis technique used for the analysis of data taken with the Whipple Observatory 10 m telescope and the results for the TeV J2032+4130 campaign. We include a short discussion of the physical mechanisms that may be responsible for the observed gamma-ray emission, based on possible association with known astrophysical objects, in particular Cygnus OB2.

The first atmospheric Cherenkov telescope of VERITAS (the Very Energetic Radiation Imaging Telescope Array System) has been in operation since February 2005. We present here a technical description of the instrument and a summary of its performance. The calibration methods are described, along with the results of Monte Carlo simulations of the telescope and comparisons between real and simulated data. The analysis of TeV $\gamma$-ray observations of the Crab Nebula, including the reconstructed energy spectrum, is shown to give results consistent with earlier measurements. The telescope is operating as expected and has met or exceeded all design specifications.

Galaxy clusters might be sources of TeV gamma rays emitted by high-energy protons and electrons accelerated by large scale structure formation shocks, galactic winds, or active galactic nuclei. Furthermore, gamma rays may be produced in dark matter particle annihilation processes at the cluster cores. We report on observations of the galaxy clusters Perseus and Abell 2029 using the 10 m Whipple Cherenkov telescope during the 2003-2004 and 2004-2005 observing seasons. We apply a two-dimensional analysis technique to scrutinize the clusters for TeV emission. In this paper we first determine flux upper limits on TeV gamma-ray emission from point sources within the clusters. Second, we derive upper limits on the extended cluster emission. We subsequently compare the flux upper limits with EGRET upper limits at 100 MeV and theoretical models. Assuming that the gamma-ray surface brightness profile mimics that of the thermal X-ray emission and that the spectrum of cluster cosmic rays extends all the way from thermal energies to multi-TeV energies with a differential spectral index of -2.1, our results imply that the cosmic ray proton energy density is less than 7.9% of the thermal energy density for the Perseus cluster.

We present a new measurement of the cosmic-ray positron fraction at energies between 5 and 15 GeV with the balloon-borne HEAT-pbar instrument in the spring of 2000. The data presented here are compatible with our previous measurements, obtained with a different instrument. The combined data from the three HEAT flights indicate a small positron flux of non-standard origin above 5 GeV. We compare the new measurement with earlier data obtained with the HEAT-e+- instrument, during the opposite epoch of the solar cycle, and conclude that our measurements do not support predictions of charge sign dependent solar modulation of the positron abundance at 5 GeV.

We present a new measurement of atmospheric muons made during an ascent of the High Energy Antimatter Telescope balloon experiment. The muon charge ratio mu+/mu- as a function of atmospheric depth in the momentum interval 0.3-0.9 GeV/c is presented. The differential mu- intensities in the 0.3-50 GeV/c range and for atmospheric depths between 4-960 g/cm^2 are also presented. We compare these results with other measurements and model predictions. We find that our charge ratio is ~1.1 for all atmospheric depths and is consistent, within errors, with other measurements and the model predictions. We find that our measured mu- intensities are also consistent with other measurements, and with the model predictions, except at shallow atmospheric depths.

A decade after the discovery of TeV gamma-rays from the blazar Mrk 421 (Punch et al. 1992), the list of TeV blazars has increased to five BL Lac objects: Mrk 421 (Punch et al. 1992; Petry et al. 1996; Piron et al. 2001), Mrk 501 (Quinn et al. 1996; Aharonian et al. 1999; Djannati-Atai et al. 1999), 1ES2344+514 (Catanese et al. 1998), H1426+428 (Horan et al. 2000, 2002; Aharonian et al. 2002; Djannati-Atai et al. 2002) and 1ES1959+650 (Nishiyama et al. 1999; Konopelko et al. 2002; Holder et al. 2002). In this paper we report results from recent observations of Mrk 421, H1426+428 and 1ES1959+650 using the Whipple Observatory 10 m telescope.

The BL Lac object H1426+428 was recently detected as a high energy gamma-ray source by the VERITAS collaboration (Horan et al. 2002). We have reanalyzed the 2001 portion of the data used in the detection in order to examine the spectrum of H1426+428 above 250 GeV. We find that the time-averaged spectrum agrees with a power law of the shape dF/dE = 10^(-7.31 +- 0.15(stat) +- 0.16(syst)) x E^(-3.50 +- 0.35(stat) +- 0.05(syst)) m^(-2)s^(-1)TeV^(-1) The statistical evidence from our data for emission above 2.5 TeV is 2.6 sigma. With 95% c.l., the integral flux of H1426+428 above 2.5 TeV is larger than 3% of the corresponding flux from the Crab Nebula. The spectrum is consistent with the (non-contemporaneous) measurement by Aharonian et al. (2002) both in shape and in normalization. Below 800 GeV, the data clearly favours a spectrum steeper than that of any other TeV Blazar observed so far indicating a difference in the processes involved either at the source or in the intervening space.

The detection of spectral variability of the gamma-ray blazar Mrk 421 at TeV energies is reported. Observations with the Whipple Observatory 10m gamma-ray telescope taken in 2000/2001 revealed exceptionally strong and long-lasting flaring activity. Flaring levels of 0.4 to 13 times that of the Crab Nebula flux provided sufficient statistics for a detailed study of the energy spectrum between 380 GeV and 8.2 TeV as a function of flux level. These spectra are well described by a power law with an exponential cutoff. There is no evidence for variation in the cutoff energy with flux, and all spectra are consistent with an average value for the cutoff energy of 4.3 TeV. The spectral index varies between 1.89 in a high flux state and 2.72 in a low state. The correlation between spectral index and flux is tight when averaging over the total 2000/2001 data set. Spectral measurements of Mrk~421 from previous years (1995/96 and 1999) by the Whipple collaboration are consistent with this flux-spectral index correlation, which suggest that this may be a constant or a long-term property of the source. If a similar flux-spectral index correlation were found for other gamma-ray blazars, this universal property could help disentangle the intrinsic emission mechanism from external absorption effects.

The observation of a small change in spectral slope, or 'knee' in the fluxes of cosmic rays near energies 10^15 eV has caused much speculation since its discovery over 40 years ago. The origin of this feature remains unknown. A small workshop to review some modern experimental measurements of this region was held at the Adler Planetarium in Chicago, USA in June 2000. This paper summarizes the results presented at this workshop and the discussion of their interpretation in the context of hadronic models of atmospheric airshowers.

We present a new measurement of the antiproton to proton abundance ratio, pbar/p, in the cosmic radiation. The HEAT-pbar instrument, a balloon borne magnet spectrometer with precise rigidity and multiple energy loss measurement capability, was flown successfully in Spring 2000, at an average atmospheric depth of 7.2 g/cm^2. A total of 71 antiprotons were identified above the vertical geomagnetic cut-off rigidity of 4.2 GV. The highest measured proton energy was 81 GeV. We find that the pbar/p abundance ratio agrees with that expected from a purely secondary origin of antiprotons produced by primary protons with a standard soft energy spectrum.

Exceptionally strong and long lasting flaring activity of the blazar Markarian 421 (Mrk 421) occurred between January and March 2001. Based on the excellent signal-to-noise ratio of the data we derive the energy spectrum between 260 GeV - 17 TeV with unprecedented statistical precision. The spectrum is not well described by a simple power law even with a curvature term. Instead the data can be described by a power law with exponential cutoff: $\rm {{dN}\over{dE}} \propto \rm E^{-2.14 \pm 0.03_{stat}} \times e^{-E/E_{0}} m^{-2} s^{-1} TeV^{-1}$ with $\rm E_{0} = 4.3 \pm 0.3_{stat} TeV$. Mrk 421 is the second $\gamma$-ray blazar that unambiguously exhibits an absorption-like feature in its spectral energy distribution at 3-6 TeV suggesting that this may be a universal phenomenon, possibly due to the extragalactic infra-red background radiation.

The accurate determination of the elemental composition of cosmic rays at high energies is expected to provide crucial clues on the origin of these particles. Previous direct measurements of composition have been limited by experiment collecting power, resulting in marginal statistics above $10^{14}$ eV, precisely the region where the ``knee'' of the cosmic-ray energy spectrum is starting to develop. In contrast, indirect measurements using extensive air showers can produce sufficient statistics in this region but generate elemental measurements which have relatively large uncertainties. Here we discuss a technique which has become possible through the use of modern ground-based Cerenkov imaging detectors. We combine a measurement of the Cerenkov light produced by the incoming cosmic-ray nucleus in the upper atmosphere with an estimate of the total nucleus energy produced by the extensive air shower initiated when the particle interacts deeper in the atmosphere. The emission regions prior to and after the first hadronic interaction can be separated by an imaging Cerenkov system with sufficient angular and temporal resolution. Monte Carlo simulations indicate an expected charge resolution of $\Delta Z/Z <5%$ for incident iron nuclei in the region of the ``knee'' of the cosmic-ray energy spectrum. This technique also has the intriguing possibility to unambiguously discover nuclei heavier than iron at energies above 10$^{14}$ eV. The identification and rejection of background produced by charged particles in ground based gamma-ray telescopes is also discussed.

The small change in the spectral slope of the overall intensity of cosmic rays near 1 PeV may be associated with the endpoint energy of supernova shock acceleration. A crucial test of this connection and other ideas of the origin of the spectral `knee' is the reliable determination of the variation of elemental composition in this region. Recent measurements at the DICE/CASA/MIA air shower installation in Dugway, Utah, USA have provided several independent air shower parameters for each event. The derivation of elemental composition from a combination of Cherenkov size, depth of shower maximum in the atmosphere, muon size and electron size at ground level and the reliability of these results are discussed. There is no evidence from these data for a large change in the mean mass of cosmic rays across the `knee'. These measurements show the cosmic rays are composed of ~ 70% protons and alpha-particles near total energies of 10PeV.

We give an overview of the current status and scientific goals of VERITAS, a proposed hexagonal array of seven 10 m aperture imaging Cherenkov telescopes. The selected site is Montosa Canyon (1390 m a.s.l.) at the Whipple Observatory, Arizona. Each telescope, of 12 m focal length, will initially be equipped with a 499 element photomultiplier camera covering a 3.5 degree field of view. A central station will initiate the readout of 500 MHz FADCs upon receipt of multiple telescope triggers. The minimum detectable flux sensitivity will be 0.5% of the Crab Nebula flux at 200 GeV. Detailed simulations of the array's performance are presented elsewhere at this meeting. VERITAS will operate primarily as a gamma-ray observatory in the 50 GeV to 50 TeV range for the study of active galaxies, supernova remnants, pulsars and gamma-ray bursts.

Cosmic rays at the Earth include a secondary component originating in collisions of primary particles with the diffuse interstellar gas. The secondary cosmic rays are relatively rare but carry important information on the Galactic propagation of the primary particles. The secondary component includes a small fraction of antimatter particles, positrons and antiprotons. In addition, positrons and antiprotons may also come from unusual sources and possibly provide insight into new physics. For instance, the annihilation of heavy supersymmetric dark matter particles within the Galactic halo could lead to positrons or antiprotons with distinctive energy signatures. With the High-Energy Antimatter Telescope (HEAT) balloon-borne instrument, we have measured the abundances of positrons and electrons at energies between 1 and 50 GeV. The data suggest that indeed a small additional antimatter component may be present that cannot be explained by a purely secondary production mechanism. Here we describe the signature of the effect and discuss its possible origin.

Observations of cosmic-ray electrons and positrons have been made with a new balloon-borne detector, HEAT (the "High-Energy Antimatter Telescope"), first flown in 1994 May from Fort Sumner, NM. We describe the instrumental approach and the data analysis procedures, and we present results from this flight. The measurement has provided a new determination of the individual energy spectra of electrons and positrons from 5 GeV to about 50 GeV, and of the combined "all-electron" intensity (e+ + e-) up to about 100 GeV. The single power-law spectral indices for electrons and positrons are alpha = 3.09 +/- 0.08 and 3.3 +/- 0.2, respectively. We find that a contribution from primary sources to the positron intensity in this energy region, if it exists, must be quite small.

The Dual Imaging Cerenkov Experiment (DICE) was designed and operated for making elemental composition measurements of cosmic rays near the knee of the spectrum at several PeV. Here we present the first results using this experiment from the measurement of the average location of the depth of shower maximum, <X_max>, in the atmosphere as a function of particle energy. The value of <X_max> near the instrument threshold of ~0.1 PeV is consistent with expectations from previous direct measurements. At higher energies there is little change in composition up to ~5 PeV. Above this energy <X_max> is deeper than expected for a constant elemental composition implying the overall elemental composition is becoming lighter above the knee region. These results disagree with the idea that cosmic rays should become on average heavier above the knee. Instead they suggest a transition to a qualitatively different population of particles above 5 PeV.

A next generation atmospheric Cherenkov observatory is described based on the Whipple Observatory $\gamma$-ray telescope. A total of nine such imaging telescopes will be deployed in an array that will permit the maximum versatility and give high sensitivity in the 50 GeV - 50 TeV band (with maximum sensitivity from 100 GeV to 10 TeV).

Two measurements of the cosmic-ray positron fraction as a function of energy have been made using the High Energy Antimatter Telescope (HEAT) balloon-borne instrument. The first flight took place from Ft. Sumner, New Mexico in 1994, and yielded results above the geomagnetic cutoff energy of 4.5 GeV. The second flight from Lynn Lake, Manitoba in 1995 permitted measurements over a larger energy interval, from 1 GeV to 50 GeV. In this letter we present results on the positron fraction based on data from the Lynn Lake flight, and compare these with the previously published results from the Ft. Sumner flight. The results confirm that the positron fraction does not increase with energy above ~10 GeV, although a small excess above purely secondary production cannot be ruled out. At low energies the positron fraction is slightly larger than that reported from measurements made in the 1960's. This effect could possibly be a consequence of charge dependence in the level of solar modulation.