Results are presented of a harmonic analysis of the large scale cosmic-ray anisotropy as observed by the Milagro observatory. We show a two-dimensional display of the sidereal anisotropy pro jections in right ascension generated by the fitting of three harmonics to 18 separate declination bands. The Milagro observatory is a water Cherenkov detector located in the Jemez mountains near Los Alamos, New Mexico. With a high duty cycle and large field-of-view, Milagro is an excellent instrument for measuring this anisotropy with high sensitivity at TeV energies. The analysis is conducted using a seven year data sample consisting of more than 95 billion events, the largest such data set in existence. We observe an anisotropy with a magnitude around 0.1% for cosmic rays with a median energy of 6 TeV. The dominant feature is a deficit region of depth (2.49 +/- 0.02 stat. +/- 0.09 sys.)x10^(-3) in the direction of the Galactic North Pole centered at 189 degrees right ascension. We observe a steady increase in the magnitude of the signal over seven years.

An analysis of 7 years of Milagro data performed on a 10-degree angular scale has found two localized regions of excess of unknown origin with greater than 12 sigma significance. Both regions are inconsistent with gamma-ray emission with high confidence. One of the regions has a different energy spectrum than the isotropic cosmic-ray flux at a level of 4.6 sigma, and it is consistent with hard spectrum protons with an exponential cutoff, with the most significant excess at ~10 TeV. Potential causes of these excesses are explored, but no compelling explanations are found.

Diffuse $\gamma$-ray emission produced by the interaction of cosmic-ray particles with matter and radiation in the Galaxy can be used to probe the distribution of cosmic rays and their sources in different regions of the Galaxy. With its large field of view and long observation time, the Milagro Gamma Ray Observatory is an ideal instrument for surveying large regions of the Northern Hemisphere sky and for detecting diffuse $\gamma$-ray emission at very high energies. Here, the spatial distribution and the flux of the diffuse $\gamma$-ray emission in the TeV energy range with a median energy of 15 TeV for Galactic longitudes between 30$^\circ$ and 110$^\circ$ and between 136$^\circ$ and 216$^\circ$ and for Galactic latitudes between -10$^\circ$ and 10$^\circ$ are determined. The measured fluxes are consistent with predictions of the GALPROP model everywhere except for the Cygnus region ($l\in[65^\circ,85^\circ]$). For the Cygnus region, the flux is twice the predicted value. This excess can be explained by the presence of active cosmic ray sources accelerating hadrons which interact with the local dense interstellar medium and produce gamma rays through pion decay.

Recent rapid localizations of short, hard gamma-ray bursts (GRBs) by the Swift and HETE satellites have led to the observation of the first afterglows and the measurement of the first redshifts from this type of burst. Detection of >100 GeV counterparts would place powerful constraints on GRB mechanisms. Seventeen short duration (< 5 s) GRBs detected by satellites occurred within the field of view of the Milagro gamma-ray observatory between 2000 January and 2006 December. We have searched the Milagro data for >100 GeV counterparts to these GRBs and find no significant emission correlated with these bursts. Due to the absorption of high-energy gamma rays by the extragalactic background light (EBL), detections are only expected for redshifts less than ~0.5. While most long duration GRBs occur at redshifts higher than 0.5, the opposite is thought to be true of short GRBs. Lack of a detected VHE signal thus allows setting meaningful fluence limits. One GRB in the sample (050509b) has a likely association with a galaxy at a redshift of 0.225, while another (051103) has been tentatively linked to the nearby galaxy M81. Fluence limits are corrected for EBL absorption, either using the known measured redshift, or computing the corresponding absorption for a redshift of 0.1 and 0.5, as well as for the case of z=0.

A survey of Galactic gamma-ray sources at a median energy of ~20 TeV has been performed using the Milagro Gamma Ray Observatory. Eight candidate sources of TeV emission are detected with pre-trials significance $>4.5\sigma$ in the region of Galactic longitude $l\in[30^\circ,220^\circ]$ and latitude $b\in[-10^\circ,10^\circ]$. Four of these sources, including the Crab nebula and the recently published MGRO J2019+37, are observed with significances $>4\sigma$ after accounting for the trials involved in searching the 3800 square degree region. All four of these sources are also coincident with EGRET sources. Two of the lower significance sources are coincident with EGRET sources and one of these sources is Geminga. The other two candidates are in the Cygnus region of the Galaxy. Several of the sources appear to be spatially extended. The fluxes of the sources at 20 TeV range from ~25% of the Crab flux to nearly as bright as the Crab.

The diffuse gamma radiation arising from the interaction of cosmic ray particles with matter and radiation in the Galaxy is one of the few probes available to study the origin of the cosmic rays. Milagro is a water Cherenkov detector that continuously views the entire overhead sky. The large field-of-view combined with the long observation time makes Milagro the most sensitive instrument available for the study of large, low surface brightness sources such as the diffuse gamma radiation arising from interactions of cosmic radiation with interstellar matter. In this paper we present spatial and flux measurements of TeV gamma-ray emission from the Cygnus Region. The TeV image shows at least one new source MGRO J2019+37 as well as correlations with the matter density in the region as would be expected from cosmic-ray proton interactions. However, the TeV gamma-ray flux as measured at ~12 TeV from the Cygnus region (after excluding MGRO J2019+37) exceeds that predicted from a conventional model of cosmic ray production and propagation. This observation indicates the existence of either hard-spectrum cosmic-ray sources and/or other sources of TeV gamma rays in the region.

The Milagro gamma-ray observatory employs a water Cherenkov detector to observe extensive air showers produced by high energy particles interacting in the Earth's atmosphere. Milagro has a wide field of view and high duty cycle, monitoring the northern sky almost continuously in the 100 GeV to 100 TeV energy range. Milagro is, thus, uniquely capable of searching for very high-energy emission from gamma-ray bursts (GRBs) during the prompt emission phase. Detection of >100 GeV counterparts would place powerful constraints on GRB mechanisms. Twenty-five satellite-triggered GRBs occurred within the field of view of Milagro between January 2000 and December 2001. We have searched for counterparts to these GRBs and found no significant emission from any of the burst positions. Due to the absorption of high-energy gamma rays by the extragalactic background light, detections are only expected to be possible for redshifts less than ~0.5. Three of the GRBs studied have measured redshifts. GRB 010921 has a redshift low enough (0.45) to allow an upper limit on the fluence to place an observational constraint on potential GRB models.

The Milagro telescope monitors the northern sky for 100 GeV to 100 TeV transient emission through continuous very high energy wide-field observations. The large effective area and ~100 GeV energy threshold of Milagro allow it to detect very high energy (VHE) gamma-ray burst emission with much higher sensitivity than previous instruments and a fluence sensitivity at VHE energies comparable to that of dedicated gamma-ray burst satellites at keV to MeV energies. Even in the absence of a positive detection, VHE observations can place important constraints on gamma-ray burst (GRB) progenitor and emission models. We present limits on the VHE flux of 40 s -- 3 h duration transients nearby to earth, as well as sensitivity distributions which have been corrected for gamma-ray absorption by extragalactic background light and cosmological effects. The sensitivity distributions suggest that the typical intrinsic VHE fluence of GRBs is similar or weaker than the keV -- MeV emission, and we demonstrate how these sensitivity distributions may be used to place observational constraints on the absolute VHE luminosity of gamma-ray bursts for any GRB emission and progenitor model.